INSTITUTE OF MATHEMATICS THE CZECH ACADEMY OF SCIENCES. Derivation of the Navier-Stokes-Poisson system with radiation for an accretion disk

Transcription

1 INSTITUTE OF MATHEMATICS THE CZECH ACADEMY OF SCIENCES Derivation of the Navier-Stokes-Poisson system with radiation for an accretion disk Bernard Ducomet Šárka Nečasová Milan Pokorný M. Angeles Rodríguez-Bellido Preprint No PRAHA 27

2

3 Derivation of the Navier Stokes Poisson system with radiation for an accretion disk Bernard Ducomet, Šárka Nečasová2, Milan Pokorný 3, M. Angeles Rodríguez Bellido 4 CEA, DAM, DIF, F-9297 Arpajon, France 2 Institute of Mathematics of the Academy of Sciences of the Czech Republic Žitná 25, 5 67 Praha, Czech Republic 3 Charles University, Faculty of Mathematics and Physics Mathematical Inst. of Charles University, Sokolovská 83, Prague 8, Czech Republic 4 Dpto. Ecuaciones Diferenciales y Análisis Numérico and IMUS, Facultad de Matemáticas, Universidad de Sevilla, C/Tarfia, s/n, 42 Sevilla, Spain Abstract We study the 3-D compressible barotropic radiation fluid dynamics system describing the motion of the compressible rotating viscous fluid with gravitation and radiation confined to a straight layer ɛ = ω (, ɛ, where ω is a 2-D domain. We show that weak solutions in the 3-D domain converge to the strong solution of the rotating 2-D Navier Stokes Poisson system with radiation in ω as ɛ for all times less than the maximal life time of the strong solution of the 2-D system when the Froude number is small (F r = O( ɛ, the rotating pure 2-D Navier Stokes system with radiation in ω as ɛ when F r = O(. Key words: Navier Stokes Poisson system, radiation, rotation, Froude number, accretion disk, weak solution, thin domain, dimension reduction. Introduction Our aim in this work is the rigorous derivation of the equations describing objects called accretion disks which are quasi planar structures observed in various places in the universe. From a naive point of view, if a massive object attracts matter distributed around it through the Newtonian gravitation in presence of a high angular momentum, this matter is not accreted isotropically around the central object but forms a thin disk around it. As the three main ingredients claimed by astrophysicists for explaining the existence of such objects are gravitation, angular momentum and viscosity (see [2] [22] [24] for detailed

4 presentations, a reasonable framework for their study seems to be a viscous self-gravitating rotating fluid system of equations. These disks are indeed three-dimensional but their size in the third dimension is usually very small, therefore they are often modeled as two-dimensional structures. Our goal in this paper is to derive rigorously the fluid equations of the disk from the equations set in a thin cylinder of thickness ɛ by passing to the limit ɛ and applying recent techniques of dimensional reduction introduced and applied in various situations by P. Bella, E. Feireisl, D. Maltese, A. Novotný and R. Vodák (see [2], [7], [27] and [28]. The mathematical model which we consider is the compressible barotropic Navier Stokes Poisson system with radiation ([8], [9], [] describing the motion of a viscous radiating fluid confined in a bounded straight layer ɛ = ω (, ɛ, where ω R 2 has smooth boundary. Moreover, as we suppose a global rotation of the system, some new terms appear due to the change of frame. Concerning gravitation a modelization difficulty appears as we consider the restriction to ɛ of the solution of the Poisson equation in R 3 : when the thickness of the cylinder tends to zero, a simple argument shows that the gravitational potential given by the Poisson equation in the whole space goes to zero. So if we want to recover the presence of gravitation at the limit, and then keep track of the physical situation, we will have to impose some scaling conditions. In fact as the limit problem will not depend on x 3, the flow is stratified and we expect that the scaling involves naturally the Froude number; see also [7]. More precisely, the system of equations giving the evolution of the mass density ϱ = ϱ(t, x and the velocity field u = u(t, x = (u, u 2, u 3, as functions of the time t (, T and the spatial coordinate x = (x, x 2, x 3 ɛ R 3, reads as follows: t ϱ div x (ϱ u =, (. t (ϱ u div x (ϱ u u x p(ϱ ϱ χ u = div x S ϱ x φ ϱ x χ x 2 S F. (.2 On the right-hand side of (.2 the radiative momentum S F appears, given by S F = (σ a σ s ςi d ς σ dν, (.3 S 2 where the unknown function I = I(t, x, ς, ν is the radiative intensity; see below for more details concerning the quantities describing the radiative effects. The gravitational body forces are represented by the force term ϱ x φ, where the potential φ obeys Poisson s equation φ = 4πG(ηϱ ( ηg in (, T ɛ. (.4 Above, G is the Newton constant and g is a given function, modelling the external gravitational effect. Solving (.4 in the whole space and supposing that ϱ is extended by outside ɛ, we have ηϱ(t, y ( ηg( y φ(t, x = G d y. (.5 R 3 x y The parameter η may take the values or : for η = self-gravitation is present and for η = gravitation acts only as an external field (some astrophysicists consider self-gravitation 2

5 of accretion disks as small compared to the external attraction by a given massive central object modelled by g, see [24]. Note that for the simplicity reasons we assume the external gravitation to be time independent. We suppose that g belongs to the regularity class such that integral (.5 converges. Moreover, since in the momentum equation the term x φ appears, we also need that R 3 K( x y ( ηϱ(t, y ( ηg( y d y <, where K( x y = x y. The effect of radiation is incorporated into the system through the radiative intensity I = I(t, x, ς, ν, depending, besides the variables t, x, on the direction vector ς S 2, where S 2 denotes the unit sphere in R 3, and the frequency ν >. The action of radiation is then expressed in term of integral average S F with respect to the variables ς and ν. The evolution of the compressible viscous barotropic flow is coupled to radiation through radiative transfer equation [4] which reads c ti ς x I = S, (.6 where c is the speed of light. The radiative source S := S a S s is the sum of an emission absorption term S a := σ a (B(ν, ϱ I and a scattering contribution S s := σ s (Ĩ I, where Ĩ := 4π I d S 2 ς σ. The radiation source S then reads ( S = σ a (B I σ s I d ς σ I. (.7 4π S 2 We further assume: Isotropy: The coefficients σ a, σ s are independent of ς. Grey hypothesis: The coefficients σ a, σ s are independent of ν. The function B = B(ν, ϱ measures the distance from equilibrium and is a barotropic equivalent of the Planck function. Furthermore, we take σ s (ϱ, σ a (ϱ c, (.8 σ a (ϱb(ν, ϱ ( B(ν, ϱ h(ν, h L (, (.9 for any ϱ. Note that relations (.8.9 represent cut-off hypotheses at large density. We need one more assumption on the radiative quantities, ϱ σ a (ϱ, ϱ σ s (ϱ, ϱ B(ϱ, ν, B(ϱ, ν c 2. (. Assumption (.9 is needed in the a priori estimate to get existence of a weak solution, assumption (. will be important later in order to get estimates of the remainder in the relative entropy inequality. 3

6 Our system is globally rotating at uniform velocity χ around the vertical direction e 3 and we denote χ = χ e 3. The Coriolis acceleration ϱ χ u and the centrifugal force term ϱ x χ x 2 is therefore present (see [5]. The pressure is a given function of density satisfying hypotheses p C([, C ((,, p( =, p (ϱ > for all ϱ >, lim ϱ p (ϱ = a > (. ϱγ for a certain γ > 3/2. The viscous stress tensor S fulfils Newton s rheological law determined by S = µ ( x u t x u 2 3 div x u I ξ div x u I, (.2 where µ > is the shear viscosity coefficient and ξ is the bulk viscosity coefficient. Finally, the system is supplemented with the initial conditions ϱ(, x = ϱ,ɛ ( x, u(, x = u,ɛ ( x, I(, x, ς, ν = Ĩ,ɛ( x, ς, ν, x ɛ, ς S 2, ν R (.3 and with the boundary conditions. Here, the situation is more complex. For the velocity, we consider the no slip boundary conditions on the boundary part ω (, ɛ (the lateral part of the domain u ω (,ɛ = (.4 and slip boundary condition on the boundary part ω {, ɛ} (the top and bottom part of the layer u n ω {,ɛ} =, [S( x u n] n ω {,ɛ} =. (.5 Let us remark that we have n = ± e 3 on ω {, ɛ}, hence the first condition in (.5 can be rewritten as u 3 = on ω {, ɛ}. (.6 We imposed the slip condition on the boundary ω {, ɛ} in order to avoid difficulties in passing to the infinitely thin limit; using the no slip boundary condition on the top and bottom part of the layer would imply that the velocity converges to zero when we let ɛ. Similar problem we meet with the radiative intensity. We consider at the lateral part of the boundary the condition { } I(t, x, ς, ν = for ( x, ς Γ ( x, ς ( x, ς ω (, ɛ S 2, ς n. (.7 Considering the same condition also on the top and bottom part of the layer (i.e., for x ω {, ɛ} would lead to a situation we try to avoid: in the limit, the radiation disappears. We therefore consider I(t, x, { ς, ν = I(t, x, ς 2( ς n n, ν } for ( x, ς Γ 2 ( x, ς ( x, ς ω {, ɛ} S 2, ς n. (.8 4

7 This boundary condition is called specular reflection. More details needed for our paper will be given later, see also [] for further comments and different possibilities. Our proof will be based on the relative entropy inequality, developed by Feireisl, Novotný and coworkers in [3] and [2]. Recall, however, that the relative entropy inequality was first introduced in the context of hyperbolic equations in the work of C. Dafermos [6], then developed by A. Mellet and A. Vasseur [9], L. Saint-Raymond [25] and finally extended to the compressible barotropic case by P. Germain [5]. Remark. The relativistic version of system (..7 has been introduced by Pomraning [23] and Mihalas and Weibel Mihalas [2] and investigated more recently in astrophysics and laser applications (in the inviscid case by Lowrie, Morel and Hittinger [6] and Buet and Desprès [3], with a special attention to asymptotic regimes. In the remaining part of this section we suitably rescale our system of equations and formulate the primitive and the target system. Section 2 contains definition of the weak solution to our system. Section 3 deals with the existence of solutions to the target system. In Section 4 we present the relative entropy inequality and state the convergence result for our thin disk model. Last Section 5 contains the proof of the convergence result.. Formal scaling analysis, primitive system and target system We rescale our problem to a fixed domain. To this aim, we introduce (x h, ɛx 3 ɛ (x h, x 3, where x h = (x, x 2 ω, x 3 (,, however, keep the notation ϱ for the density, u for the velocity and I for the radiative intensity. We further denote ɛ = ( h, ɛ x 3, div ɛ u = div h u h ɛ x 3 u 3, x h = (x, x 2, u h = (u, u 2, h = ( x, x2, div h u h = x u x2 u 2. Moreover, in order to identify the appropriate limit regime, we perform a general scaling. Since we are only interested in the behaviour of the Froude number, we set all other nondimensional numbers immediately equal to one. The continuity equation reads now the momentum equation is t ϱ div ɛ (ϱ u =, (.9 t (ϱ u div ɛ (ϱ u u ɛ p(ϱ ϱ χ u = div ɛ S( ɛ u F r 2 ϱ ɛφ ϱ ɛ χ x 2 S F, (.2 5

8 and the transport equation has the form where t I ς x I = S = σ a (B I σ s ( 4π S 2 I d ς σ I, (.2 ɛ φ(t, x = ɛ ηϱ(t, y (x y, x 2 y 2, ɛ(x 3 y 3 d y (.22 ( x h y h 2 ɛ 2 (x 3 y ( ηg( y (x y, x 2 y 2, ɛx 3 y 3 d y =: ɛηφ R 3 ( x h y h 2 (ɛx 3 y ( η Φ 2 =: Φ, 2 χ = (,,, ɛ χ x 2 = ( h χ x 2, = (x, x 2,, (.23 x 2 x 2 2 and recall We denote (cf. (.3 S F = (σ a σ s S 2 ςi d ς σ dν. (.24 ϱ(, x = ϱ,ɛ ( x, u(, x = u,ɛ ( x, I(, x, ς, ν = I,ɛ ( x, ς, ν, x, ς S 2, ν R. (.25 We now distinguish two cases with respect to the behaviour of the Froude number, namely F r and F r ɛ. In order to avoid technicalities, we directly consider either F r = ɛ or F r =. Furthermore, according to the choice of the Froude number, we have to consider the correct form of the gravitational potential, namely in the former the self-gravitation and in the latter the external gravitation force. In the latter, we could also include the self-gravitation, it would, however, disappear after the limit passage ɛ. Supposing F r = ɛ and η =, we get the primitive system t ϱ div ɛ (ϱ u =, (.26 t (ϱ u div ɛ (ϱ u u ɛ p(ϱ ϱ χ u = div ɛ S( ɛ u ϱφ ϱ ɛ χ x 2 S F (.27 ( t I ς ɛ I = σ a (B I σ s I d ς σ I. (.28 4π S 2 Next, taking F r = and η =, the primitive system reads t ϱ div ɛ (ϱ u =, (.29 t (ϱ u div ɛ (ϱ u u ɛ p(ϱ ϱ χ u = div ɛ S( ɛ u ϱφ 2 ϱ ɛ χ x 2 S F, (.3 ( t I ς ɛ I = σ a (B I σ s I d ς σ I. (.3 4π S 2 Our goal is to investigate the limit process ɛ in the systems of equations ( and (.29.3, respectively, under the assumptions that initial data [ϱ,ɛ, u,ɛ, I,ɛ ] converge in a certain sense to [r, V, J ] = [r,h, ( w,h,, J,h ]. Let us return back to the former, i.e. F r = ɛ and η =. As the target system does not depend on the vertical variable x 3, we expect that the sequence [ϱ ɛ, u ɛ, I ɛ ] of weak 6

9 solutions to ( will converge to [r, V, J] for V = [ w, ], where w = (w, w 2 and the triple [r(t, x h, w(t, x h, J(t, x h, ς, ν] solves the following 2-D rotating Navier Stokes Poisson system with radiation in the domain (, T ω t r div h (r w =, (.32 r t wr w h w h p(rr( χ w h = div h S h ( h wr h φr h ( χ x h 2 S F h, (.33 ( t J ς h J = σ a (r (B J σ s (r J d ς σ J, (.34 4π S 2 with the formula r(t, y h φ(t, x h = x h y h d y h, (.35 ω where S h ( h w = µ ( h w ( h w T ( µ div h w I h ξ divh w I h. (.36 3 Above, I h is the unit tensor in R 2 2, S F h = (σ a σ s ς h J d ς σ dν, (.37 S 2 and ( χ w h = ( w 2, χw, ( χ x h 2 = x h 2, ς h = (ς, ς 2. When F r =, we also expect that the sequence [ϱ ɛ (t, x, u ɛ (t, x, I ɛ (t, x, ς, ν] of weak solutions to (.29.3 will converge to [r, V, J], where the velocity vector V is as above, [r(t, x h, w(t, x h, J(t, x h, ς, ν] solves now the 2-D rotating Navier Stokes system with radiation and external gravitational force t r div h (r w =, (.38 r t wr w h w h p(rr( χ w h = div h S h ( h wr h φr h ( χ x h 2 S F h, (.39 ( t J ς h J = σ a (B J σ s J d ς σ J, (.4 4π S 2 where g( y φ(t, x h = d y. (.4 R 3 xh y h 2 y3 2 Observe that, through formula (.35, the gravitational contribution in the target momentum equation for F r = ɛ is the tangential gradient of a single layer potential which actually is different from the analogous quantity deriving from the solution of the 2-D Poisson equation h φ = Gr, which would lead to the well-known logarithmic expression. Finally we check, as stressed by Maltese and Novotný [7], that the bulk viscosity coefficient is modified in the limit (compare (.36 with (.2. Our aim is now to prove that solutions of ( and (.29.3 converge in a certain sense (to be precised to the unique solution of ( and (.38.4, respectively. Note also that considering the boundary conditions of the type (.7 on the whole boundary of we would get in the limit that J. Our method would yield that the solutions to ( and (.29.3, respectively, would converge to the same system as above, however, without the radiation. 7

10 2 Weak solutions of the primitive system We consider the rescaled problems ( and (.29.3, respectively, with boundary conditions u ω (, =, (2. and I(t, x, ς, ν = for ( x, ς Γ u n ω {,} =, [S( x u n] n ω {,} =, (2.2 { } ( x, ς ( x, ς ω (, S 2, ς n, (2.3 I(t, x, { ς, ν = I(t, x, ς 2( ς n n, ν } for ( x, ς Γ 2 ( x, ς ( x, ς ω {, } S 2, ς n. (2.4 We define the adapted functional space W,2, n (; R3 = { u W,2 (; R 3 : u n ω {,} =, u ω (, = }. In the weak formulation of the Navier Stokes Poisson system, equation of continuity (.26 is replaced by its weak version τ ( ϱϕ(τ, d x ϱ ɛ, ϕ(, d x = ϱ t ϕ u ɛ ϕ dt d x, (2.5 satisfied for all τ (, T ] and any test function ϕ Cc ([, T. Similarly, the momentum equation (.27 is replaced by ϱ u ϕ(τ, d x ϱ ɛ, u ɛ, ϕ(, d x = τ τ ( ϱ u t ϕ ϱ u u : ɛ ϕ ϱ( χ u ϕ p(ϱdiv ɛ ϕ d x dt (2.6 ( S : ɛ ϕ ϱ Φ j ϕ ϱ ɛ χ x 2 ϕ S F ϕ d x dt, for any ϕ Cc ([, T ; R 3 such that ϕ [,T ] ω {,} = and ϕ 3 [,T ] ω {,} =. Above, j = if η = (i.e. F r = ɛ and j = 2 if η = (i.e. F r =. The radiative transport equation is satisfied in the following sense τ S 2 Iϕ(τ, d ς σ dν d x S 2 I ς ɛ ϕ d ς σ dν d x dt τ ω {,} τ τ ω (, S 2 { ς n } 8 I t ϕ d ς σ dν d x dt S 2 I ς Nϕ d ς σ dν d x σ dt S 2 { ς n } I ς Nϕ d ς σ dν d x σ dt (2.7

11 τ = ω {,} S 2 { ς n } S 2 I,ɛ ϕ(, d ς σ dν d x I(t, x, ς 2( ς n n, ν ς Nϕ d ς σ dν d x σ dt τ S 2 Sϕ d ς σ dν d x dt for all ϕ C ([, T ], where N = ( n, n 2, ɛ n 3 with n the external normal to. Note that ς n = ±ς 3 on ω {, }. Moreover, denoting ϱ p(s H(ϱ = ϱ s 2 ds, (2.8 and E R (I = I d ς σ dν, (2.9 S 2 the energy inequality τ [ 2 ϱ u 2 H(ϱ E R (I] (τ, d x τ τ [ ϱφ j u ϱ ɛ χ x 2 u S ] F u S d ς σ dν d x dt S 2 S( ɛ u : ɛ u d x dt d x dt (2. [ ] 2 ϱ,ɛ u,ɛ 2 H(ϱ,ɛ E R (I,ɛ d x holds for a.e. τ (, T, where j =, 2 as above. Its validity is closely connected to the following result. Lemma 2. [Darrozes Guiraud] Under our assumptions, we have for a.a. (t, x (, T I ς n d ς σ dν. S 2 The proof of the lemma can be found in []. We are now in position to define weak solutions of our primitive system. Definition 2. We say that ϱ, u, I is a weak solution of problem ( and (.29.3, respectively, if ϱ, for a.a. (t, x in (, T, ϱ L (, T ; L γ (, I for a.a. (t, x, ς, ν in (, T S 2 (,, u L 2 (, T ; W,2, n (; R3, I L ((, T S 2 (, L (, T ; L ( S 2 (,, I L (, T ; L 2 ( L (S 2 (,, and if ϱ, u, I satisfy the integral identities (2.5, (2.6, (2.7 together with the total energy inequality (2. and the integral representation of the gravitational force (.22. 9

12 We have the following existence result for the primitive system Proposition 2. Assume that ω R 2 is a domain with compact boundary of class C 2ν, ν >. Suppose that the stress tensor is given by (.2 and p verifies (., the boundary conditions are given by ( and the initial data satisfy the conditions H(ϱ,ɛ L (, ϱ,ɛ, ϱ,ɛ = M ɛ >, I,ɛ ( I, I,ɛ (, ν h(ν for a certain h L (,, ( 2 ϱ,ɛ u,ɛ 2 H(ϱ,ɛ E R (I,ɛ d x <. Let γ > 3/2 if η = or γ > 2 7 if η = and let the external force g Lp (R 3 for p = if γ > 6 and p = 6γ 7γ 6 for 3 2 < γ 6. Then problems ( and (.29.3, respectively, admit at least one finite energy weak solution according to Definition 2.. More details can be found in []. Note that the different boundary conditions for the radiation intensity do not cause any troubles due to Lemma 2.. Using this result, in fact, the existence of the solution can be shown using the approach given in [9] when χ = (non rotating case and for no slip condition on. It is first easy to see that the centrifugal term can be treated in the same way as the gravitational term in [9] and that the Coriolis term may be absorbed in the energy by a Gronwall argument. Finally the slip conditions on top and bottom of the domain may be accommodated using the argument of Vodák [28]. 3 Strong solution of the target system We consider our target system ( and (.38.4, respectively, with the boundary conditions w ω = (3. and { J(t, x, ς, ν = } for ( x, ς Γ ( x, ς ( x, ς ω S 2, ς n. Let (r,, J be a given constant state with r >, and J = B(ν, r. We denote (3.2 e := r r L (ω w H (ω;r 2 E R (J E R H (ω T L 2 (ω;r 2 V L 4 (ω;r 4, (3.3 where V is the initial vorticity (recall that Vij = jwi iwj, E R = 4πc B(ν, r dν, T = (r ( µ h w (ξ 3 µ hdiv h w h p(r

13 and E := e h r L 2 (ω;r 4 h r L α (ω;r 2 (3.4 h T L 2 (ω;r 4 h J L 2 (ω;r 2 h r L α (ω;r 2, for an arbitrary fixed α such that 3 < α < 6. The following result holds Proposition 3. Let p C 2 (,. Let (r, w, J H 3 (ω; R 4, inf ω J >, inf ω r > and assume the following compatibility condition ( r h p(r r ( χ V div h S h ( h w r h φ h r h χ x 2 ω = (3.5 holds, where φ h = φ (see (.35 for η = and φ h = φ (see (.4 for η = and V = ( w,. There exist positive constants δ and Γ > depending on the data such that if E Γδ, the triple (r, w, J is the unique classical solution to the Navier Stokes Poisson system with radiation ( and (.38.4 in (, T ω for any T > such that (r, w, J C([, T ]; H 3 (ω; R 4, sup r r L (ω r/2, t t r C([, T ]; H 2 (ω, ( t w, t J C([, T ]; H (ω; R 3 L 2 (, T ; H 2 (ω; R 3. Moreover, there exists δ > such that if e δ, then sup t T ( r r 2 L 2 (ω w 2 L 2 (ω;r 2 J J 2 L 2 (ω xj 2 L 2 (ω;r 2 Γe 2, and ( sup r r L ( J J L ( Γe, t T sup E R (J E R L (ω t T 2 E R. The proof of the Proposition 3. follows from [] and [8]. Remark 3. In fact this solution (r, w, J can be defined in the whole domain = ω (, by the triple (r, V, J, where V = ( w, and all quantities are constant in x 3. Another possible strong solution can be constructed on short time intervals when no restriction on the size of the initial data is imposed. The result reads

14 Proposition 3.2 Let p C 2 (,. Let (r, w, J H 3 (ω; R 4, inf ω J >, inf ω r > and assume the following compatibility condition ( r h p(r r ( χ V div h S h ( h w r h φ h r h χ x 2 ω = (3.6 holds, where φ h = φ (see (.35 for η = and φ h = φ (see (.4 for η = and V = ( w,. There exist positive constant T depending on the data such that on (, T, there exists triple (r, w, J, the unique classical solution to the Navier Stokes Poisson system with radiation ( and (.38.4 such that (r, w, J C([, T ]; H 3 (ω; R 4, t r C([, T ]; H 2 (ω, ( t w, t J C([, T ]; H (ω; R 3 L 2 (, T ; H 2 (ω. Proof of Proposition 3.2 can be deduced from []. See [4] or [26] for a similar type of results. 4 Relative entropy inequality Let us introduce, in the spirit of [7], a relative entropy inequality which is satisfied by any weak solution (ϱ, u, I of the rotating Navier Stokes Poisson system (..6. We define the relative entropy functional E(ϱ, u, I r, V (, J = 2 ϱ u V 2 E(ϱ, r I J 2 d ς σ dν d x, (4. 2 S 2 with E(ϱ, r = H(ϱ H (r(ϱ r H(r, where (r, V, J is a triple of arbitrary smooth enough functions where only r, V are arbitrary and J satisfies the transport equation for J with the boundary condition (2.4 on ω {, } and (2.3 on ω (,, or J fulfills (2.3 on ω (, and is independent of x 3. Note that the latter case is exactly that we will need later. Then we have Lemma 4. Let all assumptions of Proposition 2. be satisfied and I L 2 ( S 2 (,. Let F r = ɛ and let (ϱ, u, I be a finite energy weak solution of system ( (then j = or F r = and let (ϱ, u, I be a finite energy weak solution of system (.29.3 (then j = 2 in the sense of Definition 2.. Then (ϱ, u, I satisfies the relative entropy inequality τ E(ϱ, u, I r, V, J(τ S( ɛ ( u V : ɛ ( u V d x dt E(ϱ, u, I r, V, J( R(ϱ, u, I, r, V, J, (4.2 2

15 where the remainder R is τ R(ϱ, u, I, r, V, J = τ τ S 2 τ ϱ ( tv u ɛv ( V u d x dt τ S( ɛv : ɛ ( V u d x dt t p(r ϱ ] r u ɛp(r p(ϱdiv ɛv [( ϱ r [ ϱ( χ u ϱ ɛ χ x 2 ϱ Φ j S F ] ( u V d x dt d x dt (4.3 [ ] σ a (ϱ(b(ϱ, ν I σ a (r(b(r, ν J (I J d ς σ dν d x dt τ σ s (r ( 4π [ ( σ s (ϱ I d ς σ I S 4π 2 S ] 2 J d ς σ J (I J d ς σ dν d x dt, S 2 for any triple (r, V, J of test functions such that r C ([, T ], r >, V C ([, T ] ; R 3, V ω (, =, and either V 3 ω {,} = or V is independent of x 3, and J satisfies the transport equation for J with the boundary condition (2.4 at ω {, } and (2.3 at ω (,, or J fulfills (2.3 at ω (, and is independent of x 3. Proof: Using ϕ = 2 V 2 as test function in (2.5 we get 2 ϱ V 2 τ (τ, d x 2 ϱ,ɛ V ( 2 (, d x = ϱ V tv u ɛv V d x dt. (4.4 Using ϕ = V as test function in (2.6 yields ϱ u V (τ, d x ϱ,ɛ u,ɛ V (, d x τ = (ϱ u tv ϱ u ( u ɛ V ϱ( χ u V p(ϱ div ɛv S( ɛ u : ɛv d x dt τ ( ϱφ j V ϱ ɛ χ x 2 V S F V d x dt. (4.5 3

16 Above, j = for η = and j = 2 for η =. Using ϕ = H (r as test function in (2.5 leads to ϱh (r(τ d x ϱ,ɛ H (r( d x = τ ( ϱ t H (r ϱ u ɛ H (r d x dt. Note that rh (r = H(r p(r, therefore r t H (r = t p(r and r ɛ H (r = ɛ p(r. Employing these identities gives ( ( (ϱ rh (r H(r (τ d x (ϱ rh (r H(r ( d x = τ ( ( ϱ t p(r ϱ r r u ɛp(r d x dt. (4.6 Taking difference between the weak formulation of the transport equation for I and J, using as test function ϕ = I J yields (I J 2 (τ, d ς σ dν d x (I J 2 (, d ς σ dν d x 2 S 2 2 S 2 τ (I J 2 ς 2 N d ς σ dν d x σ dt (4.7 2 = τ τ 2 τ ω {,} ω (, ω {,} S 2 { ς n } S 2 { ς n } S 2 { ς n } (I J 2 ς N d ς σ dν d x σ dt (I J 2 (t, x, ς 2(ς n n, ν ς N d ς σ dν d x σ dt [σ a (ϱ(b(ϱ, ν I σ a (r(b(r, ν J] (I J dν d ς σ d x dt S 2 τ ( ( σ s (ϱ I d ς σ I S 4π 2 S ( 2 σ s (r J d ς σ J (I J d ς σ dν d x dt, 4π S 2 where N = (n, n 2, ɛ n 3. Adding ( and (2. (without the part connected to the radiative transfer equation and recalling that the boundary integrals in (4.7 are nonnegative, we end up with τ E(ϱ, u, I r, V, J(τ S( ɛ ( u V : ɛ ( u V d x dt τ E(ϱ, u, I r, V, J( ϱ ( tv u ɛv ( V u d x dt τ S( ɛv : ɛ ( V u d x dt 4

17 τ τ τ which yields (4.3. S 2 [ ( ϱ t p(r ϱ ] r r u ɛp(r p(ϱdiv ɛv d x dt [ ϱ( χ u ϱ ɛ χ x 2 ϱφ j S ] F ( u V d x dt [ ] σ a (ϱ(b(ϱ, ν I σ a (r(b(r, ν J (I J d ς σ dν d x dt τ σ s (r ( 4π 4. Convergence result [ ( σ s (ϱ I d ς σ I S 4π 2 S ] 2 J d ς σ J (I J d ς σ dν d x dt, S 2 We aim at proving the following result. Theorem 4. Suppose that the pressure p satisfy hypothesis (., and that the stress tensor is given by (.2. Let r, w, J satisfy assumptions of Proposition 3. or 3.2 and let T > be the time interval of existence of the strong solution to the problem ( or (.38.4, respectively, corresponding to r, w, J. In addition to hypotheses of Proposition 2., we suppose that I L 2 ( S 2 (,, (. and either F r =, η =, γ > 3 2 and g Lp (R 3 with p = for γ > 6 and p = γ ( 3 2, 6], and for all x h ω or F r = ɛ, η = and γ 2 5. g( yy 3 R 3 ( d y = x h y h 2 y γ 7γ 6 for Let (ϱ ɛ, u ɛ, I ɛ be a sequence of weak solutions to the 3D compressible Navier-Stokes- Poisson system with radiation ( or (.29.3 with ( emanating from the initial data ϱ, u, I. Suppose that E(ϱ,ɛ, u, I r, V, J, (4.8 where V = [ w, ] and all quantities are extended constantly in x 3 to. Then ess sup t [,T ]E(ϱ ɛ, u ɛ, I ɛ, r, V, J, (4.9 u ɛ V = ( w, strongly in L 2 (, T ; W,2 (; R 3, (4. and the triple (r, w, J restricted to ω satisfies the 2D rotating Navier Stokes Poisson system with radiation ( or (.38.4, respectively, with the boundary condition ( on the time interval [, T ]. 5

18 Remark 4. From (4.9 it follows in addition to (4. ϱ ɛ r in C weak ([, T ]; L γ (, ϱ ɛ r a.e. in (, T, (4. and I ɛ J strongly in L (, T ; L 2 ( S 2 R. (4.2 Corollary 4. Suppose that the pressure p satisfy hypothesis (., and that the stress tensor is given by (.2. Assume that [ϱ ɛ,, u ɛ,, I ɛ, ], ϱ ɛ, satisfy ϱ ɛ, (x dx 3 r weakly in L (ω, (4.3 ϱ ɛ, (x u ɛ, dx 3 r w weakly in L (ω; R 2, (4.4 I ɛ, (x dx 3 J weakly in L (ω S 2 R, (4.5 where r, w, J belong to the regularity class of Propositions 3. and 3.2, and [ ] [ ] 2 ϱ,ɛ u,ɛ 2 I,ɛ 2 H(ϱ,ɛ d x 2 ϱ,ɛ r 2 J 2 H(r d x h. Let [ϱ ɛ, u ɛ, I ɛ ] be a sequence of weak solutions to the 3-D compressible Navier Stokes Poisson system with radiation (..6 emanating from the initial data [ϱ ɛ,, u ɛ,, I ɛ, ]. Then properties ( hold. 5 Proof of Theorem Preliminaries We can easily verify that S( ɛ u : ɛ u = ω (ξ 2 3 µ div ɛ u 2 µ( ɛ u 2 ɛ u : ( ɛ u T (5. for any u W,2 (; R 3. However, for any u W,2, n (; R3 we have ɛ u : ( ɛ u T d x = (div ɛ v 2 d x. Thus for any u W,2, n (; R3 S( ɛ u : ɛ u d x C u 2 W,2 (;R 3, (5.2 6

19 provided µ > and ξ. Moreover, under same assumptions S h ( h w : h w d x h C w 2 W,2 (ω;r 2 (5.3 ω for any w W,2 (ω; R 2. Moreover, note that we also have the Poincaré inequality in the form w L2 (ω;r 2 c h w L2 (ω;r 4 (5.4 for any w W,2 (ω; R 2. Due to the energy equality (2. and Korn s inequality (5.2 above, we have the following bounds for the sequence (ϱ ɛ, u ɛ, I ɛ ϱ ɛ L (,T ;L γ ( ϱ ɛ u ɛ L (,T ;L 2 (;R 3 (5.5 u ɛ L2 (,T ;W,2 (;R 3 I ɛ L (,T ;L ( S 2 R C with the constant C independent of ɛ. These estimates hold if γ 2 5 (if η = or under the assumptions on g from Theorem 4. (if η =, for any γ > 3 2. Note that the limit on γ comes from the gravitational potential, as ϱ ɛ ( y(x y, x 2 y 2, ɛ(x 3 y 3 ( d y C ϱ ɛ L ( x h y h 2 ɛ 2 (x 3 y 3 3 p ( L p ( for < p <, with C independent of ɛ. Thus T ϱ ɛφ2 u ɛ d x dt C ϱ ɛ L (,T ;L γ ( u ɛ L2 (,T ;L 6 (;R 3 Φ 2 C ϱ ɛ 2 L (,T ;L γ ( u ɛ L 2 (,T ;L 6 (;R 3 if γ 2 5. On the other hand, T ϱ ɛφ u ɛ d x dt C ϱ ɛ L (,T ;L γ ( u ɛ L2 (,T ;L 6 (;R 3 Φ C ϱ ɛ L (,T ;L γ ( u ɛ L2 (,T ;L 6 ( R 3 g Lp (R 3 with p from Theorem 4., as g( y(x y, x 2 y 2, ɛx 3 y 3 R 3 ( d y ( x h y h 2 ɛ 2 (x 3 y 3 3 6γ C g Lp (R 3 L 5γ 6 ( 6γ L (,T ;L 5γ 6 (;R 3 6γ L (,T ;L 5γ 6 (;R 3 where we used the embedding W,p L 6γ 5γ 6. Moreover, we may deduce the following estimate for the radiative intensity. Assuming that I belongs to L 2 ( S 2 R, multiplying (.6 by I we get 2 ti 2 2 ς xi 2 = σ a (B II σ s 7 ( 4π S 2 I d ς σ I I.

20 Consequently, denoting Ĩ(t, x, ν = I(t, x, ς, ν d ς σ, 4π S 2 we deduce, after integrating integrating the above expression and using Lemma 2., that τ I 2 (τ, d ς σ dν d x σ a (B I 2 S 2 S 2 d ς σ dν d x dt (5.6 2 τ ( Ĩ 2 σ s I d ς σ dν d x dt 2 S 2 S 2 I 2 d ς dν d x 4π τ σ a B 2 dν d x dt C. We now recall the necessary definitions of essential and residual sets. 5.2 Essential and residual sets For two numbers < ϱ ϱ <, the essential and residual subsets of are defined for a.e. t (, T as follows: { Oess(t ϱɛ = x } 2 ϱ ϱ ɛ(t, x 2ϱ, Ores(t ϱɛ = \Oess(t. ϱɛ (5.7 For any function h defined for a.e. (t, x (, T, we write [h] ϱɛ ess(t, x = h(t, x O ϱɛ ess(t(x, In the sequel we will choose ϱ = inf (,T r and ϱ = sup (,T r. From [7] we have [h] ϱɛ res(t, x = h(t, x O ϱɛ res(t(x. (5.8 Lemma 5. Let < a < b <. There exists a constant C = C(a, b > such that for all ϱ [, and r [a, b] where E(ϱ, r C(a, b ( O ϱ res ϱ γ O ϱ res (ϱ r 2 O ϱ ess, E(ϱ, r = H(ϱ H (r(ϱ r H(r, and ϱ = a, ϱ = b in the definition of the essential set. A consequence of this result is the lower bound E(ϱ, u, I r, V, J (5.9 ( C(ϱ, ϱ O ϱ res [ϱ γ ] O ϱ res [ϱ r] 2 Oess ϱ ϱ u V 2 (I J 2 d ς σ dν d x. S 2 8

21 5.3 Estimates of the remainder In what follows, we plan to use in the relative entropy inequality (4.2 as smooth test functions the solution to the 2-D rotating Navier Stokes Poisson system with rotation constructed in Section 3. To this aim, we slightly rearrange the terms in the remainder (4.3 in order to be able to use the validity of the 2-D system. However, we keep writing all the integrals over and assume for a moment that all functions which are independent of x 3 are constant in this variable, and the third velocity component is zero. Denoting by (ϱ ɛ, u ɛ, I ɛ the solution of the primitive system we get τ R(ϱ ɛ, u ɛ, I ɛ, r, V, J(τ τ ( = ϱ ɛ ( uɛ V ɛv ( V uɛ d x dt τ (ϱ ɛ r ( tv V ɛv ( V u ɛ d x dt τ ( r ( tv V ɛv divɛ S( ɛv ɛ p(r r (χ V r ɛ χ x 2 S F (r, J ( V u ɛ d x dt τ [ ( ϱ t p(r ϱ r r u ɛp(r p(ϱdiv ɛv ɛ p(r ( V u ɛ ] d x dt [ ϱ ɛ ( χ u ɛ r( χ V (ϱ ɛ r ɛ χ x 2 S F (ϱ ɛ, I ɛ ] S F (r, J ( V u ɛ d x dt τ ϱ ɛφj ( V u ɛ d x dt [ ] (σ a (ϱ ɛ (B(ϱ ɛ, ν I ɛ σ a (r(b(r, ν J (I ɛ J d ς σ dν d x dt S 2 τ [ ( σ s (ϱ ɛ I ɛ d ς I ɛ S 4π 2 S 2 ( ] 8 σ s (r J d ς J (I ɛ J d ς σ dν d x dt =: R j. 4π S 2 τ In what follows, we will estimates the terms R j for j =, 2,..., 8. j= 5.3. Estimate of R We have τ R = τ h w L (;R 4 ϱ ɛ u ɛ V 2 L ( dt ϱ ɛ ( u ɛ V ɛv ( V uɛ d x dt τ A(tE(ϱ ɛ, u ɛ, I ɛ r, V, J(t dt. Recall that we used estimate (5.9, that fact that V = ( w, and note that due to Section 3 we know that A = h w L (ω;r 4 L (, T. 9

22 5.3.2 Estimate of R 2 We first consider the part of the integral over the essential set and use again estimate (5.9 from Lemma 5.. τ O ϱɛ ess( (ϱ ɛ r ( tv V ɛv ( V u ɛ d x dt with τ δ t w w h w L (ω;r 2 [ϱ ɛ r] ϱɛ τ u ɛ V 2 dt C(δ L 2 (;R 2 τ ess L 2 ( u ɛ V L 2 (;R 3 B 2 (te(ϱ ɛ, u ɛ, I ɛ r, V, J dt, B = t w w h w L (ω;r 2 L2 (, T. For the residual part we consider separately the regions {ϱ < ϱ/2} and {ϱ > 2ϱ}. Then τ {ϱ<ϱ/2} (ϱ ɛ r ( tv V ɛv ( u ɛ V d x dt Finally τ ϱ τ δ τ τ ϱ τ O ϱɛ res t V V ɛ V uɛ V d x dt t w w h w L (ω;r 2 res Oϱɛ 2 u ɛ V dt L2 (;R 3 u ɛ V 2 τ dt C(δ B 2 (te(ϱ ɛ, u ɛ, I ɛ r, V, J dt. L 2 (;R 3 τ {ϱ>2ϱ} (ϱ ɛ r ( tv V ɛv ( u ɛ V d x dt O ϱɛ res ϱɛ t V V ɛ V ϱɛ u ɛ V d x dt t w w h w L (ω;r 2 [ϱ ɛ] ϱɛ res /2 L ( ϱ ɛ u ɛ V 2 /2 τ C B(tE(ϱ ɛ, u ɛ, I ɛ r, V, J dt. Note that we used ϱ ɛ d x = ϱ,ɛ d x C independently of ɛ. L ( dt dt Estimate of R 3 We use the fact that (r, w, J solves the target system. Therefore we have τ R 3 = r ɛ φ ( V u ɛ d x dt in the case when F r = ɛ, and R 3 = τ r ɛ φ ( V uɛ d x dt in the case when F r =. We will use this fact in the treatment of the term R 6. 2

23 5.3.4 Estimate of R 4 Since we have ( Therefore τ R 4 = [ ϱ ɛ r ( u ɛ V ɛp(r r t p(r = p (r t r = p (rdiv h (r w, ( t p(r = p (r(ϱ ɛ rdiv h w w h p(r ϱ ɛ r (r ϱ ɛ ( p(ϱ ɛ p(r p (r(ϱ ɛ r ] div h w d x dt = R4 R4. 2 In order to estimate the first term, we use a similar approach as in the estimate of R 2. We divide the integral into three parts: over the essential set, the set where ϱ ɛ < ϱ/2 and where ϱ ɛ > 2ϱ. Then τ R, 4 = O ϱɛ (r ϱ ess ɛ ɛp(r ( u ɛ V r d x dt with and, finally, τ h p(r r τ δ u ɛ V L (ω;r 2 [ϱ ɛ r] ϱɛ 2 dt C(δ L 2 (;R 3 C = τ R,2 4 = τ ϱ τ δ ϱ u ɛ V τ h p(r r τ R,3 4 = τ τ h p(r r h p(r r ess L2 ( τ. u ɛ V dt L 2 (;R 3 C 2 (te(ϱ ɛ, u ɛ, I ɛ r, V, J dt L 2 (, T, L (ω;r 2 {ϱ<ϱ/2} (ϱ ɛ r ɛp(r ( u ɛ V r d x dt O ϱɛ ɛ p(r res r u ɛ V d x dt Ores ϱɛ L2 (( u ɛ V dt L2 (;R 3 L (ω;r 2 2 dt C(δ L 2 (;R 3 L (ω;r 2 τ C 2 (te(ϱ ɛ, u ɛ, I ɛ r, V, J dt, {ϱ>2ϱ} (ϱ ɛ r ɛp(r ( u ɛ V r d x dt ɛ p(r O ϱɛ res ϱɛ r ϱɛ u ɛ V d x dt [ϱ ɛ ] ϱɛ res /2 ϱ ɛ u ɛ V 2 /2 L ( τ C C(tE(ϱ ɛ, u ɛ, I ɛ r, V, J dt. 2 L ( dt

24 Similarly, we will deal with R4. 2 Using Taylor formula and the regularity of the pressure, and dividing the integral over the essential and residual sets, we have τ R 2, [ 4 = (p(ϱɛ p (r(ϱ ɛ r p(r ] ϱ ɛ ess div h w d x dt with τ τ C div h w L (ω [ϱ ɛ r] ϱɛ ess 2 L 2 ( dt C D(tE(ϱ ɛ, u ɛ, I ɛ r, V, J dt Using the bound we can estimate D = div h w L ( L (, T. [(p(ϱɛ p (r(ϱ ɛ r p(r] ϱɛ R 2,2 4 = τ τ C div h w L (ω Estimate of R 5 We write R 5 = τ res ( ϱ ɛ res [ϱ γ ɛ ] ϱɛ res, [(p(ϱ ɛ p (r(ϱ ɛ r p(r] ϱɛ res div h w d x dt ( O ϱɛ τ [ϱγ res ɛ ] ϱɛ res dt C D(tE(ϱ ɛ, u ɛ, I ɛ r, V, J dt. [ ϱ ɛ χ ( u ɛ V (ϱ ɛ r( χ V (ϱ ɛ r ɛ χ x 2 (σ a (ϱ ɛ σ s (ϱ ɛ ( σ a (r σ s (r σ a (ϱ ɛ σ s (ϱ ɛ Easily, as in the estimate of R and R 2, we have with R 5 R 2 5 R 3 5 C Due to (.8 we have R 4 5 δ τ δ u ɛ V τ τ u ɛ V S 2 ς(j I ɛ d ς σ dν S 2 ςj d ς σ dν E 2 (te(ϱ ɛ, u ɛ, I ɛ r, V, J dt δ E = ( w L (ω;r 2 L 2 (, T. 2 dt C L 2 (;R 3 τ 2 dt C L 2 (;R 3 22 τ ] ( V u ɛ d x dt = τ u ɛ V 5 R j 5. j= 2 dt L 2 (;R 3 S 2 (I ɛ J 2 d ς σ dν d x dt E(ϱ ɛ, u ɛ, I ɛ r, V, J dt.

25 Similarly, using also (., we get τ R5 5 δ u ɛ V 2 dt C L 2 (;R 3 with τ F = J L (;L ((, S 2 L 2 (, T. F 2 (te(ϱ ɛ, u ɛ, I ɛ r, V, J dt Estimate of R 6 For the gravitational potential, we have to consider both cases separately. We start with the simpler one. i.e. with the case F r =. Here, only the gravitational effect of other objects than the fluid itself is considered. Recall that g(yy 3 R 3 ( d y =. x h y h 2 y3 2 3 We combine the term R 3 with R 6. Therefore we have to verify τ lim ( ( V u ɛ ϱ ɛ ( x g( y (5. ɛ R 3 [ ( x h y h, y 3 ( x h y h 2 y ( x h y h, ɛx 3 y 3 ] 3 23 ( x h y h 2 (ɛx 3 y d y d x dt =. First note that due to our assumption on the integrability of g and proceeding similarly as in the estimate of R 2 (replacing the L 2 estimate of V u ɛ by the L 6 estimate is is enough to verify that [ lim r( ( V u ɛ ɛ ( x h y h, y 3 ( x h y h 2 y ( x h y h, ɛx 3 y 3 ] 3 23 ( x h y h 2 (ɛx 3 y g( y (5. R 3 d y d x = for a.a. τ (, T. Moreover, it is not difficult to verify that (note that to get estimates independent of ɛ of the integral over R 3 is easy it remains to verify [ ( x h y h, y 3 lim g( y ɛ R 3 ( x h y h 2 y ( x h y h, ɛx 3 y 3 ] 3 23 ( d y = x h y h 2 (ɛx 3 y for all x h ω, x 3 (,, t (, T and g C c (R 3. As ( ( x h y h, y 3 lim ɛ ( x h y h 2 y ( x h y h, ɛx 3 y ( = x h y h 2 (ɛx 3 y for a.a. (x h, x 3, (y h, y 3 R 3, τ (, T, and ( x h y h, ɛx 3 y 3 ( L x h y h 2 (ɛx 3 y x h y h 2 (ɛx 3 y 3 2 loc (R 3 ɛ [, ], 23

26 the Lebesgue dominated convergence theorem yields the required identity (5.. The case of the self-gravitation (F r = ɛ is more complex. Here, we have to show that ( V [ ϱ ɛ (t, y( x h y h, ɛ(x 3 y 3 u ɛ ϱ ɛ (t, x ( x h y h 2 ɛ 2 (x 3 y 3 2 d y r(t, x r(t, y h ] 3 ɛ x h y h d y h d x δ V u ɛ 2 L 6 (;R 3 C(δ; r, V E(ϱ ɛ, u ɛ, I ɛ r, V, J H ɛ, ω (5.2 where H ɛ = o(ɛ as ɛ. The derivative of the integral over ω with respect to x 3 is indeed zero. First of all, for γ 2 5, as in (5.5, using the decomposition to the essential and the residual set and proceeding as in the estimates of the remainder above, we can show that it is enough to verify that lim r [ r(t, y h ( x h y h, ɛ 2 (x 3 y 3 V ɛ ( x h y h 2 ɛ 2 (x 3 y 3 2 d y r(t, y h ] 3 ɛ ω x h y h d y h d x = for a.a. τ (, T. Using the change of the variables to integrate over ɛ it is enough to show r(t, y h ( x h y h, ɛ(x 3 y 3 Note that lim ɛ [ ( x h y h 2 ɛ 2 (x 3 y d y ɛ ɛ ω r(t, y h x h y h d y r(t, y h ( x h y h h = v.p. ω x h y h 3 2 ω r(t, y h ] x h y h d y h =. where v.p. means the integral in the principal value sense. Thus it remains to show ɛr(t, y h (x 3 y 3 lim ɛ ( d y =, (5.3 x h y h 2 ɛ 2 (x 3 y d y h, and lim ɛ r(t, y h ( x h y h ( x h y h 2 ɛ 2 (x 3 y 3 2 d y = v.p. 3 ω r(t, y h ( x h y h x h y h 3 d y h. (5.4 We fix x ω, >, sufficiently small, and denote B ( x = { x ω; x x < } and C ( x = { x ; x h x <, < x 3 < }. We first consider (5.3. Fix δ >. Using the change of variables (from back to ɛ it is not difficult to see that there exists > such that for any < ɛ, < x 3 < it holds ɛr(t, y h (x 3 y 3 ( x y h 2 ɛ 2 (x 3 y 3 2 d y < δ 3 C ( x and for this > there exists ɛ > such that we have for any < ɛ ɛ ɛr(t, y h (x 3 y 3 ( x y h 2 ɛ 2 (x 3 y 3 2 d y < δ 3 which yields (5.3. \C ( x 24

27 C ( x x h y h x h y h 3 In order to verify (5.4, we proceed similarly. Since is a singular integral kernel in the sense of Calderón Zygmund, for a fixed x ω, < x 3 < and δ > there exists > such that r(t, y h ( x y h ( x y h 2 ɛ 2 (x 3 y 3 2 d y < δ, 3 and v.p. We fix such >. Using that B ( x r(t, y h ( x h y h x h y h 3 d y h < δ. ( x y h 2 ɛ 2 (x 3 y 3 2 as ɛ 3 x y h 3 for any y h ω, < x 3, y 3 <, except x = y h, we see that for the above fixed > there exists ɛ > such that for any < ɛ ɛ r(t, y h ( x h y h ( r(t, y h ( x h y h d y v.p. x h y h 2 ɛ 2 (x 3 y x h y h 3 d y h < δ, \C ( x hence we get (5.4. ω\b ( x Estimate of R 7 and R 8 Repeating the arguments from the estimate of R5 4 and R5, 5 using (.8. (in particular, the Lipschitz continuity of B, σ a and σ s in the density we easily verify that Conclusion R 7 R 8 C τ ( F (t E(ϱɛ, u ɛ, I ɛ r, V, J dt. Collecting all of the previous estimates, plugging them into the relative entropy inequality and taking δ small enough, we end with the inequality where K L (, T and E(ϱ ɛ, u ɛ, I ɛ r, V, J(τ h ɛ (τ τ K(tE(ϱ ɛ, u ɛ, I ɛ r, V, J dt, h ɛ (τ = E(ϱ,ɛ, u,ɛ, I,ɛ r, V, J H ɛ (τ where H ɛ (τ for ɛ in L (, T. Hence, it implies by virtue of Gronwall s lemma E(ϱ ɛ, u ɛ, I ɛ r, V, J(τ h ɛ (τ τ h ɛ (tk(te τ K(s ds t dt for a.a. τ [, T ], which establishes (4.9. Returning back to the relative entropy inequality (4.2, we verify (4. which finishes the proof of Theorem

28 Acknowledgements: B. D. is partially supported by the ANR project INFAMIE (ANR-5-CE4- Š. N. is supported by the Czech Science Foundation, grant No S and by RVO Part of this paper was written during her stay in CEA and she would like to thank Prof. Ducomet for his hospitality during her stay. M. P. was supported by the Czech Science Foundation, grant No S. M.A.R.B. was partially supported by MINECO grant MTM P (Ministerio de Economía y Competitividad, Spain with the participation of FEDER. She would also like to thank Prof. Nečasová for her hospitality during the stay in Prague. References [] G. Allaire, F. Golse, Transport et diffusion, Lecture notes, available on [2] P. Bella, E. Feireisl, A. Novotný, Dimensional reduction for compressible viscous fluids, Acta Appl. Math., 34 (24 2. [3] C. Buet and B. Després. Asymptotic analysis of fluid models for the coupling of radiation and hydrodynamics, J. Quant. Spectroscopy Rad. Transf. 85 ( [4] S. Chandrasekhar, Radiative transfer, Dover Publications, Inc., New York, 96. [5] A.R. Choudhuri, The physics of fluids and plasmas, an introduction for astrophysicists, Cambridge University Press, 998. [6] C.M. Dafermos, The second law of thermodynamics and stability, Arch. Rational Mech. Anal. 7 ( [7] B. Ducomet, M. Caggio, Š. Nečasová, M. Pokorný, The rotating Navier Stokes Fourier Poisson system on thin domains, submitted. [8] B. Ducomet, E. Feireisl, Š. Nečasová, On a model of radiation hydrodynamics, Ann. I. H. Poincaré-AN 28 ( [9] B. Ducomet, E. Feireisl, H. Petzeltová, I. Straškraba, Global in time weak solutions for compressible barotropic self-gravitating fluids, Discrete and Continuous Dynamical Systems ( [] B. Ducomet, Š. Nečasová, Global smooth solution of the Cauchy problem for a model of radiative flow, Ann. della Scuola Norm. Sup. di Pisa 4 ( [] B. Ducomet, Š. Nečasová, Non equilibrium diffusion limit in a barotropic radiative flow, in Proceedings of the International Conference on Recent Advances in PDEs and Applications, in honor of Hugo Beiro da Veiga s 7th birthday, Levico Terme, Italy, , editor(s: Vicentiu D. Radulescu, Adlia Sequeira, Vsevolod A. Solonnikov, Recent Advances in Partial Differential Equations and Applications, American Mathematical Society, Providence, 26,

29 [2] E. Feireisl, B.J. Jin, A. Novotný, Relative entropies, suitable weak solutions, and weakstrong uniqueness for the compressible Navier Stokes system, J. Math. Fluid Mech. 4 ( [3] E. Feireisl, A. Novotný, Weak Strong Uniqueness Property for the Full Navier Stokes Fourier System, Arch. Rational Mech. Anal. 24 ( [4] W. Fiszdon, W.M. Zajaczkowski, The initial boundary value problem for the flow of a barotropic viscous fluid, global in time, Applicable Analysis 5 ( [5] P. Germain, Weak-strong uniqueness for the isentropic compressible Navier Stokes system, J. Math. Fluid Mech. 3 ( [6] R.B. Lowrie, J.E. Morel, J.A. Hittinger, The coupling of radiation and hydrodynamics, The Astrophysical Journal 52 ( [7] D. Maltese, A. Novotný, Compressible Navier-Stokes equations in thin domains, J. Math. Fluid Mech. 6 ( [8] D. Matsumura, T. Nishida, Initial boundary value problems for the equations of motion of compressible viscous and heat-conductive fluid, Commun. Math. Phys. 89 ( [9] A. Mellet, A. Vasseur, Existence and uniqueness of global strong solutions for onedimensional compressible Navier Stokes equations, SIAM J. Math. Anal. 39 ( [2] B. Mihalas and B. Weibel Mihalas, Foundations of radiation hydrodynamics, Dover Publications, Dover, 984. [2] A. Montesinos Armijo, Review: Accretion disk theory, ArXiv:23.685v [astro-ph.he] 3 Mar 22. [22] G.I. Ogilvie, Accretion disks, in Fluid dynamics and dynamos in astrophysics and geophysics, A.M. Soward, C.A. Jones, D.W. Hughes, N.O. Weiss Edts., pp. 28, CRC Press, 25. [23] G.C. Pomraning, Radiation hydrodynamics, Dover Publications, New York, 25. [24] J.E. Pringle, Accretion disks in astrophysics, Ann. Rev. Astron. Astrophys. 9 ( [25] L. Saint-Raymond, Hydrodynamic limits: some improvements of the relative entropy method, Ann. I. H. Poincaré-AN 26 ( [26] A. Valli, W. Zajaczkowski, Navier Stokes equations for compressible fluids: global existence and qualitative properties of the solutions in the general case, Comm. Math. Phys. 3 ( [27] R. Vodák, Asymptotic analysis of three-dimensional Navier Stokes equations for compressible nonlinearly viscous fluids, Dynamics of PDE 5 (

30 Powered by TCPDF ( [28] R. Vodák, Asymptotic analysis of steady and nonsteady Navier Stokes equations for barotropic compressible flows, Acta Appl. Math. (