Week 5 - Mid term project -Stability analysis in the unsteady problem

Stability Analysis:

1. The Stability of a Numerical Scheme depends on how the error propagates in the equation at different time steps
2. If the error magnifies during the course of equation solving then it causes the solution to get unstable
3. If the error remains stable during the equation solving then it remains stable.
4. An important criterion which indicates at the stability of an equation is the Courant Friedrichs Lewy Condition or CFL condition which determines conditions necessary for convergence while solving certain partial differential equations numerically by finite differences methods
5. The stability of the 2D Heat Conduction equation is `c=alpha(1/(deltax^2)+1/(deltay^2))`

Explicit Form :

Number of iteration done are 1808

Time taken for solving equation is 0.018894

Discussion:

As the CFL number becomes more than 0.25 we get an unstable solution. With CFL number less than 0.25 we get a stable solution.

Implicit Form:

• Implicit form is stable for large values of dt and the value of C ie CFL number does not affect the stabilty of the solution.
• Implicit Schemes time accuracy decreases with increase in time steps.
• Implicit Schemes are unconditionally stable and it uses entire domain to calculate each time steps therefore not affected by the CFL number.