Well-balanced and Asymptotic Schemes for Friedrichs systems
|
|
- Ross Fisher
- 6 years ago
- Views:
Transcription
1 B. Després (LJLL- University Paris VI) C. Buet (CEA) Well-balanced and Asymptotic Schemes for Friedrichs systems thanks to T. Leroy (PhD CEA), X. Valentin (PhD CEA), C. Enaux (CEA), P. Lafitte (Centrale Paris) and E. Franck (Inria) B. Després (LJLL-University Paris VI) C. Buet (CEA) thanks to T. Leroy (PhD CEA), X. Valentin (PhD CEA), C. Enaux (CEA), P. Lafitte (Centrale Paris) and E. Franck (Inria) Support of CEA is acknowledged p. 1 / 9
2 Section 1 Numhyp 015 p. / 9
3 Discretize tu + f (U) = S(U,... ) with FV-P 0 -cell centered schemes U n+1 j Tentative definition : U n j t + f n j+ 1 f n j 1 x = S(U n j,... ) WB-FV scheme = exact for stationary solutions AP (asymptotic preserving) schemes = WB with stiffness. Literature on WB : it is a one case after another strategy. Can we have a global view (classification) on WB methods? Numhyp 015 p. / 9
4 Some references Cargo-LeRoux, A WB for a model of an atmosphere with gravity CRAS I 318 (1994) Greenberg-Leroux, A WB scheme for the numerical processing of source terms in hyp. eq., Audusse-Bouchut-Bristeau-Klein-Perthame, A fast and stable well-balanced scheme with hydrostatic reconstruction for shallow water flows, (004). Bouchut. Nonlinear stability of finite volume methods for hyperbolic conservation laws, and well-balanced schemes for sources, Frontiers in Mathematics series, 004. Coquel-Godlewski. Asymptotic preserving scheme for Euler system with large friction, 011. T. Muller and A. Pfeiffer, Well-balanced simulation of geophysical flows via the shallow water equations with bottom topography : consistency and numerical computations, 014. Larsen-Morel-Miller, Asymptotic solution of num. transport solutions in optical thick regimes L. Gosse and G. Toscani, An AP WB scheme for the hyp. heat equations. CRAS I, 00. S. Jin, AP schemes for multiscale kinetic and hyperbolic equations : a review, 010. Buet-Després-Franck Design of AP schemes for the hyp. heat equation on unstructured meshes, 01. Buet-Després-Franck, AP Schemes on Distorted Meshes for Friedrichs Systems with Stiff Relaxation : Application to Angular Models in Linear Transport, 014. L. Gosse, Computing Qualitatively Correct Approximations of Balance Laws..., 013. Numhyp 015 p. 3 / 9
5 Model problem : Friedrichs systems with relaxation A paradigm is the p-system with friction : linearization yields { tτ xu = 0, tu + xp(τ) = g σu. = { tτ 1 xu 1 = 0, tu 1 x ( c0(x) τ 1 ) = σu1. One gets the linear hyperbolic heat equation tu + x(a(x)u) = R(x)U, U = (c 0(x)τ 1, u 1) t, ( ) ( ) c where A = c 0(x) = A t and R = 0 (x) σ Notice the compatibility relation R + R t + A (x) 0 U t + (U, A(x)U) x = 1 ( R(x) + R(x) t + A (x)u, U ) 0. Numhyp 015 p. 4 / 9
6 S n and DOM Another example is DOM (Discrete Ordinates Method) for transfer Approximate I (x, t, µ) = ti + µ xi = σ (< I > I ). n n w i f i (x, t)δ(µ µ i ) + w i g i (x, t)δ(µ + µ i ). i=1 i=1 Normalize U = ( diag(w i ) f, diag(w i ) g) t R n, with w = ( w 1,..., w n, w 1,..., w n) R n. One gets the S n model tu + A xu = RU, A = A t, R = w w + I d, Numhyp 015 p. 5 / 9
7 Main idea Generic model problem is tu + x(a(x)u) + y (B(x)U) = R(x)U, U(t, x) R n, where x = (x, y) R, A(x) = A(x) t, B(x) = B(x) t. The vectorial space of stationary states, needed for WB analysis, is : U = {x U(x); x(a(x)u) + y (B(x)U) = R(x)U}. Main idea - consider the adjoint equation : tv + A(x) t xv + B(x) t y V = R(x) t V, V (t, x) R n. Magic property - this is just duality : t(u, V ) + x(a(x)u, V ) + y (B(x)U, V ) = 0. Numhyp 015 p. 6 / 9
8 Strategy The vectorial space of dual stationary states is V = { x V (x); A(x) t xv + B(x) t y V = R(x) t U }. Pick up V p V and define α p = (U, V p) R which is solution of a conservative equation tα p + x(a(x)u, V p) + y (B(x)U, V p) = 0 a) In 1D, assemble the system of conservation laws for α = (α p) R dim(v) = R dim(u). b) Discretize the new system of conservation laws. c) Write the new scheme for the original variable U. Numhyp 015 p. 7 / 9
9 For simplicity : A is constant with det(a) 0 So dim(v) = dim(u) = n. A t xv = R t V xv = A t R t V V (x) = e A t R t x V (0) where e A t R t x is a matrix exponential. Proposition : Define the change-of-basis matrix P(x) = e RA 1x, and the change of unknown α = P(x)U U = P 1 (x)α. The new conservative system rewrites tα + x(q(x)α) = 0 with a matrix defined through a global change of basis Q(x) = P(x)AP 1 (x) Note Q(x) is similar to A, so has the same eigenvalues. But the eigenvectors depend on x. Numhyp 015 p. 8 / 9
10 Discretization For additional simplification (once again), A is also symmetric : A = A t R n n. One has the spectral decomposition Au p = λ pu p, λ p 0 ( ) u p, U j+ 1 U j = 0, λ p > 0, ( ) u p, U j+ 1 U j+1 = 0, λ p < 0. Intermediate state L R The intermediate state is computed in function of L = j and R = j + 1. The FV discretization of tα + xβ = 0 with β = Q(x)α writes where β n j+ 1 α n+1 j α n j t + β n j+ 1 is the flux at time step t n = n t. β n j 1 x j = 0 Numhyp 015 p. 9 / 9
11 Diagonalisation of Q The right and left spectral decompositions of Q = Q(x ) = P(x )AP(x ) 1 Q are { Q(x )r p = λ pr p, r p = P(x )u p, Q t (x )s p = λ ps p, s p = P(x ) t u p. The solution at interface of the system t ( s p, β ) + λ p x ( s p, β ) = 0 yields a family of. One-state : defined by the well-posed linear system { ( s p, β β L ) = 0, λp > 0, ( s p, β β R ) = 0, λp < 0. It defines a first solver ϕ : R n R n R R n ϕ(β L, β R, x ) = β. Numhyp 015 p. 10 / 9
12 The WB scheme with the one-state solver First formulation with the global change of basis α n+1 α n j j t + ϕ(β n j,βn j+1,x j+ 1 ) ϕ(βj 1 n,βn j,x j 1 ) x j = 0 n ϕ(β 1 Second formulation Un+1 U n j j j + P(x t j ),βn j+1,x j+ 1 ) ϕ(βj 1 n,βn j,x j 1 ) x j = 0 Third purely local formulation, the one for implementation, U n+1 j U n j t U j+ + A 1 U j 1 x j + P(x j) 1 P(x j+ 1 ) I AU j+ x 1 + I P(x j) 1 P(x j 1 ) AU j j x 1 = 0. j where the flux U j+ 1 is solution of the linear system (with x ± j = x j± 1 x j ) ( ) u p, U j+ 1 e A 1 R x + j U j = 0, λ p > 0, ( ) u p, U j+ 1 e A 1 R x j+1u j+1 = 0, λ p < 0. Property : the scheme is WB. Proof : If β j = β for all j, the solution is stationary in the first formulation. Numhyp 015 p. 11 / 9
13 Schematic Structure of the one-state solver. Intermediate state N L L Two-states solver. Same principle but with upwinding of the eigenvectors { ( s L p, β β L ) = 0, λp > 0, ( s R p, β β R ) = 0, λp < 0, N R R It defines ψ : R n R n R R R n : ψ(β L, β R, x L, x R ) = β. Intermediate state N L N R L R Th. : If R + R t 0, then { } sp L { } s R λ p p>0 λ p<0 is linearly independent. Numhyp 015 p. 1 / 9
14 A 1 R = ( 0 σ 0 0 Interpretation for : t p + x u = 0, t u + x p = σu ) is nilpotent. So e A 1Rx = I + A 1 Rx. One-state solver = Steady state approximation of Huang-Liu 1986, or Jin-Levermore 1996 p j+ 1 u j+ 1 = p j + p j+1 = p j p j+1 + (1 σ x + j )u j (1 + σ x j+1 )u j+1, + (1 σ x + j )u j + (1 + σ x j+1 )u j+1. Two-state solver= Equal to the Gosse-Toscani scheme 00. p 1 σ x 1+σ x j+ 1 j+1 = (p +σ x + σ x j u j ) + + j (p +σ x + σ x j+1 u j+1 ), j j+1 j j+1 u 1 j+ 1 = (u j + u j+1 + p j p j+1 ). +σ x + σ x j j+1 - Immediate interest for non linear fluid with σ the friction parameter. Numhyp 015 p. 13 / 9
15 Section 3 Numhyp 015 p. 14 / 9
16 #1 : P n method, Valentin (PhD-CEA)+Enaux+Laffitte This is a challenging problem with a vast literature in neutron propagation. Machorro 007 : consider µ r u + 1 µ u = σu + q with r u(r, µ) = n i=0 u i(r)p i (µ) in the domain Solve tu + A r U + 1 GU = σu + b. r Valentin-Enaux-Lafitte propose to use Magnus series to compute the local matrix exponentials. Numhyp 015 p. 14 / 9
17 P n results : Machorro 007 Discontinuity of coefficients. flux dip, oscillations, higher cost with DG. Numhyp 015 p. 15 / 9
18 P n results : Valentin 015 No flux dip, no oscillations, reduced cost. Numhyp 015 p. 16 / 9
19 # : Transport-Doppler effect, Leroy (PhD-CEA)+Buet+D. Solve : tρ = κ ν νρ + σ(ν)(b(ν) ρ), v, t > 0 with κ = xu R. 3 The SWB (Spectrally Well Balanced) scheme writes where tρ j = ρ(ν ν, ρ ) = ρ e κ 3 ν ν σ(ν j ) ρ(ν j ν j+1, ρ j+1 ) ρ j 1 M(ν j+1, ν j ) ν j σ(s) s ds + 3 ν B(s)σ(s) s κ ν s τ is the stationary solution at ν j and the additional term M(ν j+1, ν) = e 3 κ ν j+1 ν j σ(s) s ds gives the correct limit κ 0. σ(τ) dτds τ Numhyp 015 p. 17 / 9
20 Upwing/WB/SWB WB and SWB correct WB non correct for κ = 0 Numhyp 015 p. 18 / 9
21 #3 : System with a zero eigenvalue This is a serious issue, related to resonance in hyperbolic theory. Consider P 1 model coupled a linear temperature equation tp + xu = τ(t p), tu + xp = σu, tt = τ(p T ), A = A t = and R = R t = τ 0 τ 0 σ 0 τ 0 τ Assume σ, τ > 0 : the solutions of the adjoint stationary equation satisfy T = p, so xû = 0 and x p = σû. Therefore. dim(v) = dim(u) = < 3 with basis 1 σx V 1 = 0 1 and V = 0 σx. Numhyp 015 p. 19 / 9
22 One abstract solution= shift the spectrum Consider progressive solutions of the adjoint equation { tv + A xv = RV, tv + ξ xv = 0, ξ R, where ξ R is arbitrary. It yields the shifted problem A ξ xv = RV, A ξ = A ξi. So det(a ξ ) 0. Probably non correct in the regime ξ 0, so useless. Numhyp 015 p. 0 / 9
23 Ad-hoc integration in time Set α 1 = (U, V 1) = p + T, α = (U, V ) = σx(p + T ) + u and α 3 = T tα 1 + x( σxα ( 1 + α ) = 0, ) tα + x (1 σ x )α 1 + σxα α 3 = 0, tα 3 = τ(α 1 α 3). The last equation is more an ODE (e τt α 3) = τe τt α 1 t α 3(x, t) = e τt τe τs α 1(x, s)ds + e τt α 3(x, 0). 0 Explicit Euler approximation over t reads α 3(x, t) 1 ( 1 e τ t) α 1(x, 0) + e τ t α 3(x, 0). It yields { tα 1 + x( σxα 1 + α ) = 0, tα + x ( (1 σ x )α 1 + σxα 1 ( 1 e τ t ) α 1 e τ t α 3(0) ) = 0. Numhyp 015 p. 1 / 9
24 Section 4 Numhyp 015 p. / 9
25 #4 : D Assume σ > 0 is constant tp + xu + y v = 0, tu + xp = σu, tv + y p = σv, Here det(a) = 0, det(b) = 0 and AB BA. The adjoint stationary states ( p, û, v) are xû + y v = 0, x p = σû, y p = σ v. So (û, v) = 1 p and p = 0. Therefore dim(v) = dim(u) = σ V = {( p, û, v); p is harmonic and (û, v) = 1σ } p. Numhyp 015 p. / 9
26 New formulation Set V n = V { p is an harmonic polynomial of degree n}. p V 1 is equivalent to p = a + bx + cy. V 1 yields three test functions 1 V 1 = 0 0, V = σx 1 0, V 3 = Set α i = (U, V i ) for i = 1,, 3. The new system is α 1 m 1 m t α + x α 1 + σxm 1 + y σxm α 3 σym σym σy 0 1 = 0 with m 1 = σxα 1 + α and m = σyα 1 + α 3. The original variable is recovered as U = P(x) 1 α where P(x) = σx 1 0 is non singular, and P(x)P(y) = P(y)P(x). σy 0 1 Numhyp 015 p. 3 / 9
27 Convergence Performed on the problem with an additional small parameter { tp + 1 ε u = 0, tu + 1 ε p = σ ε u, where σ > 0 is given and ε > 0 is a small parameter. In the limit ε 0 +, one gets asymptotically the parabolic equation tp 1 p = 0. σ Theorem (D on general grids) : The implicit edge-based-fv WB scheme (= Gosse-Toscani if 1D on uniform mesh) is convergent p ε h, t p ε, u ε h, t u ε L ([0,T ] Ω) C(h t 1 ), uniformly with respect to small parameter ε (0, 1] (numerics show convergence order between 1 and ). Sketch of the proof : linear estimates but of hyperbolic-parabolic type, balance between dissipativity and source terms, general grid, control of the constants,.... Numhyp 015 p. 4 / 9
28 Non AP feature of the standard edge based VF scheme element size L1 Norm L Norm h^ Error for the VF4 scheme Numhyp 015 p. 5 / 9
29 Meshes and P 1 (E. Franck) Initial data is Dirac mass at center Kershaw mesh Random mesh non WB-AP WB-FV edge, non AP WB-AP corner Mesh imprint vastly reduced with WB-AP corner-based-flux FV scheme. Numhyp 015 p. 6 / 9
30 Convergence order Analytical solution in the torus : p = f + ε tf and u = ε f with σ σ tf + ε σ t f 1 f = 0, σ f (t, x, y) = α(t) cos(lπx)) cos(lπy). α (t) + ε σ α (t) + L π σ α(t) = 0, α(t) = λ λ λ 1 e λ 1t λ 1 λ λ 1 e λ t : 1: 0.5: 0.5: 0: order one L Error e N In tests : ε = 0.01(40h) τ for τ {0, 1, 1, 1, }. 4 1 experimental convergence. Numhyp 015 p. 7 / 9
31 #5 : Lagrangian fluid dynamics, comparison standard FV versus WB-FV Initial stage of Rayleigh-Taylor instability : source term = gravity. D P1 standard FV spurious velocities WB-FV vortex is correct Numhyp 015 p. 8 / 9
32 Conclusions Unification of WB based on the dual equation essentially ideas : duality and linear, essentially families of WB. General principle : any WB-FV is a standard FV. AMC 015, Buet-D : and references therein. Validated for Friedrichs systems of large size inspired by radiation/neutron inspired problems. Any solver can be used : high order, even FEM is possible, or central schemes,... The change-of-basis matrix P(x) seems the good object to manipulate, particularly in multid. Ongoing : convergence estimates, zero eigenvalue problem. Open problems : - Comparison with other methods for transport equations. - Convergence for AP with this method : e A 1 R x ε - Extension to the full non linear case. Numhyp 015 p. 9 / 9
Model adaptation in hierarchies of hyperbolic systems
Model adaptation in hierarchies of hyperbolic systems Nicolas Seguin Laboratoire J.-L. Lions, UPMC Paris 6, France February 15th, 2012 DFG-CNRS Workshop Nicolas Seguin (LJLL, UPMC) 1 / 29 Outline of the
More informationAngular momentum preserving CFD on general grids
B. Després LJLL-Paris VI Thanks CEA and ANR Chrome Angular momentum preserving CFD on general grids collaboration Emmanuel Labourasse (CEA) B. Després LJLL-Paris VI Thanks CEA and ANR Chrome collaboration
More informationAll-regime Lagrangian-Remap numerical schemes for the gas dynamics equations. Applications to the large friction and low Mach regimes
All-regime Lagrangian-Remap numerical schemes for the gas dynamics equations. Applications to the large friction and low Mach regimes Christophe Chalons LMV, Université de Versailles Saint-Quentin-en-Yvelines
More informationAll-regime Lagrangian-Remap numerical schemes for the gas dynamics equations. Applications to the large friction and low Mach coefficients
All-regime Lagrangian-Remap numerical schemes for the gas dynamics equations. Applications to the large friction and low Mach coefficients Christophe Chalons LMV, Université de Versailles Saint-Quentin-en-Yvelines
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 informationAdmissibility and asymptotic-preserving scheme
Admissibility and asymptotic-preserving scheme F. Blachère 1, R. Turpault 2 1 Laboratoire de Mathématiques Jean Leray (LMJL), Université de Nantes, 2 Institut de Mathématiques de Bordeaux (IMB), Bordeaux-INP
More informationAn asymptotic preserving unified gas kinetic scheme for the grey radiative transfer equations
An asymptotic preserving unified gas kinetic scheme for the grey radiative transfer equations Institute of Applied Physics and Computational Mathematics, Beijing NUS, Singapore, March 2-6, 2015 (joint
More informationUn schéma volumes finis well-balanced pour un modèle hyperbolique de chimiotactisme
Un schéma volumes finis well-balanced pour un modèle hyperbolique de chimiotactisme Christophe Berthon, Anaïs Crestetto et Françoise Foucher LMJL, Université de Nantes ANR GEONUM Séminaire de l équipe
More informationFDM for wave equations
FDM for wave equations Consider the second order wave equation Some properties Existence & Uniqueness Wave speed finite!!! Dependence region Analytical solution in 1D Finite difference discretization Finite
More informationOn a simple model of isothermal phase transition
On a simple model of isothermal phase transition Nicolas Seguin Laboratoire Jacques-Louis Lions Université Pierre et Marie Curie Paris 6 France Micro-Macro Modelling and Simulation of Liquid-Vapour Flows
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 informationPDEs, part 1: Introduction and elliptic PDEs
PDEs, part 1: Introduction and elliptic PDEs Anna-Karin Tornberg Mathematical Models, Analysis and Simulation Fall semester, 2013 Partial di erential equations The solution depends on several variables,
More informationAlexei F. Cheviakov. University of Saskatchewan, Saskatoon, Canada. INPL seminar June 09, 2011
Direct Method of Construction of Conservation Laws for Nonlinear Differential Equations, its Relation with Noether s Theorem, Applications, and Symbolic Software Alexei F. Cheviakov University of Saskatchewan,
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 informationLarge time-step and asymptotic-preserving numerical schemes for gaz dynamics with stiff sources
Lagrange-Projection strategy for the gas-dynamics equations Relaxation for the Lagrangian step Source terms and notion of consistency in the integral sense Propertie Large time-step and asymptotic-preserving
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 informationWell-Balanced Schemes for the Euler Equations with Gravity
Well-Balanced Schemes for the Euler Equations with Gravity Alina Chertock North Carolina State University chertock@math.ncsu.edu joint work with S. Cui, A. Kurganov, S.N. Özcan and E. Tadmor supported
More informationGeneralised Summation-by-Parts Operators and Variable Coefficients
Institute Computational Mathematics Generalised Summation-by-Parts Operators and Variable Coefficients arxiv:1705.10541v [math.na] 16 Feb 018 Hendrik Ranocha 14th November 017 High-order methods for conservation
More informationModeling and simulation of bedload transport with viscous effects
Introduction Modeling and simulation of bedload transport with viscous effects E. Audusse, L. Boittin, M. Parisot, J. Sainte-Marie Project-team ANGE, Inria; CEREMA; LJLL, UPMC Université Paris VI; UMR
More informationSolving PDEs with freefem++
Solving PDEs with freefem++ Tutorials at Basque Center BCA Olivier Pironneau 1 with Frederic Hecht, LJLL-University of Paris VI 1 March 13, 2011 Do not forget That everything about freefem++ is at www.freefem.org
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 informationMath 660-Lecture 23: Gudonov s method and some theories for FVM schemes
Math 660-Lecture 3: Gudonov s method and some theories for FVM schemes 1 The idea of FVM (You can refer to Chapter 4 in the book Finite volume methods for hyperbolic problems ) Consider the box [x 1/,
More informationModelling and numerical methods for the diffusion of impurities in a gas
INERNAIONAL JOURNAL FOR NUMERICAL MEHODS IN FLUIDS Int. J. Numer. Meth. Fluids 6; : 6 [Version: /9/8 v.] Modelling and numerical methods for the diffusion of impurities in a gas E. Ferrari, L. Pareschi
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 informationModeling uncertainties with kinetic equations
B. Després (LJLL- University Paris VI) thanks to G. Poette (CEA) and D. Lucor (Orsay) B. Perthame (LJLL) E. Trélat (LJLL) Modeling uncertainties with kinetic equations B. Després (LJLL-University Paris
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 informationPalindromic Discontinuous Galerkin Method
Palindromic Discontinuous Galerkin Method David Coulette, Emmanuel Franck, Philippe Helluy, Michel Mehrenberger, Laurent Navoret To cite this version: David Coulette, Emmanuel Franck, Philippe Helluy,
More informationLucio Demeio Dipartimento di Ingegneria Industriale e delle Scienze Matematiche
Scuola di Dottorato THE WAVE EQUATION Lucio Demeio Dipartimento di Ingegneria Industriale e delle Scienze Matematiche Lucio Demeio - DIISM wave equation 1 / 44 1 The Vibrating String Equation 2 Second
More informationFinite volumes for complex applications In this paper, we study finite-volume methods for balance laws. In particular, we focus on Godunov-type centra
Semi-discrete central schemes for balance laws. Application to the Broadwell model. Alexander Kurganov * *Department of Mathematics, Tulane University, 683 St. Charles Ave., New Orleans, LA 708, USA kurganov@math.tulane.edu
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 Dept. of
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 informationNumerical Oscillations and how to avoid them
Numerical Oscillations and how to avoid them Willem Hundsdorfer Talk for CWI Scientific Meeting, based on work with Anna Mozartova (CWI, RBS) & Marc Spijker (Leiden Univ.) For details: see thesis of A.
More informationConservation Laws: Systematic Construction, Noether s Theorem, Applications, and Symbolic Computations.
Conservation Laws: Systematic Construction, Noether s Theorem, Applications, and Symbolic Computations. Alexey Shevyakov Department of Mathematics and Statistics, University of Saskatchewan, Saskatoon,
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 Survey of Computational High Frequency Wave Propagation II. Olof Runborg NADA, KTH
A Survey of Computational High Frequency Wave Propagation II Olof Runborg NADA, KTH High Frequency Wave Propagation CSCAMM, September 19-22, 2005 Numerical methods Direct methods Wave equation (time domain)
More informationHyperbolic Systems of Conservation Laws. in One Space Dimension. II - Solutions to the Cauchy problem. Alberto Bressan
Hyperbolic Systems of Conservation Laws in One Space Dimension II - Solutions to the Cauchy problem Alberto Bressan Department of Mathematics, Penn State University http://www.math.psu.edu/bressan/ 1 Global
More informationNon-Conforming Finite Element Methods for Nonmatching Grids in Three Dimensions
Non-Conforming Finite Element Methods for Nonmatching Grids in Three Dimensions Wayne McGee and Padmanabhan Seshaiyer Texas Tech University, Mathematics and Statistics (padhu@math.ttu.edu) Summary. In
More informationFinite Element methods for hyperbolic systems
Finite Element methods for hyperbolic systems Eric Sonnendrücker Max-Planck-Institut für Plasmaphysik und Zentrum Mathematik, TU München Lecture notes Wintersemester 14-15 January 19, 15 Contents 1 Introduction
More informationProjection Dynamics in Godunov-Type Schemes
JOURNAL OF COMPUTATIONAL PHYSICS 142, 412 427 (1998) ARTICLE NO. CP985923 Projection Dynamics in Godunov-Type Schemes Kun Xu and Jishan Hu Department of Mathematics, Hong Kong University of Science and
More informationUniformly accurate averaging numerical schemes for oscillatory evolution equations
Uniformly accurate averaging numerical schemes for oscillatory evolution equations Philippe Chartier University of Rennes, INRIA Joint work with M. Lemou (University of Rennes-CNRS), F. Méhats (University
More informationA kinetic scheme for transient mixed flows in non uniform closed pipes: a global manner to upwind all the source terms
A kinetic scheme for transient mixed flows in non uniform closed pipes: a global manner to upwind all the source terms C. BOURDARIAS M. ERSOY S. GERBI LAMA-Université de Savoie de Chambéry, France September
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 informationLECTURE 3: DISCRETE GRADIENT FLOWS
LECTURE 3: DISCRETE GRADIENT FLOWS Department of Mathematics and Institute for Physical Science and Technology University of Maryland, USA Tutorial: Numerical Methods for FBPs Free Boundary Problems and
More informationOptimal Prediction for Radiative Transfer: A New Perspective on Moment Closure
Optimal Prediction for Radiative Transfer: A New Perspective on Moment Closure Benjamin Seibold MIT Applied Mathematics Mar 02 nd, 2009 Collaborator Martin Frank (TU Kaiserslautern) Partial Support NSF
More informationarxiv: v2 [math.ap] 1 Jul 2011
A Godunov-type method for the shallow water equations with discontinuous topography in the resonant regime arxiv:1105.3074v2 [math.ap] 1 Jul 2011 Abstract Philippe G. efloch 1 and Mai Duc Thanh 2 1 aboratoire
More informationLecture 12: Detailed balance and Eigenfunction methods
Miranda Holmes-Cerfon Applied Stochastic Analysis, Spring 2015 Lecture 12: Detailed balance and Eigenfunction methods Readings Recommended: Pavliotis [2014] 4.5-4.7 (eigenfunction methods and reversibility),
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 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 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 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 informationDedicated to the 70th birthday of Professor Lin Qun
Journal of Computational Mathematics, Vol.4, No.3, 6, 39 5. ANTI-DIFFUSIVE FINITE DIFFERENCE WENO METHODS FOR SHALLOW WATER WITH TRANSPORT OF POLLUTANT ) Zhengfu Xu (Department of Mathematics, Pennsylvania
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 informationUnderstand the existence and uniqueness theorems and what they tell you about solutions to initial value problems.
Review Outline To review for the final, look over the following outline and look at problems from the book and on the old exam s and exam reviews to find problems about each of the following topics.. Basics
More informationNumerical Modeling of Stratified Shallow Flows: Applications to Aquatic Ecosystems
Numerical Modeling of Stratified Shallow Flows: Applications to Aquatic Ecosystems Marica Pelanti 1,2, Marie-Odile Bristeau 1 and Jacques Sainte-Marie 1,2 1 INRIA Paris-Rocquencourt, 2 EDF R&D Joint work
More informationFinite-Volume-Particle Methods for Models of Transport of Pollutant in Shallow Water
Journal of Scientific Computing ( 006) DOI: 0.007/s095-005-9060-x Finite-Volume-Particle Methods for Models of Transport of Pollutant in Shallow Water Alina Chertock, Alexander Kurganov, and Guergana Petrova
More informationLECTURE # 0 BASIC NOTATIONS AND CONCEPTS IN THE THEORY OF PARTIAL DIFFERENTIAL EQUATIONS (PDES)
LECTURE # 0 BASIC NOTATIONS AND CONCEPTS IN THE THEORY OF PARTIAL DIFFERENTIAL EQUATIONS (PDES) RAYTCHO LAZAROV 1 Notations and Basic Functional Spaces Scalar function in R d, d 1 will be denoted by u,
More informationA Well-Balanced Path-Integral f-wave Method for Hyperbolic Problems with Source Terms
A Well-Balanced Path-Integral f-wave Method for Hyperbolic Problems with Source Terms Randall J. LeVeque 1 Abstract. Systems of hyperbolic partial differential equations with source terms (balance laws)
More informationLarge time step and asymptotic preserving numerical schemes for the gas dynamics equations with source terms
Large time step and asymptotic preserving numerical schemes for the gas dynamics equations with source terms Christophe Chalons, Samuel Kokh, Mathieu Girardin To cite this version: Christophe Chalons,
More informationA Godunov-type method for the diphasic Baer-Nunziato model
A Godunov-type method for the diphasic Baer-Nunziato model C. Chalons Joint work with : A. Ambroso, P.-A. Raviart with benefit from discussions with F. Coquel (CNRS, Ecole Polytechnique/CMAP) N. Seguin
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 informationOn the Asymptotic Preserving property of the Unified Gas Kinetic Scheme for the diffusion limit of linear kinetic models
On the Asymptotic Preserving property of the Unified Gas Kinetic Scheme for the diffusion limit of linear kinetic models Luc Mieussens To cite this version: Luc Mieussens. On the Asymptotic Preserving
More informationHigh Order Accurate Runge Kutta Nodal Discontinuous Galerkin Method for Numerical Solution of Linear Convection Equation
High Order Accurate Runge Kutta Nodal Discontinuous Galerkin Method for Numerical Solution of Linear Convection Equation Faheem Ahmed, Fareed Ahmed, Yongheng Guo, Yong Yang Abstract This paper deals with
More informationENO and WENO schemes. Further topics and time Integration
ENO and WENO schemes. Further topics and time Integration Tefa Kaisara CASA Seminar 29 November, 2006 Outline 1 Short review ENO/WENO 2 Further topics Subcell resolution Other building blocks 3 Time Integration
More informationLinear Hyperbolic Systems
Linear Hyperbolic Systems 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, 2014 1 / 56 We study some basic
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 informationTwo-layer shallow water system and its applications
Proceedings of Symposia in Applied Mathematics Two-layer shallow water system and its applications Jihwan Kim and Randall J. LeVeque Abstract. The multi-layer shallow water system is derived by depth averaging
More informationMath 7824 Spring 2010 Numerical solution of partial differential equations Classroom notes and homework
Math 7824 Spring 2010 Numerical solution of partial differential equations Classroom notes and homework Jan Mandel University of Colorado Denver May 12, 2010 1/20/09: Sec. 1.1, 1.2. Hw 1 due 1/27: problems
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 informationApplication of the Kurganov Levy semi-discrete numerical scheme to hyperbolic problems with nonlinear source terms
Future Generation Computer Systems () 65 7 Application of the Kurganov Levy semi-discrete numerical scheme to hyperbolic problems with nonlinear source terms R. Naidoo a,b, S. Baboolal b, a Department
More informationNonlinear stability of semidiscrete shocks for two-sided schemes
Nonlinear stability of semidiscrete shocks for two-sided schemes Margaret Beck Boston University Joint work with Hermen Jan Hupkes, Björn Sandstede, and Kevin Zumbrun Setting: semi-discrete conservation
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 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 informationSecond Order Elliptic PDE
Second Order Elliptic PDE T. Muthukumar tmk@iitk.ac.in December 16, 2014 Contents 1 A Quick Introduction to PDE 1 2 Classification of Second Order PDE 3 3 Linear Second Order Elliptic Operators 4 4 Periodic
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 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 information23 Elements of analytic ODE theory. Bessel s functions
23 Elements of analytic ODE theory. Bessel s functions Recall I am changing the variables) that we need to solve the so-called Bessel s equation 23. Elements of analytic ODE theory Let x 2 u + xu + x 2
More informationAn asymptotic-preserving micro-macro scheme for Vlasov-BGK-like equations in the diffusion scaling
An asymptotic-preserving micro-macro scheme for Vlasov-BGK-like equations in the diffusion scaling Anaïs Crestetto 1, Nicolas Crouseilles 2 and Mohammed Lemou 3 Saint-Malo 13 December 2016 1 Université
More information256 Summary. D n f(x j ) = f j+n f j n 2n x. j n=1. α m n = 2( 1) n (m!) 2 (m n)!(m + n)!. PPW = 2π k x 2 N + 1. i=0?d i,j. N/2} N + 1-dim.
56 Summary High order FD Finite-order finite differences: Points per Wavelength: Number of passes: D n f(x j ) = f j+n f j n n x df xj = m α m dx n D n f j j n= α m n = ( ) n (m!) (m n)!(m + n)!. PPW =
More informationA Hamiltonian Numerical Scheme for Large Scale Geophysical Fluid Systems
A Hamiltonian Numerical Scheme for Large Scale Geophysical Fluid Systems Bob Peeters Joint work with Onno Bokhove & Jason Frank TW, University of Twente, Enschede CWI, Amsterdam PhD-TW colloquium, 9th
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 informationThe RAMSES code and related techniques 4. Source terms
The RAMSES code and related techniques 4. Source terms Outline - Optically thin radiative hydrodynamics - Relaxation towards the diffusion limit - Hydrodynamics with gravity source term - Relaxation towards
More informationAnti-diffusive finite difference WENO methods for shallow water with. transport of pollutant
Anti-diffusive finite difference WENO methods for shallow water with transport of pollutant Zhengfu Xu 1 and Chi-Wang Shu 2 Dedicated to Professor Qun Lin on the occasion of his 70th birthday Abstract
More informationNumerical Simulation of Rarefied Gases using Hyperbolic Moment Equations in Partially-Conservative Form
Numerical Simulation of Rarefied Gases using Hyperbolic Moment Equations in Partially-Conservative Form Julian Koellermeier, Manuel Torrilhon May 18th, 2017 FU Berlin J. Koellermeier 1 / 52 Partially-Conservative
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 informationSpace-time Discontinuous Galerkin Methods for Compressible Flows
Space-time Discontinuous Galerkin Methods for Compressible Flows Jaap van der Vegt Numerical Analysis and Computational Mechanics Group Department of Applied Mathematics University of Twente Joint Work
More informationFullSWOF a software for overland flow simulation
FullSWOF a software for overland flow simulation O. Delestre, S. Cordier, F. Darboux, M. Du, F. James, C. Laguerre, C. Lucas & O. Planchon LJA Dieudonn e & Polytech Nice Sophia SIMHYDRO 2012 Sophia-Antipolis
More informationBasic Principles of Weak Galerkin Finite Element Methods for PDEs
Basic Principles of Weak Galerkin Finite Element Methods for PDEs Junping Wang Computational Mathematics Division of Mathematical Sciences National Science Foundation Arlington, VA 22230 Polytopal Element
More informationPeriodic oscillations in the Gross-Pitaevskii equation with a parabolic potential
Periodic oscillations in the Gross-Pitaevskii equation with a parabolic potential Dmitry Pelinovsky 1 and Panos Kevrekidis 2 1 Department of Mathematics, McMaster University, Hamilton, Ontario, Canada
More informationHierarchical Reconstruction with up to Second Degree Remainder for Solving Nonlinear Conservation Laws
Hierarchical Reconstruction with up to Second Degree Remainder for Solving Nonlinear Conservation Laws Dedicated to Todd F. Dupont on the occasion of his 65th birthday Yingjie Liu, Chi-Wang Shu and Zhiliang
More informationHIGH ORDER FINITE VOLUME SCHEMES BASED ON RECONSTRUCTION OF STATES FOR SOLVING HYPERBOLIC SYSTEMS WITH NONCONSERVATIVE PRODUCTS
MATHEMATICS OF COMPUTATION Volume 75, Number 255, July 26, Pages 113 1134 S 25-5718(6)1851-5 Article electronically published on March 21, 26 HIGH ORDER FINITE VOLUME SCHEMES BASED ON RECONSTRUCTION OF
More informationElectromagnetic Modeling and Simulation
Electromagnetic Modeling and Simulation Erin Bela and Erik Hortsch Department of Mathematics Research Experiences for Undergraduates April 7, 2011 Bela and Hortsch (OSU) EM REU 2010 1 / 45 Maxwell s Equations
More informationUne décomposition micro-macro particulaire pour des équations de type Boltzmann-BGK en régime de diffusion
Une décomposition micro-macro particulaire pour des équations de type Boltzmann-BGK en régime de diffusion Anaïs Crestetto 1, Nicolas Crouseilles 2 et Mohammed Lemou 3 La Tremblade, Congrès SMAI 2017 5
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 informationAn Accurate Deterministic Projection Method for Hyperbolic Systems with Stiff Source Term
An Accurate Deterministic Projection Method for Hyperbolic Systems with Stiff Source Term Alexander Kurganov Department of Mathematics, Tulane University, 683 Saint Charles Avenue, New Orleans, LA 78,
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 informationWELL-BALANCED SCHEMES TO CAPTURE NON-EXPLICIT STEADY STATES. PART 1: RIPA MODEL
WELL-BALANCED SCHEMES TO CAPTURE NON-EXPLICIT STEADY STATES. PART 1: RIPA MODEL VIVIEN DESVEAUX, MARKUS ZENK, CHRISTOPHE BERTHON, AND CHRISTIAN KLINGENBERG Abstract. The present paper concerns the derivation
More informationHIGH-ORDER ACCURATE METHODS BASED ON DIFFERENCE POTENTIALS FOR 2D PARABOLIC INTERFACE MODELS
HIGH-ORDER ACCURATE METHODS BASED ON DIFFERENCE POTENTIALS FOR 2D PARABOLIC INTERFACE MODELS JASON ALBRIGHT, YEKATERINA EPSHTEYN, AND QING XIA Abstract. Highly-accurate numerical methods that can efficiently
More informationLecture 12: Detailed balance and Eigenfunction methods
Lecture 12: Detailed balance and Eigenfunction methods Readings Recommended: Pavliotis [2014] 4.5-4.7 (eigenfunction methods and reversibility), 4.2-4.4 (explicit examples of eigenfunction methods) Gardiner
More informationApplication of Nodal Discontinuous Glaerkin Methods in Acoustic Wave Modeling
1 Application of Nodal Discontinuous Glaerkin Methods in Acoustic Wave Modeling Xin Wang ABSTRACT This work will explore the discontinuous Galerkin finite element method (DG-FEM) for solving acoustic wave
More informationFinite Difference and Finite Element Methods
Finite Difference and Finite Element Methods Georgy Gimel farb COMPSCI 369 Computational Science 1 / 39 1 Finite Differences Difference Equations 3 Finite Difference Methods: Euler FDMs 4 Finite Element
More information