Questions tagged [finite-difference]
Referring to the discretization of derivatives by Finite differences, and its applications to numerical solutions of partial differential equations.
791
questions
1
vote
0
answers
60
views
Can the standard multigrid performance be used for time-dependent PDEs?
Consider a time dependent pde(i.e u(x,t)).I know when only space-coarsening is used the standard multigrid performance can be applied but what if instead we use only time-coarsening?Can we apply the ...
- 471
1
vote
1
answer
173
views
Simulating advection - diffusion problem in a network of 1D pipe
I'm interested in solving the following advection-diffusion system in a 1D network of pipes.
$$\frac{\partial C}{\partial t} = D\frac{\partial ^2 C}{\partial x^2} - v\frac{\partial C}{\partial x}$$
...
- 447
1
vote
0
answers
85
views
Solving a system of PDEs with no-flux boundary conditions (finite difference discretization)
I am interested in solving a system of linear PDEs with the finite difference method and I'm having trouble to solve the no-flux boundary condition correctly.
\begin{align}
\frac{\partial n}{\partial ...
- 11
1
vote
0
answers
83
views
Solving diffusion equation using finite difference method
I am solving an 1-dimensional diffusion equation with Neumann boundary condition at outlet and constant concentration, C, at the inlet. In the end, I want to observe how the concentration diffuses ...
- 447
3
votes
3
answers
2k
views
Lattice Boltzmann methods vs Navier stokes/ other eulerian methods for *water* simulation
Note, there is already a question here, however the answers don't answer the original question, let alone specific considerations when dealing with nearly in-compressible fluids (water). Another ...
- 151
2
votes
1
answer
88
views
Stability of a finite-difference scheme for the reaction-diffusion equation
I currently need to solve numerically the following reaction-diffusion equation:
$$\partial_tu=\partial^2_xu+u-u^2$$
For this purpose, I use the following numerical scheme (Crank-Nicolson??):
$$ \...
- 129
1
vote
1
answer
75
views
Should the derivative of an array be calculated array by array or element by element in CFD codes?
I am making my own finite difference computational magnetohydrodynamic code in Fortran 90. Looking at other codes they appear to calculate for example their $x$-derivatives, bb of their variables, e.g....
- 219
2
votes
1
answer
293
views
1-D boundary value problem: How implement mixed boundary conditions using a FD method?
I have been given a convection-diffusion ODE modeling the steady state temperature of a pipe (through which flows a fluid) as
$$-\frac{d}{dz}\left(\kappa \frac{dT}{dz} \right)+v\rho C\frac{dT}{dz}=Q(...
- 207
2
votes
1
answer
135
views
Confusion about Zabusky and Kruskal's stepper for the KdV equation
In Zabusky and Kruskal's paper about solitons, they derive the following update for the Korteweg de Vries equation (their footnote 6):
\begin{align*}
u_{i}^{j+1} = u_{i}^{j-1} - \frac{1}{3} \frac{k}{...
- 2,031
1
vote
0
answers
36
views
Finite difference/element method : time step and spatial resolution close to a finite singularity
I'm using the finite element method (FEM), but my question is quite a global question. It's related to this question but it is not the same.
Let's assume we have this equation : $$\partial_t c - u\...
- 171
3
votes
0
answers
98
views
numerical instabilities in Fluid Dynamics, Finite Element Method
I'm looking for references to understand where the numerical instabilities come from in hydrodynamics in general, and notably when the Péclet number: $Pe>1$. I'm using the finite element method.
...
- 171
1
vote
0
answers
46
views
Predictor-Corrector vs. Deferred Difference Corrections
I want to use the Numerov method but keep higher-order terms from the Taylor expansions. In the literature, I found the term "Deferred Difference Corrections" for the procedure of first solving the ...
- 11
3
votes
1
answer
232
views
Why do many people use FDM method to solve Stokes equations, i.e., saddle point matrix?
For numerical methods of the Stokes equations, with appropriate boundary:
$$-\nabla^{2} \vec{u}+\nabla p=\overrightarrow{0}$$
$$\nabla \cdot \vec{u}=0$$
one may use FDM (finite difference method) ...
- 971
0
votes
1
answer
1k
views
Implementing structured grid boundary conditions using NumPy arrays?
I am making a toy code in Python to solve the advection equation $$u_t + cu_x = 0$$
with, for example, periodic boundary conditions.
Background information
The numerical grid is specified like this:
...
2
votes
0
answers
231
views
Implementation of boundary conditions for 1D Euler equations
I'm trying to solve 1D Euler equations with gravity in spherical coordinates using a finite-difference TVD MacCormack method on a non-uniform grid of $N$ components, following the method provided in ...
2
votes
0
answers
46
views
Verification on pressure predictor method for CFD code
I have developed a python code for a lid-drive cavity model. However, my results are not converging. The algorithm of my code looks like this:
Euler Momentum Equation looks like this:
$$\frac{u^{n+1}...
- 123
1
vote
0
answers
144
views
FDM discretization of equation on the boundary
In order to simulate the following equation using FDM
$$u_t(t,x)-u_{xx}(t,x)=0, \quad (t,x) \in (0,1)\times (0,1)$$
$$(u_t(t,x)-u_{x}(t,x))\rvert_{x=0}=0, \quad t \in (0,1)$$
$$(u_t(t,x)+u_{x}(t,x))\...
- 111
1
vote
0
answers
96
views
WENO scheme on curvilinear coordinates
I've been developing a curvilinear FVM code. So far I've implemented the PPM scheme and am looking into adding WENO schemes. So far I've been discretizing the grid metrics using a second-order central....
- 11
0
votes
1
answer
652
views
Need an example Legendre-Gauss-Radau pseudospectral differentiation matrix or Matlab code
I'm trying to implement various kinds of pseudospectral methods for direct optimization in Matlab using IPOPT. I've got some working Legendre-Gauss-Lobatto code, but would like to use the flipped ...
- 133
2
votes
1
answer
364
views
Solving nonlinear PDE with finite difference based on Newton-Krylov
I am now working on solving MHD equations with finite difference method, which include nonlinear equations:
$$
\frac{\partial\rho}{\partial t}+\nabla\cdot\left[\left(\rho_0+\rho\right){v}\right]-\...
- 21
2
votes
1
answer
108
views
What method of Finite difference is this?
I am reviewing Numerical Recipes method on solving ODEs via relaxation (Chapter 18.3 in the 3rd edition) and they chose a finite difference method I am unfamiliar with (Equation 18.3.2):
\begin{...
- 135
2
votes
1
answer
154
views
Interpolating the gradient of a cylindrically symmetric potential field expected obey the Laplace equation, especially near/across r=0
The script below tries to implement a Jacobi iterative relaxation of a potential field for an electrostatic lens.
In order to plot electric field lines and calculate trajectories for charged particles,...
- 872
1
vote
0
answers
47
views
Dealing with spurious oscillations in particle tracking methods
I work on modelling high intensity discharge xenon-filled lamps.
The model governing the discharge is quite complex and sadly includes fluid dynamics.
After some time, I managed to implement a finite-...
- 131
1
vote
1
answer
108
views
Finite difference methods
I am currently applying the finite difference method to the solution of the diffusion equation.
I think that a problem has occurred, and is as follows, my explicit method is the most accurate when ...
- 33
0
votes
0
answers
54
views
Neumann boundary conditions on arbitrary surface for finite difference diffusion
I am facing the following problem, formulated in practical terms: I have a region $\Omega$ in two or three dimensions, represented as a binary mask, and an initial density $u_0$ within that region ...
- 101
1
vote
2
answers
271
views
Second derivative using Fornberg finite difference method
I have some discrete data, non-equispaced in $x$, $y=f(x)$.
I want to use a numerical finite difference method to calculate the second derivatives of $y$, at some point.
I am using the Fornberg ...
- 201
2
votes
0
answers
82
views
How to determine the order of convergence of the Euler-Maruyama method?
This question is originally posted in Quant.StackExchange but has been unanswered for some time so I ask in here.
To make this simple let us consider the Geometric Brownian Motions (GBM).
My ...
- 255
3
votes
2
answers
120
views
Finite difference for a highly nonlinear equation - The wind within the forest
Based on the Navier-Stokes equations and a few parameterizations, the horizontal steady-state wind $u(z)$ within a forest of height H satisfies:
$$
a\left(\frac{du}{dz}\right)^2 + b\frac{du}{dz} \...
- 33
1
vote
1
answer
546
views
Dirichlet boundary conditions in the 1D Heat Equation
Please consider the assignment I have uploaded on the picture.
I am confused about the functions $g_L(t)$, $g_R(t)$ and $\eta(x)$. What are they and how do I find them... My question:
Is it possbile ...
- 255
3
votes
2
answers
206
views
Damped Harmonic Oscillation. Efficient algorithm to find the parameters resulting in threshold oscillation amplitude
Let's assume, that we have damped harmonic oscillation of a body in the form of a cone, immersed in a liquid.
Equilibrium condition of the body is:
$$m\overrightarrow{a} = \overrightarrow{F_\text{...
0
votes
1
answer
301
views
Numerically solving a non-linear PDE
I have this non-linear partial differential equation.
$$
\frac{\partial C}{\partial t}=\left(\frac{\partial C}{\partial x}\right)^2+C\frac{\partial^2 C}{\partial x^2}
$$
I want to use the finite ...
- 33
0
votes
0
answers
156
views
Double mach reflection at a inclined wedge
I am running into a strange problem when solving the 2D compressible Euler equations on a inclined wedge. To elaborate, my top boundary condition seems to emitting some type of instability. I have ...
- 75
0
votes
1
answer
103
views
Numerical Stability of a Generalized Spatial Discretization Scheme
After reading the matrix stability chapter (10) of Hirsch [1], I decided to dive in the reference list of the chapter. One of the papers [2], which is cited as reference shows an very interesting ...
- 77
2
votes
2
answers
110
views
Stability of Crank-Nicolson for $u_t = iu_{xx}+2iu$
I want to use the Crank-Nicolson scheme to solve the equation
$$u_t = iu_{xx}+2iu$$
Here's the analysis: Suppose we make a grid, with $k = dt$ and $h = dx$, the usual notation, and also $u_j^n = u(...
- 139
1
vote
0
answers
40
views
Boundary conditions for a Non-linear Schrödinger equation using an extended crank nicolson scheme
I try to solve numerically the following PDE for $E(r, z)$ with a cylindrical symmetrie (i. e. $E(r, z) = E(-r, z)$).
$\frac{\partial E}{\partial z} = \frac{i}{2k} \Delta E + \mathcal{N}(E)$
Where $...
- 11
5
votes
0
answers
85
views
Numerical methods for the continuity equation with Sobolev vector field
Consider the continuity equation
$$
\partial_t \rho(x,t) + \operatorname{div} (b(x,t) \rho(x,t)) = 0, \qquad t \in [0,T], \quad x \in \mathbb R^N,
$$
with $b \in L^1((0,T), W^{1,p}(\mathbb R^N))$.
...
- 51
1
vote
1
answer
151
views
How to obtain linear tridiagonal system from PDE
I'm trying to re-solve the governing equations in hydraulic fracturing modeling as instructed step by step in a paper.
After (A-9), the author stated that by substituting A-6, A-8 and A-9 into ...
- 11
1
vote
2
answers
936
views
CFL equation for non-linear equation
I am trying to solve numerically (obviously) inviscid Burgers' equation with the finite difference method. The equation is the following:
$$ \displaystyle \partial_t u + u \, \partial_x u = 0 $$
...
user30503
0
votes
1
answer
219
views
Good C, C++ library for efficient grid search / tuples, ideally with bindings to Eigen
I have a $q$-dimensional grid, known at run, not compile-time, that has $50$ points in each direction and hence $50^3$ combinations that I would like to first build and then call a function with each ...
- 183
3
votes
0
answers
65
views
Use of non-typical values of $\theta$ in theta-methods
The theta-method is a popular solution for solving time-transient PDEs (or ODEs), which consists of solving the general equation for each time step:
$$
\frac{u^{n+1} - u^{n}}{\Delta t} + (\theta f(u^{...
- 153
1
vote
0
answers
449
views
High-accuracy numerical differentiation
I have a $200 \times 200$ matrix representing the values taken by a function over an equally spaced grid. I would like to perform derivatives on it.
I am interested in its gradient (i.e. its ...
- 203
0
votes
1
answer
109
views
Error for the finite differences scheme -- Advection equation
Consider the advection equation (1D in space)
$$
\frac{\partial u}{\partial t} + V\, \frac{\partial u}{\partial x}=0
$$
and we solve it numerically on $[0,1]\times [0,1]\ni (t,x)$ using a forward ...
- 259
1
vote
1
answer
270
views
Approximation Error in a Finite Difference Approximation of the Square of Derivative
First Part: (First-order derivative)
Assuming $f$ is an infinitely differential function everywhere, the Taylor series of $f(x + h)$ at $x$ is
\begin{align}\tag{1}
f(x + h) = f(x) + hf'(x) + \frac{1}...
- 95
1
vote
0
answers
72
views
Oscillations when solving parabolic heat equation with FTCS
I'm wondering if someone could help me out, or point me in a direction of how I can understand the following oscillations that occur when I solve the Porous Medium Equation $$u_t = u_{xx}^{m+1}$$ ...
- 111
3
votes
0
answers
100
views
Numerical solution to N-dimensional diffusion on simplex?
Assume I have a system of at least (but generally only) $N+1$ points in an $N$-dimensional space ($N > 3$ is possible). At each of these points $x_i, i=1,...,N+1$ I know an initial potential/...
- 203
0
votes
2
answers
510
views
Non-linear Boundary Value Problem. How to compute the Jacobian?
Consider a Boundary Value Problem:
$$
\delta u''+u(u'-1) =0 \Leftrightarrow u''=\frac{-u(u'-1)}{\delta}=:f(t,u',u), \\
u(0)=a, u(1)=b
$$
$\delta,a,b$ are known parameters. I want to implement Newton'...
- 146
1
vote
3
answers
992
views
How to make a less diffusive code to solve 2D advection equation?
I would like to solve the following differential equation numerically in 2D,
$$\frac{\partial z^-}{\partial t}+(\vec{B}\cdot\vec{\nabla})z^-=0,$$
see Wikipedia if you are curious about what the ...
- 219
0
votes
1
answer
149
views
Actual global error vs theoretical global error: How to combine theory with practice
I have implemented an Adams Bashforth 4 method to solve an Initial Value Problem for an ODE and I am testing it against the test equation:
$y'=\lambda y$ with $y(0)=1$ with the exact solution: $y(t)=...
- 255
0
votes
1
answer
1k
views
Discretization Neumann boundary condition
I'm currently working with the following Poisson equation with mixed boundary conditions, including a Neumann boundary condition.
$$\Delta u = f\\ u(x,0) = g_1(x), 0<x<1\\u(0,y) = g_2(y), 0<...
- 216
2
votes
1
answer
270
views
Numerical integration of Fokker-Planck equation allowing for negative drift?
The Fokker-Planck equation (a.k.a Kolmogorov forward equation or Smoluchowski equation) describes the evolution of a probability density function and numerical integration of the FPE should conserve ...
- 123