Convergence of the MAC scheme for incompressible flows
|
|
- Kelly Johnson
- 5 years ago
- Views:
Transcription
1 Convergence of the MAC scheme for incompressible flows T. Galloue t?, R. Herbin?, J.-C. Latche??, K. Mallem????????? Aix-Marseille Universite I.R.S.N. Cadarache Universite d Annaba
2 Calanque de Cortiou Marseille main sewer
3 Incompressible Variable density flows Bounded domain (0, T ), R d, T > 0. Non homogeneous fluid : density ρ is not constant Incompressible fluid : div ū = 0, Incompressible variable density Navier-Stokes equations: t ρ + div( ρū) = 0, t ρ ū + div( ρ ū ū) ū + p = 0, div ū = 0, u = 0, u t=0 = u 0, ρ t=0 = ρ 0 ρ min > 0, (VDNS)
4 Variable density NS equations: a priori estimates maximum principle for ρ divū = 0 and t ρ + div( ρū) = 0 t ρ + ū ρ = 0 (transport equation) 0 < min ρ 0 ρ max ρ 0. kinetic energy balance: Multiply momentum balance eq by ū and use mass balance (twice) t ( 1 2 ρ ū 2 ) + div( 1 2 ρ ū 2 u) ū ū + p ū = 0. L ((0, T ); L 2 () d ) and L 2 ((0, T ); H 1 0 ()d ) bound: 1 t ρ(., t) ū(., t) 2 dx ū 2 dx dt = 1 ρ 0 ū 0 2 dx, t (0, T ). 2
5 Variable density NS equations: weak solutions Weak Solutions: ρ 0 L () such that ρ 0 > 0 and u 0 L 2 () d. A weak solution is a pair (ρ, u) such that: ρ L ((0, T ) ) and ρ > 0 a.e. in (0, T ). ū L ((0, T ); L 2 () d ) L 2 ((0, T ); H 1 0 ()d ) et divū = 0 a.e. in (0, T ). For any φ in C c ( [0, T )), T ( ρ t φ + ρū φ) dx dt = ρ 0 φ(., 0) dx. 0 for any v in C c ( [0, T ))d such that divv = 0, T ( ) ρū t v ( ρū ū) : v + ū : v dx dt = ρ 0 ū 0 v(., 0) dx. 0 Existence of a weak solution Simon 90. Discretization Liu Walkington Discontinuous Galerkin 2007 Latché Saleh Rannacher Turek 2016
6 The standard Marker-And-Cell scheme (Harlow Welsh 65) structured grids, pressure at the cell-center and normal velocity components at mid-edges, staggered meshes. Advantages and drawback + minimal number of unknowns, + no pressure stabilization needed, + simplicity and robustness of the scheme, needs local refinement for complex geometries. Convergence analysis: Shin Strickwerda 96, 97, stability, inf-sup condidion Nicolaides Wu 96, (ω, ψ) formulation Blanc 05 irregular rectangles, finite volume approach
7 The MAC grid (or Arakawa C) Primal mesh (pressure) T : K : Primal cell. E set of faces of T. E (i) ; i = 1,..., d : set of faces orthogonal to the i-th component of e i. Dual mesh (velocity) : Dσ dual cell associate to the face σ E (i). D σ K D σ D σ ρ = K T ρ K χ K, u (i) = σ E (i) u σχ Dσ
8 From the continuous to the discrete equations: time-implicit discretization Continuous equations ρ L ((0, T ) ), ū L ((0, T ); L 2 () d ) L 2 ((0, T ); H 1 0 ()d ), div ū = 0, t ρ + div( ρū) = 0, t ( ρ) ū + div( ρ ū ū) ū + p = 0, Discrete equations ρ (n+1), p (n+1) L T, u (n+1) H E, div T u (n+1) = 0, ð t ρ (n+1) T + div T ((ρu) (n+1) ) = 0 ð t (ρu) (n+1) + c E ((ρu) (n+1) )u (n+1) E u (n+1) + E p (n+1) = 0, Discrete time derivative: ð t (φ) (n+1) = 1 δt (φ(, t n+1) φ(, t n)), φ H T or φ H E. Discrete operators: div T u, div T ρu, c E ((ρu))u, E u, E p...
9 The time-implicit MAC scheme Discrete equations ρ, p piecewise constant on the primal mesh, u piecewise constant on the dual mesh ρ (n+1), p (n+1) L T, u (n+1) H E, div T u (n+1) = 0, ð t ρ (n+1) T + div T ((ρu) (n+1) ) = 0 ð t (ρu) (n+1) + c E ((ρu) (n+1) )u (n+1) E u (n+1) + E p (n+1) = 0, Discrete operators div T u, div T (ρu) : For K T : div K u (n+1) = 1 K 1 δt (ρ(n+1) K σ E(K ) ρ (n) K ) + 1 K σ u K,σ = 0, with u K,σ = ±u σ, σ E(K ) with F K,σ = σ ρ up σ u K,σ, F (n+1) K,σ = 0, ρ up σ : upwind approximation of ρ on σ = ρ(n+1) > 0. K
10 The time-implicit MAC scheme Discrete equations, ρ, p piecewise constant on the primal mesh, and u piecewise constant on the dual mesh ρ (n+1), p (n+1) L T, u (n+1) H E, div T u (n+1) = 0, ð t ρ (n+1) T + div T ((ρu) (n+1) ) = 0 ð t (ρu) (n+1) + c E ((ρu) (n+1) )u (n+1) E u (n+1) + E p (n+1) = 0, For σ E : 1 δt (ρ(n+1) σ u σ ρ (n) σ u(n) σ ) + 1 F (n+1) σ,ɛ D σ u(n+1) ɛ ( u) (n+1) σ ɛ Ē(Dσ ) ρ σ and F σ,ɛ chosen later. centered approximation for u ɛ, with: ( u) σ = two point flux FV approximation on the velocity meshes, E p (n+1) = ( p) σ, ( p) σ = σ D (p L p K )n K,σ, σ E σ int K D σ ɛ = σ σ + ( p) (n+1) σ σ σ = K L ɛ == σ σ σ = 0, σ E int. L
11 Discrete duality, coercivity, inf-sup property Discrete divergence-gradient duality q div T v + E q v = 0. Coercivity d u 1,E = D σ E (i) u (i) u σ = E u L 2 () d C u L 2 () d i=1 σ E (i) Inf-sup property q L T, v H E ; q = div T v and v E C p L 2 () d (discrete version of Necas Lemma, Strickwerda 90, Blanc 05)
12 Discrete kinetic energy inequality, choice of ρ σ and Discrete equations If a mass balance holds on the dual cells: D σ 2δt (ρ σ ρ σ ) + σ E F σ,ɛ = 0, and ρ σ > 0, then 1 ( ρ σuσ 2 ρ 2δt σ(uσ) 2) D σ Proof Choice of ρ σ and F σ,ɛ ɛ Ẽ(D σ) ɛ=σ σ F σ,ɛu σu σ ( u) σu σ+(ðp) σu σ 0 Herbin, Kheriji, Latché 2013 τ ɛ τ F σ,ɛ = D σ ρ σ = D K,σ ρ K + D L,σ ρ L, 1 2 (F K,τ + F K,τ ) if ɛ σ, ɛ τ τ 1 2 ( F K,σ + F K,σ ) if ɛ K = [σ, σ ] K ɛ K σ = K L D σ L
13 Convergence analysis Theorem : (T m, E m), δt m sequence of meshes and time steps h m 0, δt m 0 + suitable assumptions on the mesh, There exists ρ m, u m solution to the scheme, and, up to a subsequence ρ m ρ in L p ( (0, T )), p [1, + [, u m ū in L 2 ( (0, T )) ( ρ, ū) is a weak solution to (VDNS) Sketch of proof A priori estimates on the approximate solutions ρ m, u m existence of a solution to the scheme. Compactness thanks to discrete functional analysis tools: Up to a subsequence, ρ m ρ, u m ū for some norms. Passage to the limit in the weak form of the scheme, ( ρ, ū) is a weak solution to (VDNS).
14 Convergence analysis: estimates and existence 1 - A priori estimates. Maximum principle (upwind discretization of the mass balance): min ρ 0(x) ρ max ρ 0(x). x x Kinetic energy identity + grad-div duality: u L ((0,T );L 2 () d ) + u L 2 ((0,T );H 1 E,0 ) C Estimate of the momentum balance convection form: C E (ρu)v w C ηt ρ L () u L 4 () d v L 4 () d w 1,E,0 C ηt ρ L () u 1,E,0 v 1,E,0 w 1,E, Existence of a solution thanks to a topological degree argument.
15 Easy case : homogeneous fluid, constant density Estimate on u m in L 2 (0, T ; H 1 E,0 ) Estimate on ð t u m N 1 ð t u m L 1 (0,T ;E ) = δt ð t u E C where v E = max E E E ϕ E v ϕ dx ; n=0 E ϕ 1,E,0 1 Discrete Aubin-Simon theorem = compactness on u m in L 2 ((0, T ), L 2 ()) Discrete Aubin-Simon theorem 1 p < +, B Banach space, B m B, dim B m < +, Xm and Ym norms on B m s.t. : if w m Xm is bounded then: (i) w B; w m w in B up to a subsequence, (ii) If w m w in B, and w m Ym 0 then w = 0. Let T > 0 and (u m) m N be a sequence of L p ((0, T ), B), piecewise constant on the time intervals, such that 1. ( T 0 um p Xm dt) m N is bounded, 2. ( T 0 ðum ym dt) m N is bounded. Then there exists u L p ((0, T ), B) such that, up to a subsequence, u m u in L p ((0, T ), B). u m ū L 2 ((0, T ), L 2 ()) + (u m) m N bounded in L 2 ((0, T ), H 1 E,0 ) ū L 2 ((0, T ), H 1 0 ()).
16 Variable density, estimate on the translates No estimate on ð t u m Estimate on the time translates (continuous case: Boyer Fabrie 13 chapter 2); T τ u m(x, t + τ) u m(x, t) 2 ρ max dx dt C ηt,t ( u m 3 0 ρ L 2 (H min E,0 ) + 1) τ + δt Kolmogorov compactness theorem compactness of u m in L 2 ((0, T ), L 2 ()). τ + δt : No estimate on time derivative. Will not generalize to the compressible case because of ρ min
17 Variable density, convergence Convergence proof, sketch : u m ū in L 2 (0, T ; L 2 () d ), ρ m ρ weakly in L 2 (0, T ; L 2 ()) T T ρ m u m ψ ρ ū ψ for any ψ L ( (0, T )) d ( ) 0 0 In particular for ψ = mϕ, T T div m(ρ m u m) ϕ dxdt = ρ m u m mϕ dxdt ( ρ ū) ϕ dxdt 0 0 (thanks to weak BV estimate) ð t ρ m t ρ, ð t (ρ m u m) t (ρu), divu m divū div mū = 0 (linear operators) ρ mu m bounded in L 2 (0, T ; L 2 () d ) = ρ mu m converges weakly in L 2 (0, T ; L 2 () d ). Hence by ( ), ρ mu m ρū weakly in L 2 (0, T ; L 2 () d ) and since u m u in L 2 (0, T ; L 2 () d ) ρ mu m u m ϕ ρū u ϕ
18 Variable density case, convergence, alternate proof Compactness of ρ m in L 2 (H 1 ) thanks to Aubin-Simon with B = H 1 () X = L 2 () Y = W 1,1 () X = L 2 () compactly embedded in B = H 1 () t ρ m bounded in L 1 (0, T, W 1,1 ()) ρ m bounded in L 2 (0, T, L 2 ()) ρ m compact in L 2 (0, T, H 1 ()) u m bounded in L 2 (0, T, H 1 0 ()) hence (weak strong product) ρ mu m ρū weakly in L 1. Same ideas for the convergence of ρ mu m u m... Thierry Gallouët, 2017
19 Recent and on going work Locally refined grids : with Chénier, Eymard, Gallouët for modified MAC, Latché, Piar, Saleh for RT grids. Stability and weak consistency for the full NS equations (with D. Grapsas, W. Kheriji, J.-C. Latché) Convergence for steady-state barotropic NS (γ 3) (with J.-C. Latché, T. Gallouët and D. Maltese) Error estimates for barotropic NS, adaptation of strong weak uniqueness (with D. Maltese, and A. Novotny). Stability and weak consistency for the Euler equations with higher order schemes and entropy consistency (with T. Gallouët, J.-C. Latché, N. Therme) Low Mach limit for barotropic NS (with J.-C. Latché, K. Saleh). Convergence of the fully discrete pressure correction scheme. TODO Convergence for the time-dependent barotropic compressible Navier-Stokes. TODO
On a class of numerical schemes. for compressible flows
On a class of numerical schemes for compressible flows R. Herbin, with T. Gallouët, J.-C. Latché L. Gastaldo, D. Grapsas, W. Kheriji, T.T. N Guyen, N. Therme, C. Zaza. Aix-Marseille Université I.R.S.N.
More informationarxiv: v1 [math.na] 17 Jan 2017
CONVERGENCE OF THE MARKER-AND-CELL SCHEME FOR THE INCOMPRESSIBLE NAVIER-STOKES EQUATIONS ON NON-UNIFORM GRIDS T. Gallouët 1, R. Herbin 2, J.-C. Latché 3 and K. Mallem 4 arxiv:1701.04553v1 [math.na] 17
More informationFinite difference MAC scheme for compressible Navier-Stokes equations
Finite difference MAC scheme for compressible Navier-Stokes equations Bangwei She with R. Hošek, H. Mizerova Workshop on Mathematical Fluid Dynamics, Bad Boll May. 07 May. 11, 2018 Compressible barotropic
More informationStudy of a numerical scheme for miscible two-phase flow in porous media. R. Eymard, Université Paris-Est
Study of a numerical scheme for miscible two-phase flow in porous media R. Eymard, Université Paris-Est in collaboration with T. Gallouët and R. Herbin, Aix-Marseille Université C. Guichard, Université
More informationConsistency analysis of a 1D Finite Volume scheme for barotropic Euler models
Consistency analysis of a 1D Finite Volume scheme for barotropic Euler models F Berthelin 1,3, T Goudon 1,3, and S Mineaud,3 1 Inria, Sophia Antipolis Méditerranée Research Centre, Proect COFFEE Inria,
More informationarxiv: v1 [math.na] 16 Oct 2007
Mathematical Modelling and Numerical Analysis Modélisation Mathématique et Analyse Numérique Will be set by the publisher arxiv:070.2987v [math.na] 6 Oct 2007 AN UNCONDITIONNALLY STABLE PRESSURE CORRECTION
More informationEntropy estimates for a class of schemes for the euler equations
Entropy estimates for a class of schemes for the euler equations Thierry Gallouët, Raphaele Herbin, J.-C Latché, N Therme To cite this version: Thierry Gallouët, Raphaele Herbin, J.-C Latché, N Therme.
More informationVector and scalar penalty-projection methods
Numerical Flow Models for Controlled Fusion - April 2007 Vector and scalar penalty-projection methods for incompressible and variable density flows Philippe Angot Université de Provence, LATP - Marseille
More informationVariable Exponents Spaces and Their Applications to Fluid Dynamics
Variable Exponents Spaces and Their Applications to Fluid Dynamics Martin Rapp TU Darmstadt November 7, 213 Martin Rapp (TU Darmstadt) Variable Exponent Spaces November 7, 213 1 / 14 Overview 1 Variable
More informationOn two fractional step finite volume and finite element schemes for reactive low Mach number flows
Fourth International Symposium on Finite Volumes for Complex Applications - Problems and Perspectives - July 4-8, 2005 / Marrakech, Morocco On two fractional step finite volume and finite element schemes
More informationMixed Finite Element Methods. Douglas N. Arnold, University of Minnesota The 41st Woudschoten Conference 5 October 2016
Mixed Finite Element Methods Douglas N. Arnold, University of Minnesota The 41st Woudschoten Conference 5 October 2016 Linear elasticity Given the load f : Ω R n, find the displacement u : Ω R n and the
More informationIntroduction to Finite Volume projection methods. On Interfaces with non-zero mass flux
Introduction to Finite Volume projection methods On Interfaces with non-zero mass flux Rupert Klein Mathematik & Informatik, Freie Universität Berlin Summerschool SPP 1506 Darmstadt, July 09, 2010 Introduction
More informationCONVERGENCE THEORY. G. ALLAIRE CMAP, Ecole Polytechnique. 1. Maximum principle. 2. Oscillating test function. 3. Two-scale convergence
1 CONVERGENCE THEOR G. ALLAIRE CMAP, Ecole Polytechnique 1. Maximum principle 2. Oscillating test function 3. Two-scale convergence 4. Application to homogenization 5. General theory H-convergence) 6.
More information.u= 0 ρ( u t +(u. )u)= ρ g p+.[µ( u+ t u)]
THETIS is a numerical simulation tool developed by University of Bordeaux. It is a versatile code to solve different problems: fluid flows, heat transfers, scalar transports or porous mediums. The potential
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 informationOn the stability of colocated clustered finite volume simplicial discretizations for the 2D Stokes problem
On the stability of colocated clustered finite volume simplicial discretizations for the 2D Stokes problem Robert Eymard, Raphaele Herbin, Jean-Claude Latché, Bruno Piar To cite this version: Robert Eymard,
More informationAMS subject classifications. Primary, 65N15, 65N30, 76D07; Secondary, 35B45, 35J50
A SIMPLE FINITE ELEMENT METHOD FOR THE STOKES EQUATIONS LIN MU AND XIU YE Abstract. The goal of this paper is to introduce a simple finite element method to solve the Stokes equations. This method is in
More informationThe incompressible Navier-Stokes equations in vacuum
The incompressible, Université Paris-Est Créteil Piotr Bogus law Mucha, Warsaw University Journées Jeunes EDPistes 218, Institut Elie Cartan, Université de Lorraine March 23th, 218 Incompressible Navier-Stokes
More informationThe CG1-DG2 method for conservation laws
for conservation laws Melanie Bittl 1, Dmitri Kuzmin 1, Roland Becker 2 MoST 2014, Germany 1 Dortmund University of Technology, Germany, 2 University of Pau, France CG1-DG2 Method - Motivation hp-adaptivity
More informationKonstantinos Chrysafinos 1 and L. Steven Hou Introduction
Mathematical Modelling and Numerical Analysis Modélisation Mathématique et Analyse Numérique Will be set by the publisher ANALYSIS AND APPROXIMATIONS OF THE EVOLUTIONARY STOKES EQUATIONS WITH INHOMOGENEOUS
More informationPart IV: Numerical schemes for the phase-filed model
Part IV: Numerical schemes for the phase-filed model Jie Shen Department of Mathematics Purdue University IMS, Singapore July 29-3, 29 The complete set of governing equations Find u, p, (φ, ξ) such that
More informationOn Pressure Stabilization Method and Projection Method for Unsteady Navier-Stokes Equations 1
On Pressure Stabilization Method and Projection Method for Unsteady Navier-Stokes Equations 1 Jie Shen Department of Mathematics, Penn State University University Park, PA 1682 Abstract. We present some
More informationWeak-Strong Uniqueness of the Navier-Stokes-Smoluchowski System
Weak-Strong Uniqueness of the Navier-Stokes-Smoluchowski System Joshua Ballew University of Maryland College Park Applied PDE RIT March 4, 2013 Outline Description of the Model Relative Entropy Weakly
More informationEXISTENCE AND REGULARITY RESULTS FOR SOME NONLINEAR PARABOLIC EQUATIONS
EXISTECE AD REGULARITY RESULTS FOR SOME OLIEAR PARABOLIC EUATIOS Lucio BOCCARDO 1 Andrea DALL AGLIO 2 Thierry GALLOUËT3 Luigi ORSIA 1 Abstract We prove summability results for the solutions of nonlinear
More informationControl of Interface Evolution in Multi-Phase Fluid Flows
Control of Interface Evolution in Multi-Phase Fluid Flows Markus Klein Department of Mathematics University of Tübingen Workshop on Numerical Methods for Optimal Control and Inverse Problems Garching,
More informationDiscontinuous Galerkin Methods
Discontinuous Galerkin Methods Joachim Schöberl May 20, 206 Discontinuous Galerkin (DG) methods approximate the solution with piecewise functions (polynomials), which are discontinuous across element interfaces.
More information1 The Stokes System. ρ + (ρv) = ρ g(x), and the conservation of momentum has the form. ρ v (λ 1 + µ 1 ) ( v) µ 1 v + p = ρ f(x) in Ω.
1 The Stokes System The motion of a (possibly compressible) homogeneous fluid is described by its density ρ(x, t), pressure p(x, t) and velocity v(x, t). Assume that the fluid is barotropic, i.e., the
More informationConvergence analysis of a finite volume method for the Stokes system using non-conforming arguments
IMA Journal of Numerical Analysis (2005) 25, 523 548 doi:10.1093/imanum/dri007 Advance Access publication on February 7, 2005 Convergence analysis of a finite volume method for the Stokes system using
More informationConvergence of Finite Volumes schemes for an elliptic-hyperbolic system with boundary conditions
Convergence of Finite Volumes schemes for an elliptic-hyperbolic system with boundary conditions Marie Hélène Vignal UMPA, E.N.S. Lyon 46 Allée d Italie 69364 Lyon, Cedex 07, France abstract. We are interested
More informationA second order scheme on staggered grids for detonation in gases
A second order scheme on staggered grids for detonation in gases Chady Zaza IRSN PSN-RES/SA2I/LIE August 1, 2016 HYP 2016, RWTH Aachen Joint work with Raphaèle Herbin, Jean-Claude Latché and Nicolas Therme
More informationLecture Note III: Least-Squares Method
Lecture Note III: Least-Squares Method Zhiqiang Cai October 4, 004 In this chapter, we shall present least-squares methods for second-order scalar partial differential equations, elastic equations of solids,
More informationA Least-Squares Finite Element Approximation for the Compressible Stokes Equations
A Least-Squares Finite Element Approximation for the Compressible Stokes Equations Zhiqiang Cai, 1 Xiu Ye 1 Department of Mathematics, Purdue University, 1395 Mathematical Science Building, West Lafayette,
More informationConvergence analysis of a locally stabilized collocated finite volume scheme for incompressible flows
Convergence analysis of a locally stabilized collocated finite volume scheme for incompressible flows Robert Eymard,, Raphaele Herbin, Jean-Claude Latché, Bruno Piar To cite this version: Robert Eymard,,
More informationAn adaptive projection scheme for all-mach flows applied to shock diffraction
An adaptive projection scheme for all-mach flows applied to shock diffraction CEA DEN/DM2S/STMF/LMEC AMU I2M UMR 7373/AA Group January 20, 2015 All-Mach solver? Approximate Riemann solvers for M 1 (Guillard-Viozat,
More informationSTOCHASTIC NAVIER-STOKES EQUATIONS FOR COMPRESSIBLE FLUIDS
STOCHASTIC NAVIER-STOKES EUATIONS FOR COMPRESSIBLE FLUIDS DOMINIC BREIT AND MARTINA HOFMANOVÁ Abstract. We study the Navier-Stokes equations governing the motion of an isentropic compressible fluid in
More informationAA210A Fundamentals of Compressible Flow. Chapter 1 - Introduction to fluid flow
AA210A Fundamentals of Compressible Flow Chapter 1 - Introduction to fluid flow 1 1.2 Conservation of mass Mass flux in the x-direction [ ρu ] = M L 3 L T = M L 2 T Momentum per unit volume Mass per unit
More informationThe gradient discretisation method
The gradient discretisation method Jérôme Droniou, Robert Eymard, Thierry Gallouët, Cindy Guichard, Raphaele Herbin To cite this version: Jérôme Droniou, Robert Eymard, Thierry Gallouët, Cindy Guichard,
More informationApproximation of fluid-structure interaction problems with Lagrange multiplier
Approximation of fluid-structure interaction problems with Lagrange multiplier Daniele Boffi Dipartimento di Matematica F. Casorati, Università di Pavia http://www-dimat.unipv.it/boffi May 30, 2016 Outline
More informationVariational and Topological methods : Theory, Applications, Numerical Simulations, and Open Problems 6-9 June 2012, Northern Arizona University
Variational and Topological methods : Theory, Applications, Numerical Simulations, and Open Problems 6-9 June 22, Northern Arizona University Some methods using monotonicity for solving quasilinear parabolic
More informationMixed Finite Elements Method
Mixed Finite Elements Method A. Ratnani 34, E. Sonnendrücker 34 3 Max-Planck Institut für Plasmaphysik, Garching, Germany 4 Technische Universität München, Garching, Germany Contents Introduction 2. Notations.....................................
More informationTHE STOKES SYSTEM R.E. SHOWALTER
THE STOKES SYSTEM R.E. SHOWALTER Contents 1. Stokes System 1 Stokes System 2 2. The Weak Solution of the Stokes System 3 3. The Strong Solution 4 4. The Normal Trace 6 5. The Mixed Problem 7 6. The Navier-Stokes
More informationOn a Suitable Weak Solution of the Navier Stokes Equation with the Generalized Impermeability Boundary Conditions
Proceedings of the 3rd IASME/WSEAS Int. Conf. on FLUID DYNAMICS & AERODYNAMICS, Corfu, Greece, August -, 5 pp36-41 On a Suitable Weak Solution of the Navier Stokes Equation with the Generalized Impermeability
More informationWeek 6 Notes, Math 865, Tanveer
Week 6 Notes, Math 865, Tanveer. Energy Methods for Euler and Navier-Stokes Equation We will consider this week basic energy estimates. These are estimates on the L 2 spatial norms of the solution u(x,
More informationRelative entropies, suitable weak solutions, and weak-strong uniqueness for the compressible Navier-Stokes system
Relative entropies, suitable weak solutions, and weak-strong uniqueness for the compressible Navier-Stokes system Institute of Mathematics, Academy of Sciences of the Czech Republic, Prague joint work
More informationZdzislaw Brzeźniak. Department of Mathematics University of York. joint works with Mauro Mariani (Roma 1) and Gaurav Dhariwal (York)
Navier-Stokes equations with constrained L 2 energy of the solution Zdzislaw Brzeźniak Department of Mathematics University of York joint works with Mauro Mariani (Roma 1) and Gaurav Dhariwal (York) LMS
More informationMeasure-valued - strong uniqueness for hyperbolic systems
Measure-valued - strong uniqueness for hyperbolic systems joint work with Tomasz Debiec, Eduard Feireisl, Ondřej Kreml, Agnieszka Świerczewska-Gwiazda and Emil Wiedemann Institute of Mathematics Polish
More informationFrom Boltzmann Equations to Gas Dynamics: From DiPerna-Lions to Leray
From Boltzmann Equations to Gas Dynamics: From DiPerna-Lions to Leray C. David Levermore Department of Mathematics and Institute for Physical Science and Technology University of Maryland, College Park
More informationON LIQUID CRYSTAL FLOWS WITH FREE-SLIP BOUNDARY CONDITIONS. Chun Liu and Jie Shen
DISCRETE AND CONTINUOUS Website: http://aimsciences.org DYNAMICAL SYSTEMS Volume 7, Number2, April2001 pp. 307 318 ON LIQUID CRYSTAL FLOWS WITH FREE-SLIP BOUNDARY CONDITIONS Chun Liu and Jie Shen Department
More informationIntroduction to discrete functional analysis techniques for the numerical study of diffusion equations with irregular data
Introduction to discrete functional analysis techniques for the numerical study of diffusion equations with irregular data arxiv:1505.02567v2 [math.na] 24 Feb 2016 J. Droniou 1 February 25, 2016 Abstract
More informationWELL POSEDNESS OF PROBLEMS I
Finite Element Method 85 WELL POSEDNESS OF PROBLEMS I Consider the following generic problem Lu = f, where L : X Y, u X, f Y and X, Y are two Banach spaces We say that the above problem is well-posed (according
More informationDiplomarbeit. Ausgefu hrt am Institut fu r. Analysis und Scientific Computing. der Technischen Universita t Wien. unter der Anleitung von
Diplomarbeit Pressure Robust Discretizations for Navier Stokes Equations: Divergence-free Reconstruction for Taylor-Hood Elements and High Order Hybrid Discontinuous Galerkin Methods Ausgefu hrt am Institut
More informationNumerical Methods for the Navier-Stokes equations
Arnold Reusken Numerical Methods for the Navier-Stokes equations January 6, 212 Chair for Numerical Mathematics RWTH Aachen Contents 1 The Navier-Stokes equations.............................................
More informationHomogenization of the compressible Navier Stokes equations in domains with very tiny holes
Homogenization of the compressible Navier Stokes equations in domains with very tiny holes Yong Lu Sebastian Schwarzacher Abstract We consider the homogenization problem of the compressible Navier Stokes
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 informationConvection and total variation flow
Convection and total variation flow R. Eymard joint work wit F. Boucut and D. Doyen LAMA, Université Paris-Est marc, 2nd, 2016 Flow wit cange-of-state simplified model for Bingam fluid + Navier-Stokes
More informationAn Upwinding Cell-Centered Method with Piecewise Constant Velocity over Covolumes
An Upwinding Cell-Centered Method with Piecewise Constant Velocity over Covolumes So Hsiang Chou, 1 Panayot S. Vassilevski 2 1 Department of Mathematics and Statistics, Bowling Green State University,
More informationFinite Volume Methods Schemes and Analysis course at the University of Wroclaw
Finite Volume Methods Schemes and Analysis course at the University of Wroclaw Robert Eymard 1, Thierry Gallouët 2 and Raphaèle Herbin 3 April, 2008 1 Université Paris-Est Marne-la-Vallée 2 Université
More informationAn Equal-order DG Method for the Incompressible Navier-Stokes Equations
An Equal-order DG Method for the Incompressible Navier-Stokes Equations Bernardo Cockburn Guido anschat Dominik Schötzau Journal of Scientific Computing, vol. 40, pp. 188 10, 009 Abstract We introduce
More informationFormulation of the problem
TOPICAL PROBLEMS OF FLUID MECHANICS DOI: https://doi.org/.43/tpfm.27. NOTE ON THE PROBLEM OF DISSIPATIVE MEASURE-VALUED SOLUTIONS TO THE COMPRESSIBLE NON-NEWTONIAN SYSTEM H. Al Baba, 2, M. Caggio, B. Ducomet
More informationINSTITUTE OF MATHEMATICS THE CZECH ACADEMY OF SCIENCES
INSTITUTE OF MATHEMATICS THE CZECH ACADEMY OF SCIENCES Asymptotic preserving error estimates for numerical solutions of compressible Navier-Stokes equations in the low Mach number regime Eduard Feireisl
More informationGoal. Robust A Posteriori Error Estimates for Stabilized Finite Element Discretizations of Non-Stationary Convection-Diffusion Problems.
Robust A Posteriori Error Estimates for Stabilized Finite Element s of Non-Stationary Convection-Diffusion Problems L. Tobiska and R. Verfürth Universität Magdeburg Ruhr-Universität Bochum www.ruhr-uni-bochum.de/num
More informationNumerical Analysis of Higher Order Discontinuous Galerkin Finite Element Methods
Numerical Analysis of Higher Order Discontinuous Galerkin Finite Element Methods Contents Ralf Hartmann Institute of Aerodynamics and Flow Technology DLR (German Aerospace Center) Lilienthalplatz 7, 3808
More informationSimulating Drag Crisis for a Sphere Using Skin Friction Boundary Conditions
Simulating Drag Crisis for a Sphere Using Skin Friction Boundary Conditions Johan Hoffman May 14, 2006 Abstract In this paper we use a General Galerkin (G2) method to simulate drag crisis for a sphere,
More informationEuler Equations: local existence
Euler Equations: local existence Mat 529, Lesson 2. 1 Active scalars formulation We start with a lemma. Lemma 1. Assume that w is a magnetization variable, i.e. t w + u w + ( u) w = 0. If u = Pw then u
More informationUNIFORM PRECONDITIONERS FOR A PARAMETER DEPENDENT SADDLE POINT PROBLEM WITH APPLICATION TO GENERALIZED STOKES INTERFACE EQUATIONS
UNIFORM PRECONDITIONERS FOR A PARAMETER DEPENDENT SADDLE POINT PROBLEM WITH APPLICATION TO GENERALIZED STOKES INTERFACE EQUATIONS MAXIM A. OLSHANSKII, JÖRG PETERS, AND ARNOLD REUSKEN Abstract. We consider
More informationZdzislaw Brzeźniak. Department of Mathematics University of York. joint works with Mauro Mariani (Roma 1) and Gaurav Dhariwal (York)
Navier-Stokes equations with constrained L 2 energy of the solution Zdzislaw Brzeźniak Department of Mathematics University of York joint works with Mauro Mariani (Roma 1) and Gaurav Dhariwal (York) Stochastic
More informationGLOBAL WEAK SOLUTION TO THE MAGNETOHYDRODYNAMIC AND VLASOV EQUATIONS
GLOBAL WEAK SOLUTION TO THE MAGNETOHYDRODYNAMIC AND VLASOV EQUATIONS ROBIN MING CHEN, JILONG HU, AND DEHUA WANG Abstract. An initial-boundary value problem for the fluid-particle system of the inhomogeneous
More informationSimulation of free surface fluids in incompressible dynamique
Simulation of free surface fluids in incompressible dynamique Dena Kazerani INRIA Paris Supervised by Pascal Frey at Laboratoire Jacques-Louis Lions-UPMC Workshop on Numerical Modeling of Liquid-Vapor
More informationSTUDY OF DISCRETE DUALITY FINITE VOLUME SCHEMES FOR THE PEACEMAN MODEL
SUDY OF DISCREE DUALIY FINIE VOLUME SCHEMES FOR HE PEACEMAN MODEL C. CHAINAIS-HILLAIRE, S. KRELL, AND A. MOUON Abstract. In this paper, we are interested in the finite volume approximation of a system
More informationSingular Limits of the Klein-Gordon Equation
Singular Limits of the Klein-Gordon Equation Abstract Chi-Kun Lin 1 and Kung-Chien Wu 2 Department of Applied Mathematics National Chiao Tung University Hsinchu 31, TAIWAN We establish the singular limits
More informationError estimates for the discretization of the velocity tracking problem
Numerische Mathematik manuscript No. (will be inserted by the editor) Error estimates for the discretization of the velocity tracking problem Eduardo Casas 1, Konstantinos Chrysafinos 2 1 Departamento
More informationLeast-Squares Finite Element Methods
Pavel В. Bochev Max D. Gunzburger Least-Squares Finite Element Methods Spri ringer Contents Part I Survey of Variational Principles and Associated Finite Element Methods 1 Classical Variational Methods
More informationSome remarks on grad-div stabilization of incompressible flow simulations
Some remarks on grad-div stabilization of incompressible flow simulations Gert Lube Institute for Numerical and Applied Mathematics Georg-August-University Göttingen M. Stynes Workshop Numerical Analysis
More informationOn the domain dependence of solutions to the compressible Navier-Stokes equations of an isothermal fluid
Nečas Center for Mathematical Modeling On the domain dependence of solutions to the compressible Navier-Stokes equations of an isothermal fluid Nikola Hlaváčová Preprint no. 211-2 Research Team 1 Mathematical
More informationINSTITUTE of MATHEMATICS. ACADEMY of SCIENCES of the CZECH REPUBLIC
INSTITUTEofMATHEMATICS Academy of Sciences Czech Republic INSTITUTE of MATHEMATICS ACADEMY of SCIENCES of the CZECH REPUBLIC Error estimates for a numerical method for the compressible Navier-Stokes system
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 informationA NOTE ON THE LADYŽENSKAJA-BABUŠKA-BREZZI CONDITION
A NOTE ON THE LADYŽENSKAJA-BABUŠKA-BREZZI CONDITION JOHNNY GUZMÁN, ABNER J. SALGADO, AND FRANCISCO-JAVIER SAYAS Abstract. The analysis of finite-element-like Galerkin discretization techniques for the
More informationPARABOLIC CONTROL PROBLEMS IN SPACE-TIME MEASURE SPACES
ESAIM: Control, Optimisation and Calculus of Variations URL: http://www.emath.fr/cocv/ Will be set by the publisher PARABOLIC CONTROL PROBLEMS IN SPACE-TIME MEASURE SPACES Eduardo Casas and Karl Kunisch
More informationSpatial discretization scheme for incompressible viscous flows
Spatial discretization scheme for incompressible viscous flows N. Kumar Supervisors: J.H.M. ten Thije Boonkkamp and B. Koren CASA-day 2015 1/29 Challenges in CFD Accuracy a primary concern with all CFD
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 informationA primer on Numerical methods for elasticity
A primer on Numerical methods for elasticity Douglas N. Arnold, University of Minnesota Complex materials: Mathematical models and numerical methods Oslo, June 10 12, 2015 One has to resort to the indignity
More informationAsymptotic preserving error estimates for numerical solutions of compressible Navier-Stokes equations in the low Mach number regime
Asymptotic preserving error estimates for numerical solutions of compressible Navier-Stokes equations in the low Mach number regime Eduard Feireisl Mária Lukáčová-Medviďová Šárka Nečasová Antonín Novotný
More informationANALYSIS AND CONVERGENCE OF A COVOLUME METHOD FOR THE GENERALIZED STOKES PROBLEM
MATHEMATICS OF COMPUTATION Volume 66, Number 217, January 1997, Pages 85 104 S 0025-5718(97)00792-8 ANALYSIS AND CONVERGENCE OF A COVOLUME METHOD FOR THE GENERALIZED STOKES PROBLEM S. H. CHOU Abstract.
More information1. Introduction. The Stokes problem seeks unknown functions u and p satisfying
A DISCRETE DIVERGENCE FREE WEAK GALERKIN FINITE ELEMENT METHOD FOR THE STOKES EQUATIONS LIN MU, JUNPING WANG, AND XIU YE Abstract. A discrete divergence free weak Galerkin finite element method is developed
More informationTRANSPORT IN POROUS MEDIA
1 TRANSPORT IN POROUS MEDIA G. ALLAIRE CMAP, Ecole Polytechnique 1. Introduction 2. Main result in an unbounded domain 3. Asymptotic expansions with drift 4. Two-scale convergence with drift 5. The case
More informationFinite volumes schemes preserving the low Mach number limit for the Euler system
Finite volumes schemes preserving the low Mach number limit for the Euler system M.-H. Vignal Low Velocity Flows, Paris, Nov. 205 Giacomo Dimarco, Univ. de Ferrara, Italie Raphael Loubere, IMT, CNRS, France
More informationEXISTENCE OF SOLUTIONS TO THE CAHN-HILLIARD/ALLEN-CAHN EQUATION WITH DEGENERATE MOBILITY
Electronic Journal of Differential Equations, Vol. 216 216), No. 329, pp. 1 22. ISSN: 172-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu EXISTENCE OF SOLUTIONS TO THE CAHN-HILLIARD/ALLEN-CAHN
More informationConservation of Mass. Computational Fluid Dynamics. The Equations Governing Fluid Motion
http://www.nd.edu/~gtryggva/cfd-course/ http://www.nd.edu/~gtryggva/cfd-course/ Computational Fluid Dynamics Lecture 4 January 30, 2017 The Equations Governing Fluid Motion Grétar Tryggvason Outline Derivation
More informationA Bound-Preserving Fourth Order Compact Finite Difference Scheme for Scalar Convection Diffusion Equations
A Bound-Preserving Fourth Order Compact Finite Difference Scheme for Scalar Convection Diffusion Equations Hao Li Math Dept, Purdue Univeristy Ocean University of China, December, 2017 Joint work with
More informationSECOND-ORDER FULLY DISCRETIZED PROJECTION METHOD FOR INCOMPRESSIBLE NAVIER-STOKES EQUATIONS
Tenth MSU Conference on Differential Equations and Computational Simulations. Electronic Journal of Differential Equations, Conference 3 (06), pp. 9 0. ISSN: 07-669. URL: http://ejde.math.txstate.edu or
More informationOn quasi-richards equation and finite volume approximation of two-phase flow with unlimited air mobility
On quasi-richards equation and finite volume approximation of two-phase flow with unlimited air mobility B. Andreianov 1, R. Eymard 2, M. Ghilani 3,4 and N. Marhraoui 4 1 Université de Franche-Comte Besançon,
More informationChapter 2. General concepts. 2.1 The Navier-Stokes equations
Chapter 2 General concepts 2.1 The Navier-Stokes equations The Navier-Stokes equations model the fluid mechanics. This set of differential equations describes the motion of a fluid. In the present work
More informationOn some singular limits for an atmosphere flow
On some singular limits for an atmosphere flow Donatella Donatelli Dipartimento di Ingegneria e Scienze dell Informazione e Matematica Università degli Studi dell Aquila 67100 L Aquila, Italy donatella.donatelli@univaq.it
More informationFINITE ELEMENT APPROXIMATION OF STOKES-LIKE SYSTEMS WITH IMPLICIT CONSTITUTIVE RELATION
Proceedings of ALGORITMY pp. 9 3 FINITE ELEMENT APPROXIMATION OF STOKES-LIKE SYSTEMS WITH IMPLICIT CONSTITUTIVE RELATION JAN STEBEL Abstract. The paper deals with the numerical simulations of steady flows
More informationAdaptive methods for control problems with finite-dimensional control space
Adaptive methods for control problems with finite-dimensional control space Saheed Akindeinde and Daniel Wachsmuth Johann Radon Institute for Computational and Applied Mathematics (RICAM) Austrian Academy
More informationarxiv: v1 [math.na] 30 Jan 2018
FINITE ELEMENT CONVERGENCE FOR THE TIME-DEPENDENT JOULE HEATING PROBLEM WITH MIXED BOUNDARY CONDITIONS arxiv:181.115v1 [math.na] 3 Jan 218 MAX JENSEN 1, AXEL MÅLQVIST 2, AND ANNA PERSSON 2 Abstract. We
More informationA Two-Grid Stabilization Method for Solving the Steady-State Navier-Stokes Equations
A Two-Grid Stabilization Method for Solving the Steady-State Navier-Stokes Equations Songul Kaya and Béatrice Rivière Abstract We formulate a subgrid eddy viscosity method for solving the steady-state
More informationProof of the existence (by a contradiction)
November 6, 2013 ν v + (v )v + p = f in, divv = 0 in, v = h on, v velocity of the fluid, p -pressure. R n, n = 2, 3, multi-connected domain: (NS) S 2 S 1 Incompressibility of the fluid (divv = 0) implies
More informationGradient Schemes for an Obstacle Problem
Gradient Schemes for an Obstacle Problem Y. Alnashri and J. Droniou Abstract The aim of this work is to adapt the gradient schemes, discretisations of weak variational formulations using independent approximations
More informationExact controllability of the superlinear heat equation
Exact controllability of the superlinear heat equation Youjun Xu 1,2, Zhenhai Liu 1 1 School of Mathematical Sciences and Computing Technology, Central South University, Changsha, Hunan 410075, P R China
More information