A Godunov-type method for the diphasic Baer-Nunziato model
|
|
- Rolf Pitts
- 5 years ago
- Views:
Transcription
1 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 (University Paris 6/LJLL) Strasbourg, April, /35
2 THE 7-EQUATION MODEL In 1D and dimensionless form, the model reads α k α k + u I t x =Θ(p k p l ), t (α k k)+ x (α k ku k )=0, t (α k ku k )+ x (α k( ku 2 k + p k)) p I α k x =α k kg Λ(u k u l ), t (α k ke k )+ x (α k( ke k + p k )u k ) p I u I α k x =α k kgu k p I Θ(p k p l ) u I Λ(u k u l ) withα 1 +α 2 = 1 u I, p I : interfacial velocity and pressure (to be precised) We note that the system is of the general nonconservative form U t + x F(U)+B(U) U x = S(U) 2/35
3 THE 7-EQUATION MODEL In 1D and dimensionless form, the model reads α k α k + u I t x =Θ(p k p l ), t (α k k)+ x (α k ku k )=0, t (α k ku k )+ x (α k( ku 2 k + p k)) p I α k x =α k kg Λ(u k u l ), t (α k ke k )+ x (α k( ke k + p k )u k ) p I u I α k x =α k ku k g p I Θ(p k p l ) u I Λ(u k u l ) We assume that the drag force and pressure relaxation coefficients are given by Θ= θ(u) ǫ 2 for a small parameterǫ. Then we have Λ= λ(u) ǫ 2 u 1 u 2 p 2 p 1 =O(ǫ 2 ), u 2 u 1 =O(ǫ) 3/35
4 ASYMPTOTIC ANALYSIS Following the Chapman-Enskog method, assume that { pr = p 1 p 2 = 0+ǫp 1 r+o(ǫ 2 ) u r = u 1 u 2 = 0+ǫu 1 r +O(ǫ2 ) and set ρ=α 1 ρ 1 +α 2 ρ 2 ρu=α 1 ρ 1 u 1 +α 2 ρ 2 u 2 ρe=α 1 ρ 1 e 1 +α 2 ρ 2 e 2 ρy=α 2 ρ 2 A first-order approximation w.r.t.ǫ of the 7-equation model is given by the following differential drift-flux model t ρ+ x ρu=0 t ρy+ x (ρyu+ρy(1 Y)u r )=0 t ρu+ x (ρu 2 + p+ρy(1 Y)u 2 r )=ρg t ρe+ x (ρeu+pu+ρy(1 Y)u 2 r u)=ρgu with the following Darcy-like differential relation for the relative velocity u r u r =ǫ 2ρY(ρ ρy) ( 1 1 ) xp λ(u) ρ 1 ρ 2 ρ 4/35
5 ASYMPTOTIC ANALYSIS The solutions of the 7-equation model behave asǫ 0 like the ones of the drift-flux model with u r u r =ǫ 2ρY(ρ ρy) ( 1 1 ) xp λ(u) ρ 1 ρ 2 ρ Note that in the case of permanent flows (i.e. long-time limit solutions in the Lagrangian frame) x p= ρg so that the drift law reads u r u r = ǫ 2ρY(ρ ρy) ( 1 1 )g λ(u) ρ 1 ρ 2 See Ambroso-Chalons-Coquel-Galié-Godlewski-Raviart-Seguin See Dellacherie See Guillard-Murrone 5/35
6 MATHEMATICAL PROPERTIES Eigenvalues of the Jacobian matrix F (U)+B(U) are always real and given by where c k is the sound speed of phase k u I u k u k ± c k k=1, 2 The system is weakly hyperbolic : resonance occurs if u I = u k ± c k Eigenvalues u k ± c k are GNL and eigenvalues u k are LD Eigenvalue u I if LD provided that χα 1 1 u I =βu 1 + (1 β)u 2, β= χα 1 1+ (1 χ)α 2 2 whereχis a constant such thatχ [0, 1], for instanceχ=0, 1, 1 2 See for instance Gallouet-Hérard-Seguin 6/35
7 MATHEMATICAL PROPERTIES Regarding the interfacial pressure, we set p I =µp 1 + (1 µ)p 2, µ=µ(u) [0, 1] The choice ofµis related to entropy considerations Note s k = s k ( k,ε k ) the specific entropy and T k the temperature of the phase k Letφ=φ(s) be a decreasing function. We set η=η(u)= 2 α k kφ(s k ), q=q(u)= k=1 The following entropy inequality holds true if η t + q x 0 (1 β)t 2 µ=µ(β)= βt 1 + (1 β)t 2 The pair (η, q) is a mathematical entropy pair ifφis convex 2 α k kφ(s k )u k. k=1 7/35
8 NUMERICAL DIFFICULTIES The Riemann problem is difficult to solve due to the large number of waves (7) and intermediate states (6) the presence of the nonconservative products 1 p I x α k and p I u I x α k the possibly strong nonlinearities of the pressure laws p k t u 1 u I u 2 u k c k u k + c k U L U R x FIG.: General structure of a Riemann problem 1 These are defined thanks to the Riemann invariants 8/35
9 OBJECTIVES Our objective is to propose an Approximate Riemann Solver : with the simplest possible structure (with a predefined waves position, without underlying nonlinear algorithms) with an accurate treatment of the nonconservative products able to deal with any equation of state and any choice (u I, p I ) leading to an asymptotic-preserving Godunov-type scheme (upwinding the sources...) 9/35
10 HOW TO GET THE ASYMPTOTIC-PRESERVING PROPERTY? How to get the asymptotic-preserving (AP) property? In a recent work by Chalons, Coquel, Godlewski, Raviart and Seguin (M3AS, 2010), devoted to the Euler equations with friction and gravity t ρ+ x ρu=0 t ρu+ x (ρu 2 + p)=ρ(g 1 ǫ φ(u)) t ρe+ x (ρeu+pu)=ρ(gu 1 ǫ ψ(u)), we proved that a good stretagy to get the AP property was to consider approximate Riemann solvers that are consistent with the integral form of the full model, i.e. including both the convective part and the sources. 10/35
11 CONSISTENCY OF SIMPLE APPROXIMATE RIEMANN SOLVERS We consider the nonconservative system with sources U t + x F(U)+B(U) U x = S(U) A simple Approximate Riemann Solver has the following form x U 1 = U L, t<σ 1, W ( x t ; U L, U R )= U k, σ k 1 < x t<σ k, k= 2,.., m, x U m+1 = U R, t>σ m. whereσ k =σ k (U L, U R ), 1 k m and the intermediate states may depend on =( x, t) Following Gallice (2002), we remind how to design a consistent and simple Approximate Riemann Solver with initial data { UL if x<0, U(x, 0)=U 0 (x)= U R if x>0 11/35
12 CONSISTENCY IN THE INTEGRAL SENSE Definition. The simple Approximate Riemann solver is said to be consistent with the integral form if there exists a matrix B (U L, U R ) and a vector S (U L, U R ) such that on the first hand lim U L, U R U 0 on the other hand B (U L, U R )=B(U), lim U L, U R U 0 F+B (U L, U R ) U x S (U L, U R )= S (U L, U R )=S(U) m σ k (U k+1 U k ) under the usual CFL condition and where we have set U=U R U L and F=F(U R ) F(U L ). k=1 In the case of a usual system of conservation laws, we recover the classical relation m F= σ k (U k+1 U k ) k=1 12/35
13 SO... Our ojective is to propose an Approximate Riemann Solver for α k α k + u I t x = 0, t (α k k)+ x (α k ku k )=0, t (α k ku k )+ x (α k( ku 2 k + p k)) p I α k x =α k kg Λ(u k u l ), t (α k ke k )+ x (α k( ke k + p k )u k ) p I u I α k x =α k kgu k u I Λ(u k u l ) with the property of consistency with the integral form of the full model 13/35
14 HOWTODEALWITHGENERALEQUATIONSOFSTATE? How to deal with general equations of state? Like in a HLLC-approach, we linearize the system replacingα k p k withπ k : α k α k + u I t x = 0 α k k t + α k ku k x = 0 α k ku k t α k ke k t + x (α k ku 2 k +Π k) p I α k x =α k kg K(u k u l ) + x ((α k ke k +Π k )u k ) p I u I α k x =α k kgu k u I K(u k u l ) Note that if we setπ k,l =α k,l p k,l andπ k,r =α k,r p k,r, the consistency relations F+B (U L, U R ) U x S (U L, U R )= m σ k (U k+1 U k ) k=1 are identical the system under consideration being linearized or not 14/35
15 HOW TO GET THE SIMPLEST POSSIBLE STRUCTURE? How to get the simplest possible structure? We focus on the following structure which consists in : "averaging" the waves u 1, u 2 and u I into a single coupling-wave ũ I approximating the characteristic speeds u k ± c k by u k ± a k /ρ k t u I u k,l a k,l /ρ k,l u k,r + a k,r /ρ k,r U L U R x FIG.: General structure of an approximate Riemann solution 15/35
16 THEN... We determine the intermediate states by : imposing Rankine-Hugoniot jump conditions of the linearized system across the discontinuities of the Approximate Riemann Solver imposing the validity of the consistency relations We then get : explicit formulas for the intermediate states an accurate treatment of the non conservative products (the stationary coupling waves are exactly computed) the propertiesα k (0, 1) and (α k ρ k )>0 the conservativity of (α k ρ k ), (α 1 ρ 1 u 1 +α 2 ρ 2 u 2 ) and (α 1 ρ 1 E 1 +α 2 ρ 2 E 2 ) 16/35
17 NUMERICAL RESULTS We considered several 1D test cases to assess the strategy : Riemann problems : isolated coupling waves general Riemann problems including comparisons with other schemes available in the literature (HLL-scheme proposed by Saurel-Abgrall, the scheme proposed by Andrianov-Saurel-Warnacke, VFRoe-scheme by Gallouët-Hérard-Seguin, Godunov-scheme by Schwendeman-Wahle-Kapila) More typical two-phase flows : Ransom faucet test case sedimentation desequilibrium in velocities bubbly column and asymptotic behavior 17/35
18 TEST 1A Test 1a : isolated coupling wave (Gallouët-Hérard-Seguin) We choose and u I = α 1ρ 1 u 1 +α 2 ρ 2 u 2, χ=0.5 α 1 ρ 1 +α 2 ρ 2 α 1,L = 0.9, (ρ 1, u 1, p 1 ) L = (1, 100, 10 5 ) (ρ 2, u 2, p 2 ) L = (1, 100, 10 5 ) α 1,R = 0.5, (ρ 1, u 1, p 1 ) R = (0.125, 100, 10 5 ) (ρ 2, u 2, p 2 ) R = (0.125, 100, 10 5 ) 18/35
19 TEST 1A 1 Godunov-type VFRoe exact Godunov-type VFRoe exact 1 Godunov-type VFRoe exact FIG.: x versusα 1,ρ 1,ρ 2 19/35
20 TEST 1A 110 Godunov-type VFRoe exact 110 Godunov-type VFRoe exact Godunov-type VFRoe exact Godunov-type VFRoe exact FIG.: x versus u 1, u 2, p 1, p 2 20/35
21 TEST 1B Test 1b : isolated coupling wave (Andrianov-Saurel-Warnacke) We chooseχ=1so that and u I = u 1 α 1,L = 0.8, (ρ 1, u 1, p 1 ) L = (2, 0.3, 5) (ρ 2, u 2, p 2 ) L = (1, 2, 1) α 1,R = 0.3, (ρ 1, u 1, p 1 ) R = (2, 0.3, ) (ρ 2, u 2, p 2 ) R = (0.1941, , 0.1) 21/35
22 TEST 1B pts exact pts exact pts exact FIG.: x versusα 1,ρ 1, u 1 22/35
23 TEST 2 Test 2 : a more general Riemann problem taken from the paper by Schwendeman-Wahle-Kapila We chooseχ=1so that and u I = u 1 α 1,L = 0.8, (ρ 1, u 1, p 1 ) L = (1, 0, 1) (ρ 2, u 2, p 2 ) L = (0.2, 0, 0.3) α 1,R = 0.3, (ρ 1, u 1, p 1 ) R = (1, 0, 1) (ρ 2, u 2, p 2 ) L = (1, 0, 1). We compare our results with the ones given by the schemes proposed by Saurel-Abgrall, Andrianov-Saurel-Warnacke, Schwendeman-Wahle-Kapila 23/35
24 TEST pts exact pts exact pts exact pts exact FIG.:α 1 (x), p 1 (α 1 ),ρ 1 (x),ρ 2 (x) 24/35
25 TEST (a) 1.1 (b) α p x α 1.06 (c) 1 (d) ρ ρ x x FIG.:α 1 (x), p 1 (α 1 ),ρ 1 (x),ρ 2 (x) (by courtesy of Schwendeman-Wahle-Kapila) caption : G HLL (blue), G ASW (green) et G 1 (red) 25/35
26 TEST 2 Riemann invariants associated with u I are K 0 = u I = u 1 K 1 = (1 α 1 )ρ 2 (u 2 u I ) K 2 = (1 α 1 )ρ 2 (u 2 u I ) 2 +α 1 p 1 + (1 α 1 )p 2 K 3 = γp 2 (γ 1)ρ (u 2 u I ) 2 K 4 = p 2 ρ γ 2 Forα 1 varying fromα 1,L toα 1,R, these equations provide a parametrization of the coupling wave u I. This is compared with the "boundary layer" induced by the diffusion of the numerical solution. Constants K i, i=0,..., 4 are calculated using the exact value of the right state associated withα 1 =α 1,R similarly as in the paper Schwendeman-Wahle-Kapila. We observe a very good agreement 26/35
27 TEST 3 : INFLUENCE OF χ chi=0 chi=0.5 chi= chi=0 chi=0.5 chi=1 1.2 chi=0 chi=0.5 chi= FIG.: x versusα 1,ρ 1,ρ 2 27/35
28 TEST 3 : INFLUENCE OF χ chi=0 chi=0.5 chi=1 0 chi=0 chi=0.5 chi= chi=0 chi=0.5 chi=1 1.1 chi=0 chi=0.5 chi= FIG.: x versus u 1, u 2, p 1, p 2 28/35
29 INFLUENCE OF THE SOURCES External forces : g=9.81 Pressure relaxation : Drag force : θ(u)= 1 τ p α 1 (1 α 1 ) p 1 + p 2 (p 1 p 2 ) λ(u)= 1 8 C Da int ρ 1 u 1 u 2 with C D = 0.5, a int = 3(1 α 1), r b = r b 29/35
30 TEST 4 : RANSOM FAUCET TEST CASE VFRoe Godunov-type 15 VFRoe Godunov-type VFRoe Godunov-type VFRoe Godunov-type FIG.: Leftα 2, right u 1. Top 2000 points, bottom points 30/35
31 TEST 5 :SEDIMENTATION t=0.2 t=0.4 t=0.6 t=0.8 t=1.0 t=1.2 t= FIG.:α 1 withτ p = s and 400 cells 31/35
32 TEST 5 :SEDIMENTATION t=0.2 t=0.4 t=0.6 t=0.8 t=1.0 t=1.2 t= /35
33 TEST 7 :BUBBLYCOLUMN Evolution of a two-phase flow in a vartical bubbly column of 1m-length Boundary conditions : inlet, we impose u 1 = 5m/s, u 2 = 15m/s,α 1 = 0.97,ρ 1 = 10 3 kg/m 3,ρ 2 = 1kg/m 3 and the mixture pressure p= Pa outlet, we impose the mixture pressure p= Pa Asymptotic behavior : x p= ρg, p=α 1 p 1 +α 2 p 2, ρ=α 1 ρ 1 +α 2 ρ 2, et u 2 r= α 1(1 ρ 2 ρ 1 )g 33/35
34 TEST 7 :BUBBLYCOLUMN 30 u_r^2 asymptotic relative velocity 0.2 u_r^2 asymptotic relative velocity e+07 p ^7 - rho g (x-1) e e e e e e e e FIG.: Stationary solutions of Test 7. x versus u 2 r, u 2 r (close-up), p 34/35
35 OPEN PROBLEMS Open problems are : the positivity of the internal energies the AP property (seems to be ok from the numerical experiments) an entropy inequality (related to the numerical choice of p I?) 35/35
Large 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 informationModel 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 informationHervé Guillard, INRIA Projet Smash, B.P. 93, Sophia-Antipolis Cedex, France,
TRAVELING WAVE ANALYSIS OF TWO-PHASE DISSIPATIVE MODELS Hervé Guillard, INRIA Projet Smash, B.P. 93, 06902 Sophia-Antipolis Cedex, France, Herve.Guillard@sophia.inria.fr Joint work with : Mathieu Labois,
More informationAnnalisa Ambroso 1, Christophe Chalons 2, Frédéric Coquel 3 and Thomas. Introduction
Mathematical Modelling and Numerical Analysis Modélisation Mathématique et Analyse Numérique Will be set by the publisher RELAXATION AND NUMERICAL APPROXIMATION OF A TWO-FLUID TWO-PRESSURE DIPHASIC MODEL
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 informationOn the Modeling and Simulation of a Laser-Induced Caviation Bubble
1 On the Modeling and Simulation of a Laser-Induced Caviation Bubble Ali Zein Supervisor: Prof. Dr. Gerald Warnecke Co-supervisor: Dr. Maren Hantke Institut für Analysis und Numerik Otto-von-Guericke-Universität
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 informationWELL-BALANCED ROE-TYPE NUMERICAL SCHEME FOR A MODEL OF TWO-PHASE COMPRESSIBLE FLOWS
J. Korean Math. Soc. 51 2014), No. 1, pp. 163 187 http://dx.doi.org/10.4134/jkms.2014.51.1.163 WELL-BALANCED ROE-TYPE NUMERICAL SCHEME FOR A MODEL OF TWO-PHASE COMPRESSIBLE FLOWS Mai Duc Thanh Abstract.
More informationModeling of two-phase flow > Direct steam generation in solar thermal power plants
Modeling of two-phase flow > Direct steam generation in solar thermal power plants Pascal Richter RWTH Aachen University EGRIN Pirirac-sur-Mer June 3, 2015 Novatec Solar GmbH Direct steam generation Solar
More informationNumerical investigation of shock wave dense particles cloud interaction
25 th ICDERS August 2 7, 215 Leeds, UK Numerical investigation of shock wave dense particles cloud interaction Pavel S. Utkin Moscow Institute of Physics and Technology Dolgoprudny, Moscow region, Russia
More informationTHE numerical simulation of the creation and evolution
Proceedings of the World Congress on Engineering Vol III WCE, July 4-6,, London, U.K. Numerical Simulation of Compressible Two-phase Flows Using an Eulerian Type Reduced Model A. Ballil, Member, IAENG,
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 informationA SPLITTING METHOD FOR THE ISENTROPIC BAER-NUNZIATO TWO-PHASE FLOW MODEL
ESAIM: PROCEEDINGS, Vol.?, 01, 1-10 Editors: Will be set by the publisher A SPLITTING METHOD FOR THE ISENTROPIC BAER-NUNZIATO TWO-PHASE FLOW MODEL Frédéric Coquel 1, Jean-Marc Hérard and Khaled Saleh,
More informationModélisation et simulation d écoulements diphasiques avec contrainte granulaire
Modélisation et simulation d écoulements diphasiques avec contrainte granulaire T. Gallouët, P. Helluy, J.-M. Hérard, J. Nussbaum LATP Marseille IRMA Strasbourg EDF Chatou ISL Saint-Louis Context Flow
More information0.3.4 Burgers Equation and Nonlinear Wave
16 CONTENTS Solution to step (discontinuity) initial condition u(x, 0) = ul if X < 0 u r if X > 0, (80) u(x, t) = u L + (u L u R ) ( 1 1 π X 4νt e Y 2 dy ) (81) 0.3.4 Burgers Equation and Nonlinear Wave
More informationA Numerical Study of Compressible Two-Phase Flows Shock and Expansion Tube Problems
A Numerical Study of Compressible Two-Phase Flows Shock and Expansion Tube Problems Dia Zeidan,a) and Eric Goncalves 2 School of Basic Sciences and Humanities, German Jordanian University, Amman, Jordan
More informationFinite Volume simulation of cavitating flows
Finite Volume simulation of cavitating flows Thomas Barberon, Philippe Helluy ISITV, Laboratoire MNC, BP 56, 83162 La Valette, France Abstract We propose a numerical method adapted to the modelling of
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 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 informationA Study of Detonation Evolution and Structure for a Model of Compressible Two-Phase Reactive Flow
A Study of Detonation Evolution and Structure for a Model of Compressible Two-Phase Reactive Flow D. W. Schwendeman, C. W. Wahle and A. K. Kapila Department of Mathematical Sciences, Rensselaer Polytechnic
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 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 informationComputations of non-reacting and reacting two-fluid interfaces
Computations of non-reacting and reacting two-fluid interfaces Kunkun Tang Alberto Beccantini Christophe Corre Lab. of Thermal Hydraulics & Fluid Mechanics (LATF), CEA Saclay Lab. of Geophysical & Industrial
More informationCouplage de modèles en thermohydraulique. Troisième rapport de synthèse
Couplage de modèles en thermohydraulique. Troisième rapport de synthèse Annalisa Ambroso, Benjamin Boutin, Christophe Chalons, Frédéric Coquel, Thomas Galié Edwige Godlewski, Frédéric Lagoutière, Pierre-Arnaud
More informationRecent advances in numerical methods for compressible two-phase flow with heat & mass transfers
Recent advances in numerical methods for compressible two-phase flow with heat & mass transfers Keh-Ming Shyue Institute of Applied Mathematical Sciences National Taiwan University Joint work with Marica
More informationHyperbolic Systems of Conservation Laws. in One Space Dimension. I - Basic concepts. Alberto Bressan. Department of Mathematics, Penn State University
Hyperbolic Systems of Conservation Laws in One Space Dimension I - Basic concepts Alberto Bressan Department of Mathematics, Penn State University http://www.math.psu.edu/bressan/ 1 The Scalar Conservation
More informationMathematical modeling of impact of two metal plates using two-fluid approach
Journal of Physics: Conference Series PAPER OPEN ACCESS Mathematical modeling of impact of two metal plates using two-fluid approach To cite this article: P S Utkin and S V Fortova 018 J. Phys.: Conf.
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 informationCommissariat à l Energie Atomique - Saclay Department of Nuclear Energy Fluid Modeling and Simulation. MIT, November 3rd, 2008
ON THE NUMERICAL SIMULATION OF GENERAL MULTI-FIELD TWO-PHASE FLOW MODELS Anela KUMBARO Commissariat à l Energie Atomique - Saclay Department of Nuclear Energy Fluid Modeling and Simulation MIT, November
More informationEFFECT OF LIQUID PHASE COMPRESSIBILITY ON MODELING OF GAS-LIQUID TWO-PHASE FLOWS USING TWO-FLUID MODEL
EFFECT OF LIQUID PHASE COMPRESSIBILITY ON MODELING OF GAS-LIQUID TWO-PHASE FLOWS USING TWO-FLUID MODEL Vahid SHOKRI 1*,Kazem ESMAEILI 2 1,2 Department of Mechanical Engineering, Sari Branch, Islamic Azad
More informationA Hybrid Algorithm for the Baer-Nunziato Model Using the Riemann Invariants
A Hybrid Algorithm for the Baer-Nunziato Model Using the Riemann Invariants Smadar Karni Gerardo Hernández-Dueñas Abstract The paper considers the Baer-Nunziato model for two-phase flow in porous media,
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 informationOn 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 informationCapSel Roe Roe solver.
CapSel Roe - 01 Roe solver keppens@rijnh.nl modern high resolution, shock-capturing schemes for Euler capitalize on known solution of the Riemann problem originally developed by Godunov always use conservative
More informationEulerian Two-Phase Flow CFD Simulation Using a Compressible and Equilibrium Eight- Equation Model. Y. Wang 1 and R. D. Reitz
ILASS Americas 27th Annual Conference on Liquid Atomization and Spray Systems, Raleigh, NC, May 2015 Eulerian Two-Phase Flow CFD Simulation Using a Compressible and Equilibrium Eight- Equation Model Y.
More informationA Very Brief Introduction to Conservation Laws
A Very Brief Introduction to Wen Shen Department of Mathematics, Penn State University Summer REU Tutorial, May 2013 Summer REU Tutorial, May 2013 1 / The derivation of conservation laws A conservation
More informationCapSel Euler The Euler equations. conservation laws for 1D dynamics of compressible gas. = 0 m t + (m v + p) x
CapSel Euler - 01 The Euler equations keppens@rijnh.nl conservation laws for 1D dynamics of compressible gas ρ t + (ρ v) x = 0 m t + (m v + p) x = 0 e t + (e v + p v) x = 0 vector of conserved quantities
More informationThe RAMSES code and related techniques 2- MHD solvers
The RAMSES code and related techniques 2- MHD solvers Outline - The ideal MHD equations - Godunov method for 1D MHD equations - Ideal MHD in multiple dimensions - Cell-centered variables: divergence B
More informationThe Riemann problem and a high-resolution Godunov method for a model of compressible two-phase flow
The Riemann problem and a high-resolution Godunov method for a model of compressible two-phase flow D.W. Schwendeman, C.W. Wahle 2, and A.K. Kapila Department of Mathematical Sciences, Rensselaer Polytechnic
More informationx a(x) H(U), , H(U) =
EXACT RIEMANN SOLUTIONS TO COMPRESSIBLE EULER EQUATIONS IN DUCTS WITH DISCONTINUOUS CROSS SECTION EE HAN, MAREN HANTKE, AND GERALD WARNECKE Abstract. We determine completely the exact Riemann solutions
More informationA 5-wave relaxation solver for the shallow water MHD system
A 5-wave relaxation solver for the shallow water MHD system François Bouchut, Xavier Lhébrard Abstract The shallow water magnetohydrodynamic system describes the thin layer evolution of the solar tachocline.
More informationComparison of Approximate Riemann Solvers
Comparison of Approximate Riemann Solvers Charlotte Kong May 0 Department of Mathematics University of Reading Supervisor: Dr P Sweby A dissertation submitted in partial fulfilment of the requirement for
More informationRelaxation approximation of the Euler equations
Relaxation approximation of the Euler equations Christophe Chalons, Jean-François Coulombel Laboratoire Jacques-Louis Lions and Université Paris 7 Boîte courrier 87, 2 place Jussieu, 75252 PARIS CEDEX
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 informationNumerical methods for FCI Part IV Multi-temperature fluid models Hele-Shaw models
B. Després+ X. Blanc LJLL-Paris VI+CEA Thanks to same collegues as before plus C. Buet, H. Egly and R. Sentis methods for FCI Part IV Multi-temperature fluid s s B. Després+ X. Blanc LJLL-Paris VI+CEA
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 informationPositivity-preserving cell-centered Lagrangian schemes for multi-material compressible flows: From first-order to high-orders
Positivity-preserving cell-centered Lagrangian schemes for multi-material compressible flows: From first-order to high-orders François Vilar a,, Chi-Wang Shu a, Pierre-Henri Maire b a Division of Applied
More informationSCHEDULE. Gonca Aki (Berlin) An incompressible diffuse flow with phase transition
SCHEDULE TUESDAY AFTERNOON 2:00pm Sergey Gavrilyuk (Marseille) Diffuse interface model for compressible fluid - compressible elastic-plastic solid interaction 2:45pm Gonca Aki (Berlin) An incompressible
More informationc 1999 Society for Industrial and Applied Mathematics
SIAM J. SCI. COMPUT. Vol. 21, No. 3, pp. 1115 1145 c 1999 Society for Industrial and Applied Mathematics A SIMPLE METHOD FOR COMPRESSIBLE MULTIFLUID FLOWS RICHARD SAUREL AND RÉMI ABGRALL Abstract. A simple
More informationInternational Engineering Research Journal
Special Edition PGCON-MECH-7 Development of high resolution methods for solving D Euler equation Ms.Dipti A. Bendale, Dr.Prof. Jayant H. Bhangale and Dr.Prof. Milind P. Ray ϯ Mechanical Department, SavitribaiPhule
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 informationHyperbolic Systems of Conservation Laws. I - Basic Concepts
Hyperbolic Systems of Conservation Laws I - Basic Concepts Alberto Bressan Mathematics Department, Penn State University Alberto Bressan (Penn State) Hyperbolic Systems of Conservation Laws 1 / 27 The
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 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 informationA Fluid-Mixture Type Algorithm for Compressible Multicomponent Flow with Mie Grüneisen Equation of State
Journal of Computational Physics 171, 678 707 (2001 doi:10.1006/jcph.2001.6801, available online at http://www.idealibrary.com on A Fluid-Mixture Type Algorithm for Compressible Multicomponent Flow with
More informationEntropy stable schemes for degenerate convection-diffusion equations
Entropy stable schemes for degenerate convection-diffusion equations Silvia Jerez 1 Carlos Parés 2 ModCompShock, Paris 6-8 Decmber 216 1 CIMAT, Guanajuato, Mexico. 2 University of Malaga, Málaga, Spain.
More informationWEAKLY IMPLICIT NUMERICAL SCHEMES FOR A TWO-FLUID MODEL
SIAM J. SCI. COMPUT. Vol. 6, No. 5, pp. 1449 1484 c 005 Society for Industrial and Applied Mathematics WEAKLY IMPLICIT NUMERICAL SCHEMES FOR A TWO-FLUID MODEL STEINAR EVJE AND TORE FLÅTTEN Abstract. The
More informationCOMPUTATIONAL METHODS AND ALGORITHMS Vol II - Computational Methods for Compressible Flow Problems - Remi Abgrall
COMPUTATIONAL METHODS FOR COMPRESSIBLE FLOW PROBLEMS Rémi Abgrall Université Bordeaux I, Talence Cedex, France Keywords: computational methods, numerical schemes, 1- D problems, multidimensional problems,
More informationA Simple Model for Cavitation with Non-condensable Gases
A Simple Model for Cavitation with Non-condensable Gases Mathieu Bachmann, Siegfried Müller, Philippe Helluy, Hélène Mathis To cite this version: Mathieu Bachmann, Siegfried Müller, Philippe Helluy, Hélène
More informationEntropy-based moment closure for kinetic equations: Riemann problem and invariant regions
Entropy-based moment closure for kinetic equations: Riemann problem and invariant regions Jean-François Coulombel and Thierry Goudon CNRS & Université Lille, Laboratoire Paul Painlevé, UMR CNRS 854 Cité
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 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 informationVarious lecture notes for
Various lecture notes for 18311. R. R. Rosales (MIT, Math. Dept., 2-337) April 12, 2013 Abstract Notes, both complete and/or incomplete, for MIT s 18.311 (Principles of Applied Mathematics). These notes
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 informationA Finite Volume Code for 1D Gas Dynamics
A Finite Volume Code for 1D Gas Dynamics Michael Lavell Department of Applied Mathematics and Statistics 1 Introduction A finite volume code is constructed to solve conservative systems, such as Euler
More informationSHOCK WAVES FOR RADIATIVE HYPERBOLIC ELLIPTIC SYSTEMS
SHOCK WAVES FOR RADIATIVE HYPERBOLIC ELLIPTIC SYSTEMS CORRADO LATTANZIO, CORRADO MASCIA, AND DENIS SERRE Abstract. The present paper deals with the following hyperbolic elliptic coupled system, modelling
More informationForce analysis of underwater object with supercavitation evolution
Indian Journal of Geo-Marine Sciences Vol. 42(8), December 2013, pp. 957-963 Force analysis of underwater object with supercavitation evolution B C Khoo 1,2,3* & J G Zheng 1,3 1 Department of Mechanical
More informationLecture Notes on Hyperbolic Conservation Laws
Lecture Notes on Hyperbolic Conservation Laws Alberto Bressan Department of Mathematics, Penn State University, University Park, Pa. 16802, USA. bressan@math.psu.edu May 21, 2009 Abstract These notes provide
More informationLecture 5.7 Compressible Euler Equations
Lecture 5.7 Compressible Euler Equations Nomenclature Density u, v, w Velocity components p E t H u, v, w e S=c v ln p - c M Pressure Total energy/unit volume Total enthalpy Conserved variables Internal
More informationNon-linear Wave Propagation and Non-Equilibrium Thermodynamics - Part 3
Non-linear Wave Propagation and Non-Equilibrium Thermodynamics - Part 3 Tommaso Ruggeri Department of Mathematics and Research Center of Applied Mathematics University of Bologna January 21, 2017 ommaso
More informationStochastic representation of random positive-definite tensor-valued properties: application to 3D anisotropic permeability random fields
Sous la co-tutelle de : LABORATOIRE DE MODÉLISATION ET SIMULATION MULTI ÉCHELLE CNRS UPEC UNIVERSITÉ PARIS-EST CRÉTEIL UPEM UNIVERSITÉ PARIS-EST MARNE-LA-VALLÉE Stochastic representation of random positive-definite
More informationChapter 1. Introduction
Chapter 1 Introduction Many astrophysical scenarios are modeled using the field equations of fluid dynamics. Fluids are generally challenging systems to describe analytically, as they form a nonlinear
More informationA General Technique for Eliminating Spurious Oscillations in Conservative Schemes for Multiphase and Multispecies Euler Equations
A General Technique for Eliminating Spurious Oscillations in Conservative Schemes for Multiphase and Multispecies Euler Equations Ronald P. Fedkiw Xu-Dong Liu Stanley Osher September, 2 Abstract Standard
More informationThe Riemann problem. The Riemann problem Rarefaction waves and shock waves
The Riemann problem Rarefaction waves and shock waves 1. An illuminating example A Heaviside function as initial datum Solving the Riemann problem for the Hopf equation consists in describing the solutions
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 informationCOMPRESSIBLE TWO-PHASE FLOWS BY CENTRAL AND UPWIND SCHEMES
ESAIM: M2AN Vol. 38, N o 3, 24, pp. 477 493 DOI:.5/m2an:2424 ESAIM: Mathematical Modelling and Numerical Analysis COMPRESSIBLE TWO-PHASE FLOWS BY CENTRAL AND UPWIND SCHEMES Smadar Karni, Eduard Kirr 2,
More informationNumerical simulations of gas-particle flows with combustion
Numerical simulations of gas-particle flows with combustion Julien NUSSBAUM French-german Research Institut of Saint-Louis 5, rue du Général Cassagnou BP 34, 68301 SAINT LOUIS CEDEX Philippe HELLUY Laboratoire
More informationAn Added-Mass Partition Algorithm for Fluid-Structure Interactions of Compressible Fluids and Nonlinear Solids
An Added-Mass Partition Algorithm for Fluid-Structure Interactions of Compressible Fluids and Nonlinear Solids J. W. Banks a,,2, W. D. Henshaw a,,3,, A. K. Kapila a,4, D. W. Schwendeman a,,3,4 a Department
More informationSolving the Payne-Whitham traffic flow model as a hyperbolic system of conservation laws with relaxation
Solving the Payne-Whitham traffic flow model as a hyperbolic system of conservation laws with relaxation W.L. Jin and H.M. Zhang August 3 Abstract: In this paper we study the Payne-Whitham (PW) model as
More informationA high order adaptive finite element method for solving nonlinear hyperbolic conservation laws
A high order adaptive finite element method for solving nonlinear hyperbolic conservation laws Zhengfu Xu, Jinchao Xu and Chi-Wang Shu 0th April 010 Abstract In this note, we apply the h-adaptive streamline
More informationInfo. No lecture on Thursday in a week (March 17) PSet back tonight
Lecture 0 8.086 Info No lecture on Thursday in a week (March 7) PSet back tonight Nonlinear transport & conservation laws What if transport becomes nonlinear? Remember: Nonlinear transport A first attempt
More informationNUMERICAL SOLUTION OF HYPERBOLIC PARTIAL DIFFERENTIAL EQUATIONS
NUMERICAL SOLUTION OF HYPERBOLIC PARTIAL DIFFERENTIAL EQUATIONS JOHN A. TRANGENSTEIN Department of Mathematics, Duke University Durham, NC 27708-0320 Ш CAMBRIDGE ЩР UNIVERSITY PRESS Contents 1 Introduction
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 informationComputational Astrophysics
16 th Chris Engelbrecht Summer School, January 2005 3: 1 Computational Astrophysics Lecture 3: Magnetic fields Paul Ricker University of Illinois at Urbana-Champaign National Center for Supercomputing
More informationA wave propagation method for compressible multicomponent problems on quadrilateral grids
A wave propagation method for compressible multicomponent problems on quadrilateral grids K.-M. Shyue Department of Mathematics, National Taiwan University, Taipei 16, Taiwan Abstract We describe a simple
More informationNotes: Outline. Shock formation. Notes: Notes: Shocks in traffic flow
Outline Scalar nonlinear conservation laws Traffic flow Shocks and rarefaction waves Burgers equation Rankine-Hugoniot conditions Importance of conservation form Weak solutions Reading: Chapter, 2 R.J.
More informationAMath 574 February 11, 2011
AMath 574 February 11, 2011 Today: Entropy conditions and functions Lax-Wendroff theorem Wednesday February 23: Nonlinear systems Reading: Chapter 13 R.J. LeVeque, University of Washington AMath 574, February
More informationNumerical Methods for Conservation Laws WPI, January 2006 C. Ringhofer C2 b 2
Numerical Methods for Conservation Laws WPI, January 2006 C. Ringhofer ringhofer@asu.edu, C2 b 2 2 h2 x u http://math.la.asu.edu/ chris Last update: Jan 24, 2006 1 LITERATURE 1. Numerical Methods for Conservation
More informationClosure conditions for non-equilibrium multi-component models
Closure conditions for non-equilibrium multi-component models Siegfried Müller, Maren Hante 2 and Pascal Richter 3 Preprint No. 44 December 204 Key words: multi-component flows, entropy, relaxation, phase
More informationNONCLASSICAL SHOCK WAVES OF CONSERVATION LAWS: FLUX FUNCTION HAVING TWO INFLECTION POINTS
Electronic Journal of Differential Equations, Vol. 2006(2006), No. 149, pp. 1 18. ISSN: 1072-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu ftp ejde.math.txstate.edu (login: ftp) NONCLASSICAL
More informationACCURATE NUMERICAL SCHEMES FOR APPROXIMATING INITIAL-BOUNDARY VALUE PROBLEMS FOR SYSTEMS OF CONSERVATION LAWS.
ACCURATE NUMERICAL SCHEMES FOR APPROXIMATING INITIAL-BOUNDARY VALUE PROBLEMS FOR SYSTEMS OF CONSERVATION LAWS. SIDDHARTHA MISHRA AND LAURA V. SPINOLO Abstract. Solutions of initial-boundary value problems
More informationSome approximate Godunov schemes to compute. shallow-water equations with topography. Universite de Provence. 39 rue Joliot Curie.
Some approximate Godunov schemes to compute shallow-water equations with topography Thierry Gallouet, Jean-Marc Herard ;, Nicolas Seguin ; Laboratoire d'analyse Topologie et Probabilites - UMR 663 Centre
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 informationSteinar Evje 1 and Tore Flåtten 1
ESAIM: M2AN Vol. 39, N o 2, 25, pp. 253 273 DOI: 1.151/m2an:2511 ESAIM: Mathematical Modelling Numerical Analysis HYBRID CENTRAL-UPWIND SCHEMES FOR NUMERICAL RESOLUTION OF TWO-PHASE FLOWS Steinar Evje
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 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 informationFlux corrected transport solver for solving 1D multiphase equations for drilling applications
Computational Methods in Multiphase Flow VII 169 Flu corrected transport solver for solving 1D multiphase equations for drilling applications P. A. Varadarajan & P. S. Hammond Schlumberger Gould Research,
More informationCoupling algorithms for hyperbolic systems
Coupling algorithms for hyperbolic systems Edwige Godlewski Laboratoire Jacques-Louis Lions Université Pierre et Marie Curie (Paris) Cemracs, July 2011 Collaborations with Frédéric Coquel ( CMAP, Cnrs),
More informationK. Ambika and R. Radha
Indian J. Pure Appl. Math., 473: 501-521, September 2016 c Indian National Science Academy DOI: 10.1007/s13226-016-0200-9 RIEMANN PROBLEM IN NON-IDEAL GAS DYNAMICS K. Ambika and R. Radha School of Mathematics
More information