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

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 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

Page1 / 5

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

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