* Your assessment is very important for improving the workof artificial intelligence, which forms the content of this project
Download 2002
Routhian mechanics wikipedia , lookup
Computational chemistry wikipedia , lookup
Computational electromagnetics wikipedia , lookup
Perturbation theory wikipedia , lookup
Post-quantum cryptography wikipedia , lookup
Computational fluid dynamics wikipedia , lookup
Multidimensional empirical mode decomposition wikipedia , lookup
Relativistic quantum mechanics wikipedia , lookup
Written Preliminary Exam COMPUTATIONAL MECHANICS Spring 2002 Consider the diffusion equation, ∂φ ∂ 2φ =ν 2 , ∂t ∂x on 0 ≤ x ≤ 1 , subject to φ (0) = 0 and φ (1) = 1 . We solve this equation on a uniform spatial mesh, with spatial step ∆x and time step ∆t , using the following two schemes. φ n +1 − φ in φ n +1 − 2φ in +1 + φ in−+11 = ν i +1 ( I ) i 2 ∆ t x ∆ φ in +1 − φ in ( II ) ∆t ν φ in++11 − 2φ in +1 + φ in−+11 ν φ in+1 − 2φ in + φ in−1 = + 2 2 2 2 ∆ x ∆ x (a) Derive the modified equation for scheme II, determine its spatial and temporal order of accuracy, and establish consistency of the scheme. (b) By inspection, determine the spatial and temporal order of accuracy of scheme I (you do not need to derive the modified equation). (c) Perform von Neumann stability analyses for schemes I and II, and derive the amplification factors. (d) Experience shows that for large ∆t scheme I converges to the steady state solution while scheme II oscillates from one iteration to the next (φ in +1 = −φ in ) . Prove this. (e) The scheme φ in +1 − φ in ∆t is stable if φ in+1 − 2φ in + φ in−1 = ν ∆x 2 ν∆t 1 ≤ . Let’s say you‘ve written a code that uses this scheme. You specify ∆x 2 2 ∆x and ∆t as input. Your code blows up. What would you try to do to fix the problem and why?