Contents of lecture 2b. Lectures 2a & 2b. Physical vs. computational coordinates [2] Physical vs. computational coordinates [1]
|
|
- Ella Neal
- 5 years ago
- Views:
Transcription
1 Contents of lecture b Lectures a & b P. A. Zegeling Mathematical Institute Utrecht University The Netherlands Parameter-free non-singular moving grids in D: Theory & properties Application to resistive Magneto-hydrodynamics Application to the D Euler equations March 7, 9 Extension to 3D and application to blow-up models Physical vs. computational coordinates [] Physical vs. computational coordinates [] The D adaptive grid can also be seen as an approximation of a coordinate transformation between computational coordinates U steep solution V mild solution (ξ, η) T Ω c := [, ] [, ] Φ (with a uniform grid partitioning) and physical coordinates (x, y) T Ω p R Y η (with a non-uniform adaptive grid) non uniform mesh X uniform mesh ξ
2 Winslow s method The adaptive grid seen as a system of springs In a variational setting, a grid-energy functional (à la Winslow) can be defined as E = ( ) T x ω x + T y ω y dξdη, Ω c X(i,j) ω (i,j) where = ( ξ, η )T and ω > is a monitor function. Minimizing the energy E yields the Euler-Lagrange equations: F (ω x) = (ω y) = F F on Ω c = [, ] [, ] with BCs: x ξ= = x L, y η= = y L, x ξ= = x R, y η= = y U, x n ξ= = x n ξ= = y n η= = y n η= =. F Regularity of the transformation in D A few additional properties of the D grid Theorem [Clément, Hagmeijer & Sweers, 96]: Let ω c >, ω C, (Ω c ) and ω ξ, ω η C γ ( Ω c ), for some γ (, ). unique solution (x, y) C ( Ω c ), which is a bijection from Ω c into itself. Moreover, the acobian satisfies: = x ξ y η x η y ξ >. The three main ingredients of proof are: Carleman-Hartman-Winter Theorem (3D??) ordan Curve Theorem (3D??) Maximum principle for elliptic PDEs (3D ok!) Equidistribution in D: x ξ }{{} ω = cst Winslow in D: (x ξ ) (y η ) (x η ) (y ξ ) = ω = cst (remember: = x ξ y η x η y ξ ) the transformation behind Winslow s method is not a harmonic mapping, but it is related to it a counterexample can be given for the 3D (harmonic) case, for which the transformation looses its regularity (Liao et al 94) several components in the proof of Theorem 3 can not be applied in 3D either...
3 The monitor function ω Smoothness of the new monitor function Arclength-type (AL-) monitor: ω = + α u u α is a (problem-dependent) adaptivity -parameter which controls the amount of adaptivity. New monitor (Beckett & Mackenzie ): ω = α(t) + u m, with α(t) = u m dξdη Ω c m = : better scaling and more adaptivity than for m = Define (Huang ) ω = + γ u ( γ) RR Ωc u dξdη Ω γ = c ω dξdη Ω c ω dξdη with γ [, ) For γ =, we have ω = + RR u Ωc u 5% of the grid dξdη points is concentrated in regions of high spatial derivatives, since Ω c ω dξdη the total # of grid points Ω c ω dξdη the # of grid points in the steep layer = a smoother distributed grid than for the AL-monitor (with constant α) and α(t) is automatically computed! Smoothing of the AL-monitor function Resistive Magneto-HydroDynamics With the new monitor, application of a filter or smoother to the grid or monitor values is not necessary. Normally, smoother transitions in a general non-uniform grid can be obtained (and are needed!) by working with the smoothed value S(ω i+,j+ ) = 4 ω i+,j+ + 8 (ω i+ 3,j+ + ω i,j+ + ω i+,j+ 3 +ω i+,j ) + 6 (ω i,j + ω i,j+ 3 + ω i+ 3,j + ω i+ 3,j+ 3 ) In the numerical experiments we denote this with filter on or filter off (working merely with ω i+,j+ values i.e. S(ω) = ω). (ρv) ρ + (ρv) = + (ρvv BB) + p tot = e + (ve + vp tot BB v) = η m ( B) B + (vb Bv) = η m B p tot = p + B, p = (γ )(e ρv B )
4 Applications of models Kinematic flux expulsion in D [] B 8πµ ρ v B = (v B) + η m B with v = z =, B = (B, B, ), B = B B B = η m B + v y v B y + B v y B v y B = η m B v x + v B x B v x + B v x Kinematic flux expulsion in D [] Kinematic flux expulsion in D [3] B := A, where A := vector potential A = v ( A) η m ( A) B 3 = ( A = A ) : Four-cell convection (Weiss 66): v (x, y) = sin(πx) cos(πy) v (x, y) = cos(πx) sin(πy) A 3 A 3 = v x v A 3 y + η m A 3 < η m, A 3 t= = x where v (x, y) and v (x, y) are given and satisfy A 3 x= =, A 3 x= =, A 3 y= = A 3 y= v x + v y =
5 Global solution behaviour [] Global solution behaviour [] < η m write formal asymptotic expansion in η m : A 3 (x, y, t) = a 3,j (x, y, t) ηm j j= Substitute expansion into PDE and check first-order term (setting η m = ); this gives hyperbolic PDE: a 3, = v a 3, initial solution stays constant on sub-characteristics given by: As v =, the only critical points in ODE system are center points or saddle points ( field amplification ). At some point of time, the solution a 3, (x, y, t) can not satisfy the boundary conditions of the original PDE model: boundary (and internal) layers are formed of width O( η m ) ( magnetic flux concentrates at edges of convective cells ). For t the solution reaches a non-trivial steady-state, in which the diffusion and convection terms settle down to an equilibrium. (ẋ, ẏ) T = v Global solution behaviour [3] The algorithm for the D model S S S C C advection-diffusion-reaction PDE decoupling of physical and grid PDEs grid PDEs: system of heat eq s with artificial time S S S fixed time steps t C C Implicit-explicit time integration: -SBDF S S S freezing of non-linear terms in PDEs BiCGstab + diagonal preconditioning
6 Transformation of the PDE model [] Transformation of the PDE model [] A 3 A 3 = v x v A 3 y + η m A 3 ξ = ξ(x, y, t), η = η(x, y, t), θ = t Using the chain rule of differentiation: We can also find that A 3,t = A 3,θ θ t + A 3,ξ ξ t + A 3,η η t. Using a similar relation for η t gives us A 3,t = A 3,θ + A 3,ξ (x ηy θ x θ y η ) + A 3,η (x θy ξ x ξ y θ ). Since ξ x = y η, ξ y = x η, η y = x ξ and η x = y ξ, we find for the first-order spatial derivative terms: A 3,x = (A 3,ξy η A 3,η y ξ ) Recall: ξ t = x θ ξ x y θ ξ y = (x θy η y θ x η ) and A 3,y = (A 3,ηx ξ A 3,ξ x η ). = x ξ y η x η y ξ Transformation of the PDE model [3] Numerical results; four-cell convection [] PDE in computational coordinates becomes: Weiss model, Four cell convection; velocity field (u(x,y),v(x,y)) T 3 A 3,θ + A 3,ξ (x η y θ x θ y η ) + A 3,η (x θ y ξ x ξ y θ ) 5 = A 3,ξ ( v y η + v x η ) + A 3,η (v y ξ v x ξ )+ η m [( x η+y η A 3,ξ) ξ ( x ξx η +y ξ y η A 3,η ) ξ Y β(t) 5 ( x ξx η +y ξ y η A 3,ξ ) η + ( x ξ +y ξ A 3,η) η ] = x ξ y η x η y ξ X # of TIMESTEPS
7 D case D Euler Numerical results; four-cell convection [] D case D Euler y x y x y D case x D Euler The D Euler equations ρ ρu + ρv x E ρu ρu + p + y ρuv u(e + p) y Numerical results; four-cell convection [3] D case x D Euler A typical non-uniform finite volume cell Ai+,j+ ρv ρuv ρv + p = v (E + p) ρ : density (ρu, ρv )T : momentum vector E : total energy p : pressure = (γ )(E ρ u +v ) Write hyperbolic PDE system as: Q F (Q) F (Q) + + = G(x, y, Q) x y n3 i+,j+ A n i+/,j+/ i,j i,j+ n4 i+,j n
8 Discretization of the grid PDEs A conservative solution-updating method Given a non-uniform partitioning {A i+,j+ } i,j of Ω p, where A i+,j+ is a quadrangle with four vertices x i+k,j+l, k, l (denote numerical approximations to x = x(ξ, η) by x i,j = x(ξ i, η j )) Subdivide Ω c = {(ξ, η) ξ, η } uniformly: (ξ i, η j ) ξ i = i ξ, η j = j η; i I ξ +, j I η + Discretize the elliptic system of grid PDEs by second-order central finite differences Apply a Gauß-Seidel iteration method to the resulting system of algebraic equations Having computed the new grid, the solution values Q have to be updated on this grid by interpolation (Tang et al 3 use a conservative interpolation method) Simple linear interpolation is not good enough Denote with x i,j & x i,j the coordinates of old and new grid points (x i,j moves to x i,j after Gauß-Seidel iterations) Using a perturbation technique and assuming small grid speeds, then the following mass-conservation is satisfied: i,j Ãi+,j+ Q i+,j+ where A is the area of cell A = i,j A i+,j+ Q i+,j+ Finite volume discretization on non-uniform grids [] Finite volume discretization on non-uniform grids [] Q + F (Q) + F (Q) = G(x, y, Q) x y Integration over the finite control volume A i+,j+ : Q dxdy + A i+,j+ 4 F n l (Q) (x,y) sl ds = s l l= A i+,j+ G dxdy Q n+ i+,j+ = Q n i+,j+ t A i+,j+ { F (Q n n i+,j )+F (Q n n i+ 3,j+ ) 4 } +F (Q n n 3 i+,j+ )+F (Q n 3 n 4 i,j+ )+ F + (Q n n l i+,j+ ) + t G n i+,j+ l= The time step size t is determined every time step by with F n l (Q) = F n l x + F n l y and F n l Discretization = F + n l + F n l t = min( x, y) CFL, max λ where λ are the eigenvalues of the acobian matrix F Q
9 D case D Euler Decoupling of the grids and physical PDEs D case D Euler D Riemann problem: shock waves Step Partition Ωc uniformly; give initial partitioning of Ωp ; compute initial grid values by a cell average of control volume Ai+,j+ based on initial data Q(x, y, ). In a loop over the time steps, update grid and solution and evaluate the PDE: Step a Move grid xi,j to x i,j by solving the discretized grid PDEs using Gauß-Seidel iterations. Step b Compute Q i+,j+ based on conservative interpolation. Repeat step a and step b for a fixed number of Gauß-Seidel iterations. Step 3 Evaluate the physical PDEs by the finite volume method on the grid x i,j to obtain Qn+ at t n+. i+,j+ The first test example is a two-dimensional Riemann problem (config. 4 in Lax & Liu 98) and has the following initial data: (.,.,.,.) if x >, y >, (65, 939,., 5) if x <, y >, (ρ, u, v, p)t= = (., 939, 939,.) if x <, y <, (65,., 939, 5) if x >, y <. They correspond to a left forward shock, right backward shock, upper backward shock and a lower forward shock, resp. The spatial domain Ωp is [, ] [, ] and Tend = 5. Step 4 Repeat steps a, b and 3 until t = Tend is reached. D case D Euler Comparison: AL vs. new monitor vs. uniform grid RUN # I II III IV V VI VII VIII IX monitor AL-monitor AL-monitor AL-monitor AL-monitor AL-monitor new monitor new monitor uniform (4 4) uniform (6 6) α m - filter on on on on off off off - runtime h 5m h 47m h 3m 5h 48m 3h 4m h57m h 9m h 5m 3h 5m Table: A shock wave model (configuration 4 from Lax & Liu 98) D case D Euler Numerical results for the AL-monitor left: α =., middle: α =, right: α =
10 D case D Euler Numerical results for the new monitor D case D Euler The adaptive grid: a close-up near (, 5) left: AL-monitor (α = & filter off); middle: AL-monitor (α = & filter on); right: new monitor with m = (no filter needed) left: uniform grid 4 and 6, middle: m =, right: m = D case D Euler The Rayleigh-Taylor instability model? This instability happens on an interface between fluids with different densities when an acceleration is directed from the heavy fluid to the light fluid? It has a fingering nature, with bubbles of light fluid rising into the ambient heavy fluid and spikes of heavy fluid falling into the light fluid? Here we have non-homogeneous BCs and an extra source term in the PDE-system: G = (,, ρg, ρvg )T, where g is the acceleration due to gravity D case D Euler Numerical results for Rayleigh-Taylor instability [] Results at t = and t =.
11 Extension of Winslow to 3D Extension of Winslow to 3D Minimization of grid-energy yields in 3D: with (ω x) = (ω y) = (ω z) =, (x, y, z) T [, ] 3 ω = α(t) + u, α(t) = Ω c u dξdηdζ Minimization of grid-energy yields in 3D: with (ω x) = x τ (ω y) = y τ (ω z) = z τ, (x, y, z)t [, ] 3 ω = α(t) + u, α(t) = Ω c u dξdηdζ Non-singular transformation? Theory? Non-singular transformation? Theory? Transformation of a blow-up PDE models (intro) [] The PDE u = u + up transforms via (x, y, z, t) (ξ, η, ζ, θ); with t = θ to: u θ + [u ξ( x θ (y η z ζ y ζ z η ) y θ (x ζ z η x η z ζ ) z θ (x η y ζ x ζ y η )) + u η ( x θ (y ζ z ξ y ξ z ζ ) Application areas: y θ (x ξ z ζ x ζ z ξ ) z θ (x ζ y ξ x ξ y ζ )) + u ζ ( x θ (y ξ z η y η z ξ ) y θ (x η z ξ x ξ z η ) z θ (x ξ y η x η y ξ ))] = [ (y η z ζ y ζ z η ) + (xζ z η x η z ζ ) + (x η y ζ x ζ y η )! u ξ + ξ «(yη z ζ y ζ z η )(y ζ z ξ y ξ z ζ ) + (x ζ z η x η z ζ )(x ξ z ζ x ζ z ξ ) + (x η y ζ x ζ y η )(x ζ y ξ x ξ y ζ ) u η + ξ «(yη z ζ y ζ z η )(y ξ z η y η z ξ ) + (x ζ z η x η z ζ )(x η z ξ x ξ z η ) + (x η y ζ x ζ y η )(x ξ y η x η y ξ ) u ζ + ξ «(yζ z ξ y ξ z ζ )(y η z ζ y ζ z η ) + (x ξ z ζ x ζ z ξ )(x ζ z η x η z ζ ) + (x ζ y ξ x ξ y ζ )(x η y ζ x ζ y η ) u ξ +... η Combustion models & chemical reaction dynamics Population dynamics: motion of colonies of micro-organisms Plasma physics: wave motion in fluids and electromagnetic fields with = z ζ (x ξ y η x η y ξ ) z η (x ξ y ζ x ζ y ξ ) + z ξ (x η y ζ x ζ y η )
12 models (intro) [] Consider the ODE: Exact solution: { u = up, p > u() = u u(t) = [(p )(T t)] p, T = u p (p ) u(t) models (intro) [3] Consider the PDE: { u = u + up, p > u Ω =, u(x, ) = u (x) Kaplan 963: if u smooth and large enough, then the solution u is regular for every t < T, but lim u(, t) L t T = + u t= t=t solution blows up at time T t models (intro) [4] Adaptive time steps with a Sundman transformation [] Another PDE example: { u = u u + u p, x (, π), t > x u(x, ) = u (x), x (, π), u(, t) = u(π, t) =, t > Let f = π u(x, t) sin(x) dx, then ḟ = π ( u x u p ) sin(x) dx f Via Hölder s inequality: ḟ f + f p p If f () > p p, then f (t) in finite time By Cauchy-Schwartz: f u L (,π) sin(x) L (,π) Explicit Euler with fixed t for { u = up u() = : u n+ = u n + t (u n ) p Numerical solution exists for all t n = n t (i.e. no blow up...) u(t) exact solution num. solution Thus f implies u L (,π) u leaves L (, π) in finite time u t= t=t t
13 Adaptive time steps with a Sundman transformation [] Adaptive time steps with a Sundman transformation [3] Explicit Euler + central FD s for u t = u xx + u p : In terms of scaling invariance and self-similarity: with t n = u n+ j u n j t n θ u n p = un j+ un j + u n j ( x) and constant θ ( a Sundman time transformation t(θ) = θ + (u n j ) p u p for sufficiently small x (+ stab. restriction on t n ), the numerical solution blows up at T x and lim x T x = T (Abia et al, ) ) Consider the ODE u = u, i.e. p =, then using a fictive computational time variable θ gives rise to a new ODE system with { du dθ = u dt dθ = u that is invariant under the scaling t λt, u λ u and for which the numerical solution u n uniformly approximates the true self-similar solution u(t) = t of the original ODE for θ. (Budd, Piggott, Leimkuhler et al) The numerical algorithm in 3D in 3D [] decoupling of blow-up and grid PDEs grid PDEs: system of heat equations with artificial time central finite differences on non-uniform grid for freezing of non-linear terms in PDEs in each time step Case : u(x, y, z, ) = exp( 3((x ) + (y ) + (z ) ) (p = 3, 3 / 3 /4 3 grid; with or without Sundman t(θ)) implicit Euler for and explicit for reaction term BiCGstab + ILU-preconditioning for underlying linear systems variable t using Sundman-transformation Maximum value of u as a function of time t.
14 in 3D [] in 3D [3] Case : Case (cntd.): u(x, y, z, ) = sin(πx) sin(πy) sin(πz) (p = 3, 3 / 3 /4 3 grid; with or without Sundman t(θ)) Adaptive grids and numerical solutions at two points of time just before blow-up (note the scale of the solution in each of the plots) Maximum value of u as a function of time t. (blow-up time in our experiments T from Liang & Lin, 5:.749) in 3D [4] Summary Case (cntd.): D adaptive moving grids: theory and applications! D simple domains : some theory but many applications : lack of theory and still only a few applications Adaptive grid and numerical solution just before blow-up (right plot: close-up of grid and solution)
Basic Aspects of Discretization
Basic Aspects of Discretization Solution Methods Singularity Methods Panel method and VLM Simple, very powerful, can be used on PC Nonlinear flow effects were excluded Direct numerical Methods (Field Methods)
More informationNumerical Solutions to Partial Differential Equations
Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University Numerical Methods for Partial Differential Equations Finite Difference Methods
More informationIntroduction to multiscale modeling and simulation. Explicit methods for ODEs : forward Euler. y n+1 = y n + tf(y n ) dy dt = f(y), y(0) = y 0
Introduction to multiscale modeling and simulation Lecture 5 Numerical methods for ODEs, SDEs and PDEs The need for multiscale methods Two generic frameworks for multiscale computation Explicit methods
More informationChapter Two: Numerical Methods for Elliptic PDEs. 1 Finite Difference Methods for Elliptic PDEs
Chapter Two: Numerical Methods for Elliptic PDEs Finite Difference Methods for Elliptic PDEs.. Finite difference scheme. We consider a simple example u := subject to Dirichlet boundary conditions ( ) u
More informationTutorial 2. Introduction to numerical schemes
236861 Numerical Geometry of Images Tutorial 2 Introduction to numerical schemes c 2012 Classifying PDEs Looking at the PDE Au xx + 2Bu xy + Cu yy + Du x + Eu y + Fu +.. = 0, and its discriminant, B 2
More informationGame Physics. Game and Media Technology Master Program - Utrecht University. Dr. Nicolas Pronost
Game and Media Technology Master Program - Utrecht University Dr. Nicolas Pronost Soft body physics Soft bodies In reality, objects are not purely rigid for some it is a good approximation but if you hit
More informationNumerical Solutions to Partial Differential Equations
Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University Introduction to Hyperbolic Equations The Hyperbolic Equations n-d 1st Order Linear
More informationA Space-Time Expansion Discontinuous Galerkin Scheme with Local Time-Stepping for the Ideal and Viscous MHD Equations
A Space-Time Expansion Discontinuous Galerkin Scheme with Local Time-Stepping for the Ideal and Viscous MHD Equations Ch. Altmann, G. Gassner, F. Lörcher, C.-D. Munz Numerical Flow Models for Controlled
More informationRelevant self-assessment exercises: [LIST SELF-ASSESSMENT EXERCISES HERE]
Chapter 6 Finite Volume Methods In the previous chapter we have discussed finite difference methods for the discretization of PDEs. In developing finite difference methods we started from the differential
More informationScientific Computing: An Introductory Survey
Scientific Computing: An Introductory Survey Chapter 11 Partial Differential Equations Prof. Michael T. Heath Department of Computer Science University of Illinois at Urbana-Champaign Copyright c 2002.
More informationChapter 3. Finite Difference Methods for Hyperbolic Equations Introduction Linear convection 1-D wave equation
Chapter 3. Finite Difference Methods for Hyperbolic Equations 3.1. Introduction Most hyperbolic problems involve the transport of fluid properties. In the equations of motion, the term describing the transport
More informationFinite Difference Methods for Boundary Value Problems
Finite Difference Methods for Boundary Value Problems October 2, 2013 () Finite Differences October 2, 2013 1 / 52 Goals Learn steps to approximate BVPs using the Finite Difference Method Start with two-point
More informationNumerical Analysis of Differential Equations Numerical Solution of Parabolic Equations
Numerical Analysis of Differential Equations 215 6 Numerical Solution of Parabolic Equations 6 Numerical Solution of Parabolic Equations TU Bergakademie Freiberg, SS 2012 Numerical Analysis of Differential
More informationChapter 10 Exercises
Chapter 10 Exercises From: Finite Difference Methods for Ordinary and Partial Differential Equations by R. J. LeVeque, SIAM, 2007. http://www.amath.washington.edu/ rl/fdmbook Exercise 10.1 (One-sided and
More informationPartitioned Methods for Multifield Problems
C Partitioned Methods for Multifield Problems Joachim Rang, 6.7.2016 6.7.2016 Joachim Rang Partitioned Methods for Multifield Problems Seite 1 C One-dimensional piston problem fixed wall Fluid flexible
More informationAdditive Manufacturing Module 8
Additive Manufacturing Module 8 Spring 2015 Wenchao Zhou zhouw@uark.edu (479) 575-7250 The Department of Mechanical Engineering University of Arkansas, Fayetteville 1 Evaluating design https://www.youtube.com/watch?v=p
More informationNotes: Outline. Diffusive flux. Notes: Notes: Advection-diffusion
Outline This lecture Diffusion and advection-diffusion Riemann problem for advection Diagonalization of hyperbolic system, reduction to advection equations Characteristics and Riemann problem for acoustics
More informationPDE Solvers for Fluid Flow
PDE Solvers for Fluid Flow issues and algorithms for the Streaming Supercomputer Eran Guendelman February 5, 2002 Topics Equations for incompressible fluid flow 3 model PDEs: Hyperbolic, Elliptic, Parabolic
More informationFinite difference methods for the diffusion equation
Finite difference methods for the diffusion equation D150, Tillämpade numeriska metoder II Olof Runborg May 0, 003 These notes summarize a part of the material in Chapter 13 of Iserles. They are based
More informationÉquation de Burgers avec particule ponctuelle
Équation de Burgers avec particule ponctuelle Nicolas Seguin Laboratoire J.-L. Lions, UPMC Paris 6, France 7 juin 2010 En collaboration avec B. Andreianov, F. Lagoutière et T. Takahashi Nicolas Seguin
More informationLast time: Diffusion - Numerical scheme (FD) Heat equation is dissipative, so why not try Forward Euler:
Lecture 7 18.086 Last time: Diffusion - Numerical scheme (FD) Heat equation is dissipative, so why not try Forward Euler: U j,n+1 t U j,n = U j+1,n 2U j,n + U j 1,n x 2 Expected accuracy: O(Δt) in time,
More informationIntroduction to numerical schemes
236861 Numerical Geometry of Images Tutorial 2 Introduction to numerical schemes Heat equation The simple parabolic PDE with the initial values u t = K 2 u 2 x u(0, x) = u 0 (x) and some boundary conditions
More informationInverse Lax-Wendroff Procedure for Numerical Boundary Conditions of. Conservation Laws 1. Abstract
Inverse Lax-Wendroff Procedure for Numerical Boundary Conditions of Conservation Laws Sirui Tan and Chi-Wang Shu 3 Abstract We develop a high order finite difference numerical boundary condition for solving
More informationHOMEWORK 4: Numerical solution of PDEs for Mathmatical Models, Analysis and Simulation, Fall 2011 Report due Mon Nov 12, Maximum score 6.0 pts.
HOMEWORK 4: Numerical solution of PDEs for Mathmatical Models, Analysis and Simulation, Fall 2011 Report due Mon Nov 12, 2012. Maximum score 6.0 pts. Read Strang s (old) book, Sections 3.1-3.4 and 6.4-6.5.
More information2.29 Numerical Fluid Mechanics Spring 2015 Lecture 13
REVIEW Lecture 12: Spring 2015 Lecture 13 Grid-Refinement and Error estimation Estimation of the order of convergence and of the discretization error Richardson s extrapolation and Iterative improvements
More information0.2. CONSERVATION LAW FOR FLUID 9
0.2. CONSERVATION LAW FOR FLUID 9 Consider x-component of Eq. (26), we have D(ρu) + ρu( v) dv t = ρg x dv t S pi ds, (27) where ρg x is the x-component of the bodily force, and the surface integral is
More informationx n+1 = x n f(x n) f (x n ), n 0.
1. Nonlinear Equations Given scalar equation, f(x) = 0, (a) Describe I) Newtons Method, II) Secant Method for approximating the solution. (b) State sufficient conditions for Newton and Secant to converge.
More information2.29 Numerical Fluid Mechanics Spring 2015 Lecture 9
Spring 2015 Lecture 9 REVIEW Lecture 8: Direct Methods for solving (linear) algebraic equations Gauss Elimination LU decomposition/factorization Error Analysis for Linear Systems and Condition Numbers
More informationHyperbolic Systems of Conservation Laws. in One Space Dimension. I - Basic concepts. Alberto Bressan. Department of Mathematics, Penn State University
Hyperbolic Systems of Conservation Laws in One Space Dimension I - Basic concepts Alberto Bressan Department of Mathematics, Penn State University http://www.math.psu.edu/bressan/ 1 The Scalar Conservation
More informationNumerical Solutions to Partial Differential Equations
Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University The Implicit Schemes for the Model Problem The Crank-Nicolson scheme and θ-scheme
More informationFinite Volume for Fusion Simulations
Finite Volume for Fusion Simulations Elise Estibals, Hervé Guillard, Afeintou Sangam To cite this version: Elise Estibals, Hervé Guillard, Afeintou Sangam. Finite Volume for Fusion Simulations. Jorek Meeting
More informationFinite Volume Method
Finite Volume Method An Introduction Praveen. C CTFD Division National Aerospace Laboratories Bangalore 560 037 email: praveen@cfdlab.net April 7, 2006 Praveen. C (CTFD, NAL) FVM CMMACS 1 / 65 Outline
More informationPartial Differential Equations
Partial Differential Equations Introduction Deng Li Discretization Methods Chunfang Chen, Danny Thorne, Adam Zornes CS521 Feb.,7, 2006 What do You Stand For? A PDE is a Partial Differential Equation This
More informationAM 205: lecture 14. Last time: Boundary value problems Today: Numerical solution of PDEs
AM 205: lecture 14 Last time: Boundary value problems Today: Numerical solution of PDEs ODE BVPs A more general approach is to formulate a coupled system of equations for the BVP based on a finite difference
More informationConsider a volume Ω enclosing a mass M and bounded by a surface δω. d dt. q n ds. The Work done by the body on the surroundings is.
The Energy Balance Consider a volume Ω enclosing a mass M and bounded by a surface δω. δω At a point x, the density is ρ, the local velocity is v, and the local Energy density is U. U v The rate of change
More informationApplied Math Qualifying Exam 11 October Instructions: Work 2 out of 3 problems in each of the 3 parts for a total of 6 problems.
Printed Name: Signature: Applied Math Qualifying Exam 11 October 2014 Instructions: Work 2 out of 3 problems in each of the 3 parts for a total of 6 problems. 2 Part 1 (1) Let Ω be an open subset of R
More informationDifferential equations, comprehensive exam topics and sample questions
Differential equations, comprehensive exam topics and sample questions Topics covered ODE s: Chapters -5, 7, from Elementary Differential Equations by Edwards and Penney, 6th edition.. Exact solutions
More informationFinite Volume Schemes: an introduction
Finite Volume Schemes: an introduction First lecture Annamaria Mazzia Dipartimento di Metodi e Modelli Matematici per le Scienze Applicate Università di Padova mazzia@dmsa.unipd.it Scuola di dottorato
More informationGeometric PDE and The Magic of Maximum Principles. Alexandrov s Theorem
Geometric PDE and The Magic of Maximum Principles Alexandrov s Theorem John McCuan December 12, 2013 Outline 1. Young s PDE 2. Geometric Interpretation 3. Maximum and Comparison Principles 4. Alexandrov
More informationLecture 5: Single Step Methods
Lecture 5: Single Step Methods J.K. Ryan@tudelft.nl WI3097TU Delft Institute of Applied Mathematics Delft University of Technology 1 October 2012 () Single Step Methods 1 October 2012 1 / 44 Outline 1
More informationDiffusion / Parabolic Equations. PHY 688: Numerical Methods for (Astro)Physics
Diffusion / Parabolic Equations Summary of PDEs (so far...) Hyperbolic Think: advection Real, finite speed(s) at which information propagates carries changes in the solution Second-order explicit methods
More informationFDM for parabolic equations
FDM for parabolic equations Consider the heat equation where Well-posed problem Existence & Uniqueness Mass & Energy decreasing FDM for parabolic equations CNFD Crank-Nicolson + 2 nd order finite difference
More informationRecapitulation: Questions on Chaps. 1 and 2 #A
Recapitulation: Questions on Chaps. 1 and 2 #A Chapter 1. Introduction What is the importance of plasma physics? How are plasmas confined in the laboratory and in nature? Why are plasmas important in astrophysics?
More informationThe one-dimensional equations for the fluid dynamics of a gas can be written in conservation form as follows:
Topic 7 Fluid Dynamics Lecture The Riemann Problem and Shock Tube Problem A simple one dimensional model of a gas was introduced by G.A. Sod, J. Computational Physics 7, 1 (1978), to test various algorithms
More informationPDEs, part 3: Hyperbolic PDEs
PDEs, part 3: Hyperbolic PDEs Anna-Karin Tornberg Mathematical Models, Analysis and Simulation Fall semester, 2011 Hyperbolic equations (Sections 6.4 and 6.5 of Strang). Consider the model problem (the
More informationBasics of Discretization Methods
Basics of Discretization Methods In the finite difference approach, the continuous problem domain is discretized, so that the dependent variables are considered to exist only at discrete points. Derivatives
More informationApplying Asymptotic Approximations to the Full Two-Fluid Plasma System to Study Reduced Fluid Models
0-0 Applying Asymptotic Approximations to the Full Two-Fluid Plasma System to Study Reduced Fluid Models B. Srinivasan, U. Shumlak Aerospace and Energetics Research Program, University of Washington, Seattle,
More informationFinite Differences for Differential Equations 28 PART II. Finite Difference Methods for Differential Equations
Finite Differences for Differential Equations 28 PART II Finite Difference Methods for Differential Equations Finite Differences for Differential Equations 29 BOUNDARY VALUE PROBLEMS (I) Solving a TWO
More informationSoft Bodies. Good approximation for hard ones. approximation breaks when objects break, or deform. Generalization: soft (deformable) bodies
Soft-Body Physics Soft Bodies Realistic objects are not purely rigid. Good approximation for hard ones. approximation breaks when objects break, or deform. Generalization: soft (deformable) bodies Deformed
More informationAn Introduction to the Discontinuous Galerkin Method
An Introduction to the Discontinuous Galerkin Method Krzysztof J. Fidkowski Aerospace Computational Design Lab Massachusetts Institute of Technology March 16, 2005 Computational Prototyping Group Seminar
More informationAdvection / Hyperbolic PDEs. PHY 604: Computational Methods in Physics and Astrophysics II
Advection / Hyperbolic PDEs Notes In addition to the slides and code examples, my notes on PDEs with the finite-volume method are up online: https://github.com/open-astrophysics-bookshelf/numerical_exercises
More informationTime stepping methods
Time stepping methods ATHENS course: Introduction into Finite Elements Delft Institute of Applied Mathematics, TU Delft Matthias Möller (m.moller@tudelft.nl) 19 November 2014 M. Möller (DIAM@TUDelft) Time
More informationLecture 4: Numerical solution of ordinary differential equations
Lecture 4: Numerical solution of ordinary differential equations Department of Mathematics, ETH Zürich General explicit one-step method: Consistency; Stability; Convergence. High-order methods: Taylor
More informationPH.D. PRELIMINARY EXAMINATION MATHEMATICS
UNIVERSITY OF CALIFORNIA, BERKELEY Dept. of Civil and Environmental Engineering FALL SEMESTER 2014 Structural Engineering, Mechanics and Materials NAME PH.D. PRELIMINARY EXAMINATION MATHEMATICS Problem
More informationA Second Order Discontinuous Galerkin Fast Sweeping Method for Eikonal Equations
A Second Order Discontinuous Galerkin Fast Sweeping Method for Eikonal Equations Fengyan Li, Chi-Wang Shu, Yong-Tao Zhang 3, and Hongkai Zhao 4 ABSTRACT In this paper, we construct a second order fast
More informationFinite difference method for elliptic problems: I
Finite difference method for elliptic problems: I Praveen. C praveen@math.tifrbng.res.in Tata Institute of Fundamental Research Center for Applicable Mathematics Bangalore 560065 http://math.tifrbng.res.in/~praveen
More informationu n 2 4 u n 36 u n 1, n 1.
Exercise 1 Let (u n ) be the sequence defined by Set v n = u n 1 x+ u n and f (x) = 4 x. 1. Solve the equations f (x) = 1 and f (x) =. u 0 = 0, n Z +, u n+1 = u n + 4 u n.. Prove that if u n < 1, then
More informationCharacteristics for IBVP. Notes: Notes: Periodic boundary conditions. Boundary conditions. Notes: In x t plane for the case u > 0: Solution:
AMath 574 January 3, 20 Today: Boundary conditions Multi-dimensional Wednesday and Friday: More multi-dimensional Reading: Chapters 8, 9, 20 R.J. LeVeque, University of Washington AMath 574, January 3,
More information12 The Heat equation in one spatial dimension: Simple explicit method and Stability analysis
ATH 337, by T. Lakoba, University of Vermont 113 12 The Heat equation in one spatial dimension: Simple explicit method and Stability analysis 12.1 Formulation of the IBVP and the minimax property of its
More informationIntegrodifferential Hyperbolic Equations and its Application for 2-D Rotational Fluid Flows
Integrodifferential Hyperbolic Equations and its Application for 2-D Rotational Fluid Flows Alexander Chesnokov Lavrentyev Institute of Hydrodynamics Novosibirsk, Russia chesnokov@hydro.nsc.ru July 14,
More informationNumerical Methods for Engineers and Scientists
Numerical Methods for Engineers and Scientists Second Edition Revised and Expanded Joe D. Hoffman Department of Mechanical Engineering Purdue University West Lafayette, Indiana m MARCEL D E К К E R MARCEL
More informationNumerical Solution Techniques in Mechanical and Aerospace Engineering
Numerical Solution Techniques in Mechanical and Aerospace Engineering Chunlei Liang LECTURE 9 Finite Volume method II 9.1. Outline of Lecture Conservation property of Finite Volume method Apply FVM to
More informationNumerical Methods for PDEs
Numerical Methods for PDEs Problems 1. Numerical Differentiation. Find the best approximation to the second drivative d 2 f(x)/dx 2 at x = x you can of a function f(x) using (a) the Taylor series approach
More informationNumerical Methods for Conservation Laws WPI, January 2006 C. Ringhofer C2 b 2
Numerical Methods for Conservation Laws WPI, January 2006 C. Ringhofer ringhofer@asu.edu, C2 b 2 2 h2 x u http://math.la.asu.edu/ chris Last update: Jan 24, 2006 1 LITERATURE 1. Numerical Methods for Conservation
More informationMath 127C, Spring 2006 Final Exam Solutions. x 2 ), g(y 1, y 2 ) = ( y 1 y 2, y1 2 + y2) 2. (g f) (0) = g (f(0))f (0).
Math 27C, Spring 26 Final Exam Solutions. Define f : R 2 R 2 and g : R 2 R 2 by f(x, x 2 (sin x 2 x, e x x 2, g(y, y 2 ( y y 2, y 2 + y2 2. Use the chain rule to compute the matrix of (g f (,. By the chain
More informationStability of Feedback Solutions for Infinite Horizon Noncooperative Differential Games
Stability of Feedback Solutions for Infinite Horizon Noncooperative Differential Games Alberto Bressan ) and Khai T. Nguyen ) *) Department of Mathematics, Penn State University **) Department of Mathematics,
More informationGoverning Equations of Fluid Dynamics
Chapter 3 Governing Equations of Fluid Dynamics The starting point of any numerical simulation are the governing equations of the physics of the problem to be solved. In this chapter, we first present
More informationA Framework for Analyzing and Constructing Hierarchical-Type A Posteriori Error Estimators
A Framework for Analyzing and Constructing Hierarchical-Type A Posteriori Error Estimators Jeff Ovall University of Kentucky Mathematics www.math.uky.edu/ jovall jovall@ms.uky.edu Kentucky Applied and
More informationApplied Mathematics 205. Unit III: Numerical Calculus. Lecturer: Dr. David Knezevic
Applied Mathematics 205 Unit III: Numerical Calculus Lecturer: Dr. David Knezevic Unit III: Numerical Calculus Chapter III.3: Boundary Value Problems and PDEs 2 / 96 ODE Boundary Value Problems 3 / 96
More information2.2. Methods for Obtaining FD Expressions. There are several methods, and we will look at a few:
.. Methods for Obtaining FD Expressions There are several methods, and we will look at a few: ) Taylor series expansion the most common, but purely mathematical. ) Polynomial fitting or interpolation the
More informationME Computational Fluid Mechanics Lecture 5
ME - 733 Computational Fluid Mechanics Lecture 5 Dr./ Ahmed Nagib Elmekawy Dec. 20, 2018 Elliptic PDEs: Finite Difference Formulation Using central difference formulation, the so called five-point formula
More informationFluid Animation. Christopher Batty November 17, 2011
Fluid Animation Christopher Batty November 17, 2011 What distinguishes fluids? What distinguishes fluids? No preferred shape Always flows when force is applied Deforms to fit its container Internal forces
More informationGetting started: CFD notation
PDE of p-th order Getting started: CFD notation f ( u,x, t, u x 1,..., u x n, u, 2 u x 1 x 2,..., p u p ) = 0 scalar unknowns u = u(x, t), x R n, t R, n = 1,2,3 vector unknowns v = v(x, t), v R m, m =
More informationSymmetry Methods for Differential and Difference Equations. Peter Hydon
Lecture 2: How to find Lie symmetries Symmetry Methods for Differential and Difference Equations Peter Hydon University of Kent Outline 1 Reduction of order for ODEs and O Es 2 The infinitesimal generator
More informationBlock-Structured Adaptive Mesh Refinement
Block-Structured Adaptive Mesh Refinement Lecture 2 Incompressible Navier-Stokes Equations Fractional Step Scheme 1-D AMR for classical PDE s hyperbolic elliptic parabolic Accuracy considerations Bell
More information1 Assignment 1: Nonlinear dynamics (due September
Assignment : Nonlinear dynamics (due September 4, 28). Consider the ordinary differential equation du/dt = cos(u). Sketch the equilibria and indicate by arrows the increase or decrease of the solutions.
More informationUNIVERSITY of LIMERICK OLLSCOIL LUIMNIGH
UNIVERSITY of LIMERICK OLLSCOIL LUIMNIGH Faculty of Science and Engineering Department of Mathematics and Statistics END OF SEMESTER ASSESSMENT PAPER MODULE CODE: MA4006 SEMESTER: Spring 2011 MODULE TITLE:
More informationWeighted Essentially Non-Oscillatory limiters for Runge-Kutta Discontinuous Galerkin Methods
Weighted Essentially Non-Oscillatory limiters for Runge-Kutta Discontinuous Galerkin Methods Jianxian Qiu School of Mathematical Science Xiamen University jxqiu@xmu.edu.cn http://ccam.xmu.edu.cn/teacher/jxqiu
More informationQualifying Examination
Summer 24 Day. Monday, September 5, 24 You have three hours to complete this exam. Work all problems. Start each problem on a All problems are 2 points. Please email any electronic files in support of
More informationIntro to Research Computing with Python: Partial Differential Equations
Intro to Research Computing with Python: Partial Differential Equations Erik Spence SciNet HPC Consortium 28 November 2013 Erik Spence (SciNet HPC Consortium) PDEs 28 November 2013 1 / 23 Today s class
More informationFinite Difference Methods for
CE 601: Numerical Methods Lecture 33 Finite Difference Methods for PDEs Course Coordinator: Course Coordinator: Dr. Suresh A. Kartha, Associate Professor, Department of Civil Engineering, IIT Guwahati.
More informationIntroduction to PDEs and Numerical Methods Tutorial 5. Finite difference methods equilibrium equation and iterative solvers
Platzhalter für Bild, Bild auf Titelfolie hinter das Logo einsetzen Introduction to PDEs and Numerical Methods Tutorial 5. Finite difference methods equilibrium equation and iterative solvers Dr. Noemi
More informationWaves in a Shock Tube
Waves in a Shock Tube Ivan Christov c February 5, 005 Abstract. This paper discusses linear-wave solutions and simple-wave solutions to the Navier Stokes equations for an inviscid and compressible fluid
More informationDiscontinuous Galerkin methods for fractional diffusion problems
Discontinuous Galerkin methods for fractional diffusion problems Bill McLean Kassem Mustapha School of Maths and Stats, University of NSW KFUPM, Dhahran Leipzig, 7 October, 2010 Outline Sub-diffusion Equation
More informationIdentification Methods for Structural Systems. Prof. Dr. Eleni Chatzi Lecture March, 2016
Prof. Dr. Eleni Chatzi Lecture 4-09. March, 2016 Fundamentals Overview Multiple DOF Systems State-space Formulation Eigenvalue Analysis The Mode Superposition Method The effect of Damping on Structural
More informationWell-balanced DG scheme for Euler equations with gravity
Well-balanced DG scheme for Euler equations with gravity Praveen Chandrashekar praveen@tifrbng.res.in Center for Applicable Mathematics Tata Institute of Fundamental Research Bangalore 560065 Higher Order
More informationA recovery-assisted DG code for the compressible Navier-Stokes equations
A recovery-assisted DG code for the compressible Navier-Stokes equations January 6 th, 217 5 th International Workshop on High-Order CFD Methods Kissimmee, Florida Philip E. Johnson & Eric Johnsen Scientific
More informationLecture 9 Approximations of Laplace s Equation, Finite Element Method. Mathématiques appliquées (MATH0504-1) B. Dewals, C.
Lecture 9 Approximations of Laplace s Equation, Finite Element Method Mathématiques appliquées (MATH54-1) B. Dewals, C. Geuzaine V1.2 23/11/218 1 Learning objectives of this lecture Apply the finite difference
More informationThe Lattice Boltzmann method for hyperbolic systems. Benjamin Graille. October 19, 2016
The Lattice Boltzmann method for hyperbolic systems Benjamin Graille October 19, 2016 Framework The Lattice Boltzmann method 1 Description of the lattice Boltzmann method Link with the kinetic theory Classical
More informationChapter 3 Second Order Linear Equations
Partial Differential Equations (Math 3303) A Ë@ Õæ Aë áöß @. X. @ 2015-2014 ú GA JË@ É Ë@ Chapter 3 Second Order Linear Equations Second-order partial differential equations for an known function u(x,
More informationAlberto Bressan. Department of Mathematics, Penn State University
Non-cooperative Differential Games A Homotopy Approach Alberto Bressan Department of Mathematics, Penn State University 1 Differential Games d dt x(t) = G(x(t), u 1(t), u 2 (t)), x(0) = y, u i (t) U i
More informationMath background. Physics. Simulation. Related phenomena. Frontiers in graphics. Rigid fluids
Fluid dynamics Math background Physics Simulation Related phenomena Frontiers in graphics Rigid fluids Fields Domain Ω R2 Scalar field f :Ω R Vector field f : Ω R2 Types of derivatives Derivatives measure
More informationNUMERICAL METHODS FOR ENGINEERING APPLICATION
NUMERICAL METHODS FOR ENGINEERING APPLICATION Second Edition JOEL H. FERZIGER A Wiley-Interscience Publication JOHN WILEY & SONS, INC. New York / Chichester / Weinheim / Brisbane / Singapore / Toronto
More informationImplementation of Implicit Solution Techniques for Non-equilibrium Hypersonic Flows
Short Training Program Report Implementation of Implicit Solution Techniques for Non-equilibrium Hypersonic Flows Julian Koellermeier RWTH Aachen University Supervisor: Advisor: Prof. Thierry Magin von
More informationMagnetized Target Fusion: Insights from Mathematical Modelling
: Insights from Mathematical Modelling, UBC, PhD Candidate - Applied Math Brian Wetton, UBC, Supervisor Michael Ward, UBC, Committee Member Sandra Barsky, General Fusion, Committee Member November 28,
More information1. Differential Equations (ODE and PDE)
1. Differential Equations (ODE and PDE) 1.1. Ordinary Differential Equations (ODE) So far we have dealt with Ordinary Differential Equations (ODE): involve derivatives with respect to only one variable
More informationThe Shallow Water Equations
If you have not already done so, you are strongly encouraged to read the companion file on the non-divergent barotropic vorticity equation, before proceeding to this shallow water case. We do not repeat
More informationImplicit kinetic relaxation schemes. Application to the plasma physic
Implicit kinetic relaxation schemes. Application to the plasma physic D. Coulette 5, E. Franck 12, P. Helluy 12, C. Courtes 2, L. Navoret 2, L. Mendoza 2, F. Drui 2 ABPDE II, Lille, August 2018 1 Inria
More informationChapter 6. Finite Element Method. Literature: (tiny selection from an enormous number of publications)
Chapter 6 Finite Element Method Literature: (tiny selection from an enormous number of publications) K.J. Bathe, Finite Element procedures, 2nd edition, Pearson 2014 (1043 pages, comprehensive). Available
More informationImprovements of Unsteady Simulations for Compressible Navier Stokes Based on a RK/Implicit Smoother Scheme
Improvements of Unsteady Simulations for Compressible Navier Stokes Based on a RK/Implicit Smoother Scheme Oren Peles and Eli Turkel Department of Applied Mathematics, Tel-Aviv University In memoriam of
More information