Questions tagged [implicit-methods]
Implicit methods are timestepping methods that use an inversion at every timestep. This allows for much better stability properties than explicit methods, though it comes with a serious speed penalty in some cases. Examples of implicit methods include Backward Euler and Crank-Nicholson.
63
questions
0
votes
1
answer
53
views
Local truncation error of given implicit 1-step scheme
I'm given the 1-step implicit scheme $$y_{n+1} = y_n + \frac{h}{6}[4f(t_n, y_n) + 2f(t_{n+1}, y_{n+1}) + hf'(t_n, y_n)],$$
where $y'(t) = f(t, y)$, and I'm seeking the scheme's local truncation error. ...
- 1
0
votes
0
answers
37
views
First-order modified Patankar–Euler scheme (MPE)
Is the first-order Modified Patankar–Euler scheme (MPE) an implicit or explicit method?
Is there an open-source code implementing the MPE scheme for a system of ODEs?
2
votes
1
answer
155
views
Does this second-order implicit Runge-Kutta method have a name?
I am studying the time-integration of the following paper,
Young, L. C. (1981). A finite-element method for reservoir simulation. Society of Petroleum Engineers Journal, 21(01), 115-128.
A copy (PDF)...
- 544
1
vote
1
answer
112
views
Do Explicit Methods Always Require an Analytical Solution
Following some comments from another question I wanted to ask: does an explicit method always require some sort of analytical function/solution?
Let's take Euler for example. You have a function $f$ ...
- 33
1
vote
0
answers
88
views
Accuracy of the Crank-Nicolson method for non-linear, inhomogeneous heat equation
I am currently coding a solution to the following PDE:
$\frac{\partial T }{\partial t} =\frac{\partial}{\partial \theta}(A(\theta ,\phi )\frac{\partial T }{\partial \theta}) +\frac{\partial }{\partial ...
- 11
1
vote
1
answer
56
views
How can I derive a second order implicit method for this coupled ODE update?
I appreciate this might be an easy question, but I've managed to get myself quite thoroughly confused
I'm trying to solve a system of physics equations that look as follows
$$
\frac{\partial \mathbf{E}...
- 11
8
votes
2
answers
524
views
Are stiffness and instability equivalent?
To the best of my knowledge, stiffness of ordinary differential equations is difficult to capture but can be roughly described as problems where explicit methods don't work while implicit ones do. ...
- 231
2
votes
0
answers
43
views
Calculating the species mass consumption from implicit reaction-term in diffusion-reaction equation
The 1D diffusion equation with a chemical source term has the following form:
$$\frac{\partial Y}{\partial t} = D \frac{\partial^2 Y}{\partial x^2} - k Y,$$
where $Y$ is the molar concentration of the ...
- 21
4
votes
1
answer
257
views
Stiff ODE solver in the web browser
I'd like to make a web application that lets people play with solving ODE systems, changing parameters with sliders etc. but instead of doing the computations on the server side, solving the equations ...
- 433
3
votes
1
answer
476
views
For implicit schemes, is there any general result that says numerical diffusion increases with smaller timesteps (for CFL<1) as in explicit schemes?
For the first-order explicit upwind scheme, it can be easily shown that, if one keeps the same grid size and progressively decreases the time step below the max allowed one (i.e. below CFL~1) the ...
- 435
1
vote
1
answer
71
views
Implicit methods for variable coefficients based on equations of state
For example I have an equation that goes something like
$
\partial_t \rho = -\nabla\cdot (\rho u) + \nabla \cdot(D(\rho, T) \nabla \rho) + \rho_s
$
($\rho, \rho_s, u, T$ are coupled with a few other ...
- 159
0
votes
1
answer
225
views
How to fill matrix entries for two-dimensional implicit finite-difference for the general case
If I have derived a finite-difference formula for a 2D problem, for example something like:
$af_{i,j}+bf_{i-1,j}+cf_{i,j-1}+df_{i-1,j-1}=g_{i,j}$
where f is the unknown function on a grid and ...
- 111
6
votes
1
answer
957
views
Understanding butcher tableau when it comes to implicit methods
I've been learning about butcher tables and am having some difficulty understanding how to translate them when it comes to implicit methods. Specifically, I'm looking at backwards Euler:
\begin{array}...
- 183
2
votes
0
answers
66
views
Can Taylor methods be used effectively on stiff ODEs?
Cleve Moler has stated that "all numerical methods for stiff odes are implicit." However, I don't know whether this statement is a mathematical fact, or an simply an observation. Moreover, many ...
- 2,145
2
votes
1
answer
952
views
GMRES vs Newton-GMRES for Solving nonlinear PDE's
Often when numerically solving nonlinear PDE's using method of lines approach with an implicit integrator a system of nonlinear equations have to be solved.
To be more specific, let's say we have ...
- 95
1
vote
1
answer
122
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
1
vote
1
answer
297
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 ...
- 11
5
votes
1
answer
176
views
Mass matrix and BDF time integration
I have a system of nonlinear equations on the general form:
\begin{align}
\mathbf{M}(\bar{y})\dot{\bar{y}} =\bar{f}(\bar{y},t)
\end{align}
Where $\mathbf{M}(\bar{y})$ is a matrix and $\bar{f}$ is a ...
- 95
3
votes
1
answer
916
views
ODE: should Euler implicit be more accurate than Euler explicit for a given computational step?
I am aware than Euler explicit is conditionally stable, and Euler implicit is unconditionally stable. And I am aware that it is probably pointless to use Euler implicit with a small computational step ...
- 435
1
vote
1
answer
160
views
Stability Criteria for Numerical Solution of Windkessel Ordinary Differential Equation
I'm trying to solve this equation (Windkessel equation) numerically as:
$$C \frac{d P}{d t} + \frac{P}{R} = Q(t)$$
Where $C$ is compliance, $R$ is resistance, $P$ is pressure, and $Q(t)$ is a known ...
- 173
5
votes
1
answer
947
views
Solving a PDE implicitly by iteration in python
Connected to this question here on Computational Science, I've posted a follow-up question on how to solve a PDE using an implicit scheme like Crank-Nicholson in general in this question on SO.
But I ...
- 265
5
votes
1
answer
427
views
Solving PDE implicitly or explicitly depending on stiffness
I've got a system of several PDEs for a multitude of parts which represent real hydraulic parts like pipes or thermal energy storages. Each of these parts may have an arbitrary number of nodes and/or ...
- 265
2
votes
1
answer
90
views
Implicit solution to Sylvester equation
Suppose a matrix $M\in\mathbb{R}^{n\times n}$ is defined as the solution to a Sylvester equation $$AM+MB=C,$$ for some fixed (known) matrices $A,B,C$. In the regime where $n$ is large, we may with ...
- 1,341
-2
votes
1
answer
382
views
Adding Non-Linear source term to 2d Implicit MATLAB code
I'm running out of time for this code so any help would be greatly appreciated. I am currently coding the 2D heat/diffusion equation in matlab but i'm having trouble adding in the source term. my ...
0
votes
1
answer
190
views
From explicit to implicit SSP Runge-Kutta time discretization for DG
In Hesthaven book (Nodal Discontinuous Galerkin Methods) he uses SSP Runge-Kutta time method which is explicit.
How can I change the explicit RK to an implicit one?
- 13
1
vote
1
answer
236
views
Implicit integration for FLIP?
I have problem with volume loss in FLIP simulation.
Unfortunately it's necessary to obey the CFL condition when using explicit integration methods (RK2 in my case) to advance particle positions using ...
- 175
0
votes
1
answer
240
views
Slight Modification to Backward Euler Stiff ODE Solver
I am trying to implement a stiff ODE solver that uses step doubling with the backward Euler method but I want to make a modification and I am unsure how to incorporate it.
I am trying to solve $y'=f(...
- 103
2
votes
2
answers
120
views
How to support or contradict a hypothesis on unconditional stability using numerical optimization
The main motivation behind my next question is that I think I derived a higher order numerical scheme for linear advection equation that is unconditionally stable using Von Neumann stability analysis.
...
- 1,226
4
votes
2
answers
282
views
Textbook/Manual on Implicit FEM Methods
I've recently been interested in learning about implicit finite element methods. I've found this post, but I'd like to learn more about them, specifically about how simulations like these from LS-DYNA ...
- 383
4
votes
1
answer
299
views
Solve an ODE with positivity-preserving property unconditionally
I have an ODE for a scalar function $u=u(t)$ of the form:
$$
\frac{du}{dt}=L(u).
$$
Here the function $L=L(u)$ satisfies:
$$
L(0)=0, \quad L'(u)\le0.
$$
Then it is easy to see that the solution $u=u(t)...
- 246
1
vote
0
answers
65
views
How to solve bring my implicit equation to closed form? [duplicate]
I have a simulated system as shown in the following image:
$L_0$ is attached to two other bodies $L_1$ and $L_2$. Furthermore, body $L_3$ is also in the simulation (it is attached to $L_2$ but that ...
- 121
-1
votes
1
answer
137
views
Why is this method for simulating a system of springs and masses unstable?
I have a computer simulation system of bodies connected by springs, so their movement is governed by:
$x_{n+1} = x_n-\Delta tk(x_n-r)$
Where $r$ is the idea distance between every two bodies, and $\...
- 121
2
votes
0
answers
63
views
Resolving a stiff hyperbolic problem with Neumann boundary conditions
I am trying to numerically resolve the equation for an Euler-Bernoulli beam that is inextensible, unshearable, and subject to planar deformations:
$$\rho I(s) \frac{\partial^2 \theta}{\partial t^2}(s,...
- 187
3
votes
1
answer
274
views
LU-SGS and boundary conditions
I am trying to understand how boundary conditions are implemented when one uses the nonlinear LU-SGS algorithm for Euler equations. Most papers describe the Gauss-Seidel sweep over mesh cells, but do ...
- 405
5
votes
2
answers
2k
views
Why is Crank-Nicolson considered implicit in time?
From Wikipedia:
Explicit methods calculate the state of a system at a later time from
the state of the system at the current time, while implicit methods
find a solution by solving an equation ...
- 153
1
vote
1
answer
408
views
Implicit ODE solver with discontinuous derivatives
I want to implement an implicit ODE solver, but don't know what to do when the differential equations (DEs) have discontinuities of the form:
More common type:
$$\lim_{x\rightarrow 0}\frac{\sin(x)}{x}...
- 113
2
votes
1
answer
276
views
Not getting correct numerical solution for Advection-Diffusion-Reaction eqn
Objective: I am trying to numerically solve $C(x,y,t)$ from the following advection-diffusion-reaction equation in 2D space (x,y) and time. I will be testing my numerical solution with an approximate ...
user5510
1
vote
1
answer
574
views
ADR equation implicit solution: Penta-diagonal matrix for a 2D $N\times N$ system
Objective: I am trying to simulate the following advection-diffusion-reaction equation in 2D space (x,y) and time.
$$\begin{align}
\text{ADR Equation: }\frac{\partial C}{\partial t} + \nabla\left(v.C ...
user5510
1
vote
0
answers
137
views
Implicit time integrator for Chebyshev collocation method for linear hyperbolic system
I want to solve linear hyperbolic system using Chebyshev collocation method. As this method puts severe constraint on the time step for the explicit time integration, I decided to switch to implicit ...
- 508
0
votes
1
answer
109
views
Spectral Coefficients of Implicit Finite Difference Solution
This is something I've been trying to figure out for a long time, and all I have is vague numerical results. I'm trying to answer the following question analytically:
Suppose I have a time dependent ...
- 245
7
votes
2
answers
378
views
Initializing implicit linear multistep methods
A sixth order backward differentiation formula (BDF) need six (five plus initial value) previous solutions to get started. How I can get these previous solutions? I need a method accurate to sixth ...
- 223
1
vote
1
answer
299
views
CFD: Doubt with time convergence in advection fully implicit upwind scheme
I'm trying to solve an advection - convection problem using an implicit upwind scheme - you can see here the finite difference discretization used.
I start the model (built from scratch on Scilab) ...
- 11
1
vote
2
answers
8k
views
Fluent time step size (transient implicit solver)
I am studying fluid-structure interaction across integrated circuits packaging by using Ansys Fluent (I'm a beginner). Currently I had an issue about the time step size.
When I set time step size of ...
- 11
4
votes
1
answer
288
views
Non-conservative implementation implicit Euler
In Matlab R2013a I have implemented the Implicit Euler (time) integration scheme. To find the $x^{n+1}$ value I use fixed point iterations:
$x^{n+1} = \Delta t f(x^{n+1}) + x^n$
To test this, I use ...
- 78
2
votes
1
answer
2k
views
How can I implement the implicit Euler method for a small nonlinear system of ODEs?
I am trying to solve a system of coupled ODEs:
$$
\begin{align}
\frac{dn_A}{dt} & = e\left[j(t) - f\, θ_H\sinh\left(\frac{g\,n_A}{T}\right)\right] \\
\frac{dθ_H}{dt} & = a\left[bP\,(...
1
vote
1
answer
172
views
Implicit finite difference scheme for LWR-v PDE
I want to discretize the LWR(Lighthill-Whitham-Richards)-v partial differential equation (as shown in 1)
$$
\frac{\partial u}{\partial t} + \frac{\partial R(u)}{\partial x} = 0 \quad \text{i.e.,} \...
- 11
5
votes
1
answer
5k
views
How to add reaction and source terms to a diffusion PDE solver written with MATLAB's pdepe function?
I have the following system of equations which I'm trying to solve using Matlab's pdepe solver.
The 1-D spherical heat diffusion equation with heat generation (...
- 325
3
votes
1
answer
1k
views
Implicit heat diffusion with kinetic reactions
I am using the implicit finite difference method to discretize the 1-D transient heat diffusion equation for solid spherical and cylindrical shapes:
$$
\frac{1}{\alpha}\frac{\partial T}{\partial t} = ...
- 325
9
votes
1
answer
10k
views
Matlab solution for implicit finite difference heat equation with kinetic reactions
I am trying to model heat conduction within a wood cylinder using implicit finite difference methods. The general heat equation that I'm using for cylindrical and spherical shapes is:
Where p is the ...
- 325
16
votes
2
answers
3k
views
Is it possible to solve nonlinear PDEs without using Newton-Raphson iteration?
I am trying to understand some results and would appreciate some general comments on tackling nonlinear problems.
Fisher's equation (a nonlinear reaction-diffusion PDE),
$$
u_t = du_{xx} + \beta u ...
- 5,389