Download 2002

Survey
yes no Was this document useful for you?
   Thank you for your participation!

* Your assessment is very important for improving the workof artificial intelligence, which forms the content of this project

Document related concepts

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

Spatial analysis wikipedia , lookup

Commitment scheme wikipedia , lookup

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