Consistency analysis of a 1D Finite Volume scheme for barotropic Euler models
|
|
- Wendy Sanders
- 5 years ago
- Views:
Transcription
1 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, Sophia Antipolis Méditerranée Research Centre, Proect CASTOR 3 Univ Nice Sophia Antipolis, CNRS, LJAD, UMR 7351, France Abstract This wor is concerned with the consistency study of a 1D (staggered inetic) finite volume scheme for barotropic Euler models We prove a Lax-Wendroff-lie statement: the limit of a converging (and uniformly bounded) sequence of stepwise constant functions defined from the scheme is a wea entropic-solution of the system of conservation laws 1 Introduction The model This wor is concerned with the consistency study of a (staggered inetic) Finite Volume (FV) scheme for barotropic Euler models t ρ + x (ρv ) =, (1) t (ρv ) + x (ρv + p(ρ)) = () The unnowns are the density ρ and the velocity V The pressure ( ρ p(ρ) ) is assumed to be C (, ) ) with p(ρ) >, p (ρ) >, p (ρ), ρ > Thus, the sound speed c : ρ p (ρ) is well defined and is an increasing function We consider the problem (1)-() on the bounded domain (, L), T with the boundary conditions V (, t) = = V (L, t), t > and the initial conditions ρ(x, ) = ρ (x), V (x, ) = V (x), x (, L) with ρ, V L (, L) Let Φ : ρ > Φ(ρ) such that ρφ (ρ) Φ(ρ) = p(ρ), ρ > The quantity S = 1 ρ V +Φ(ρ) is an entropy of the system: entropy solutions to (1)-() are required to satisfy: for any ϕ Cc ( ) (, L), T ) such that ϕ, FlorentBerthelin@unicefr thierrygoudon@inriafr mineaud@unicefr S t ϕ + ( S + p(ρ) ) V x ϕ (x, t) dx dt S (x)ϕ(x, ) dx (3) 1
2 The meshes We consider a set of J + 1 points = x 1 < x < < x J < x J+1 = L The x are the edges of the so-called primal mesh T We set x +1/ =x +1 x The centers of the primal cells, x +1/ = (x + x +1 )/ for {1,, J}, realize the dual mesh T We set x = (x 1/ + x +1/ )/ for {,, J 1} and x = size(t ) = max x +1/ The adaptive time step is t and we set t=max t The scheme We analyze the scheme introduced in 1 It wors on staggered grids: the densities, ρ +1/, {1,, J}, are evaluated at centers whereas the velocities, V, {1,, J + 1}, are evaluated at edges We set, for {1,, J} and i {,, J} ρ +1/ = 1 x+1 ρ (x) dx, Vi = 1 xi+1/ V (x) dx (4) x +1/ x x i x i 1/ The density is first updated with a FV approximation on the primal mesh ( x +1/ ρ +1/ +1/) ρ + t ( F+1 F ) =, {1,, J} (5) Then, the velocity is updated with a FV approximation on the dual mesh: ( x ρ V ρ V ) ( ) + t G+1/ G 1/ + π/ +1/ π/ 1/ =, (6) for {,, J}, while V 1 = V J+1 = The density on the edges ρ is defined by x ρ = x +1/ ρ +1/ + x 1/ρ 1/, {,, J} The definition of the fluxes relies on the inetic framewor We refer the reader to 1 for details Let us introduce the two following functions F + and F F ± (ρ, V ) = ρ ξ 1I c(ρ) ξ V c(ρ) d ξ ξ We adopt the following formulas for mass fluxes: F 1 = F J+1 =, F = F + (ρ 1/, V ) + F (ρ +1/, V ), {,, J}, (7) and, for momentum fluxes: G3/ = V F (ρ 5/, V ), G J+1/ = V J F + (ρ J 1/, V J ), G+1/ =V + V ( F + (ρ 1/, V ) + F + (ρ +1/, V +1) ) +1 ( F (ρ +1/, V ) + F (ρ +3/, V +1) ), {,, J 1} The discrete pressure gradient combines a space centered scheme and time semi implicit discretization, namely it uses π / +1/ = ρ +1/ Φ (ρ +1/ ) Φ(ρ +1/ ) (8)
3 Properties of the scheme The analysis is driven by the shapes of the functions F ±, see 1, Lemma 3 Here, we shall use the following properties (i) Smoothness: (ρ, V ) (, ) R F ± (ρ, V ) are of class C 1, (ii) Consistency: F + (ρ, V ) + F (ρ, V ) = ρv, V R, ρ Under CFL conditions, see 1, the scheme preserves the positivity of the discrete density and discrete inetic and internal energies evolution equations hold Lemma 11 Let N N Assume min i ( ρ i+1/ ) > For all {,, N 1}, there exists V >, which depends only on the state (ρ, V ), such that if then, min i ( ρ i+1/ ) >, {,, N} and = = x +1/ t e +1/ e +1/ (9) t min ( x+1/ )V 1, (1) D C, with D = 1 4 x ρ ( V V ), (11) + G +1 G + π / +1/ V+1 V D t, (1) x t EK, E K, + Γ +1/ Γ 1/ + π / +1/ π/ 1/ V + D, (13) t where EK, = 1 ( ) ρ V and e +1/ = Φ(ρ +1/ ) are the inetic and internal energies The fluxes are defined by G 1 = G J+1 = and G = Φ(ρ 1/ )V x ( ) 1/ Φ t ρ 1/ Φ(ρ 1/ ), {,, J}, Γ +1/ = 1 V V+1 F + F (V V+1) F, 1/ =ρ 1/ t ( F (ρ +1/ x, V 1/ ρ ) F (ρ 1/, V + F, +1, {1,, J}, ( V ) ρ 1/ V ), and F, 1 = F,, J+1 =, F = F + (ρ 1/, V ) F (ρ +1/, V ), {,, J} The function Φ is a C extension of the function Φ (see 1, Section 43) Results As in, we prove a Lax-Wendroff-lie statement: the limit of a converging (and uniformly bounded) sequence of stepwise constant functions defined from the scheme is a wea entropic-solution of the system of conservation laws Consistency analysis Notation Assuming that = t = T, we define the reconstructions (i =, 1) ρ (i) = = =1 ρ +i +1/ χ/ +1/, π = = =1 3 π / +1/ χ/ +1/, V = = = V χ /,
4 where χ / = χ x 1/, x +1/ t, t, χ / +1/ = χ x, x +1 t, t We also introduce the following discrete norms ρ,t = max max N 1 J ρ +1/, V,T = max N N ρ 1;BV,T = t J N ρ +1/ ρ 1/, V 1;BV,T = t J = = = =1 ρ BV;1,T = x +1/ ρ +1/ ρ +1/ =1 = max V, J V +1 V, For ϕ Cc ( ) (, L), T ), we set ϕ +1/ = ϕ(x +1/, t ) and ϕ = ϕ(x, t ) The interpolate ϕ T of ϕ on the primal mesh and its discrete derivatives are defined by ϕ T (, ) = ð t ϕ T = ϕ +1/ χ1/ +1/ (, ), ϕ =1 ϕ = =1 +1/ ϕ +1/ t (, t) = T χ / +1/, ð xϕ T = ϕ i+1/ χ/ +1/ = =1 = = ϕ +1/ ϕ x (, t), t >, 1/ χ / Similarly, the interpolate ϕ T of ϕ on T and its discrete derivatives are given by ϕ T (, ) = ð t ϕ T = = = = ϕ χ 1/ (, ), ϕ T (, t) = ϕ t ϕ χ /, ð xϕ T = ϕ χ / (, t), t >, = = ϕ +1 ϕ x = =1 +1/ χ / +1/ Assumptions Let ( T m )m 1 be a sequence of meshes s t size(t m) and a familly of time steps ( t ) m, m 1 verifying t m and (1) Assume that there exists N m N s t N m 1 = t m = T The scheme defines (ρ () m, V m ) m 1 Suppose that ρ () m,t + V m,t C, ρ () m 1;BV,T + V m 1;BV,T C BV (14) holds and, in the case ( ρ p (ρ)) ρ L 1 loc (, ), ( ρ () ) 1,T m there exists ( ρ, V ) L ((, T ) (, L)) such that C We assume that ( ρ () m, V m ) ( ρ, V ) in L r ((, T ) (, L)), 1 r < (15) Main results The uniform bounds imply that there exists constants such that sup A (ρ, V ) C A, with A = F ±, ρ F ± and V F ±, ρ, V C sup B(ρ) C B, with B = Φ, Φ, and Φ ρ C +4(C +C F ± ) 4
5 Note also that we have Φ ( ρ +1/) C Φ,ρ ρ +1/,, Furthermore, we show that ρ () m BV;1,T C by using (5) which allows to dominate x +1/ ρ +1/ ρ +1/ by t C ρf ± ( ρ +1/ ρ 1/ + ρ +3/ ρ +1/ ) + C V F ± V +1 V Consequently, ρ (1) m ρ and π m p( ρ) in L r ((, T ) (, L)); with (4) and since ρ, V L (, L), we get ρ () m (, ) ρ and V m (, ) V in L r ((, L)), 1 r < Finally, in the sequel, when a function ϕ Cc ( ) (, L), T ) is given, we assume that t m and x m are sufficiently small so that ϕ(x, ), x, x 3/ x J+1/, L and ϕ(, t), t t, t N Moreover, since ϕ is smooth, ϕ Tm, ϕ T m ϕ, ð t ϕ Tm, ð t ϕ T m t ϕ, and ð x ϕ Tm, ð xϕ T m x ϕ, in L r ((, T ) (, L)), 1 r Theorem 1 Assume (14) and (15) Then, ( ρ, V ) satisfies (1)-() in the distribution sense in ( Cc ( )), (, L), T ) that is ρ t ϕ + ρ V x ϕ (x, t) dx dt ρ (x)ϕ(x, ) dx =, (16) T ρ V t ϕ + ( ρ V + p( ρ) ) x ϕ Moreover, ( ρ, V ) satisfies the entropy inequality (3) (x, t) dx dt ρ (x)v (x)ϕ(x, ) dx= (17) Proof Let ϕ C c ( (, L), T ) ) For the sae of simplicity, the index m is dropped Mass balance We multiply (5) by ϕ 1 J to obtain +1/ x +1/ (ρ +1/ ρ +1/ )ϕ +1/ = =1 } {{ } :=T 1 and sum the results for N 1 and + t ( F +1 F ) ϕ +1/ = } =1 {{ } :=T For T 1, since ϕ N +1/ =, a discrete integration by part wrt time yields Noting that T 1 = = =1 t x+1 ρ () x t x +1/ ρ ( +1/ ϕ +1/ ϕ +1/ for {,, N 1}, {1,, J}, we get T 1 = = ) J x +1/ ρ +1/ ϕ +1/ =1 ð t ϕ T dx dt = x +1/ ρ ( +1/ ϕ +1/ ) ϕ +1/ ρ () ð t ϕ T dx dt ρ () (x, )ϕ T (x, ) dx For T, by integrating by part wrt space, we readily obtain T = t = J F = 5 ( ϕ +1/ ϕ 1/)
6 Bearing in mind that x = x 1/ + x +1/, we then combine the two following expressions of mass fluxes (see (9)-(ii)) to write F F = = ρ ±1/ V x 1/ x R,± ρ 1/ V with R,± + x +1/ x = F ± (ρ +1/, V ) F ± (ρ 1/, V ) ρ +1/ V x 1/ + R, x This expression of the mass fluxes leads to T = T,1 T, with T,1 = T, = = = t t = 1 x 1/ ρ 1/ + x +1/ρ +1/ = 1 x 1/ R, x +1/ R,+ We now observe that, for {,, N 1}, {,, J}, and t x t ρ () x 1/ t x+1/ ρ () t x V ð x ϕ T dx dt = t x 1/ V ð x ϕ T dx dt = t x +1/ Summing these equalities yields T xj+1/ T,1 = ρ () V ð x ϕ T dx dt = x 3/ since ð x ϕ T (x, ) for x, x 3/ x J+1/, L V V ρ 1/ V ρ +1/ V x +1/ x ϕ +1/ ϕ 1/ x, ϕ +1/ ϕ 1/ x ϕ +1/ ϕ x 1/ ϕ +1/ ϕ x 1/ ρ () V ð x ϕ T dx dt, R,+ With (15), we pass to the limit in T 1 and T,1 We prove that ( ρ, V ) satisfies the mass conservation equation (16) by showing that T, since T, C ρf ± xϕ L V,T ρ 1;BV,T x x Momentum balance We multiply (6) by ϕ and sum for N 1 and J to obtain T 3 +T 4 +T 5 = with ( T 3 = x ρ V ρ V ) ϕ, T 4 = = t = = = ( G +1/ G 1/ ) ϕ, T 5 = For T 3, integrating by part wrt time yields T 3 = = = = x ρ V ( ϕ ϕ t ) = =, ( / π +1/ ) π/ 1/ ϕ x ρ V ϕ 6
7 Next, we observe that t x+1/ ρ () t x 1/ V ð t ϕ T dx dt = V ϕ ϕ t = x ρ V t x+1/ ρ () t x 1/ ( ϕ ϕ ) Summing these equalities for {,, N 1}, {,, J} yields T 3 = = T xj+1/ ρ () x 3/ V ð t ϕ T dx dt ρ () V ð t ϕ T dx dt xj+1/ ρ () x 3/ dx dt (x, )V (x, )ϕ T (x, ) dx, ρ () (x, )V (x, )ϕ T (x, ) dx For T 4 and T 5, we first integrate by part wrt space and obtain T 4 = t = =1 G+1/ ( ϕ +1 ϕ ), T5 = We then use the following expression of the momentum flux G+1/ = 1 (V ) ρ ( ) +1/ + V +1 + Q +1/, Q +1/ = 1 V R,+ + 1 V +1 R, +1 to write T 4 = T 4,1 T 4, with T 4,1 = = t J π / ( +1/ ϕ +1 ) ϕ = =1 ( 1 V +1 V ) F + (ρ +1/, V +1 ) F (ρ +1/, V ) t =1 1 (V ρ ) ( +1/ + V ) (ϕ +1 T 4, = t J Q ( +1/ ϕ +1 ) ϕ = =1 +1 ) ϕ Next, we observe that t ( ) x+1 ρ () ( ) ð V ( ) t V x ϕ T dx dt = t ρ + V +1 ( +1/ ϕ +1 x ) ϕ and, consequently, summing for {,, N 1}, {1,, J}, Similarly, for T 5, we get T 4,1 = ρ () ( ) ð V x ϕ T dx dt T 5 = π ð xϕ T dx dt With (15), we pass to the limit in T 3, T 4,1 and T 5 We obtain that ( ρ, V) satisfies the momentum balance equation (17) by showing that T 4, Indeed, we have ( ) T 4, V,T x ϕ L C ρf ± ρ 1;BV,T + C F ± V 1;BV,T x x 7
8 Entropy inequality We now assume that ϕ Kinetic energy We multiply (13) by t ϕ and sum for N 1 and J We obtain to get T 6 +T 7 + T 8 with T 6 = x EK, E K, ϕ J, T 7 = t Γ +1/ Γ 1/ ϕ, = = T 8 = = t = Integrating by part wrt time yields T 6 = x EK, ϕ ϕ = = = = = π / +1/ π/ 1/ V ϕ + x EK,ϕ =1 1 ( ) ð L ρ() 1 V t ϕ T dx dt ρ() For T 7, we write Γ +1/ = 1 ) 4 +1/( ρ V 3 ( ) + V S +1/ = V V +1 = = D ϕ (x, ) ( V (x, ) ) ϕt (x, ) dx 4 S +1/ where R, +1 R,+ + (V+1 V ) F, + F, +1 ρ +1/ (V + V+1) Integration by part wrt space leads to T 7 = T 7,1 T 7, with T 7,1 = 1 ( ) 3ð ρ() V x ϕ T dx dt, T 7, = 1 4 Finally T 7, x since it is dominated by C ρf ± = t J S+1/ =1 ϕ +1 ϕ x x ϕ L V,T V,T ρ 1;BV,T + ( ) C F ± + V,T ρ,t V 1;BV,T Internal energy Multiply (1) by t ϕ to get T 9 + T 1 + T 11 with T 9 = = =1 T 11 = +1/ x +1/ e +1/ e +1/ ϕ = t =1 π / ( +1/ V and sum for N 1 and 1 J +1/, T 1 = +1 V Owing to integration by part wrt time, we get T 9 = = = =1 x +1/ Φ(ρ +1/ )( ϕ +1/ ϕ +1/ Φ ( ρ () ) ðt ϕ T dx dt = t ) ϕ +1/ =1 = =1 G +1 G ϕ +1/, D ϕ +1/ ) J x +1/ Φ(ρ +1/ )ϕ +1/, =1 Φ ( ρ () (x, ) ) ϕ T (x, ) dx 8
9 For T 1, we rewrite the flux as follows G = 1 x x 1/ Φ(ρ 1/ ) + x +1/Φ(ρ +1/ V ) + U 1, + U, + U 3,, U1, = ( e 1/ V V ), U, = x +1/ e x +1/ e 1/ V, U3, = x 1/ Φ( ρ t 1/) Φ(ρ 1/ ) It leads to T 1 = T 1, T 1,1 T 1, T 1,3 with T 1, = Φ ( ρ () ) V ð x ϕ T dx dt, T 1,i = = t = The term T 1,1 can be bounded as follows T 1,1 C Φ,ρ x ϕ t Ui, ( ϕ +1/ 1/) ϕ, i = 1,, 3 = = x ρ 1/ V V (18) Since a min(a, b)+ b a, we get ρ 1/ ρ + ρ +1/ ρ 1/ This leads to T 1,1 C Φ,ρ x ϕ ( L t J ) x ρ V V + V,T ρ 1;BV,T x = } = {{ } :=T Writing ρ =ρ (ρ ρ ) and using the Cauchy-Schwarz inequality yields ) 1/ ( T (T L ρ,t t = = D ) 1/+ V,T ρ BV;1,T t t 1/ It finally leads to T 1,1 t 1 + x The term T 1, can be bounded as follows T 1, C Φ V,T x ϕ L ρ 1;BV,T x x We now turn to T 1,3 We remar that 1/ ρ 1/ ρ t ( x 1/ C ρf ± ρ +1/ ρ 1/ + ρ 1/ Hence, using the same bound as for T 1,1 yields ( ( ) T 1,3 C Φ x ϕ L C ρf ± + V,T ) ρ 1;BV,T x + T V V ) t 1/ + x Pressure terms It remains to get the limit of T 8 + T 11 = T 1, T 1,1 T 1, T 1,3 9
10 with T 1, = π V ð xϕ T dx dt, T 1,1 = = = D ( ϕ +1/ ) ϕ, T 1, = 1 t J π / ( +1/ V V + V+1 V )( +1 ϕ +1 ) ϕ, = =1 T 1,3 = 1 = t =1 We bound T 1,1 and T 1,3 as follows T 1,1 ( xϕ L = = Note that π / +1/ D π / ( +1/ V ) x x, +1 V )( ϕ +1/ ϕ ( C Φ + C Φ,ρ ) ρ +1/ It readily leads to +1 ) ϕ T 1,3 C π 4 ( ) ( ) xxϕ L V 1;BV,T x x T 1, ( C Φ + C Φ,ρ ) x ϕ L T t 1 With (15), we pass to the limit in T 6, T 7,1, T 9, T 1, and T 1, We arrive at (3) since the other terms tend to References 1 F Berthelin, T Goudon, and S Mineaud Kinetic schemes on staggered grids for barotropic Euler models: entropy-stability analysis 14 in revision R Herbin, J-C Latché, and T T Nguyen Explicit staggered schemes for the compressible Euler equations ESAIM:Proc, 4:83 1, 13 1
Convergence of the MAC scheme for incompressible flows
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 Calanque de Cortiou
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 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 informationHigh Order Semi-Lagrangian WENO scheme for Vlasov Equations
High Order WENO scheme for Equations Department of Mathematical and Computer Science Colorado School of Mines joint work w/ Andrew Christlieb Supported by AFOSR. Computational Mathematics Seminar, UC Boulder
More informationFluid Dynamics. Part 2. Massimo Ricotti. University of Maryland. Fluid Dynamics p.1/17
Fluid Dynamics p.1/17 Fluid Dynamics Part 2 Massimo Ricotti ricotti@astro.umd.edu University of Maryland Fluid Dynamics p.2/17 Schemes Based on Flux-conservative Form By their very nature, the fluid equations
More informationScalar conservation laws with moving density constraints arising in traffic flow modeling
Scalar conservation laws with moving density constraints arising in traffic flow modeling Maria Laura Delle Monache Email: maria-laura.delle monache@inria.fr. Joint work with Paola Goatin 14th International
More informationA FV Scheme for Maxwell s equations
A FV Scheme for Maxwell s equations Convergence Analysis on unstructured meshes Stephanie Lohrengel * Malika Remaki ** *Laboratoire J.A. Dieudonné (UMR CNRS 6621), Université de Nice Sophia Antipolis,
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 informationSemi-Lagrangian Formulations for Linear Advection Equations and Applications to Kinetic Equations
Semi-Lagrangian Formulations for Linear Advection and Applications to Kinetic Department of Mathematical and Computer Science Colorado School of Mines joint work w/ Chi-Wang Shu Supported by NSF and AFOSR.
More informationSome Notes on Elliptic Regularity
Some Notes on Elliptic Regularity Here we consider first the existence of weak solutions to elliptic problems of the form: { Lu = f, u = 0, (1) and then we consider the regularity of such solutions. The
More informationMicro-macro methods for Boltzmann-BGK-like equations in the diffusion scaling
Micro-macro methods for Boltzmann-BGK-like equations in the diffusion scaling Anaïs Crestetto 1, Nicolas Crouseilles 2, Giacomo Dimarco 3 et Mohammed Lemou 4 Saint-Malo, 14 décembre 2017 1 Université de
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 informationRelaxation methods and finite element schemes for the equations of visco-elastodynamics. Chiara Simeoni
Relaxation methods and finite element schemes for the equations of visco-elastodynamics Chiara Simeoni Department of Information Engineering, Computer Science and Mathematics University of L Aquila (Italy)
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 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 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 informationRegularizations of general singular integral operators
Regularizations of general singular integral operators Texas A&M University March 19th, 2011 This talk is based on joint work with Sergei Treil. This work is accepted by Revista Matematica Iberoamericano
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 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 informationExistence Theory for the Isentropic Euler Equations
Arch. Rational Mech. Anal. 166 23 81 98 Digital Object Identifier DOI 1.17/s25-2-229-2 Existence Theory for the Isentropic Euler Equations Gui-Qiang Chen & Philippe G. LeFloch Communicated by C. M. Dafermos
More informationAnomalous transport of particles in Plasma physics
Anomalous transport of particles in Plasma physics L. Cesbron a, A. Mellet b,1, K. Trivisa b, a École Normale Supérieure de Cachan Campus de Ker Lann 35170 Bruz rance. b Department of Mathematics, University
More informationExercises - Chapter 1 - Chapter 2 (Correction)
Université de Nice Sophia-Antipolis Master MathMods - Finite Elements - 28/29 Exercises - Chapter 1 - Chapter 2 Correction) Exercise 1. a) Let I =], l[, l R. Show that Cl) >, u C Ī) Cl) u H 1 I), u DĪ).
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 informationLois de conservations scalaires hyperboliques stochastiques : existence, unicité et approximation numérique de la solution entropique
Lois de conservations scalaires hyperboliques stochastiques : existence, unicité et approximation numérique de la solution entropique Université Aix-Marseille Travail en collaboration avec C.Bauzet, V.Castel
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 informationat time t, in dimension d. The index i varies in a countable set I. We call configuration the family, denoted generically by Φ: U (x i (t) x j (t))
Notations In this chapter we investigate infinite systems of interacting particles subject to Newtonian dynamics Each particle is characterized by its position an velocity x i t, v i t R d R d at time
More informationSTOKES PROBLEM WITH SEVERAL TYPES OF BOUNDARY CONDITIONS IN AN EXTERIOR DOMAIN
Electronic Journal of Differential Equations, Vol. 2013 2013, No. 196, pp. 1 28. ISSN: 1072-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu ftp ejde.math.txstate.edu STOKES PROBLEM
More informationLARGE-TIME ASYMPTOTICS, VANISHING VISCOSITY AND NUMERICS FOR 1-D SCALAR CONSERVATION LAWS
LAGE-TIME ASYMPTOTICS, VANISHING VISCOSITY AND NUMEICS FO 1-D SCALA CONSEVATION LAWS L. I. IGNAT, A. POZO, E. ZUAZUA Abstract. In this paper we analyze the large time asymptotic behavior of the discrete
More informationFinite volumes for complex applications In this paper, we study finite-volume methods for balance laws. In particular, we focus on Godunov-type centra
Semi-discrete central schemes for balance laws. Application to the Broadwell model. Alexander Kurganov * *Department of Mathematics, Tulane University, 683 St. Charles Ave., New Orleans, LA 708, USA kurganov@math.tulane.edu
More 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 information2D compressible vortex sheets. Paolo Secchi
2D compressible vortex sheets Paolo Secchi Department of Mathematics Brescia University Joint work with J.F. Coulombel EVEQ 2008, International Summer School on Evolution Equations, Prague, Czech Republic,
More informationDidier Bresch, Pierre Emmanuel Jabin
Global Existence of Wea Solutions for Compresssible Navier Stoes Equations: Thermodynamically unstable pressure and anisotropic viscous stress tensor Didier Bresch, Pierre Emmanuel Jabin Abstract We prove
More informationNumerical schemes for short wave long wave interaction equations
Numerical schemes for short wave long wave interaction equations Paulo Amorim Mário Figueira CMAF - Université de Lisbonne LJLL - Séminaire Fluides Compréssibles, 29 novembre 21 Paulo Amorim (CMAF - U.
More informationLocal null controllability of the N-dimensional Navier-Stokes system with N-1 scalar controls in an arbitrary control domain
Local null controllability of the N-dimensional Navier-Stokes system with N-1 scalar controls in an arbitrary control domain Nicolás Carreño Université Pierre et Marie Curie-Paris 6 UMR 7598 Laboratoire
More informationFinite Volume Schemes: an introduction
Finite Volume Schemes: an introduction First lecture Annamaria Mazzia Dipartimento di Metodi e Modelli Matematici per le Scienze Applicate Università di Padova mazzia@dmsa.unipd.it Scuola di dottorato
More informationNUMERICAL ANALYSIS OF THE CAHN-HILLIARD EQUATION AND APPROXIMATION FOR THE HELE-SHAW PROBLEM, PART II: ERROR ANALYSIS AND CONVERGENCE OF THE INTERFACE
NUMERICAL ANALYSIS OF THE CAHN-HILLIARD EQUATION AND APPROXIMATION FOR THE HELE-SHAW PROBLEM, PART II: ERROR ANALYSIS AND CONVERGENCE OF THE INTERFACE XIAOBING FENG AND ANDREAS PROHL Abstract. In this
More informationDETERMINATION OF THE BLOW-UP RATE FOR THE SEMILINEAR WAVE EQUATION
DETERMINATION OF THE LOW-UP RATE FOR THE SEMILINEAR WAVE EQUATION y FRANK MERLE and HATEM ZAAG Abstract. In this paper, we find the optimal blow-up rate for the semilinear wave equation with a power nonlinearity.
More informationFinite Elements. Colin Cotter. February 22, Colin Cotter FEM
Finite Elements February 22, 2019 In the previous sections, we introduced the concept of finite element spaces, which contain certain functions defined on a domain. Finite element spaces are examples of
More informationLecture Notes on PDEs
Lecture Notes on PDEs Alberto Bressan February 26, 2012 1 Elliptic equations Let IR n be a bounded open set Given measurable functions a ij, b i, c : IR, consider the linear, second order differential
More informationGas Dynamics: Basic Equations, Waves and Shocks
Astrophysical Dynamics, VT 010 Gas Dynamics: Basic Equations, Waves and Shocks Susanne Höfner Susanne.Hoefner@fysast.uu.se Astrophysical Dynamics, VT 010 Gas Dynamics: Basic Equations, Waves and Shocks
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 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 informationSolution Sheet 3. Solution Consider. with the metric. We also define a subset. and thus for any x, y X 0
Solution Sheet Throughout this sheet denotes a domain of R n with sufficiently smooth boundary. 1. Let 1 p
More informationFinite volume method for conservation laws V
Finite volume method for conservation laws V Schemes satisfying entropy condition Praveen. C praveen@math.tifrbng.res.in Tata Institute of Fundamental Research Center for Applicable Mathematics Bangalore
More informationNonlinear stability of compressible vortex sheets in two space dimensions
of compressible vortex sheets in two space dimensions J.-F. Coulombel (Lille) P. Secchi (Brescia) CNRS, and Team SIMPAF of INRIA Futurs Evolution Equations 2006, Mons, August 29th Plan 1 2 3 Related problems
More informationConservation law equations : problem set
Conservation law equations : problem set Luis Silvestre For Isaac Neal and Elia Portnoy in the 2018 summer bootcamp 1 Method of characteristics For the problems in this section, assume that the solutions
More informationModeling and Numerical Approximation of Traffic Flow Problems
Modeling and Numerical Approximation of Traffic Flow Problems Prof. Dr. Ansgar Jüngel Universität Mainz Lecture Notes (preliminary version) Winter Contents Introduction Mathematical theory for scalar conservation
More informationOrder of Convergence of Second Order Schemes Based on the Minmod Limiter
Order of Convergence of Second Order Schemes Based on the Minmod Limiter Boan Popov and Ognian Trifonov July 5, 005 Abstract Many second order accurate non-oscillatory schemes are based on the Minmod limiter,
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 informationTHE L 2 -HODGE THEORY AND REPRESENTATION ON R n
THE L 2 -HODGE THEORY AND REPRESENTATION ON R n BAISHENG YAN Abstract. We present an elementary L 2 -Hodge theory on whole R n based on the minimization principle of the calculus of variations and some
More informationThe reconstruction problem for the Euler-Poisson system in cosmology
The reconstruction problem for the Euler-Poisson system in cosmology Grégoire LOEPER February, 008 arxiv:math/0306430v [math.ap] 7 Apr 005 Abstract The motion of a continuum of matter subject to gravitational
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 informationHere we used the multiindex notation:
Mathematics Department Stanford University Math 51H Distributions Distributions first arose in solving partial differential equations by duality arguments; a later related benefit was that one could always
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 On weak solutions to a diffuse interface model of a binary mixture of compressible
More informationUniversität des Saarlandes. Fachrichtung 6.1 Mathematik
Universität des Saarlandes U N I V E R S I T A S S A R A V I E N I S S Fachrichtung 6.1 Mathematik Preprint Nr. 253 Monotonicity of Implicit Finite Difference Methods for Hyperbolic Conservation Laws Michael
More informationNumerical methods for conservation laws with a stochastically driven flux
Numerical methods for conservation laws with a stochastically driven flux Håkon Hoel, Kenneth Karlsen, Nils Henrik Risebro, Erlend Briseid Storrøsten Department of Mathematics, University of Oslo, Norway
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 informationA general well-balanced finite volume scheme for Euler equations with gravity
A general well-balanced finite volume scheme for Euler equations with gravity Jonas P. Berberich, Praveen Chandrashekar, Christian Klingenberg Abstract We present a second order well-balanced Godunov-type
More information18.657: Mathematics of Machine Learning
8.657: Mathematics of Machine Learning Lecturer: Philippe Rigollet Lecture 3 Scribe: Mina Karzand Oct., 05 Previously, we analyzed the convergence of the projected gradient descent algorithm. We proved
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 informationAA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS
AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS 1 / 59 AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS The Finite Volume Method These slides are partially based on the recommended textbook: Culbert B.
More informationA unifying model for fluid flow and elastic solid deformation: a novel approach for fluid-structure interaction and wave propagation
A unifying model for fluid flow and elastic solid deformation: a novel approach for fluid-structure interaction and wave propagation S. Bordère a and J.-P. Caltagirone b a. CNRS, Univ. Bordeaux, ICMCB,
More informationDOUBLE-DIFFUSION MODELS FROM A HIGHLY-HETEROGENEOUS MEDIUM
DOUBLE-DIFFUSION MODELS FROM A HIGHL-HETEROGENEOUS MEDIUM R.E. SHOWALTER AND D.B. VISARRAGA Abstract. A distributed microstructure model is obtained by homogenization from an exact micro-model with continuous
More informationStochastic homogenization 1
Stochastic homogenization 1 Tuomo Kuusi University of Oulu August 13-17, 2018 Jyväskylä Summer School 1 Course material: S. Armstrong & T. Kuusi & J.-C. Mourrat : Quantitative stochastic homogenization
More informationNumerical Solutions to Partial Differential Equations
Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University Nonconformity and the Consistency Error First Strang Lemma Abstract Error Estimate
More informationModeling and simulation of bedload transport with viscous effects
Introduction Modeling and simulation of bedload transport with viscous effects E. Audusse, L. Boittin, M. Parisot, J. Sainte-Marie Project-team ANGE, Inria; CEREMA; LJLL, UPMC Université Paris VI; UMR
More informationNumerical explorations of a forward-backward diffusion equation
Numerical explorations of a forward-backward diffusion equation Newton Institute KIT Programme Pauline Lafitte 1 C. Mascia 2 1 SIMPAF - INRIA & U. Lille 1, France 2 Univ. La Sapienza, Roma, Italy September
More informationPARABOLIC LIMIT AND STABILITY OF THE VLASOV FOKKER PLANCK SYSTEM
Mathematical Models and Methods in Applied Sciences Vol. 10, No. 7 (2000 1027 1045 c World Scientific Publishing Company PARABOLIC LIMIT AND STABILITY OF THE VLASOV FOKKER PLANCK SYSTEM F. POUPAUD Laboratoire
More informationFree energy estimates for the two-dimensional Keller-Segel model
Free energy estimates for the two-dimensional Keller-Segel model dolbeaul@ceremade.dauphine.fr CEREMADE CNRS & Université Paris-Dauphine in collaboration with A. Blanchet (CERMICS, ENPC & Ceremade) & B.
More informationSobolev spaces. May 18
Sobolev spaces May 18 2015 1 Weak derivatives The purpose of these notes is to give a very basic introduction to Sobolev spaces. More extensive treatments can e.g. be found in the classical references
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 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 informationAveraging lemmas. P-E Jabin (University of Nice)
Averaging lemmas P-E Jabin (University of Nice) Presentation of the course Kinetic equations : Transport equation in phase space, i.e. on f (x, v) of x and v t f + v x f = g, t 0, x, v R d. As for hyperbolic
More informationFractal Conservation Laws: Global Smooth Solutions and Vanishing Regularization
Progress in Nonlinear Differential Equations and Their Applications, Vol. 63, 217 224 c 2005 Birkhäuser Verlag Basel/Switzerland Fractal Conservation Laws: Global Smooth Solutions and Vanishing Regularization
More informationFATIGUE EFFECTS IN ELASTIC MATERIALS WITH VARIATIONAL DAMAGE MODELS: A VANISHING VISCOSITY APPROACH
FATIGUE EFFECTS IN ELASTIC MATERIALS WITH VARIATIONAL DAMAGE MODELS: A VANISHING VISCOSITY APPROACH ROBERTO ALESSI, VITO CRISMALE, GIANLUCA ORLANDO Abstract. We study the existence of quasistatic evolutions
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 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 informationarxiv: v1 [math.ap] 29 May 2018
Non-uniqueness of admissible weak solution to the Riemann problem for the full Euler system in D arxiv:805.354v [math.ap] 9 May 08 Hind Al Baba Christian Klingenberg Ondřej Kreml Václav Mácha Simon Markfelder
More informationWEYL S LEMMA, ONE OF MANY. Daniel W. Stroock
WEYL S LEMMA, ONE OF MANY Daniel W Stroock Abstract This note is a brief, and somewhat biased, account of the evolution of what people working in PDE s call Weyl s Lemma about the regularity of solutions
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 informationAnswers to Problem Set Number 04 for MIT (Spring 2008)
Answers to Problem Set Number 04 for 18.311 MIT (Spring 008) Rodolfo R. Rosales (MIT, Math. Dept., room -337, Cambridge, MA 0139). March 17, 008. Course TA: Timothy Nguyen, MIT, Dept. of Mathematics, Cambridge,
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 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 informationApplied Math Qualifying Exam 11 October Instructions: Work 2 out of 3 problems in each of the 3 parts for a total of 6 problems.
Printed Name: Signature: Applied Math Qualifying Exam 11 October 2014 Instructions: Work 2 out of 3 problems in each of the 3 parts for a total of 6 problems. 2 Part 1 (1) Let Ω be an open subset of R
More informationAnalysis of the M 1 model: well posedness and diffusion asymptotics
Analysis of the M model: well posedness and diffusion asymptotics Thierry GOUDON Chunjin LIN June 26, 22 Abstract This paper is devoted to the analysis of the M model which arises in radiative transfer
More informationENO and WENO schemes. Further topics and time Integration
ENO and WENO schemes. Further topics and time Integration Tefa Kaisara CASA Seminar 29 November, 2006 Outline 1 Short review ENO/WENO 2 Further topics Subcell resolution Other building blocks 3 Time Integration
More informationMixed exterior Laplace s problem
Mixed exterior Laplace s problem Chérif Amrouche, Florian Bonzom Laboratoire de mathématiques appliquées, CNRS UMR 5142, Université de Pau et des Pays de l Adour, IPRA, Avenue de l Université, 64000 Pau
More informationConservation laws and some applications to traffic flows
Conservation laws and some applications to traffic flows Khai T. Nguyen Department of Mathematics, Penn State University ktn2@psu.edu 46th Annual John H. Barrett Memorial Lectures May 16 18, 2016 Khai
More informationarxiv: v2 [math.na] 7 Jul 2017
A POSERIORI ANALYSIS FOR DYNAMIC MODEL ADAPAION IN CONVECION DOMINAED PROBLEMS JAN GIESSELMANN AND RISAN PRYER Abstract. In this wor we present an a posteriori error indicator for approximation schemes
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 informationOscillating waves and optimal smoothing effect for one-dimensional nonlinear scalar conservation laws
arxiv:1302.1345v1 [math.ap] 6 Feb 2013 Oscillating waves and optimal smoothing effect for one-dimensional nonlinear scalar conservation laws Pierre Castelli and Stéphane Junca February 7, 2013 Pierre Castelli
More informationGALERKIN TIME STEPPING METHODS FOR NONLINEAR PARABOLIC EQUATIONS
GALERKIN TIME STEPPING METHODS FOR NONLINEAR PARABOLIC EQUATIONS GEORGIOS AKRIVIS AND CHARALAMBOS MAKRIDAKIS Abstract. We consider discontinuous as well as continuous Galerkin methods for the time discretization
More informationNumerical Solutions to Partial Differential Equations
Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University The Residual and Error of Finite Element Solutions Mixed BVP of Poisson Equation
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 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 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 informationGlobal Maxwellians over All Space and Their Relation to Conserved Quantites of Classical Kinetic Equations
Global Maxwellians over All Space and Their Relation to Conserved Quantites of Classical Kinetic Equations C. David Levermore Department of Mathematics and Institute for Physical Science and Technology
More informationand finally, any second order divergence form elliptic operator
Supporting Information: Mathematical proofs Preliminaries Let be an arbitrary bounded open set in R n and let L be any elliptic differential operator associated to a symmetric positive bilinear form B
More informationHamburger Beiträge zur Angewandten Mathematik
Hamburger Beiträge zur Angewandten Mathematik Numerical analysis of a control and state constrained elliptic control problem with piecewise constant control approximations Klaus Deckelnick and Michael
More informationIntroduction and Preliminaries
Chapter 1 Introduction and Preliminaries This chapter serves two purposes. The first purpose is to prepare the readers for the more systematic development in later chapters of methods of real analysis
More information