Von Neumann Analysis: Refresh Numerical Science | Any Help Appreciated

Join the discussion
Ask a follow-up here, or get your own question answered by working scientists, mathematicians and engineers — people, not an autocomplete.
Real named experts · corrections over time · the nuance an AI answer skips
1 reply · 2K views
liquidFuzz
Messages
107
Reaction score
6
I'm trying to refresh some numerical science stuff. Von Neumann analysis, if I take a slimmed down equation, convection. [tex]\frac{∂u}{∂t}+a ∇ u =0[/tex] If I'm using Euler forward, [tex]\frac{u^{n+1}-u^n}{\Delta t}+\frac{a}{2h} \left( u_{j+1}^n -u_{j-1}^n \right) =0[/tex] For [itex]\hat{u}^n = G^n\hat{u}^0[/itex] a growth factor [itex]|G|\leq1[/itex] is sufficient. This gives, [tex]\frac{u^{n+1}-u^n}{\Delta t} + \beta h \left( \frac{u_{j+1}^n-2u_j^n + u_{j-1}^n}{h^2} \right)[/tex] With the smallest possible diffusion but still stable, [itex]\beta \geq \frac{a}{2}[/itex] becomes [itex]\beta = \frac{a}{2}[/itex]. If this is plugged into the equation, [tex]\frac{u^{n+1}-u^n}{\Delta t} + \frac{a}{2h} \left( {u_{j+1}^n-2u_j^n + u_{j-1}^n} \right)[/tex] Now, according to some old notes I have the result is, [tex]\frac{u^{n+1}-u^n}{\Delta t} + \frac{1}{h} \left( {au_j^n -a u_{j-1}^n} \right)[/tex] Well this is somewhat vague for me. Specially the last step where the defraction term is finalised.

Any pointers of help would be appreciated, specially that last step.
 
Physics news on Phys.org
I faintly remember that sometimes you could simplify things using a Taylor series going forward in j and a taylor series going backwards in j, then subtracting (or adding) one from the other. At least something along those lines. Considering the j-1 term remained, I'd consider doing the taylor series over the j+1 term (eq1) and the j term (eq2) and then multiplying eq2 by a factor of 2 and then adding (or subtracting) eq1 from eq2 and see what happens. You can probably ignore higher order terms which are negligible (or your only looking at a first order scheme).
 
Last edited: