Day16-slides - 1 Atms 502 Numerical Fluid Dynamics Tue Oct...

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 Document Right 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 Tue., Oct. 17, 2006 10/18/06 Atms 502 - Fall 2006 - Jewett 2 Multistage methods (1) • Earlier schemes were single-stage The terms on the right side - the spatial derivatives - were evaluated once Thus 2-time-level and 3-time-level schemes can be single-stage • Multistage methods Try to improve accuracy by evaluating right- hand-side F terms at additional times between n ∆ t and (n+1) ∆ t. 10/18/06 Atms 502 - Fall 2006 - Jewett 3 • 2-stage schemes (amplitude behavior) Multistage methods (2) " 2 = 1 + (1 # 2 $% )( & ’ t ) 2 + $ 2 % 2 ( & ’ t ) 4 Strong damping β =1 α =1 Matsuno Weakly unstable β =1 α =1/2 Midpoint Weakly unstable β =1/2 α =1 Heun RK RK RK: weak instability/2 nd order αβ >1/2: less accuracy but damped … damping can be useful! 10/18/06 Atms 502 - Fall 2006 - Jewett 4 Adams-Bashforth • 2nd-order: Weakly unstable (P), but with damped computational mode • 3rd-order: Highly efficient 1 comp. mode limits ∆ t Higher-order methods (1) " n + 1 = " n + ( h /2) 3 F ( " n ) # F ( " n # 1 ) ( ) " n + 1 = " n + ( h /12) 23 F ( " n ) # 16 F ( " n # 1 ) + 5 F ( " n # 2 ) $ % & & & ’ ( ) ) ) 2 10/18/06 Atms 502 - Fall 2006 - Jewett 5 Runge-Kutta • 3rd-order: • RK3 / RK4 methods have two free parameters; can choose to minimize storage Higher-order methods (2) q 1 = hF ( " n ) q 2 = hF ( " 1 ) # (5/9) q 1 q 3 = hF ( " 2 ) # (153/128) q 2 " 1 = " n + (1/3) q 1 " 2 = " 1 + (15/16) q 2 " n + 1 = " 2 + (8/15) q 3 RK3 and RK4 RK2 10/18/06 Atms 502 - Fall 2006 - Jewett 6 Time differencing: summary • Lists: Order of scheme Formula Storage factor Number of full arrays needed (if explicit) • Lists: Efficiency factor Largest stable step ÷ by # evaluations of F Amplitude, phase error Max stable ω ∆ t Durran Table 2.1, pp 68-69 10/18/06 Atms 502 - Fall 2006 - Jewett 7 Method discussion (1) • For smooth solutions:...
View Full Document

{[ snackBarMessage ]}

Page1 / 5

Day16-slides - 1 Atms 502 Numerical Fluid Dynamics Tue Oct...

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