Day25-slides - 1 Atms 502 Numerical Fluid Dynamics Thu Nov...

Info iconThis preview shows pages 1–3. Sign up to view the full content.

View Full Document Right Arrow Icon

Info iconThis preview has intentionally blurred sections. Sign up to view the full version.

View Full DocumentRight Arrow Icon
This is the end of the preview. Sign up to access the rest of the document.

Unformatted text preview: 1 Atms 502 Numerical Fluid Dynamics Thu., Nov. 17, 2006 12/6/06 Atms 502 - Fall 2006 - Jewett 2 Diffusion • We considered several approaches. Add up terms, as in computer pgm #5, e.g. u3=u1+(advection)+(diffusion)+.. Compute tendencies (contribution) from each, and add, to get A (n+1) Compute one process separately (e.g. advection), update , and use this result in next process This is known as process splitting 12/6/06 Atms 502 - Fall 2006 - Jewett 3 • 1-step filter to, e.g., remove 2 ∆ x in one pass. • Add damping term to PDE " = (1 # S ) " j + S 2 " j + 1 + " j # 1 ( ) Diffusion: from last time d " j dt = # " j + 1 $ 2 " j + " $ 1 ( ) 12/6/06 Atms 502 - Fall 2006 - Jewett 4 Diffusion • Amplification factor Given: The exact amplification factor is: And this tells us… : smallest waves damped most. " ( x , t ) = " e ikx e # Mk 2 t A e = " ( t + # t ) " ( t ) = e $ Mk 2 # t ; if r = M # t # x 2 , A e = e $ Mk 2 # t = e $ r # x 2 k 2 = e $ r % 2 A e = e " r # 2 2 12/6/06 Atms 502 - Fall 2006 - Jewett 5 Diffusion • Stability example - FTCS Amplification factor: Damps 2 ∆ x most strongly Stable for 0< ν ≤ (1/2) But ν >(1/4) has λ <0 for 2 ∆ x: sign flips each step....
View Full Document

This note was uploaded on 10/06/2009 for the course ATMOSPHERI 502 taught by Professor Jewett during the Fall '09 term at University of Illinois at Urbana–Champaign.

Page1 / 5

Day25-slides - 1 Atms 502 Numerical Fluid Dynamics Thu Nov...

This preview shows document pages 1 - 3. Sign up to view the full document.

View Full Document Right Arrow Icon
Ask a homework question - tutors are online