Questions tagged [crank-nicolson]
For questions about the Crank-Nicolson method, an approach for discretizing and solving partial differential equations.
53
questions
2
votes
1
answer
278
views
Asking advice for implementation of Conservative Finite Difference Scheme for numerically solving Gross-Pitaevskii equation
I am trying to numerically solve the Gross-Pitaevskii equation for an impurity coupled with a one-dimensional weakly-interacting bosonic bath, given by (in dimensionless units):
\begin{align}
i \frac{\...
1
vote
0
answers
42
views
How can I apply a mixed boundary condition to a multi-material heat transfer problem using Crank-Nicolson?
I am working on a mixed material model for a melting material and need to enforce both a Dirichlet and Neumann type condition at the interface. Subject to an external surface heat flux at the top of ...
1
vote
0
answers
66
views
Deriving order of accuracy and interpreting a given discretization scheme when underlying method ( finite difference/volume) not known
If a spatial grid is given with time levels like this:
to solve the following model problem
Now consider the following discretization schemes:
Scheme 1
Scheme 2
Usually, to determine order of ...
0
votes
1
answer
182
views
Using Crank-Nicolson to solve Non-Linear Schrödinger equation in Python
I aim to solve the (non-linear) Schrodinger equation using the Crank-Nicolson method in Python. Here are my two functions.
...
1
vote
1
answer
304
views
Why does scipy Conjugate Gradient solver fail to converge for non-steady heat equation using Crank-Nicolson method
Could someone please explain why my implementation of the Crank-Nicolson method applied to the non-steady heat equation won't converge? There shouldn't be any nonlinear aspects to my implementation ...
1
vote
1
answer
366
views
Crank Nicolson Method with closed boundary conditions
I want to simulate 1D diffusion with a constant diffusion coefficient using the Crank-Nicolson method.
$$\frac{\partial u (x,t)}{\partial t} = D \frac{\partial^2 u(x,t)}{\partial x^2}.$$
I take an ...
0
votes
2
answers
417
views
Can the Crank-Nicolson Method Be used to Solve The Schrodinger Equation with a Time Varying Potential?
I have been following an excellent article about how to use the Crank-Nicolson method to solve the Schrodinger equation. In the article, it starts with a $V(x, y, t)$ but the potential seems to become ...
1
vote
0
answers
207
views
Solving PDE on a non-uniform grid with Crank-Nicolson scheme
I am solving a 1D diffusion-type equation with the finite-difference Crank-Nicolson (CN) scheme, and I need to densify the spatial grid around the central point. One could change the spatial variable ...
1
vote
0
answers
125
views
Crank-Nicolson vs Spectral Methods for the TDSE
The time-dependent Schroedinger equation (TDSE) depends linearly on the system's initial state $\vert \psi(0) \rangle$, such that the solution can be generally written as
$$ \vert \psi(t) \rangle = \...
2
votes
0
answers
173
views
Error in implementation of Crank-Nicolson method applied to 1D TDSE?
Some context, I've posted this question on physics SE and stack overflow. The former had nothing to offer, the latter had a great commenter that agreed with the phase looking off being one of the ...
0
votes
0
answers
89
views
Transparent Boundary Conditions for Finite Difference ADI PR 2D TDSE solution
I want to put (non-dirichlet) boundary conditions inside the code I wrote to solve the 2dim TDSE using the alternating direction implicit Peaceman - Rachford method.
$$
(1 + iB\Delta t/2 ) \psi^{n+1/2}...
4
votes
1
answer
207
views
Method to linearize highly nonlinear partial differential equation
I have a set of coupled pdes which I want to solve using finite-difference, of which one is nonlinear. The three linear pdes for quantities $T_f$, $T_s$ and $c$ are convection-diffusion-reaction-like ...
2
votes
1
answer
487
views
Solving Schrodinger Equation with finite element and Crank-Nicolson?
I have asked this in Mathematic section, but received no reply.
Please let me ask here to see if threr is any difference.
The Schrodinger equation without potential has the following form:
$$\...
2
votes
0
answers
145
views
Advection diffusion equation using Crank-Nicolson with total flux and Diriclet BCs
I am trying to model the 1D advection-diffusion equation:
$${\partial c \over \partial t} = D_c{\partial^2 c \over \partial x^2} -u{\partial c \over \partial x}.$$
With Robin boundary conditions that ...
3
votes
1
answer
148
views
Derivation of a parabolic PDE using Alternating Direction Implicit method
I have a very simple question concerning Alternating Direction Implicit (ADI) Scheme.
If I have an equation of the form:
\begin{equation*}
\frac{df(x,y,t)}{dt} = \nabla^2 f(x,y,t) + f(x,y,t)
\end{...
5
votes
0
answers
428
views
Stability of Crank-Nicholson for advection diffusion equation for spatial discretization other than finite differences second-order centered
Crank Nicholson is a time discretization method (see 4th equation here). From what I see around, you can use different space discretization, such as Finite elements. But for the linear advection-...
2
votes
0
answers
86
views
Linearising Nonlinear Coupled Partial Differential Equations - Alfvénic Diffusion
I am trying to solve the following coupled partial differential equations with a finite difference scheme:
$$\partial_tf+v\partial_zf+\partial_z\frac{1}{W}\partial_zf=0$$
$$\partial_tW+v\partial_zW-\...
3
votes
1
answer
136
views
Maintain unitary time evolution for a nonlinear ODE
I want to solve a nonlinear ODE of matrix $A(t)$
$$\mathrm{i}\dot A = A(t)M(t),\:\mathrm{with}\: M(t)=A^\dagger(t)H(t)A(t)$$ where $H(t)$ and hence $M(t)$ are Hermitian. Therefore, I presume the time ...
0
votes
0
answers
101
views
Crank-Nicolson solution of parabolic PDE with Newumann boundary conditions
I am trying to solve the non-linear parabolic PDE in $c(t,r)$
$$c_t=\frac{1}{r}(rDc_r-\alpha r^2 c)_r$$
with initial condition $c(0,r)=f(r)$
and boundary conditions $c_r(t,r_1)=\alpha r_1c_1/D$ and $...
0
votes
1
answer
187
views
FDM on nonlinear PDEs
I'm working with a 2D Navier Stokes PDE in the unstabilized version - the equation is a linear equation of the type $\frac{∂u}{∂t} = F(u,t)$.
In order to perform time discretization with FDM (finite ...
1
vote
0
answers
609
views
Time-dependent Schrodinger equation implementation in FEniCS
For our Bachelors thesis we're trying to solve the Schrodinger equation $i\partial_tu = -\nabla^2u+Vu$ in FEniCS. Given the domain $[-5, 5]^2$ with an initial value of $u_0(x, y)=e^{(-2(x^2+y^2))}$ ...
1
vote
1
answer
323
views
Crank-Nicholson for diffusion-advection vs diffusion equation
Let's consider the following 1D diffusion equation:
$\frac{\partial u}{\partial t} = xk \frac{\partial}{\partial x}(\frac{1}{x}\frac{\partial u}{\partial x})$
where we assume that the diffusion ...
6
votes
2
answers
389
views
Is the diffusion equation with Neumann and Dirichlet BCs well-posed?
I am considering the following diffusion equation:
$$\frac{\partial f}{\partial t} = \frac{\partial}{\partial x}[D(x,t)\frac{\partial f}{\partial x}]$$
over a grid ...
1
vote
1
answer
128
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 ...
2
votes
2
answers
122
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(...
1
vote
0
answers
53
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 $...
3
votes
0
answers
88
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^{...
1
vote
1
answer
313
views
Finite difference - Explicit / Implicit / Crank Nicolson - Does the implicit method require the least memory?
Examine a dynamic 2D heat equation $\dot{u} = \Delta u$ with zero boundary temperature. A standard finite difference approach is used on a rectangle using a $n\times n$ grid. For the resulting linear ...
5
votes
1
answer
779
views
Finite Differencing schemes for Convection-Diffusion equation
I'm using the Convection(/advection)-Diffusion(-Reaction) equation to calculate the temperature over time in different hydraulic parts like a pipe or a heat exchanger.
The flow/convection is always 1D,...
2
votes
1
answer
4k
views
Applying Neumann boundaries to Crank-Nicolson solution in python
Consider the heat equation
$$u_t = \kappa u_{xx}$$
with boundary conditions of
$$u(x,0)=0\\
u(0,t)=100\\
u(l,t)=0$$
Numerical analysis by pyton can be done with
...