Generalized Fourier Analyses of Semi-Discretizations of the Advection-Diffusion Equation
|
|
- Coleen Morris
- 5 years ago
- Views:
Transcription
1 Generalized Fourier Analyses of Semi-Discretizations of the Advection-Diffusion Equation Sponsored by: ASCI Algorithms - Advanced Spatial Discretization Project Mark A. Christon, homas E. Voth Computational Physics R&D (931) Mario J. Martinez Multiphase ransport Processes Department (9114) August 5, 00
2 Overview Background & Motivation Primer on Fourier Analysis 1D and D Results Summary & Conclusions
3 Choosing the best numerical method can be difficult given the plethora of methods available Developers embarking on a new code effort are faced with an array of choices for numerical methods Unstructured vs. structured grids Mesh-full vs. mesh-free Finite element vs. finite volume Local conservation vs. global conservation Use of high-order discretizations p vs. h-refinement Formulation differences, e.g., Galerkin (weighted residual) vs. aylor-series, can make many aspects of side-by-side comparisons difficult his work constitutes a first-step in a multi-methods comparison intended to identify strengths and weaknesses in the context of advection-diffusion processes.
4 Understanding the behavior of numerical methods in a common framework is the basis for comparison Galerkin FEM Control Volume FEM o o u = cos(37 ), v= sin(37 ) Convected Gaussian cone, from Gresho and Sani, pg. 4, Wiley, 1998
5 Multi-methods comparisons attempt to identify strengths and weaknesses for a broad range of attributes Numerical performance is a broad term, and includes truncation error, consistency, stability, convergence, etc. Some desirable attributes: Consistency: discretization recovers the PDE as x 0, t 0 Stability: n+ 1 n u u, e.g., stable for all t and insensitive to roundoff Convergence: numerical solutions approach the PDE solution as x For linear problems, automatic given stability & consistency via Lax s equivalence theorem Conservation: 0 th -moment should be globally conserved (at a minimum) Non-dispersive: all wavelengths propagate at the true advective speed Non-dissipative: there should be no numerical diffusion Shape-preserving: positivity, monotonicity, lack of spurious oscillations Compact operators: local stencils for the difference operators Computationally efficient: reasonable memory and CPU requirements See Baptista, Adams, Gresho, Benchmarks for the transport equation: the convection-diffusion forum and beyond, QSACOM, V47, , t 0
6 A number of global error measures may be used to assess numerical performance 1 L error (integral measure of error) L error (integral measure of squared error) L error (maximum local) ˆ dω ˆ dω ( ˆ) ˆ dω dω max ˆ max ˆ Convergence may not be observed in all norms on non-smooth data. ˆ exact solution Amplitude error or error in the peak value Maximum negative value Measure of phase error Measure of spreading ( ˆ ) ˆ max max max max, neg ˆ max ( M Mˆ ) Mˆ, M = x dω dω ( x M ) dω ( x Mˆ ) d ˆ Ω emperature Amplitude error Spreading Phase error x 3.00
7 he generalized Fourier analysis places all methods on an equal footing -- regardless of the formulation Generalized Fourier analysis considers the spectral behavior of semi-discretizations of the advection-diffusion equation Inputs + ci = α t wavenumber : k direction : θ aspect ratio : γ ransfer Function Semi Discretization M + A() c + K = 0 c ( 1 φ ) M = φm + M u = ccos ϑ, v= csinϑ k = kcos θ, k = ksinθ x y l Outputs Phase speed : c c v g v y : x g Group velocity, c c Discrete diffusivity : α α Artificial diffusivity : 1 art, P e art P c x e = α art Short wavelength behavior Grid anisotropy in terms of γθ, Asymptotic truncation error, O Grid resolution requirements p ( x )
8 he Fourier analysis places all operators on a regular grid configuration with periodic conditions FDM FEM CVFEM Fundamental sol n for A-D: mn, () t = Aexp ik( m xcosθ + n ysinθ) i ωt k αt Substitution into semi-discrete equation yields: ω, α, c = ω k, etc. For a non-dispersive medium, waves propagate at the true velocity, e.g., in 1-D, c=ω k Spatial discretization results in dispersive behavior, i.e., waves propagate at a velocity that is wavelength dependent, c ( ) = c k riangular elements require consideration of multiple regular grid configurations
9 he group velocity describes the propagation of wave packets in a dispersive medium t Group velocity is defined as: vg k[ kc] ( = kω) Wave packets consist of short-wavelength signals modulating a slow-moving envelope κ ( x, t) = aκ exp[ ικ ( x vg ( k))] t ( ) exp[ ιk( x c( k) t) ] Energy, x, t dx, contained in the wave packet moves at the group velocity Energy may not propagate with the flow, e.g., positive phase with negative group c v g -0.5 κ k x oscillations, which may be stationary, appear to move when they modulate a long-wavelength envelope 0.5 Envelope moves at the group velocity X
10 he second-order upwind FDM provides a simple prototype for the generalized Neumann analysis Complete 1-D nodal equation for M + A() c + K = 0 c m + { m 4m 1+ 3m} x α { 1 1} 0 m m + m+ = x Substitute fundamental solution, calculate phase, group, etc. m () t = Aexp ikm x i ωt k αt { c [ ( ) ( )]} 4sin k x sin k x i ω + αk = i + x α x m m 1 c x [ cos( k x) ] + 3+ cos( k x) 4cos( k x) [ ] Decomposition of the advection operator is automatic 1 1 Askew = A A, A sym = A+ A x m m + 1 m + Skew-symmetric (non-dissipative) part of advection Symbol, Aˆ Symmetric (dissipative) part of advection
11 Non-dimensional parameters characterize phase and group speed, discrete and artificial diffusivity Imaginary part of symbol, Im( Aˆ ), yields the phase speed c 1 = [ 4sin( k x) sin( k x) ] c k x Real part of symbol, Re( Aˆ ), yields the discrete and artificial diffusivity α = α + α art α 1 = [ cos( k x) ] α k x 1 α art 1 = = + art c x P k x e [ 3 cos( k x) 4cos( k x) ] Characterizing the discretization: Re( Aˆ ( k)) = 0 for all k, discretization is neutrally dissipative ( ˆ ) ( ˆ ) Re A( k) < 0 for some k, discretization is dissipative Re A( k) > 0 for some k, discretization is unstable Necessary & Sufficient: A() c is skew-symmetric M is symmetric, e.g., Galerkin FEM
12 ~c/c Phase Speed Group Speed Artificial Diffusivity x/λ k x π = x λ SOU Skew-Symmetric Advection A 1( ) skew = A A SOU Full Advection Operator ~v gx /c t = 3.75 (u t/ x) x x/λ ~ 3 x t = 3.75 (u t/ x) x 1/Pe art t = 7.5 (u t/ x) x t = 7.5 (u t/ x) ~ 3 x x
13 here is a direct relationship between order of accuracy and the flatness of phase (and diffusivity) near k x = 0 runcation error: dmn, 1 E = M A( c ) mn, dt p+ 1 p E = C x H. O. p 1 x + + Phase error: c~ p 1 = C( k ι x) + H. O. c E and phase (and diffusivity) error yield order accuracy, p: p ( x ) as 0 E = O x c~ p 1 = O ( k x) c as ( ) k x 0 c~ /c c ~ / c CD-FDM 4 CD-FDM FEM-M c k x/ π ( = x/ λ) FEM-M c
14 Analysis of non-linear methods yields an operating range for phase speed and artificial diffusivity Godunov-type central scheme (Kurganov & admor, 000) using 1.5 minmod slope limiter φ m+ 1 m = m x m x ( ) = ( ( )) φ θ min 0, max 1, θ, θ m m m m = m m+ 1 m 1 m c~ /c 0.5 FOU/CD φm 1 = 0, φm = 0 φm 1 = 0, φm = 1 φm 1 = 1, φm = 0 φm 1 = 1, φm = 1 φm 1 = 0, φm = 0 0 φm 1 = 0, φm = 1 φm 1 = 1, φm = 0 φm 1 = 1, φm = 1 1/Pe art FOU φm 1 = 0, φm = 0 φm 1 = 0, φm = 1 φm 1 = 1, φm = 0 m- m-1 m m+1 0
15 A broad cross-section of methods have been considered in this study First-order upwind (FOU) Second-order upwind (SOU) Second-order central differences (Centered FDM) hird-order upwind (OU) Fromm s method (semi-discrete version) QUICK (Quadratic Upwind Interpolation with Correct Kinematics) Least squares reconstruction (LSR(0), LSR(-1)) Galerkin finite element (FEM) Streamline-upwind Petrov-Galerkin (FEM-SUPG) Control-volume finite element (CVFEM) Streamline-upwind control-volume finite element (CVFEM-SUCV) β SUPG Opt. β = 1 15 = opt SUCV Opt. 1
16 FEM with consistent mass delivers the best overall phase speed with superconvergent asymptotic behavior Fromm's Method E λ/ x for 5% 1% error error c/c c/c Centered FDM QUICK OU CVFEM-M c FEM-M c CVFEM SUCV (β = ½) FEM SUPG (β opt ) SOU FOU FEM-M l CVFEM-M l SOU OU QUICK Fromm s FEM-M c FEM-SUPG (β opt ) FEM-SUPG (β= ½) CVFEM-M c CVFEM-SUCV (β opt ) Ο( x 4 ) Ο( x 4 ) Ο( x 6 ) Ο( x 4 ) CVFEM-SUCV (β= ½)
17 1 0 v /c g g v /c FEM with consistent mass delivers the best overall group speed with superconvergent asymptotic behavior QUICK OU Centered FDM CVFEM-M c FEM-M c FEM SUPG (β opt ) CVFEM SUCV (β = ½) SOU Fromm's -6 FOU FEM-M l CVFEM-M l SOU OU QUICK Fromm s FEM-M c FEM-SUPG (β opt ) FEM-SUPG (β= ½) CVFEM-M c Method CVFEM-SUCV (β opt ) CVFEM-SUCV (β= ½) AR Ο( x 4 ) Ο( x 4 ) Ο( x 6 ) Ο( x 4 ) 5% error λ/ x for 1% error
18 180 wo-dimensional spatial discretization introduces wavelength AND direction dependent behavior c ~ / θ c c ~ / ~c/c c SOU , 90 o.5, 67.5 o x/λ 45 o y x Grid Definitions θ θ x/λ ~c/c
19 Integrated anisotropy and error metrics used to allow objective comparison of two-dimensional results Anisotropy metric: Error metric: σ ε k ( c~ ( θ, k) c~ ( k) ) dθ dk = k θ ( c ( θ, k) c) dθ dk = θ ~ c ~ / c c ~ / c c ~ / c CVFEM-SUCV (β = ½) small anisotropy (7.3e-) small error (7.1e-) 70 FOU large anisotropy (1.7e-1) large error (3.3e-1) 70 LSR(-1) small anisotropy (8.8e-) large error (.e-1)
20 Integrated anisotropy and phase error metrics suggest FEM-SUPG and CVFEM-SUCV schemes are best c ~ / c FOU SOU OU QUICK Method Anisotropy, s 1.7e-1 1.4e-1 1.4e-1 1.5e-1 Error, e 3.3e-1.e-1.e-1.4e-1 Fromm s 1.4e-1 1.8e FOU c ~ / c LSR(0) LSR(-1) FEM-M c FEM-SUPG (β opt ) FEM-SUPG (β= ½) CVFEM-M c CVFEM-SUCV (β opt ) 9.1e- 8.8e- 1.1e-1 7.9e- 7.8e- 1.3e-1 1.1e-1.1e-1.e-1 1.4e-1 7.3e- 1.3e-1.1e-1 1.5e-1 70 CVFEM-SUCV (β= ½) CVFEM-SUCV (β = ½) 7.3e- 7.1e-
21 Unlike the physical problem, the discrete problem introduces a diffusivity that varies by wavelength Consider 1-D transient diffusion in a bar: 0.5 initial = 0 = 0 Analytical sol n is: x [ ] [ exp k t] ( x t) B sin( k x), = n n α = n 1 Numerical sol n. is similar, but with wavelength dependent diffusivity: i [ ] [ exp ~ α k t] () t B sin( k x ) = = n 1 n n i n n n α ~ x 1.5 α 0.5 x / λ
22 SUPG and SUCV formulations with consistent mass and β = β opt yield best discrete diffusivity results 1.50 α~ / α FEM - M c CVFEM - M c FEM/CVFEM - M l Method CD-FDM FEM- M L CVFEM-M L E λ/ x for 5% 1% error error FEM-M c FEM - SUPG β opt CVFEM - SUCV β opt FEM-SUPG (β opt ) FEM-SUPG (β= ½) CVFEM-M c α ~ /α 0 CVFEM-SUCV (β opt ) CVFEM-SUCV (β=½) CVFEM - SUCV β=1/ FEM - SUPG β=1/ Note: β opt is NO the optimal choice for SUCV phase 0 x /λ
23 FEM-M c and CVFEM-M c demonstrate best anisotropy and error respectively for discrete diffusivity ~ α /α Method CD-FDM FEM-M c Anisotropy, s 1.7e-1.8e- Error, e 1.8e-1 1.8e-1 FEM-SUPG (β opt ) 8.6e- 9.8e- FEM-SUPG (β=½) 1.7e-1 3.1e FEM-M c CVFEM-M c 5.3e- 5.3e- ~ α /α CVFEM-SUCV (β opt ) CVFEM-SUCV (β=½) 5.1e- 7.3e- 1.1e-1 3.e CVFEM-SUCV (β = ½)
24 Artificial diffusivity should increase with wave number suggesting poor behavior for FOU /Pe art FOU artificial diffusivity acts at ALL wavelengths Second-Order Upwind (FOU) First-Order Upwind (FOU) Fromm's Method c/c CVFEM-M c FEM-M c CVFEM SUCV (β = ½) FEM SUPG (β opt ) hird-order Upwind (FOU) QUICK 1.4 1/Pe art 1. FEM - SUPG β = 1/ CVFEM - SUPG β=1/ FEM - SUPG β opt CVFEM - SUPG β opt x /λ v g / c CVFEM-M c FEM-M c FEM SUPG (β opt ) CVFEM SUCV (β = ½) -6
25 Quadratic energy damping characteristics suggest that the FEM and CVFEM SU variants are good choices Q = Constant mode isn t damped { t t M M } t= τ + x/c t= τ t M t= τ - QUICK - FEM - SUPG β = 1/ hird-order Upwind (OU) - - FEM - SUPG β opt Q - Fromm's Method Second-Order Upwind (SOU) Q - CVFEM - SUCV β = 1/ CVFEM - SUCV β opt - - First-Order Upwind (FOU) - x /λ Finite Difference schemes - x /λ FEM-SUPG and CVFEM-SUCV schemes Signal not damped in x/c time scale
26 In the discrete world, Summary and Conclusions Waves don t propagate at the advective velocity, and don t always propagate in the direction of the wave vector Information doesn t diffuse at the continuum rate Grid bias will be present Fourier analysis can quantify these errors Some clear losers and not so obvious winners Clear losers are FOU and SOU SOU provides the WORS phase and group, requiring the most resolution for a fixed accuracy level FOU, lumped mass FEM and CVFEM also perform relatively poorly Choice of a winner is more difficult Upwind (SU) variants of FEM and CVFEM perform well with careful choice of stabilization parameter FEM variants w. consistent mass yield superconvergent phase and group with minimum mesh resolution requirements OU also exhibits superconvergent behavior Analysis of non-linear methods is possible, but results are not as sharp as for linear methods
27 A number of diagnostic metrics can be used to assess the numerical performance For + ci = 0, with ic = 0 with no-flux or periodic BC s t Global conservation Quadratic ( energy ) conservation Ability to preserve peaks Ability to preserve gradients Ability to preserve curvature Ω Ω Ω Ω Ω dω 4 dω dω i dω { } d Ω
28 Phase speed is the projection of the fluid velocity in the wave direction -- the apparent velocity Scalar advection: + u + v = 0 t c = [ uv, ] t x y Use a fundamental sol n y ( x, t) = Aexp[ ι( kxx+ kyy) ιωt] Solve for the circular frequency k t ω = ku x + kv y =kc and phase speed t c kc/ k For a non-dispersive medium, waves propagate at the true velocity, e.g., in 1-D, c = u Spatial discretization results in dispersive behavior, i.e., waves don t propagate at the true velocity and c u Waves propagate at a velocity that is wavelength dependent (Adapted from Gresho s notes, aiwan course, 1989) k t crest trough = [ k, k ] c x y x
29 wo-dimensional analysis reveals angular dependence (grid-bias) in phase, group, and diffusivities θ 1 Phase Speed x/λ ~c/c θ 1 Group Velocity (x-dir) v~ gx /c θ 18.0 Discrete Diffusivity α/α θ 1 Artificial Diffusivity 1/P art e
30 Results Outline Phase and Group Speed One-Dimensional phase and group results wo-dimensional phase results Discrete Diffusivity One-Dimensional results wo-dimensional results Artificial Diffusivity One-Dimensional results wo-dimensional results
31 Grid aspect ratio modifies the anisotropic phase behavior of the two-dimensional discretizations FEM-SUPG (β opt ) c ~ / γ = 1 c ~ / c c Method FOU SOU OU QUICK Fromm s LSR(0) FEM-M c FEM-SUPG (β opt ) FEM-SUPG (β= ½) CVFEM-M c CVFEM-SUCV (β opt ) Anisotropy, σ γ = 1 γ = 1/ 1.7e-1 1.8e-1 1.4e-1 1.e-1 1.4e-1 1.4e-1 1.5e-1 1.5e-1 1.4e-1 1.e-1 9.1e- 1.e-1 1.1e-1 9.4e- 7.9e- 5.6e- 7.8e- 6.5e- 1.3e-1 1.3e-1 1.1e-1 9.4e- 70 γ = 1/ CVFEM-SUCV (β= ½) 7.3e- 4.9e-
32 Artificial diffusivity should increase with wave number suggesting poor behavior for FOU /Pe art 1.4 First-Order Upwind (FOU) Second-Order Upwind (FOU) Fromm's Method hird-order Upwind (FOU) QUICK FOU SOU OU QUICK Fromm s Method FEM-SUPG (β opt ) FEM-SUPG (β=½).e. Ο(1) 1. CVFEM-SUPG (β opt ) 1/Pe art FEM - SUPG β = 1/ CVFEM - SUPG β=1/ CVFEM-SUPG (β=½) FEM - SUPG β opt FOU artificial diffusivity acts at ALL wavelengths CVFEM - SUPG β opt x /λ
33 he spectral dependence of the artificial diffusivity does not completely explain its damping effects Given the fundamental solution: x, t = ( x)exp( α k ( ) ) 0 art t Damping of depends on α art and k: x, t t = α k ( x)exp( ( ) α k t) art Similar dependency for Q= t M: ( ) ( Q α k exp α k t) t = art 0 art art Q 0 k = π t = 0 t =.01/α art Q(t)/Q(t=0) k = 4π x tα art
34 Least-Squares Gradient Reconstruction schemes minimize angular dependence of artificial diffusivity art 1/ Pe 0 FOU SOU OU Method Anisotropy 7.8e-.1e-1.1e-1 Diffusivity 6.9e-1 5.1e-1 1.7e-1 QUICK.1e-1.5e LSR(0) art 1/ Pe Fromm s LSR(0) LSR(-1) FEM-SUPG (β opt ) FEM-SUPG (β=½).1e-1 4.5e- 4.8e- 3.8e-1 4.5e-1.6e-1 3.1e-1 6.1e-1 1.7e-1.4e CVFEM-SUCV (β opt ) 3.4e-1 1.6e-1 CVFEM-SUCV (β=½) 4.1e-1.3e-1 70 FEM-SUPG (β = ½)
PDEs, 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 informationNumerical Solutions for Hyperbolic Systems of Conservation Laws: from Godunov Method to Adaptive Mesh Refinement
Numerical Solutions for Hyperbolic Systems of Conservation Laws: from Godunov Method to Adaptive Mesh Refinement Romain Teyssier CEA Saclay Romain Teyssier 1 Outline - Euler equations, MHD, waves, hyperbolic
More information3.4. Monotonicity of Advection Schemes
3.4. Monotonicity of Advection Schemes 3.4.1. Concept of Monotonicity When numerical schemes are used to advect a monotonic function, e.g., a monotonically decreasing function of x, the numerical solutions
More informationProf. dr. A. Achterberg, Astronomical Dept., IMAPP, Radboud Universiteit
Prof. dr. A. Achterberg, Astronomical Dept., IMAPP, Radboud Universiteit Central concepts: Phase velocity: velocity with which surfaces of constant phase move Group velocity: velocity with which slow
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 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 information2.3. Quantitative Properties of Finite Difference Schemes. Reading: Tannehill et al. Sections and
3 Quantitative Properties of Finite Difference Schemes 31 Consistency, Convergence and Stability of FD schemes Reading: Tannehill et al Sections 333 and 334 Three important properties of FD schemes: Consistency
More information5. FVM discretization and Solution Procedure
5. FVM discretization and Solution Procedure 1. The fluid domain is divided into a finite number of control volumes (cells of a computational grid). 2. Integral form of the conservation equations are discretized
More informationIn Proc. of the V European Conf. on Computational Fluid Dynamics (ECFD), Preprint
V European Conference on Computational Fluid Dynamics ECCOMAS CFD 2010 J. C. F. Pereira and A. Sequeira (Eds) Lisbon, Portugal, 14 17 June 2010 THE HIGH ORDER FINITE ELEMENT METHOD FOR STEADY CONVECTION-DIFFUSION-REACTION
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 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 informationNUMERICAL SOLUTION OF CONVECTION DIFFUSION EQUATIONS USING UPWINDING TECHNIQUES SATISFYING THE DISCRETE MAXIMUM PRINCIPLE
Proceedings of the Czech Japanese Seminar in Applied Mathematics 2005 Kuju Training Center, Oita, Japan, September 15-18, 2005 pp. 69 76 NUMERICAL SOLUTION OF CONVECTION DIFFUSION EQUATIONS USING UPWINDING
More informationTime-dependent variational forms
Time-dependent variational forms Hans Petter Langtangen 1,2 1 Center for Biomedical Computing, Simula Research Laboratory 2 Department of Informatics, University of Oslo Oct 30, 2015 PRELIMINARY VERSION
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 informationSuperposition of electromagnetic waves
Superposition of electromagnetic waves February 9, So far we have looked at properties of monochromatic plane waves. A more complete picture is found by looking at superpositions of many frequencies. Many
More informationOne Dimensional Convection: Interpolation Models for CFD
One Dimensional Convection: Interpolation Models for CFD ME 448/548 Notes Gerald Recktenwald Portland State University Department of Mechanical Engineering gerry@pdx.edu ME 448/548: 1D Convection-Diffusion
More informationFEM-Level Set Techniques for Multiphase Flow --- Some recent results
FEM-Level Set Techniques for Multiphase Flow --- Some recent results ENUMATH09, Uppsala Stefan Turek, Otto Mierka, Dmitri Kuzmin, Shuren Hysing Institut für Angewandte Mathematik, TU Dortmund http://www.mathematik.tu-dortmund.de/ls3
More informationRelevant self-assessment exercises: [LIST SELF-ASSESSMENT EXERCISES HERE]
Chapter 5 Fourier Analysis of Finite Difference Methods In this lecture, we determine the stability of PDE discretizations using Fourier analysis. First, we begin with Fourier analysis of PDE s, and then
More informationStabilization Techniques for Finite Element Analysis of Convection-Diffusion Problems
INTERNATIONAL CENTER FOR NUMERICAL METHODS IN ENGINEERING Stabilization Techniques for Finite Element Analysis of Convection-Diffusion Problems E. Oñate M. Manzan Publication CIMNE Nº-183, February 000
More informationAlgebraic flux correction and its application to convection-dominated flow. Matthias Möller
Algebraic flux correction and its application to convection-dominated flow Matthias Möller matthias.moeller@math.uni-dortmund.de Institute of Applied Mathematics (LS III) University of Dortmund, Germany
More informationOn finite element methods for 3D time dependent convection diffusion reaction equations with small diffusion
On finite element methods for 3D time dependent convection diffusion reaction equations with small diffusion Volker John and Ellen Schmeyer FR 6.1 Mathematik, Universität des Saarlandes, Postfach 15 11
More informationElectromagnetic wave propagation. ELEC 041-Modeling and design of electromagnetic systems
Electromagnetic wave propagation ELEC 041-Modeling and design of electromagnetic systems EM wave propagation In general, open problems with a computation domain extending (in theory) to infinity not bounded
More informationHydraulic Modelling for Drilling Automation
Hydraulic Modelling for Drilling Automation CASA Day Harshit Bansal April 19, 2017 Where innovation starts Team: Supervisors at TU/e : W.H.A. Schilders, N. van de Wouw, B. Koren, L. Iapichino Collaborators:
More informationDesign of optimal Runge-Kutta methods
Design of optimal Runge-Kutta methods David I. Ketcheson King Abdullah University of Science & Technology (KAUST) D. Ketcheson (KAUST) 1 / 36 Acknowledgments Some parts of this are joint work with: Aron
More informationChapter 17. Finite Volume Method The partial differential equation
Chapter 7 Finite Volume Method. This chapter focusses on introducing finite volume method for the solution of partial differential equations. These methods have gained wide-spread acceptance in recent
More informationPartial differential equations
Partial differential equations Many problems in science involve the evolution of quantities not only in time but also in space (this is the most common situation)! We will call partial differential equation
More informationTutorial School on Fluid Dynamics: Aspects of Turbulence Session I: Refresher Material Instructor: James Wallace
Tutorial School on Fluid Dynamics: Aspects of Turbulence Session I: Refresher Material Instructor: James Wallace Adapted from Publisher: John S. Wiley & Sons 2002 Center for Scientific Computation and
More information7 Hyperbolic Differential Equations
Numerical Analysis of Differential Equations 243 7 Hyperbolic Differential Equations While parabolic equations model diffusion processes, hyperbolic equations model wave propagation and transport phenomena.
More informationDivergence Formulation of Source Term
Preprint accepted for publication in Journal of Computational Physics, 2012 http://dx.doi.org/10.1016/j.jcp.2012.05.032 Divergence Formulation of Source Term Hiroaki Nishikawa National Institute of Aerospace,
More informationAdvanced numerical methods for nonlinear advectiondiffusion-reaction. Peter Frolkovič, University of Heidelberg
Advanced numerical methods for nonlinear advectiondiffusion-reaction equations Peter Frolkovič, University of Heidelberg Content Motivation and background R 3 T Numerical modelling advection advection
More informationSpatial and Modal Superconvergence of the Discontinuous Galerkin Method for Linear Equations
Spatial and Modal Superconvergence of the Discontinuous Galerkin Method for Linear Equations N. Chalmers and L. Krivodonova March 5, 014 Abstract We apply the discontinuous Galerkin finite element method
More informationAMS 529: Finite Element Methods: Fundamentals, Applications, and New Trends
AMS 529: Finite Element Methods: Fundamentals, Applications, and New Trends Lecture 25: Introduction to Discontinuous Galerkin Methods Xiangmin Jiao SUNY Stony Brook Xiangmin Jiao Finite Element Methods
More informationECE 695 Numerical Simulations Lecture 12: Thermomechanical FEM. Prof. Peter Bermel February 6, 2017
ECE 695 Numerical Simulations Lecture 12: Thermomechanical FEM Prof. Peter Bermel February 6, 2017 Outline Static Equilibrium Dynamic Equilibrium Thermal transport mechanisms Thermal transport modeling
More informationStabilization and Acceleration of Algebraic Multigrid Method
Stabilization and Acceleration of Algebraic Multigrid Method Recursive Projection Algorithm A. Jemcov J.P. Maruszewski Fluent Inc. October 24, 2006 Outline 1 Need for Algorithm Stabilization and Acceleration
More informationStable Semi-Discrete Schemes for the 2D Incompressible Euler Equations
Stable Semi-Discrete Schemes for the D Incompressible Euler Equations Doron Levy Department of Mathematics Stanford University dlevy@math.stanford.edu http://math.stanford.edu/ dlevy Maryland, May 004
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 informationHp-Adaptivity on Anisotropic Meshes for Hybridized Discontinuous Galerkin Scheme
Hp-Adaptivity on Anisotropic Meshes for Hybridized Discontinuous Galerkin Scheme Aravind Balan, Michael Woopen and Georg May AICES Graduate School, RWTH Aachen University, Germany 22nd AIAA Computational
More informationDue Tuesday, November 23 nd, 12:00 midnight
Due Tuesday, November 23 nd, 12:00 midnight This challenging but very rewarding homework is considering the finite element analysis of advection-diffusion and incompressible fluid flow problems. Problem
More informationFinite difference method for solving Advection-Diffusion Problem in 1D
Finite difference method for solving Advection-Diffusion Problem in 1D Author : Osei K. Tweneboah MATH 5370: Final Project Outline 1 Advection-Diffusion Problem Stationary Advection-Diffusion Problem in
More informationWaves Primer II Basic wave form. Waves are: 1. Important 2. Easy 3. Neat
Waves Primer II Waves are: 1. Important 2. Easy 3. Neat 23.1 Basic wave form The wave form that you will almost always use is a = Re[A c e i(kx+ly ωt) ] (23.1) Where a variable A c waveamplitude(complexinthiscase)
More information3.3. Phase and Amplitude Errors of 1-D Advection Equation
3.3. Phase and Amplitude Errors of 1-D Advection Equation Reading: Duran section 2.4.2. Tannehill et al section 4.1.2. The following example F.D. solutions of a 1D advection equation show errors in both
More informationNote the diverse scales of eddy motion and self-similar appearance at different lengthscales of the turbulence in this water jet. Only eddies of size
L Note the diverse scales of eddy motion and self-similar appearance at different lengthscales of the turbulence in this water jet. Only eddies of size 0.01L or smaller are subject to substantial viscous
More informationOn the design of higher-order FEM satisfying the discrete maximum principle
On the design of higher-order FEM satisfying the discrete maximum principle Dmitri Kuzmin Institute of Applied Mathematics (LS III), University of Dortmund Vogelpothsweg 87, D-44227, Dortmund, Germany
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 informationOn discontinuity capturing methods for convection diffusion equations
On discontinuity capturing methods for convection diffusion equations Volker John 1 and Petr Knobloch 2 1 Universität des Saarlandes, Fachbereich 6.1 Mathematik, Postfach 15 11 50, 66041 Saarbrücken, Germany,
More informationThe RAMSES code and related techniques I. Hydro solvers
The RAMSES code and related techniques I. Hydro solvers Outline - The Euler equations - Systems of conservation laws - The Riemann problem - The Godunov Method - Riemann solvers - 2D Godunov schemes -
More informationA Central Compact-Reconstruction WENO Method for Hyperbolic Conservation Laws
A Central Compact-Reconstruction WENO Method for Hyperbolic Conservation Laws Kilian Cooley 1 Prof. James Baeder 2 1 Department of Mathematics, University of Maryland - College Park 2 Department of Aerospace
More informationSemi-Lagrangian Formulations for Linear Advection Equations and Applications to Kinetic Equations
Semi-Lagrangian Formulations for Linear Advection and Applications to Kinetic Department of Mathematical and Computer Science Colorado School of Mines joint work w/ Chi-Wang Shu Supported by NSF and AFOSR.
More informationLecture Notes on Numerical Schemes for Flow and Transport Problems
Lecture Notes on Numerical Schemes for Flow and Transport Problems by Sri Redeki Pudaprasetya sr pudap@math.itb.ac.id Department of Mathematics Faculty of Mathematics and Natural Sciences Bandung Institute
More informationSingularly Perturbed Partial Differential Equations
WDS'9 Proceedings of Contributed Papers, Part I, 54 59, 29. ISN 978-8-7378--9 MTFYZPRESS Singularly Perturbed Partial Differential Equations J. Lamač Charles University, Faculty of Mathematics and Physics,
More informationNew DPG techniques for designing numerical schemes
New DPG techniques for designing numerical schemes Jay Gopalakrishnan University of Florida Collaborator: Leszek Demkowicz October 2009 Massachusetts Institute of Technology, Boston Thanks: NSF Jay Gopalakrishnan
More informationPart 1. The diffusion equation
Differential Equations FMNN10 Graded Project #3 c G Söderlind 2016 2017 Published 2017-11-27. Instruction in computer lab 2017-11-30/2017-12-06/07. Project due date: Monday 2017-12-11 at 12:00:00. Goals.
More informationScientific Computing I
Scientific Computing I Module 8: An Introduction to Finite Element Methods Tobias Neckel Winter 2013/2014 Module 8: An Introduction to Finite Element Methods, Winter 2013/2014 1 Part I: Introduction to
More informationConstruction of very high order Residual Distribution Schemes for steady problems
Construction of very high order Residual Distribution Schemes for steady problems Rémi Abgrall, Mario Ricchiuto, Cédric Tavé, Nadège Villedieu and Herman Deconinck Mathématiques Appliquées de Bordeaux,
More informationSolution Methods. Steady State Diffusion Equation. Lecture 04
Solution Methods Steady State Diffusion Equation Lecture 04 1 Solution methods Focus on finite volume method. Background of finite volume method. Discretization example. General solution method. Convergence.
More informationCFD in Heat Transfer Equipment Professor Bengt Sunden Division of Heat Transfer Department of Energy Sciences Lund University
CFD in Heat Transfer Equipment Professor Bengt Sunden Division of Heat Transfer Department of Energy Sciences Lund University email: bengt.sunden@energy.lth.se CFD? CFD = Computational Fluid Dynamics;
More informationLet s consider nonrelativistic electrons. A given electron follows Newton s law. m v = ee. (2)
Plasma Processes Initial questions: We see all objects through a medium, which could be interplanetary, interstellar, or intergalactic. How does this medium affect photons? What information can we obtain?
More informationLecture 2. Turbulent Flow
Lecture 2. Turbulent Flow Note the diverse scales of eddy motion and self-similar appearance at different lengthscales of this turbulent water jet. If L is the size of the largest eddies, only very small
More informationRECENT DEVELOPMENTS IN VERY HIGH ORDER RESIDUAL DISTRIBUTION SCHEMES
RECENT DEVELOPMENTS IN VERY HIGH ORDER RESIDUAL DISTRIBUTION SCHEMES FOR INVISCID AND VISCOUS PROBLEMS. R. Abgrall Team Bacchus INRIA Bordeaux Sud Ouest and UniversitédeBordeauxTalence,France Roscoff,
More informationForced Oscillations in a Linear System Problems
Forced Oscillations in a Linear System Problems Summary of the Principal Formulas The differential equation of forced oscillations for the kinematic excitation: ϕ + 2γ ϕ + ω 2 0ϕ = ω 2 0φ 0 sin ωt. Steady-state
More informationRECENT DEVELOPMENTS IN COMPUTATIONAL REACTOR ANALYSIS
RECENT DEVELOPMENTS IN COMPUTATIONAL REACTOR ANALYSIS Dean Wang April 30, 2015 24.505 Nuclear Reactor Physics Outline 2 Introduction and Background Coupled T-H/Neutronics Safety Analysis Numerical schemes
More informationAdvanced numerical methods for transport and reaction in porous media. Peter Frolkovič University of Heidelberg
Advanced numerical methods for transport and reaction in porous media Peter Frolkovič University of Heidelberg Content R 3 T a software package for numerical simulation of radioactive contaminant transport
More informationFourier analysis for discontinuous Galerkin and related methods. Abstract
Fourier analysis for discontinuous Galerkin and related methods Mengping Zhang and Chi-Wang Shu Abstract In this paper we review a series of recent work on using a Fourier analysis technique to study the
More informationSolar Physics & Space Plasma Research Center (SP 2 RC) MHD Waves
MHD Waves Robertus vfs Robertus@sheffield.ac.uk SP RC, School of Mathematics & Statistics, The (UK) What are MHD waves? How do we communicate in MHD? MHD is kind! MHD waves are propagating perturbations
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 informationLecture Notes on Numerical Schemes for Flow and Transport Problems
Lecture Notes on Numerical Schemes for Flow and Transport Problems by Sri Redeki Pudaprasetya sr pudap@math.itb.ac.id Department of Mathematics Faculty of Mathematics and Natural Sciences Bandung Institute
More informationWaves, the Wave Equation, and Phase Velocity. We ll start with optics. The one-dimensional wave equation. What is a wave? Optional optics texts: f(x)
We ll start with optics Optional optics texts: Waves, the Wave Equation, and Phase Velocity What is a wave? f(x) f(x-) f(x-) f(x-3) Eugene Hecht, Optics, 4th ed. J.F. James, A Student's Guide to Fourier
More informationOn the relationship of local projection stabilization to other stabilized methods for one-dimensional advection-diffusion equations
On the relationship of local projection stabilization to other stabilized methods for one-dimensional advection-diffusion equations Lutz Tobiska Institut für Analysis und Numerik Otto-von-Guericke-Universität
More informationWave propagation in an inhomogeneous plasma
DRAFT Wave propagation in an inhomogeneous plasma Felix I. Parra Rudolf Peierls Centre for Theoretical Physics, University of Oxford, Oxford OX NP, UK This version is of 7 February 208. Introduction In
More informationDimensionality influence on energy, enstrophy and passive scalar transport.
Dimensionality influence on energy, enstrophy and passive scalar transport. M. Iovieno, L. Ducasse, S. Di Savino, L. Gallana, D. Tordella 1 The advection of a passive substance by a turbulent flow is important
More information4.3 Eigen-Analysis of Spectral Derivative Matrices 203
. Eigen-Analysis of Spectral Derivative Matrices 8 N = 8 N = 6 6 6 8 8 8 8 N = 6 6 N = 6 6 6 (Round off error) 6 8 8 6 8 6 6 Fig... Legendre collocation first-derivative eigenvalues computed with 6-bit
More informationThree remarks on anisotropic finite elements
Three remarks on anisotropic finite elements Thomas Apel Universität der Bundeswehr München Workshop Numerical Analysis for Singularly Perturbed Problems dedicated to the 60th birthday of Martin Stynes
More informationLATTICE BOLTZMANN METHOD AND DIFFUSION IN MATERIALS WITH LARGE DIFFUSIVITY RATIOS
THERMAL SCIENCE: Year 27, Vol. 2, No. 3, pp. 73-82 73 LATTICE BOLTZMANN METHOD AND DIFFUSION IN MATERIALS WITH LARGE DIFFUSIVITY RATIOS by Edouard WALTHER *, Rachid BENNACER, and Caroline DE SA LMT, ENS
More informationA Petrov Galerkin formulation for the alpha-interpolation of FEM and FDM stencils. Applications to the Helmholtz equation
INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN ENGINEERING Int. J. Numer. Meth. Engng ; : 5 Published online in Wiley InterScience (www.interscience.wiley.com). DOI:./nme A Petrov Galerkin formulation
More informationPulse propagation in random waveguides with turning points and application to imaging
Pulse propagation in random waveguides with turning points and application to imaging Liliana Borcea Mathematics, University of Michigan, Ann Arbor and Josselin Garnier Centre de Mathématiques Appliquées,
More information1 Upwind scheme for advection equation with variable. 2 Modified equations: numerical dissipation and dispersion
1 Upwind sceme for advection equation wit variable coefficient Consider te equation u t + a(x)u x Applying te upwind sceme, we ave u n 1 = a (un u n 1), a 0 u n 1 = a (un +1 u n ) a < 0. CFL condition
More informationOPTIMIZED DIFFERENCE SCHEMES FOR MULTIDIMENSIONAL HYPERBOLIC PARTIAL DIFFERENTIAL EQUATIONS
Seventh Mississippi State - UAB Conference on Differential Equations and Computational Simulations, Electronic Journal of Differential Equations, Conf. 17 (009, pp. 13 5. ISSN: 107-6691. URL: http://ejde.math.txstate.edu
More informationAn explicit time-domain finite-element method for room acoustics simulation
An explicit time-domain finite-element method for room acoustics simulation Takeshi OKUZONO 1 ; Toru OTSURU 2 ; Kimihiro SAKAGAMI 3 1 Kobe University, JAPAN 2 Oita University, JAPAN 3 Kobe University,
More informationFourier transforms, Generalised functions and Greens functions
Fourier transforms, Generalised functions and Greens functions T. Johnson 2015-01-23 Electromagnetic Processes In Dispersive Media, Lecture 2 - T. Johnson 1 Motivation A big part of this course concerns
More informationNumerical methods Revised March 2001
Revised March 00 By R. W. Riddaway (revised by M. Hortal) Table of contents. Some introductory ideas. Introduction. Classification of PDE's.3 Existence and uniqueness.4 Discretization.5 Convergence, consistency
More informationThe two-dimensional streamline upwind scheme for the convection reaction equation
INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS Int. J. Numer. Meth. Fluids 2001; 35: 575 591 The two-dimensional streamline upwind scheme for the convection reaction equation Tony W. H. Sheu*,1
More informationUniversität des Saarlandes. Fachrichtung 6.1 Mathematik
Universität des Saarlandes U N I V E R S I T A S S A R A V I E N I S S Fachrichtung 6.1 Mathematik Preprint Nr. 219 On finite element methods for 3D time dependent convection diffusion reaction equations
More informationBasics on Numerical Methods for Hyperbolic Equations
Basics on Numerical Methods for Hyperbolic Equations Professor Dr. E F Toro Laboratory of Applied Mathematics University of Trento, Italy eleuterio.toro@unitn.it http://www.ing.unitn.it/toro October 8,
More informationCourant and all that. Consistency, Convergence Stability Numerical Dispersion Computational grids and numerical anisotropy
Consistency, Convergence Stability Numerical Dispersion Computational grids and numerical anisotropy c dx ε 1 he goal of this lecture is to understand how to find suitable parameters (i.e., grid spacing
More informationAProofoftheStabilityoftheSpectral Difference Method For All Orders of Accuracy
AProofoftheStabilityoftheSpectral Difference Method For All Orders of Accuracy Antony Jameson 1 1 Thomas V. Jones Professor of Engineering Department of Aeronautics and Astronautics Stanford University
More informationReview of Fundamental Equations Supplementary notes on Section 1.2 and 1.3
Review of Fundamental Equations Supplementary notes on Section. and.3 Introduction of the velocity potential: irrotational motion: ω = u = identity in the vector analysis: ϕ u = ϕ Basic conservation principles:
More informationLecture 17: Initial value problems
Lecture 17: Initial value problems Let s start with initial value problems, and consider numerical solution to the simplest PDE we can think of u/ t + c u/ x = 0 (with u a scalar) for which the solution
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 informationDissipation and Dispersion
Consider the problem with periodic boundary conditions Dissipation and Dispersion u t = au x 0 < x < 1, t > 0 u 0 = sin 40 πx u(0, t) = u(1, t) t > 0 If a > 0 then the wave is moving to the left and if
More information13.1 Ion Acoustic Soliton and Shock Wave
13 Nonlinear Waves In linear theory, the wave amplitude is assumed to be sufficiently small to ignore contributions of terms of second order and higher (ie, nonlinear terms) in wave amplitude In such a
More information1 Introduction to MATLAB
L3 - December 015 Solving PDEs numerically (Reports due Thursday Dec 3rd, carolinemuller13@gmail.com) In this project, we will see various methods for solving Partial Differential Equations (PDEs) using
More informationWeight-adjusted DG methods for elastic wave propagation in arbitrary heterogeneous media
Weight-adjusted DG methods for elastic wave propagation in arbitrary heterogeneous media Jesse Chan Department of Computational and Applied Math, Rice University ICOSAHOM 2018 July 12, 2018 Chan (CAAM)
More informationarxiv: v1 [math.na] 21 Nov 2017
High Order Finite Difference Schemes for the Heat Equation Whose Convergence Rates are Higher Than Their Truncation Errors, A. Ditkowski arxiv:7.0796v [math.na] Nov 07 Abstract Typically when a semi-discrete
More informationA Multi-Dimensional Limiter for Hybrid Grid
APCOM & ISCM 11-14 th December, 2013, Singapore A Multi-Dimensional Limiter for Hybrid Grid * H. W. Zheng ¹ 1 State Key Laboratory of High Temperature Gas Dynamics, Institute of Mechanics, Chinese Academy
More informationIntroduction to Partial Differential Equations
Introduction to Partial Differential Equations Partial differential equations arise in a number of physical problems, such as fluid flow, heat transfer, solid mechanics and biological processes. These
More informationFirst, Second, and Third Order Finite-Volume Schemes for Diffusion
First, Second, and Third Order Finite-Volume Schemes for Diffusion Hiro Nishikawa 51st AIAA Aerospace Sciences Meeting, January 10, 2013 Supported by ARO (PM: Dr. Frederick Ferguson), NASA, Software Cradle.
More informationRecent Developments in the Flux Reconstruction Method
Recent Developments in the Flux Reconstruction Method Joshua Romero, PhD candidate Kartikey Asthana, PhD candidate Jonathan Bull, Postdoctoral Scholar Antony Jameson, Professor Aerospace Computing Laboratory,
More informationPattern Formation and Chaos
Developments in Experimental Pattern Formation - Isaac Newton Institute, 2005 1 Pattern Formation and Chaos Insights from Large Scale Numerical Simulations of Rayleigh-Bénard Convection Collaborators:
More informationEINDHOVEN UNIVERSITY OF TECHNOLOGY Department of Mathematics and Computer Science. CASA-Report March2008
EINDHOVEN UNIVERSITY OF TECHNOLOGY Department of Mathematics and Computer Science CASA-Report 08-08 March2008 The complexe flux scheme for spherically symmetrie conservation laws by J.H.M. ten Thije Boonkkamp,
More informationRiemann Solvers and Numerical Methods for Fluid Dynamics
Eleuterio R Toro Riemann Solvers and Numerical Methods for Fluid Dynamics A Practical Introduction With 223 Figures Springer Table of Contents Preface V 1. The Equations of Fluid Dynamics 1 1.1 The Euler
More information