Computation of entropic measure-valued solutions for Euler equations
|
|
- Pamela Fields
- 5 years ago
- Views:
Transcription
1 Computation of entropic measure-valued solutions for Euler equations Eitan Tadmor Center for Scientific Computation and Mathematical Modeling (CSCAMM) Department of Mathematics, Institute for Physical Science & Technology University of Maryland Laboratoire Jacques-Louis Lions Universite Pierre et Marie Curie, June 7 6 Collaborators: U. Fjordholm R. Käppeli S. Mishra
2 Prologue. Perfect derivatives and conservative differences u }{{} u ν+ u ν uu }{{} ( ) uν+ u ν u ν but... u 3 u }{{} u 3 ν ( uν+ u ν b a b a u d = boundary terms of u ( u ν+ u ν ) = boundary terms (u ) {}}{ uu d = boundary terms of u telescoping sum {}}{ ( u ν u ν+ u ν u ν b a ) ) = quadratic boundary terms 4 (u4 ) {}}{ u 3 u d = quartic boundary terms of u no perfect deriv. {( }}{ u 3 ν u ν+ u ν uν 3 ) no cancellation
3 Prologue. Perfect derivatives and conservative differences ( ) Set u 3 ( u 4 ν+ uν 4 4 ( 3 u 4 ν+ uν 4 4 uν+ 3 u4 ν uν 4 u3 ν uν 3 uν 3 ( 3 u 4 ν+ u 4 u3 ν 4 ν uν+ 3 u3 ν uν+ 3 u4 ν uν 4 u3 ν uν 3 uν 3 u4 ν u 4 ν u 3 ν u 3 ν = 3 ( u 3 4 ν+ uν) ( 3 uν+ 4 u4 ν uν+ 3 u3 ν ) ) boundary terms of u ) = ) boundary terms of u 4 How do we come up with ( )? Can we conserve general multipliers η (u)u = η(u) Applications entropic solution for nonlinear conservation laws
4 Strong and weak entropic solutions Euler equations with pressure law p := (γ ) ( E ρ (u + v ) ) : ρ ρu ρv ρu t ρv + ρu + p ρuv + ρuv ρv + p = E u(e + p) v(e + p) Nonlinear conservation laws: u t + f(u) = ; u := (ρ, ρu, ρv, E) balance of conservative variables u over spatial variables = (,,...) Strong solutions pointwise values u(, t) Weak solutions give up the certainty in pointvalues Instead observe the cell averages u(, t) := d u(, t)+ dt [ t+ t τ=t f ( u(+ ), τ) dτ Compute cell averages and numerical flues: t+ t τ=t f + / / u(y, t)dy ( u( ) ], τ) dτ = d dt u ν(t) + f ν+ f ν =
5 Weak solutions and entropy stability η (u), u t + f(u) = entropy {}}{ η(u) t + perfect derivatives? { }}{ η (u), f(u) = η(u) is an entropy iff η (u), f(u) = F (u) Entropy stability: { = η(u) t + F (u) Euler equations specific entropy S = ln(pρ γ ): η(u) {}}{ ( ρs) t + F (u) {}}{ ( ρus) { = Computation impose the entropy stability requirement d dt η(u ν(t)) + F ν+ F { ν =
6 D radial shock tube problem u L, r Radial initial data u () = + ɛ(δu )() u R, > r [ ] [ (δu )() sin(π ) Perturbation of steady-state ɛ = : = (δv )() sin(π ) Convergence: Fi ɛ =. keep mesh refinement...ρ ɛ at t =.4: ] (a) 8 (b) 56 (c) (d) u u / L Figure: Approimate density of the perturbed Sod problem (ɛ =.) computed with TeCNO scheme at t =.4 Stability as ɛ...
7 y y Kelvin-Helmholtz instability u L, H < H Initial state of three layers u () = u R, H or H fied u L,R ; perturbed interface H j (, ω) = j 4 + ɛy j(, ω), j =, ɛ = : u L ( ) [ 4, 3 4 ] + u R( ) [, ] [ 3 4,] is an entropic steady state Perturbed state u (, ω) = u L ( ) [H,H ] + u R ( ) [,H ] [H,]: Y j (, ω) = a n (ω) cos(b n (ω) + nπ ) an = so that ɛy j ɛ n density, sample, t = density, sample, t = Lack of convergence: Fi ω, ɛ =. keep mesh refinement...ρ ɛ at t = : KH (a) initial 8 data: amplitude(b) perturbation 56 (left) and phase (c) 5 perturbation (right)(d) with4 ɛ =.
8 The notion of entropic weak solution is not enough (d > ) (e) L rates vs. number of gridpoints with ɛ = (f) L error with respect to ɛ, at a fied 4 mesh Figure: L differences in density ρ at time t = for the Kelvin-Helmholtz problem These are NOT numerical artifacts but... contemplate lack of uniqueness [De Lellis, Szekelyhidi, E. Chiodaroli, O. Kreml, et. al.] Just because we cannot prove that compressible flows with prescribed initial values eist doesn t mean that we cannot compute them [La 7] The question is what information is encoded in our computations?
9 Back to weak solutions ν,t = δ u(,t) (certainty of strongly convergent averages) Ron DiPerna Answer: Young measures and measure valued solutions Recall weak convergence: u n u and f(u n ) f What is the relation f and f(u )? f (, t) depends linearly and positively on f: hence f (, t) = ν,t, f in terms of measure ν = ν,t Complete description in terms of Young measure ν,t Luc Tartar Certainty: strong convergence with Dirac Mass ν,t = δ u (,t) f = f(u ) Entropy measure valued (EMV) solutions: u(, t) ν,t : u t + f(u) = t ν,t, id + ν,t, f = η(u) t + F (u) t ν,t, η + ν,t, F ν,t is the (weighted) probability of having the value ν,t, id at (, t) Measure valued solutions give up the certainty in a given solution Instead observe averages in configuration space
10 Young measures as quantifiers of uncertainty (a) t = (b) t =.5 (c) t = (d) t =.5 (e) t = Figure: The approimate PDF for density ρ ν,t at = (.5,.7). More about this later... how do we compute these Entropic Measure-Valued solutions (EMVs)?
11 Computing Entropy Measure Valued (EMV) solutions (Formally) arbitrarily high-order of accuracy; Entropy stability as a selection principle Essentially non-oscillatory in the presence of discontinuities; Convergence for linear systems; How to realize (the configuration space of) EMVs? Computationally efficient; The class of the entropy Conservative Non-Oscillatory (TeCNO) schemes which compute faithful samples required for the ensemble of EMVs b,c Fjordholm, Mishra, ET, Arbitrarily high-order accurate entropy stable schemes..., SINUM b Fjordholm, Kappeli, Mishra & ET, Construction of approimate entropy measure valued solutions for hyperbolic systems of conservation laws, FoCM (5) c FMT, On the computation of measure-valued solutions, Acta Numerica (6)
12 Entropy inequality: PDEs numerical approimations η (u), u t + f(u) = Semi-discrete approimations: d dt u ν(t) + [ ] f ν+ f ν η (u), f(u) =F (u) = η(u) t + F (u) { = = f ν+ = f(u ν p+,..., u ν+p ) Entropy conservative discretization given η find f s.t. d f f dt u ν+ ν(t) + ν =? d dt η(u ν(t)) + F ν+ F ν { = Does so that = η (u ν ), f ν+ f ν ν η(u ν (t)) { =? Fν+ F ν } ν η(u ν ())
13 Entropy variables and entropy conservative schemes Fi an entropy η(u). Set Entropy variables: v v(u) := η (u). Conveity of η( ) implies that u v is -: v ν = η (u ν )... Entropy conservation: η (u ν ), f ν+ perfect difference { }}{ vν, f f ν+ ν if and only if f?= ν Fν+ F ν perfect difference { }}{ vν+ v ν, f ν+ : Entropy conservative flu: v ν+ v ν, f = ψ(v ν+ ν+ ) ψ(v ν ) Entropy flu potential: ψ(v) := v, f(v) F (u(v)) Abbreviation entropy conservative schemes: v ν+, f = ψ ν+ ν+
14 #. Scalar eamples: f ν+ = ψ(v ν+) ψ(v ν ) v ν+ v ν Toda flow: u t + (e u ) = with Ep entropy pair: (e u ) t + (e u ) = η(u) = v(u) = e u, F (u) = eu and potential ψ(v) := vf F = v Entropy-conservative flu: f ν+ = ψ(v ν+) ψ(v ν ) v ν+ v ν = v ν+ v ν = v ν+ v ν (v ν + v ν+ ) = [ ] e uν + e u ν+ d Toda flow: dt u ν(t) = euν+(t) e uν (t) d e u ν(t) = Const. dt Inviscid Burgers: u t + ( u ) = quadratic entropy ( u ) t + ( 3 u3 ) = η(u) = u u3, v(u) = u, F (u) = 3 and potential ψ(v) := vf F = 6 u3 Entropy conservative 3 -rule : d dt u ν(t)= [ u ν+ u ] ν [ u ν u ν+ u ν ] u ν(t) =Const.
15 #. Linear equations Linear transport: How to discretize u such that u t + u = conserves the balance of quartic perfect derivatives ( u 4) t + ( u 4) =? Conservative differences of u and 4u 3 u = ( u 4) : η(u) = u 4, v = 4u 3, f (u) = u, F = u 4 and ψ(u) = 3u 4 = ψ ν+ = 3 v ν+ 4 u 3 ν+ u3 ν ( ) Indeed set u =ν 3 f ν+ 3 4 ( u 4 ν+ uν 4 uν+ 3 u3 ν ( 3 u 4 ν+ u 4 u3 ν 4 ν uν+ 3 u3 ν = 3 4 uν+ 4 u4 ν, u ( u 4 ν+ uν 4 4 u4 ν u 4 ν u 3 ν u 3 ν u4 ν u 4 ν u 3 ν u 3 ν f f ν+ ν u4 ν u 4 ν u 3 ν u 3 ν ) then uν+ 3 u3 ν ) boundary terms of u ) ( u 3 ν+ u 3 ν) ( u 4 ν+ u4 ν u 3 ν+ u3 ν = ) quartic boundary terms
16 #3. Second and higher-order accuracy The u 4 entropy conservative flu is nd-order accurate: f f ν+ ν f ν+ = 3 uν+ 4 u4 ν 4 uν+ 3 : u3 ν is nd-order approimation of u 8u 3 ν+ u f = ψ ν+ u = 3 uν+ 4 u4 ν 4 uν+ 3 3 u3 ν 4 6u ν+ u u ν+ The general case: epress the flu in the viscosity form viscosity form : f ν+ d dt u ν(t) = nd-order: = (f ν + f ν+ ) D ν+ (u ν+ u ν ) nd order centered differencing (Du ) {}}{ [ f (uν+ ) f (u ν ) ] + [{}} ]{ D u ν+ ν+ D u ν ν D := ( ( ν+ ξ ) f ( ) ) u 8 ν+ dξ u ν+ f ν+ ξ= The entropy-conservative flu is nd-order accurate Entropy conservative flues to any order of accuracy (LeFloch & Rohde )
17 #4. Systems: v ν+ v ν, f ν+ = ψ(vν+ ) ψ(v ν ) Choice of path: N linearly independent directions { r j} N Intermediate states { v j ν+ } N j= : Starting with v = v ν+ ν, and followed by ( v ν+ v ν+ v ν ) v j+ = v j + l j, v ν+ ν+ ν+ r j, j =,,..., N (v N+ = v ν+ ν+ ) [ET 3] The conservative scheme d dt u ν(t) = [ ] f f ν+ ν N ψ ( v j+ ) ( f ν+ ψ v j ) = ν+ ν+ j= l j l j, v ν+ is entropy conservative: v ν+ v ν, f ν+ = ψ(vν+ ) ψ(v ν ) j=
18 #5. Entropy conservative Euler scheme (with W.-G. Zhong) f ν+ = N ψ ( v j+ ) ( ν+ ψ v j ) ν+ l j l j, v ν+ j= Entropy function: η(u) = ρs E/e S + γ + Euler entropy variables: v(u)=η u (u)= q/θ /θ Euler entropy flu potential ψ(v) = v, f F (u) = (γ )m path in phase-space: v = v ν, v j+ = v j + l j, v ν+ r j, v 4 = v ν+ {r j } 3 : three linearly independent directions in v-space (Riemann path) j= { l j } 3 : the corresponding orthogonal system j= { m j } 3 : intermediate values of the momentum along the path j= 3 f m j+ m j = (γ ) ν+ l j, v j= ν+ lj
19 #6. Affordable entropy conservative flues for Euler eqs [Ismail & Roe & (9)] D Euler eqs u t + f + g y = : Epress an entropy conservative flu f := (f, f, f 3, f 4 ) ν+ ρ in terms of z := u z ν+ := (z ν + z ν+ ) p v and the averages z ln z ν+ := ν+ p log(z) ν+ f ν+ f ν+ = (z ) ν+ (z 4 ) ln ν+ = (z 4) ν+ (z ) ν+ + (z ) ν+ (z ) ν+ f ν+ f 3 ν+ f 4 ν+ = (z 3) ν+ (z ) f ν+ ν+ = (z ) ν+ ( γ+ γ (z ) ln ν+ f + (z ν+ ) ν+ f + (z ν+ 3 ) ν+ f 3 ν+ satisfies the compatibility relation v ν+ v ν, f ν+ = ψ(v ν+ ) ψ(v ν ) )
20 ρ p u D Euler,density,γ=.4,C v =76,κ=,entropy= ρ ln(pρ γ ), t=.5e 5, =.,T ma =..6 D Euler,velocity,γ=.4,C v =76,κ=,entropy= ρ ln(pρ γ ), t=.5e 5, =.,T ma =..9 t=. t=.5 t=..4 t=. t=.5 t= D Euler,pressure,γ=.4,C v =76,κ=,entropy= ρ ln(pρ γ ), t=.5e 5, =.,T ma =..38 D Euler, entropy conservative scheme w/ 3rd R K mtd, entropy vs time,entropy = ρ ln(pρ γ ) t=..9 t=.5 t= U (total entropy) t Entropy conservative results for Euler s Sod problem: density, velocity, pressure & entropy. spatial grids, η(u) = ρ ln (pρ ) γ
21 ρ p u.9 D Euler,density,γ=.4,C v =76,κ=,entropy= ρ ln(pρ γ ), t=.5e 5, =.5,T ma =. t=. t=.5 t=..6.4 D Euler,velocity,γ=.4,C v =76,κ=,entropy= ρ ln(pρ γ ), t=.5e 5, =.5,T ma =. t=. t=.5 t= D Euler,pressure,γ=.4,C v =76,κ=,entropy= ρ ln(pρ γ ), t=.5e 5, =.5,T ma =. t=. t=.5 t= D Euler, entropy conservative scheme w/ 3rd R K mtd, entropy vs time,entropy = ρ ln(pρ γ ) U (total entropy) t Entropy conservative density, velocity, pressure & entropy w/4 spatial grids, η(u) = ρ ln (pρ ) γ Does not enforce physical solution with numerical viscosity
22 The question of entropy stability dissipation Entropy conservation: η(u) t + F (u) = Entropy decay due to shock discontinuities (La): η(u) t + F (u) Entropy decay is balanced by perfect derivatives: η(u(, t ))d η(u(, t ))d, t > t Q. How much entropy decay is enough? Entropic measure-valued solutions: d f f dt u ν+ ν(t)+ ν = ɛ(dv ) ; d dt η(u ν(t))+ F ν+ F ν =. ɛ dictates the small scale of dissipation...
23 #7. Entropy balance in Navier-Stokes (NS) eq s A semi-discrete scheme of NS equations u t + f(u) = ɛ (D( )v ) {}} vν { η (u ν ), d dt u ν(t) + ( ) ( f f = ɛ D ν+ ν ν+ Recall f ν+ is entropy-conservative: Euler flues conserve entropy: d dt d d NS dissipate: η(u η(u ν (t)) ν (t)) = = ɛ dt dt ν ν ν d ( ρ ν S ν ) = dt ν v ν+ v ν+, f ν+ D ν = ψ ν+ η(u ν (t)) = ψ ν+ v ν+, D ν+ v ν+ v ν ) = ν viscosity heat conduction {}}{{}}{ (λ + µ) ( q ν+ ) ( ) κ ( θ ν+ ) ( ) θ ν+ ν θ ν+ ν
24 ρ p u.9 D N S,density,λ+µ=.8e 5,γ=.4,C v =76,κ=,entropy= ρ ln(pρ γ ), t=.5e 5, =.5 t=. t=.5 t=.. D N S,velocity,λ+µ=.8e 5,γ=.4,C v =76,κ=,entropy= ρ ln(pρ γ ), t=.5e 5, =.5 t=. t=.5 t= D N S,pressure,λ+µ=.8e 5,γ=.4,C v =76,κ=,entropy= ρ ln(pρ γ ), t=.5e 5, =.5 t=. t=.5 t=..38 D N S, entropy conservative scheme w/ 3rd R K mtd,entropy vs time,entropy = ρ ln(pρ γ ) U (total entropy) t Navier-Stokes equations for Sod s problem: viscosity but no heat conduction; 4 spatial grids. η(u) = ρ ln (pρ ) γ
25 ρ p u D Euler,density,γ=.4,C v =76,κ=.3,entropy= ρ ln(pρ γ ), t=.5e 5, =.5,T ma =..4 D Euler,velocity,γ=.4,C v =76,κ=.3,entropy= ρ ln(pρ γ ), t=.5e 5, =.5,T ma =..9 t=. t=.5 t=.. t=. t=.5 t= D Euler,pressure,γ=.4,C v =76,κ=.3,entropy= ρ ln(pρ γ ), t=.5e 5, =.5,T ma =..38 D Euler, entropy conservative scheme w/ 3rd R K mtd, entropy vs time,entropy = ρ ln(pρ γ ) t=..9 t=.5 t= U (total entropy) t Entropy conservative results for Sod s problem: heat conduction but no viscosity; 4 spatial grids. η(u) = ρ ln (pρ ) γ
26 ρ p u D N S,density,λ+µ=.8e 5,γ=.4,C v =76,κ=.3,entropy= ρ ln(pρ γ ), t=.5e 5, =..4 D N S,velocity,λ+µ=.8e 5,γ=.4,C v =76,κ=.3,entropy= ρ ln(pρ γ ), t=.5e 5, =..9 t=. t=.5 t=.. t=. t=.5 t= D N S,pressure,λ+µ=.8e 5,γ=.4,C v =76,κ=.3,entropy= ρ ln(pρ γ ), t=.5e 5, =..38 D N S, entropy conservative scheme w/ 3rd R K mtd,entropy vs time,entropy = ρ ln(pρ γ ) t=..9 t=.5 t= U (total entropy) t Navier-Stokes equations Sod problem: both viscosity and heat conduction; spatial grids, η(u) = ρ ln (pρ ) γ
27 ρ p u D N S,density,λ+µ=.8e 5,γ=.4,C v =76,κ=.3,entropy= ρ ln(pρ γ ), t=.5e 5, =.5. D N S,velocity,λ+µ=.8e 5,γ=.4,C v =76,κ=.3,entropy= ρ ln(pρ γ ), t=.5e 5, =.5.9 t=. t=.5 t=. t=. t=.5 t= D N S,pressure,λ+µ=.8e 5,γ=.4,C v =76,κ=.3,entropy= ρ ln(pρ γ ), t=.5e 5, =.5 t=. t=.5.9 t= D N S, entropy conservative scheme w/ 3rd R K mtd,entropy vs time,entropy = ρ ln(pρ γ ) U (total entropy) t Navier-Stokes equations Sod problem: viscosity and heat conduction; 4 spatial grids, η(u) = ρ ln (pρ ) γ
28 #8. TeCNO: Arbitrarily high-order entropy stable schemes Set numerical viscosity f ν+ f ν+ Any D ν+ is high order and f ν+ = f D ν+ ν+ (v ν+ v ν ) such that entropy stable: v ν+, D ν+ v ν+ < > will do for entropy dissipation; but... v ν+ p ; non-oscillatory? To overcome first order accuracy: D ν+ A class of arbitrarily high-order entropy-stable schemes: ( ) f ν+ := f (v) ν+ D ν+ v ν+, v ν+ := v ν+ v+ ν Entropy stability requires v ν+, D ν+ v ν+ <...? Interface values v ± : reconstructed by ENO (Harten, Engquist, Osher, Shu,..) (i) ENO reconstruction is highly accurate: D ν+ v ν+ v ν+ (ii) It is Essentially Non-Oscillatory (hence the acronym...) ; and... (iii) [FMT ] The sign property: sign v ν+ = sign v ν+ p ( )!
29 Back to computation of entropy measure valued solutions EMV solutions ν,t : t ν,t, id + ν,t, f = How to compute a realization of the Young measure ν,t? Every Young measure can be realized as a law of random variable: (Ω, Σ, P) : ν,t (E) = P(u(, t, ω) E) Initial data: Realize the law of σ = ν (,t=) by a random field u (, ω) Classical case atomic initial delta: σ = δ u () Beyond classical case uncertainty in initial measurement Compute the entropy measure valued solutions: To be well-defined, we need to make sure {u(, t, ω)} form a measurable ensemble in the right space (Ω, Σ, P) (X, B(X )) Scalar case X = L (Risebro, Schwab, Weber,...) Systems? Eample Back to D Euler realize steady KH as MC average: δ u () = Eu (, ω) M M u k (, ω), u (, ω) = a n cos(b n + nπ ) k=
30 #9. Computation of EMV solutions cont d Evolution: compute the ensemble of entropic solutions u,k (, ω) u,k (, t, ω) Realize the EMVs by propagating the law of MC simulation ν,t = M δ M u,k (,t,ω), ν,t, g(λ) = g(u,k (, t, ω)) M k= This is an approimate entropic MV solution: for all φ(, t) C φ t ν,t, id + φ ν,t, f t + φ t ν,t, id + φ ν,t, f ddt + φ( ν, φ(, ) σ ) σ, id, id d = = ν,t n R + R d ν R d and for all conve entropy pairs (η, F ) with φ C φ t ν,t, η + φ ν,t, F t + φ t ν,t, η + φ ν,t, F ddt + φ( ν, φ(, ) σ ) σ, id, id d ν,t n R + R d ν R d Entropic ν,t with weak BV-bound weak lim ν,t k = ν,t is an EMV
31 On entropic measure-valued solutions... Entropic ν,t with weak BV-bound weak lim ν,t The scalar case Wasserstein W -stability: for all ξ s in R: ν,t, u(, t) ξ d σ, u () ξ d. R d R d In particular, if σ = δ u at t = then ν,t = δ u(,t) No uniqueness for non-atomic initial data = ν,t is an EMV Systems weak-strong W -stability: if u W, classical solution ν,t, u(, t) ξ d σ, u () ξ d. R d R d There is no uniqueness but compute the observables the mean: ν,t, λ Eu (, t, ω) the variance: ν,t, λ E ( u (, t, ω) Eu ), etc.,
32 #. The computed observables (a) 8 (b) 56 (c) 5 (d) 4 Figure: Approimate sample means of the density for the Kelvin-Helmholtz problem at t = and different mesh resolutions. All results are with 4 Monte Carlo samples (a) 8 (b) 56 (c) 5 (d) 4 Figure: Approimate sample variances of the density for the Kelvin-Helmholtz problem at t = and different mesh resolutions. All results are with 4 Monte Carlo samples.
33 Comments on the computation of EMVs MC slow convergence Eu (, t, ω) Eu,M (, t, ω) < M / Key point: Each ensemble was computed with M = 4 samples...! (a) Pointwise (b) Mean (c) Variance (d) W -rate Figure: Kelvin-Helmholtz. (a) The Cauchy rates at t = for the density (y-ais) for a single sample of the Kelvin-Helmholtz problem, vs. (b)+(c) different mesh resolutions (-ais) (d) Cauchy rates in the W distance at t = for the density (y-ais) with respect to different mesh resolutions (-ais) Richtmeyer-Meshkov instability: similar computations; other methods...
34 Computing EMVs: Young measures quantify uncertainty (a) t = (b) t =.5 (c) t = (d) t =.5 (e) t = (f) t = (g) t =.5 (h) t = (i) t =.5 (j) t = Figure: The approimate PDF for density ρ at = (.5,.7) (first row) and = (.5,.8) (second row) on a grid of 4 mesh points (a) n = 8 (b) n = 56 (c) n = 5 (d) n = 4 Figure: The approimate PDF for density ρ at the points = (.5,.7) (first row) and = (.5,.8) (second row) a series of meshes
35 Entropy stability theory for difference approimations of nonlinear conservation laws..., Acta Numerica 3 FMT, ENO reconstruction and ENO interpolation are stable, FoCM FMT Arbitrarily high order accurate entropy stable ENO schemes for systems of conservation laws, SINUM FM(K b )T, Construction of approimate entropy measure valued solutions for hyperbolic systems of conservation laws, J. FoCM (5) FMT, Computation of measure-valued solutions Acta Numerica (6) Fjordholm, Lanthaler, Mishra, Statistical solutions of hyperbolic conservation laws I: Foundations, ArXiv. Fjordholm, Mishra Tadmor; b R. Kappli
36 THANK YOU
Dedicated with appreciation to Tai-Ping Liu on his 60 th birthday
ENTROPY STABLE APPROXIMATIONS OF NAVIER-STOKES EQUATIONS WITH NO ARTIFICIAL NUMERICAL VISCOSITY EITAN TADMOR AND WEIGANG ZHONG Abstract. We construct a new family of entropy stable difference schemes which
More informationVariational formulation of entropy solutions for nonlinear conservation laws
Variational formulation of entropy solutions for nonlinear conservation laws Eitan Tadmor 1 Center for Scientific Computation and Mathematical Modeling CSCAMM) Department of Mathematics, Institute for
More informationENTROPY STABLE APPROXIMATIONS OF NAVIER STOKES EQUATIONS WITH NO ARTIFICIAL NUMERICAL VISCOSITY
Journal of Hyperbolic Differential Equations Vol. 3, No. 3 (006 59 559 c World Scientific Publishing Company ENTROPY STABLE APPROXIMATIONS OF NAVIER STOKES EQUATIONS WITH NO ARTIFICIAL NUMERICAL VISCOSITY
More informationc 2012 Society for Industrial and Applied Mathematics
SIAM J. NUMER. ANAL. Vol. 5, No., pp. 544 573 c Society for Industrial and Applied Mathematics ARBITRARILY HIGH-ORDER ACCURATE ENTROPY STABLE ESSENTIALLY NONOSCILLATORY SCHEMES FOR SYSTEMS OF CONSERVATION
More informationEntropy stable schemes for compressible flows on unstructured meshes
Entropy stable schemes for compressible flows on unstructured meshes Tata Institute of Fundamental Research Center for Applicable Mathematics Bangalore deep@math.tifrbng.res.in http://math.tifrbng.res.in/
More informationNovel entropy stable schemes for 1D and 2D fluid equations
Novel entropy stable schemes for D an D flui equations Eitan Tamor an Weigang Zhong Department of Mathematics, Center for Scientific Computation an Mathematical Moeling (CSCAMM, Institute of Physical Science
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 informationOn the entropy stability of Roe-type finite volume methods
Proceedings of Symposia in Applied Mathematics On the entropy stability of Roe-type finite volume methods Mária Lukáčová - Medvid ová and Eitan Tadmor Abstract. We study the entropy stability of a class
More informationEntropy and Relative Entropy
Entropy and Relative Entropy Joshua Ballew University of Maryland October 24, 2012 Outline Hyperbolic PDEs Entropy/Entropy Flux Pairs Relative Entropy Weak-Strong Uniqueness Weak-Strong Uniqueness for
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 informationApplications of the compensated compactness method on hyperbolic conservation systems
Applications of the compensated compactness method on hyperbolic conservation systems Yunguang Lu Department of Mathematics National University of Colombia e-mail:ylu@unal.edu.co ALAMMI 2009 In this talk,
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 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 informationEntropy-stable discontinuous Galerkin nite element method with streamline diusion and shock-capturing
Entropy-stable discontinuous Galerkin nite element method with streamline diusion and shock-capturing Siddhartha Mishra Seminar for Applied Mathematics ETH Zurich Goal Find a numerical scheme for conservation
More informationA NUMERICAL STUDY FOR THE PERFORMANCE OF THE RUNGE-KUTTA FINITE DIFFERENCE METHOD BASED ON DIFFERENT NUMERICAL HAMILTONIANS
A NUMERICAL STUDY FOR THE PERFORMANCE OF THE RUNGE-KUTTA FINITE DIFFERENCE METHOD BASED ON DIFFERENT NUMERICAL HAMILTONIANS HASEENA AHMED AND HAILIANG LIU Abstract. High resolution finite difference methods
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 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 informationVelocity averaging, kinetic formulations and regularizing effects in conservation laws and related PDEs. Eitan Tadmor. University of Maryland
Velocity averaging, kinetic formulations and regularizing effects in conservation laws and related PDEs Eitan Tadmor Center for Scientific Computation and Mathematical Modeling (CSCAMM) Department of Mathematics,
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 informationFUNDAMENTALS OF LAX-WENDROFF TYPE APPROACH TO HYPERBOLIC PROBLEMS WITH DISCONTINUITIES
6th European Conference on Computational Mechanics (ECCM 6) 7th European Conference on Computational Fluid Dynamics (ECFD 7) 1115 June 018, Glasgow, UK FUNDAMENTALS OF LAX-WENDROFF TYPE APPROACH TO HYPERBOLIC
More informationOn the Dependence of Euler Equations on Physical Parameters
On the Dependence of Euler Equations on Physical Parameters Cleopatra Christoforou Department of Mathematics, University of Houston Joint Work with: Gui-Qiang Chen, Northwestern University Yongqian Zhang,
More informationFinite volume approximation of the relativistic Burgers equation on a Schwarzschild (anti-)de Sitter spacetime
Turkish Journal of Mathematics http:// journals. tubitak. gov. tr/ math/ Research Article Turk J Math 2017 41: 1027 1041 c TÜBİTAK doi:10.906/mat-1602-8 Finite volume approximation of the relativistic
More informationAffordable, entropy-consistent, Euler flux functions
Affordable, entropy-consistent, Euler flux functions (with application to the carbuncle phenomenon) Phil Roe Aerospace Engineering University 0f Michigan Ann Arbor Presented at HYP 2006 1 Entropy pairs
More informationA numerical study of SSP time integration methods for hyperbolic conservation laws
MATHEMATICAL COMMUNICATIONS 613 Math. Commun., Vol. 15, No., pp. 613-633 (010) A numerical study of SSP time integration methods for hyperbolic conservation laws Nelida Črnjarić Žic1,, Bojan Crnković 1
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 informationImproved Seventh-Order WENO Scheme
48th AIAA Aerospace Sciences Meeting Including the New Horizons Forum and Aerospace Eposition 4-7 January 2, Orlando, Florida AIAA 2-45 Improved Seventh-Order WENO Scheme Yiqing Shen Gecheng Zha Dept.
More informationThe Role of Convection and Nearly Singular Behavior of the 3D Navier-Stokes Equations
The Role of Convection and Nearly Singular Behavior of the 3D Navier-Stokes Equations Thomas Y. Hou Applied and Comput. Mathematics, Caltech PDE Conference in honor of Blake Temple, University of Michigan
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 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 informationFinite difference Hermite WENO schemes for the Hamilton-Jacobi equations 1
Finite difference Hermite WENO schemes for the Hamilton-Jacobi equations Feng Zheng, Chi-Wang Shu 3 and Jianian Qiu 4 Abstract In this paper, a new type of finite difference Hermite weighted essentially
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 informationA class of the fourth order finite volume Hermite weighted essentially non-oscillatory schemes
Science in China Series A: Mathematics Aug., 008, Vol. 51, No. 8, 1549 1560 www.scichina.com math.scichina.com www.springerlink.com A class of the fourth order finite volume Hermite weighted essentially
More informationPiecewise Smooth Solutions to the Burgers-Hilbert Equation
Piecewise Smooth Solutions to the Burgers-Hilbert Equation Alberto Bressan and Tianyou Zhang Department of Mathematics, Penn State University, University Park, Pa 68, USA e-mails: bressan@mathpsuedu, zhang
More informationarxiv: v1 [cs.na] 22 Sep 2012
Kinetic energy preserving and entropy stable finite volume schemes for compressible Euler and Navier- Stokes equations Praveen Chandrashekar TIFR Center for Applicable Mathematics, Bangalore 5665, India
More informationICES REPORT A Multilevel-WENO Technique for Solving Nonlinear Conservation Laws
ICES REPORT 7- August 7 A Multilevel-WENO Technique for Solving Nonlinear Conservation Laws by Todd Arbogast, Chieh-Sen Huang, and Xikai Zhao The Institute for Computational Engineering and Sciences The
More informationA New Fourth-Order Non-Oscillatory Central Scheme For Hyperbolic Conservation Laws
A New Fourth-Order Non-Oscillatory Central Scheme For Hyperbolic Conservation Laws A. A. I. Peer a,, A. Gopaul a, M. Z. Dauhoo a, M. Bhuruth a, a Department of Mathematics, University of Mauritius, Reduit,
More informationarxiv: v1 [math.na] 30 Jan 2018
ENTROPY CONSERVATIVE SCHEMES AND THE RECEDING FLOW PROBLEM AYOUB GOUASMI, SCOTT MURMAN, AND KARTHIK DURAISAMY arxiv:1801.1013v1 [math.na] 30 Jan 018 Abstract. This work delves into the family of entropy
More informationBlow-up or No Blow-up? the Role of Convection in 3D Incompressible Navier-Stokes Equations
Blow-up or No Blow-up? the Role of Convection in 3D Incompressible Navier-Stokes Equations Thomas Y. Hou Applied and Comput. Mathematics, Caltech Joint work with Zhen Lei; Congming Li, Ruo Li, and Guo
More informationEffects of a Saturating Dissipation in Burgers-Type Equations
Effects of a Saturating Dissipation in Burgers-Type Equations A. KURGANOV AND P. ROSENAU Tel-Aviv University Abstract We propose and study a new variant of the Burgers equation with dissipation flues that
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 informationEnergy-Preserving and Stable Approximations for the Two-Dimensional Shallow Water Equations
Energy-Preserving and Stable Approximations for the Two-Dimensional Shallow Water Equations Eitan Tadmor* and Weigang Zhong Department of Mathematics, Center for Scientific Computation and Mathematical
More informationHFVS: An Arbitrary High Order Flux Vector Splitting Method
HFVS: An Arbitrary High Order Flu Vector Splitting Method Yibing Chen, Song Jiang and Na Liu Institute of Applied Physics and Computational Mathematics, P.O. Bo 8009, Beijing 00088, P.R. China E-mail:
More informationNonlinear Regularizing Effects for Hyperbolic Conservation Laws
Nonlinear Regularizing Effects for Hyperbolic Conservation Laws Ecole Polytechnique Centre de Mathématiques Laurent Schwartz Collège de France, séminaire EDP, 4 mai 2012 Motivation Consider Cauchy problem
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 informationFoundations of Computational Mathematics: Hong Kong 2008
London Mathematical Society Lecture Note Series: 363 Foundations of Computational Mathematics: Hong Kong 008 Edited by Felipe Cucker City University of Hong Kong Allan Pinkus Technion Mike Todd Cornell
More informationFLOW AROUND A SYMMETRIC OBSTACLE
FLOW AROUND A SYMMETRIC OBSTACLE JOEL A. TROPP Abstract. In this article, we apply Schauder s fied point theorem to demonstrate the eistence of a solution to a certain integral equation. This solution
More informationNotes: Outline. Diffusive flux. Notes: Notes: Advection-diffusion
Outline This lecture Diffusion and advection-diffusion Riemann problem for advection Diagonalization of hyperbolic system, reduction to advection equations Characteristics and Riemann problem for acoustics
More informationNumerische Mathematik
Numer. Math. ) :545 563 DOI.7/s--443-7 Numerische Mathematik A minimum entropy principle of high order schemes for gas dynamics equations iangxiong Zhang Chi-Wang Shu Received: 7 July / Revised: 5 September
More informationAn Improved Non-linear Weights for Seventh-Order WENO Scheme
An Improved Non-linear Weights for Seventh-Order WENO Scheme arxiv:6.06755v [math.na] Nov 06 Samala Rathan, G Naga Raju Department of Mathematics, Visvesvaraya National Institute of Technology, Nagpur,
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 informationInverse Lax-Wendroff Procedure for Numerical Boundary Conditions of. Conservation Laws 1. Abstract
Inverse Lax-Wendroff Procedure for Numerical Boundary Conditions of Conservation Laws Sirui Tan and Chi-Wang Shu 3 Abstract We develop a high order finite difference numerical boundary condition for solving
More informationV (r,t) = i ˆ u( x, y,z,t) + ˆ j v( x, y,z,t) + k ˆ w( x, y, z,t)
IV. DIFFERENTIAL RELATIONS FOR A FLUID PARTICLE This chapter presents the development and application of the basic differential equations of fluid motion. Simplifications in the general equations and common
More information3 Generation and diffusion of vorticity
Version date: March 22, 21 1 3 Generation and diffusion of vorticity 3.1 The vorticity equation We start from Navier Stokes: u t + u u = 1 ρ p + ν 2 u 1) where we have not included a term describing a
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 informationANOMALOUS SHOCK DISPLACEMENT PROBABILITIES FOR A PERTURBED SCALAR CONSERVATION LAW
ANOMALOUS SHOCK DISPLACEMENT PROBABILITIES FOR A PERTURBED SCALAR CONSERVATION LAW JOSSELIN GARNIER, GEORGE PAPANICOLAOU, AND TZU-WEI YANG Abstract. We consider an one-dimensional conservation law with
More information0.2. CONSERVATION LAW FOR FLUID 9
0.2. CONSERVATION LAW FOR FLUID 9 Consider x-component of Eq. (26), we have D(ρu) + ρu( v) dv t = ρg x dv t S pi ds, (27) where ρg x is the x-component of the bodily force, and the surface integral is
More informationCHAPTER 4 BOUNDARY LAYER FLOW APPLICATION TO EXTERNAL FLOW
CHAPTER 4 BOUNDARY LAYER FLOW APPLICATION TO EXTERNAL FLOW 4.1 Introduction Boundary layer concept (Prandtl 1904): Eliminate selected terms in the governing equations Two key questions (1) What are the
More informationMATH 220: Problem Set 3 Solutions
MATH 220: Problem Set 3 Solutions Problem 1. Let ψ C() be given by: 0, < 1, 1 +, 1 < < 0, ψ() = 1, 0 < < 1, 0, > 1, so that it verifies ψ 0, ψ() = 0 if 1 and ψ()d = 1. Consider (ψ j ) j 1 constructed as
More informationCENTRAL DISCONTINUOUS GALERKIN METHODS ON OVERLAPPING CELLS WITH A NON-OSCILLATORY HIERARCHICAL RECONSTRUCTION
CENTRAL DISCONTINUOUS GALERKIN METHODS ON OVERLAPPING CELLS WITH A NON-OSCILLATORY HIERARCHICAL RECONSTRUCTION YINGJIE LIU, CHI-WANG SHU, EITAN TADMOR, AND MENGPING ZHANG Abstract. The central scheme of
More informationAA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS
AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS 1 / 29 AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS Hierarchy of Mathematical Models 1 / 29 AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS 2 / 29
More informationMulti-Domain Hybrid Spectral-WENO Methods for Hyperbolic Conservation Laws
Multi-Domain Hybrid Spectral-WENO Methods for Hyperbolic Conservation Laws Bruno Costa a, Wai Sun Don b, a Departamento de Matemática Aplicada, IM-UFRJ, Caia Postal 6853, Rio de Janeiro, RJ, C.E.P. 945-97,
More informationE = where γ > 1 is a constant spesific to the gas. For air, γ 1.4. Solving for p, we get. 2 ρv2 + (γ 1)E t
. The Euler equations The Euler equations are often used as a simplification of the Navier-Stokes equations as a model of the flow of a gas. In one space dimension these represent the conservation of mass,
More informationBlow-up or No Blow-up? Fluid Dynamic Perspective of the Clay Millennium Problem
Blow-up or No Blow-up? Fluid Dynamic Perspective of the Clay Millennium Problem Thomas Y. Hou Applied and Comput. Mathematics, Caltech Research was supported by NSF Theodore Y. Wu Lecture in Aerospace,
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 informationThe 2-d isentropic compressible Euler equations may have infinitely many solutions which conserve energy
The -d isentropic compressible Euler equations may have infinitely many solutions which conserve energy Simon Markfelder Christian Klingenberg September 15, 017 Dept. of Mathematics, Würzburg University,
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 informationRiemann Solvers and Numerical Methods for Fluid Dynamics
Eleuterio R Toro Riemann Solvers and Numerical Methods for Fluid Dynamics A Practical Introduction With 223 Figures Springer Table of Contents Preface V 1. The Equations of Fluid Dynamics 1 1.1 The Euler
More 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 informationThe Discontinuous Galerkin (dg) method is a popular numerical method that has received much attention
Fluid Dynamics and Co-located Conferences June 24-27, 2013, San Diego, CA 21st AIAA Computational Fluid Dynamics Conference AIAA 2013-2595 Discontinuous Galerkin method for multifluid Euler equations Marc
More informationSemi-Discrete Central-Upwind Schemes with Reduced Dissipation for Hamilton-Jacobi Equations
Semi-Discrete Central-Upwind Schemes with Reduced Dissipation for Hamilton-Jacobi Equations Steve Bryson, Aleander Kurganov, Doron Levy and Guergana Petrova Abstract We introduce a new family of Godunov-type
More informationSolving the Euler Equations!
http://www.nd.edu/~gtryggva/cfd-course/! Solving the Euler Equations! Grétar Tryggvason! Spring 0! The Euler equations for D flow:! where! Define! Ideal Gas:! ρ ρu ρu + ρu + p = 0 t x ( / ) ρe ρu E + p
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 informationNumerical Methods for Modern Traffic Flow Models. Alexander Kurganov
Numerical Methods for Modern Traffic Flow Models Alexander Kurganov Tulane University Mathematics Department www.math.tulane.edu/ kurganov joint work with Pierre Degond, Université Paul Sabatier, Toulouse
More informationAntony Jameson. AIAA 18 th Computational Fluid Dynamics Conference
Energy Estimates for Nonlinear Conservation Laws with Applications to Solutions of the Burgers Equation and One-Dimensional Viscous Flow in a Shock Tube by Central Difference Schemes Antony Jameson AIAA
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 informationPositivity-preserving high order schemes for convection dominated equations
Positivity-preserving high order schemes for convection dominated equations Chi-Wang Shu Division of Applied Mathematics Brown University Joint work with Xiangxiong Zhang; Yinhua Xia; Yulong Xing; Cheng
More information2 The incompressible Kelvin-Helmholtz instability
Hydrodynamic Instabilities References Chandrasekhar: Hydrodynamic and Hydromagnetic Instabilities Landau & Lifshitz: Fluid Mechanics Shu: Gas Dynamics 1 Introduction Instabilities are an important aspect
More informationGoverning Equations of Fluid Dynamics
Chapter 3 Governing Equations of Fluid Dynamics The starting point of any numerical simulation are the governing equations of the physics of the problem to be solved. In this chapter, we first present
More information12.1 Viscous potential flow (VPF)
1 Energy equation for irrotational theories of gas-liquid flow:: viscous potential flow (VPF), viscous potential flow with pressure correction (VCVPF), dissipation method (DM) 1.1 Viscous potential flow
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 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 informationWave operators with non-lipschitz coefficients: energy and observability estimates
Wave operators with non-lipschitz coefficients: energy and observability estimates Institut de Mathématiques de Jussieu-Paris Rive Gauche UNIVERSITÉ PARIS DIDEROT PARIS 7 JOURNÉE JEUNES CONTRÔLEURS 2014
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 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 informationEuler equations and Turbulence.
Time dependent solutions of the Navier-Stokes and Euler equations and Turbulence. Retired, Laboratoire Jacques Louis Lions, Université Pierre et Marie Luminy September 2011 Introduction The goal of this
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 informationMethods of Vanishing Viscosity for Nonlinear Conservation Laws
Methods of Vanishing Viscosity for Nonlinear Conservation Laws Gui-Qiang G. Chen Mathematical Institute, University of Oxford Gui-Qiang.Chen@maths.ox.ac.uk Cleopatra Christoforou (Cyprus) Qian Ding (Northwestern)
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 informationAnalisis para la ecuacion de Boltzmann Soluciones y Approximaciones
Analisis para la ecuacion de Boltzmann Soluciones y Approximaciones Irene M. Gamba Department of Mathematics and ICES The University of Texas at Austin Buenos Aires, June 2012 Collaborators: R. Alonso,
More informationA minimum entropy principle of high order schemes for gas dynamics. equations 1. Abstract
A minimum entropy principle of high order schemes for gas dynamics equations iangxiong Zhang and Chi-Wang Shu 3 Abstract The entropy solutions of the compressible Euler equations satisfy a minimum principle
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 informationON THE GLOBAL REGULARITY OF SUB-CRITICAL EULER-POISSON EQUATIONS WITH PRESSURE
ON THE GLOBAL REGULARITY OF SUB-CRITICAL EULER-POISSON EQUATIONS WITH PRESSURE EITAN TADMOR AND DONGMING WEI Abstract We prove that the one-dimensional Euler-Poisson system driven by the Poisson forcing
More informationThe inviscid limit to a contact discontinuity for the compressible Navier-Stokes-Fourier system using the relative entropy method
The inviscid limit to a contact discontinuity for the compressible Navier-Stokes-Fourier system using the relative entropy method Alexis Vasseur, and Yi Wang Department of Mathematics, University of Texas
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 informationDerivation of relativistic Burgers equation on de Sitter background
Derivation of relativistic Burgers equation on de Sitter background BAVER OKUTMUSTUR Middle East Technical University Department of Mathematics 06800 Ankara TURKEY baver@metu.edu.tr TUBA CEYLAN Middle
More informationRenormalized and entropy solutions of partial differential equations. Piotr Gwiazda
Renormalized and entropy solutions of partial differential equations Piotr Gwiazda Note on lecturer Professor Piotr Gwiazda is a recognized expert in the fields of partial differential equations, applied
More informationWaves and characteristics: Overview 5-1
Waves and characteristics: Overview 5-1 Chapter 5: Waves and characteristics Overview Physics and accounting: use example of sound waves to illustrate method of linearization and counting of variables
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 informationThe Scalar Conservation Law
The Scalar Conservation Law t + f() = 0 = conserved qantity, f() =fl d dt Z b a (t, ) d = Z b a t (t, ) d = Z b a f (t, ) d = f (t, a) f (t, b) = [inflow at a] [otflow at b] f((a)) f((b)) a b Alberto Bressan
More informationPOSITIVITY PROPERTY OF SECOND-ORDER FLUX-SPLITTING SCHEMES FOR THE COMPRESSIBLE EULER EQUATIONS. Cheng Wang. Jian-Guo Liu
DISCRETE AND CONTINUOUS Website: http://aimsciences.org DYNAMICAL SYSTEMS SERIES B Volume 3 Number May003 pp. 0 8 POSITIVITY PROPERTY OF SECOND-ORDER FLUX-SPLITTING SCHEMES FOR THE COMPRESSIBLE EULER EQUATIONS
More information