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.

learn more… | top users | synonyms

0
votes
1answer
45 views

Initial estimate of coupled variables while applying implicit euler to solve nonlinear system of odes

I am trying to solve a system of coupled odes having the following three equations. $$ \begin{align} \frac{dn_A}{dt} & = e\left[j(t) - f\, θ_H\sinh\left(\frac{g\,n_A}{T}\right)\right] \\ ...
1
vote
1answer
70 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.,} ...
5
votes
1answer
722 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 ...
0
votes
1answer
313 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} = ...
6
votes
1answer
1k 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 ...
9
votes
2answers
452 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 ...
6
votes
3answers
2k views

What is the difference between implicit FEM and explicit FEM?

What is the difference between explicit FEM and implicit FEM exactly? According to the post here, it seems that the only difference is whether implicit or explicit time integration is used. As I ...
8
votes
0answers
168 views

cuda and numerical methods with implicit time discretization

I am looking to port some code that resolves a set of partial differential equations (PDE) by the finite volume method in IMPLICIT form (for the time discretization). As result there is a tridiagonal ...
1
vote
1answer
139 views

CFL evolution techniques for Implicit methods

I am working on implicit schemes for Euler equations. Implicit methods allow one to use large CFL values, but is there some way to evolve CFL number from a much smaller value than desired value to ...
3
votes
2answers
250 views

Backward Euler method

Can you explain me how does the backward Euler method works? I have seen the formula and try to understand the method, but what I can't understand is why and how to use the Newton-Rapson method. Do ...
6
votes
1answer
107 views

Linearized implicit time stepping

Consider the general FD implicit time stepping scheme $\frac{x_{t+1} - x_t}{\Delta t} = f(x_{t+1})$, where $x$ is the vector variable of interest and $f$ is some function, generally non-linear. ...
12
votes
2answers
2k views

Implicit finite difference schemes for advection equation

There are numerous FD schemes for the advection equation $\frac{\partial T}{\partial t}+u\frac{\partial T}{\partial x}=0$ discuss in the web. For instance here: ...
5
votes
1answer
190 views

Does ADI/Split-operator change the stability properties of the Crank-Nicholson method?

I'm using the Crank-Nicholson method to solve the time-dependent Schrödinger equation with the split-operator method. I'm getting some weird results that are probably the result of a bug somewhere in ...
3
votes
1answer
462 views

How to obtain an implicit finite difference scheme for the wave equation?

Suppose I had the following problem: $U_{tt}=U_{xx}+U_{yy}$ in $\Omega=[0,1]\times[0,1]$ $U(x,y,0)=f(x,y)$ $U_{t}(x,y,0)=g(x,y)$ $U=0$ on $\partial \Omega$ I know that there is an explicit ...
70
votes
10answers
35k views

Recommendations for a usable, fast C++ matrix library?

Does anyone have recommendations on a usable, fast C++ matrix library? What I mean by usable is the following: Matrix objects have an intuitive interface (ex.: I can use rows and columns while ...
12
votes
2answers
182 views

What are the best practices for algorithms and implementation of multi-physics simulations?

Multi-physics simulation involves coupling multiple "physics", often with different space and/or time scales. Additionally, the single-physics codes are often written by different teams. The most ...
3
votes
1answer
135 views

Can heat distribution in an optical element irradiated by laser be oscillating?

I am modelling a heat distribution in optical element irradiated by laser. System is radially symmetric, and element is thin, i.e. heat value depends only on distance from center. Heat is received via ...
11
votes
1answer
530 views

When should implicit methods be used in the integration of hyperbolic PDEs?

Numerical methods for solving PDEs (or ODEs) fall into two broad categories: explicit and implicit methods. Implicit methods allow larger stable timesteps but require more work per step. For ...
16
votes
1answer
2k views

Why is Newton's method not converging?

I am using PETSc's nonlinear solver package SNES to solve a system of nonlinear equations obtained by discretizing a partial differential equation. How can I determine why the solver is not converging ...