Stability criteria of heat-like equation

In summary, the stability criterion for an explicit solution to the heat equation is that the time step must be less than half of the space step squared divided by the diffusivity. For an equation of the form where the diffusivity varies arbitrarily in space, the stability criterion can be determined using Von Neumann stability analysis. This involves solving for the maximum timestep allowed by considering the terms in the explicit finite differences scheme and their relative positions.
  • #1
Hypatio
151
1
I know that the criterion of stability for an explicit solution to the heat equation:

[itex]\frac{\partial T}{\partial t}=D\frac{\partial^2 T}{\partial x^2}[/itex]
is
[itex]\Delta t <\frac{1}{2}\frac{\Delta x^2}{D}[/itex]

however, what is the stability criterion for an equation of the form

[itex]\frac{\partial T}{\partial t}=D\frac{\partial^2}{\partial x^2}\left(\frac{T}{P(D)}\right)[/itex]

where P(D) indicates that P depends on the value of D, which varies arbitrarily in space.

I would like to solve this equation with explicit finite differences, so I will have a term of the form:

[itex]\frac{\partial T}{\partial t}=\frac{D+D_L}{2 \Delta x^2}(\frac{C_L}{P_L}-\frac{C}{P})+\frac{D+D_R}{2\Delta x^2}(\frac{C_R}{P_R}-\frac{C}{P})[/itex]

where subscripts L and R indicate relative position (left and right of the point of calculation).

How can I figure out the maximum timestep allowed.
 
  • #3
You want to use something call Von Neumann Stability analysis.
 

1. What is the heat-like equation?

The heat-like equation is a partial differential equation that describes the flow of heat in a given medium. It is a mathematical representation of how temperature changes over time and space.

2. What is stability criteria?

Stability criteria refers to the conditions that need to be satisfied in order for a numerical method to produce accurate and reliable solutions when solving a differential equation. In the case of the heat-like equation, stability criteria ensures that the numerical solution does not become unstable and produce meaningless results.

3. Why is stability important in the heat-like equation?

Stability is important in the heat-like equation because it ensures that the numerical solution accurately represents the physical phenomenon being modeled. Without stability, the solution may produce unrealistic or unphysical results, rendering it useless for practical applications.

4. What are some common stability criteria for the heat-like equation?

Some common stability criteria for the heat-like equation include the Courant-Friedrichs-Lewy (CFL) condition, the Von Neumann stability analysis, and the Lax equivalence theorem. These criteria involve analyzing the relationship between the time step, spatial discretization, and the physical properties of the medium being modeled.

5. How do stability criteria affect the choice of numerical methods for the heat-like equation?

The choice of numerical methods for the heat-like equation is heavily influenced by stability criteria. Some methods may be more stable than others, and therefore, are more suitable for solving the heat-like equation. Additionally, the choice of time step and spatial discretization must also satisfy stability criteria in order for the numerical solution to be accurate and reliable.

Similar threads

Replies
5
Views
1K
  • Differential Equations
Replies
3
Views
2K
Replies
4
Views
1K
  • Differential Equations
Replies
3
Views
1K
  • Differential Equations
Replies
7
Views
2K
  • Differential Equations
Replies
3
Views
1K
  • Differential Equations
Replies
1
Views
2K
  • Differential Equations
Replies
8
Views
2K
  • Differential Equations
Replies
3
Views
1K
  • Differential Equations
Replies
16
Views
901
Back
Top