(adsbygoogle = window.adsbygoogle || []).push({}); 1. The problem statement, all variables and given/known data

Perform a von Nueman stability analysis for the Crank-Nicholson FTCS scheme for the following ADE:

dc/dt + u dc/dx = D d^{2}c/dx^{2}- kc

where c is the concentration of interest, u is the advecting velocity and D is the dispersion coefficient.

We used a general forward time, centered space numerical scheme to solve the above partial differential equation ( θ = 0.5 Crank-Nicholson)

2. Relevant equations

3. The attempt at a solution

I got through most of the stability analysis including substituting in the errors for the finite difference approximation and solving for the value G:

G = [itex]\frac{1 - i 0.5 \frac{u Δt}{Δx} sin (β Δx) + \frac{D Δt}{Δx^{2}} (cos (β Δx) - 2) - 0.5k }{1+ i 0.5 \frac{u Δt}{Δx} sin (β Δx) - \frac{D Δt}{Δx^{2}} (cos (β Δx) - 2) - 0.5k }[/itex]

My question is, how do I eliminate the imaginary number term in the denominator, so that I may square G to determine whether it is less than 1 and thus stable? Do I use complex conjugates? If so, how? There are so many terms in the denominator that it just looks nightmare-ish.

**Physics Forums | Science Articles, Homework Help, Discussion**

Join Physics Forums Today!

The friendliest, high quality science and math community on the planet! Everyone who loves science is here!

The friendliest, high quality science and math community on the planet! Everyone who loves science is here!

# Homework Help: Von Nuemann Stability Analysis CN scheme

Can you offer guidance or do you also need help?

Draft saved
Draft deleted

**Physics Forums | Science Articles, Homework Help, Discussion**