Key words. Vlasov-Maxwell-Boltzmann system, dissipative structure, stability, time-decay rate
|
|
- Quentin Hopkins
- 5 years ago
- Views:
Transcription
1 DISSIPATIVE PROPERTY OF THE VLASOV-MAXWELL-BOLTZMANN SYSTEM WITH A UNIFORM IONIC BACKGROUND RENJUN DUAN Abstract. In this paper we discuss the dissipative property of near-equilibrium classical solutions to the Cauchy problem of the Vlasov-Maxwell-Boltzmann System in the whole space when the positive charged ion flow provides a spatially uniform background. The most key point of studying this coupled degenerately dissipative system here is to establish the dissipation of the electromagnetic field which turns out to be of the regularity-loss type. Precisely, for the linearized non-homogeneous system, some L energy functionals and L time-frequency functionals which are equivalent with the naturally existing ones are designed to capture the optimal dissipation rate of the system, which in turn yields the optimal L p -L q type time-decay estimates of the corresponding linearized solution operator. These results show a special feature of the one-species Vlasov-Maxwell-Boltzmann system different from the case of two-species, that is, the dissipation of the magnetic field in one-species is strictly weaker than the one in two-species. As a by-product, the global existence of solutions to the nonlinear Cauchy problem is also proved by constructing some similar energy functionals but the time-decay rates of the obtained solution still remain open. Key words. Vlasov-Maxwell-Boltzmann system, dissipative structure, stability, time-decay rate AMS subject classifications. 76X05, 76P05, 35Q0, 35B35, 35B40 1. Introduction. The Vlasov-Maxwell-Boltzmann system is an important model for plasma physics to describe the time evolution of dilute charged particles (e.g., electrons and ions in the case of two-species) under the influence of the self-consistent internally generated Lorentz forces [16]. In physical situations the ion mass is usually much larger than the electron mass so that the electrons move much faster than the ions. Thus, the ions are often described by a fixed ion background n b (x) and only the electrons move. For such simple case, the Vlasov-Maxwell-Boltzmann system takes the form of t f + x f + (E + B) f = Q(f, f), t E x B = f d, (1.1) t B + x E = 0, E = f d n b, x B = 0. Here, the unknowns are f = f(t, x, ) : (0, ) [0, ), E = E(t, x) : (0, ) and B = B(t, x) : (0, ), with f(t, x, ) standing for the number distribution function of one-species of particles (e.g., electrons) which have position x = (x 1, x, x 3 ) and velocity = ( 1,, 3 ) at time t, and E(t, x) and B(t, x) denoting the electromagnetic field in terms of the time-space variable (t, x). The initial data of the system at t = 0 is given by f(0, x, ) = f 0 (x, ), E(0, x) = E 0 (x), B(0, x) = B 0 (x). (1.) Department of Mathematics, The Chinese University of Hong Kong, Shatin, Hong Kong (rjduan@math.cuhk.edu.hk). 1
2 R.-J. DUAN Q is the bilinear Boltzmann collision operator [1] for the hard-sphere model defined by Q(f, g) = (f g fg ) ( ) ω dωd, S f = f(t, x, ), f = f(t, x, ), g = g(t, x, ), g = g(t, x, ), = [( ) ω]ω, = + [( ) ω]ω, ω S. Notice that system (1.1) in general contains physical constants such as the charge and mass of electrons and the speed of light, cf. [11]. Since our purpose in this paper is to investigate the dissipative property of solutions near global Maxwellians, those physical constants in system (1.1) have been normalized to be one for notational simplicity. Through this paper, n b (x) is assumed to be a positive constant denoting the spatially uniform density of the ionic background, and we also set n b = 1 without loss of generality Main results. We are interested in the solution to the Cauchy problem for the case when the number distribution function f(t, x, ) is near an equilibrium state M and E(t, x) and B(t, x) have small amplitudes, where M denotes the normalized Maxwellian For that, set the perturbation u by M = M() = (π) 3/ e /. f(t, x, ) = M + M 1/ u(t, x, ). Then, the Cauchy problem (1.1), (1.) can be reformulated as t u + x u + (E + B) u M 1/ E t E x B = M 1/ u d, t B + x E = 0, E = M 1/ u d, x B = 0, with initial data = Lu + Γ(u, u) + 1 Eu, (1.3) u(0, x, ) = u 0 (x, ), E(0, x) = E 0 (x), B(0, x) = B 0 (x). (1.4) Here, the linear term Lu and the nonlinear term Γ(u, u) are defined in (.1) and (.) later on. The problems to be considered are (i) whether or not any small amplitude solution [u(t), E(t), B(t)] : R + X = H N ( x ) H N ( x) H N ( x) with a properly large N for the above reformulated Cauchy problem uniquely exists for all t > 0 if initial data [u 0, E 0, B 0 ] X is sufficiently small; (ii) if so, does the solution decay in time with some explicit rate? We shall give in this paper a
3 DISSIPATIVE PROPERTY OF VLASOV-MAXWELL-BOLTZMANN SYSTEM 3 satisfactory answer to the first question and a partial answer to the second one only in the linearized level. Since these two issues have been extensively studied in different contents such as the Boltzmann equation [4, 3], the Vlasov-Poisson-Boltzmann system [7, 6] and the Vlasov-Maxwell-Boltzmann system for two-species [5], our emphasize here will be put on the study of the weakly dissipative property of the electromagnetic field and its resulting slow time-decay rate of solutions, which are even different from the case of two-species. All details for their discussions are left to the next subsection. Let us begin with the Cauchy problem on the linearized non-homogeneous Vlasov- Maxwell-Boltzmann system, in the form of t u + x u M 1/ E = Lu + h, t E x B = M 1/ u d, (1.5) t B + x E = 0, E = M 1/ u d, x B = 0, with u(0, x, ) = u 0 (x, ), E(0, x) = E 0 (x), B(0, x) = B 0 (x), (1.6) where h = h(t, x, ) denotes a given non-homogeneous source term. For simplicity, we write U = [u, E, B], U 0 = [u 0, E 0, B 0 ]. Moreover, U 0 = [u 0, E 0, B 0 ] is always supposed to satisfy the last equation of (1.5) for t = 0. The first result, concerning the naturally existing energy functional and its optimal dissipation rate, is stated as follows. Here and hereafter, ν = ν() and P are defined in (.3) and (.5), respectively; see Subsection 1.3 for more notations used in this paper. Theorem 1.1. Let N 3. Assume ν 1/ h L (HN x ) with Ph(t, x, ) = 0. Define the temporal L -energy functional EN lin (U(t)) and its dissipation rate Dlin N (U(t)) by E lin N (U(t)) u(t) L (HN x ) + [E(t), B(t)] H N x, (1.7) D lin N (U(t)) = ν 1/ {I P}u(t) L (HN x ) + xpu(t) L (HN 1 x ) + x E(t) H N x + xb(t). Hx N 3 (1.8) Then, for any smooth solution U = [u, E, B] of the Cauchy problem (1.5)-(1.6) belonging to L (HN x ) Hx N Hx N, there exists a continuous functional EN lin (U(t)) given in (3.1) such that d dt E N lin (U(t)) + λdn lin (U(t)) C ν 1/ h(t) L (HN x ) (1.9) for any t 0. Remark 1.1. The above theorem shows the precise dissipative property of the naturally existing L -energy functional u(t) L (HN x ) + [E(t), B(t)] H N x
4 4 R.-J. DUAN for a properly large N. The construction of the equivalent energy functional EN lin(u(t)) is used to capture the optimal dissipation rate DN lin (U(t)), which will also be revisited in Corollary 1.4. This is different from the case of the two-species Vlasov-Maxwell- Boltzmann system as in [5], where the dissipation rate is a little stronger due to a cancelation phenomenon between two species which was firstly observed in [17]. Specifically, if it is here supposed that h = 0 and EN lin(u 0) is finite, then not only all macroscopic quantities Pu, E, B and the highest-order derivative N x [E, B] of the electromagnetic field lose their time-space integrability, but also the same thing happens to the second-order derivative xb of the magnetic field. On the other hand, for the energy space with m-order spatial regularity for any integer m 0, its optimal dissipation rate can be described by Dm lin (U(t)) once again from Corollary 1.4. Finally, as seen from Theorem 1.5 later on, this kind of weaker dissipation property leads to some slower time-decay rates of solutions. The second result about some time-frequency functional and its optimal dissipation rate is stated as follows. Theorem 1.. Assume ν 1/ĥ L for t 0, k R3, and Ph = 0. Define the L time-frequency functional E(Û(t, k)) and its dissipation rate D(Û(t, k)) by E(Û(t, k)) û L + [Ê, ˆB], (1.10) D(Û(t, k)) = ν1/ {I P}û L + k Ê + k k 4 + (1 + k ) Ê + (1 + k ) 3 ˆB. k 1 + k [â, ˆb, ĉ] (1.11) Then, for any solution U = [u, E, B] of the Cauchy problem (1.5) and (1.6) satisfying that û L + [Ê, ˆB] is finite for t 0 and k, there indeed exists a time-frequency functional E(Û(t, k)) given in (3.30) such that t E(Û(t, k)) + λd(û(t, k)) C ν 1/ ĥ L (1.1) for any t 0 and k. Corollary 1.3. Under Theorem 1., E(Û(t, k)) further satisfies t E(Û(t, k)) + λ k 4 (1 + k ) 3 E(Û(t, k)) C ν 1/ ĥ L (1.13) for any t 0 and k. Remark 1.. Theorem 1. and Corollary 1.3 show that the magnetic field B bears the weakest dissipation property among all quantities {I P}u, a, b, c, E and B, and even the dissipation of B here is much weaker than that in the case of twospecies Vlasov-Maxwell-Boltzmann system as in [5]. This is also consistent with what has been mentioned in Remark 1.1. Corollary 1.4. Let m 0 be an integer. Assume ν 1/ m x h L x, with Ph(t, x, ) = 0. Define the L energy functional Em lin (U(t)) and its dissipation rate Dm lin (U(t)) by Em lin (U(t)) m x u(t) + m x [E(t), B(t)], (1.14) Dm lin (U(t)) = ν 1/ m x {I P}u(t) (1.15) + m x a + 1+m x x 1 [a, b, c] + 1+m x x E + +m x x 3 B.
5 DISSIPATIVE PROPERTY OF VLASOV-MAXWELL-BOLTZMANN SYSTEM 5 Then, for any smooth solution U = [u, E, B] of the Cauchy problem (1.5) and (1.6) whose m-order spatial derivative belongs to L x, L L, there indeed exists a continuous functional E lin m (U(t)) given in (3.31) such that d dt E lin m (U(t)) + λd lin m (U(t)) C ν 1/ m x h(t) (1.16) for any t 0. Remark 1.3. It is straightforward to observe that Corollary 1.4 implies Theorem 1.1 by defining E lin N (U(t)) = N m=0 E lin m (U(t)), D lin N (U(t)) = N m=0 D lin m (U(t)), and using and likewise N N x 1+m x 1 = x x 1 x m x x N 1, m=0 m=0 N x 1+m x x x N, m=0 N x +m x 3 x x N 3. m=0 Notice that the above identities and equivalent relations can be verified with respect to the frequency variable under the Fourier transform. On the other hand, it is also interesting to see that even when 0 m < N, Dm lin (U(t)) can capture the optimal dissipation rate of the naturally existing m-order energy functional Em lin (U(t)). For instance, when m = 0, the direct energy estimate on system (1.5) produces the only dissipation of the microscopic component {I P}u, which is partially contained in the optimal form D0 lin (U(t)). Furthermore, we can obtain the large-time behavior of solutions to the linearized non-homogeneous Cauchy problem. Formally, the solution to the Cauchy problem (1.5) and (1.6) is denoted by the summation of two parts, U(t) = U I (t) + U II (t), (1.17) U I (t) = A(t)U 0, U I = [u I, E I, B I ], (1.18) U II (t) = t 0 A(t s)[h(s), 0, 0]ds, U II = [u II, E II, B II ], (1.19) where A(t) is the linear solution operator for the Cauchy problem on the linearized homogeneous system corresponding to (1.5) with h = 0. Notice that U II (t) is welldefined because [h(s), 0, 0] for any 0 s t satisfies the last equation of (1.5) due to the fact that Ph(s) = 0 and hence M 1/ h(s)d = 0. For brevity, we introduce the norms H m, Zr with m 0 and r 1 given by U H m = u L (Hm x ) + [E, B] H m x, U Z r = u Zr + [E, B] L r x,
6 6 R.-J. DUAN for U = [u, E, B], and note Z = H 0. Then, the third result to describe the time-decay property of the linearized solution is stated as follows. Theorem 1.5. Let 1 p, r q, σ 0, and let m 0 be an integer. Suppose Ph = 0. Let U be defined in (1.17), (1.18) and (1.19) as the solution to the Cauchy problem (1.5)-(1.6). Then, the first part U I corresponding to the solution of the linearized homogeneous system satisfies m x U I (t) Zq C(1 + t) 3 4 ( 1 p 1 q ) m 4 U0 Zp + C(1 + t) σ m+[σ+3( 1 r 1 q )]+ x U 0 Zr, (1.0) for any t 0, and the second part U II corresponding to the solution of the linearized nonhomogeneous system with vanishing initial data satisfies m x U II (t) Z C t 0 (1 + t s) 3 ( 1 r 1 ) m ν 1/ h(s) Z p ds + C t 0 (1 + t s) σ ν 1/ m+[σ]+ x h(s) ds, (1.1) for any t 0. Here, [ ] + is defined by [σ + 3( 1 r 1 σ if σ is integer and r = q =, q )] + = [σ + 3( 1 r 1 q )] + 1 otherwise, (1.) where [ ] means the integer part of the nonnegative argument. Finally, let us go back to the Cauchy problem on the nonlinear Vlasov-Maxwell- Boltzmann system. The global existence and uniqueness of solutions are stated as follows. Theorem 1.6. Let N 4. Define L energy functional E N (U(t)) and its dissipation rate D N (U(t)) by E N (U(t)) u(t) H N x, + [E(t), B(t)] H N x, (1.3) D N (U(t)) = ν 1/ {I P}u(t) H N x, + ν 1/ x u(t) L (HN 1 x ) + x E(t) H N x + xb(t). Hx N 3 (1.4) Suppose f 0 = M+M 1/ u 0 0. There indeed exists a continuous functional E N (U(t)) given in (4.7) or (4.8) such that if initial data U 0 = [u 0, E 0, B 0 ] satisfies (1.3) 4 for t = 0 and E N (U 0 ) is sufficiently small, then the nonlinear Cauchy problem (1.3)-(1.4) admits a global solution U = [u, E, B] satisfying and f(t, x, ) = M + M 1/ u(t, x, ) 0, [u(t), E(t), B(t)] C([0, ); H N x, H N x H N x ), E N (U(t)) + λ t 0 D N (U(s))ds E N (U 0 ) for any t 0. The decay rate of the solution obtained in Theorem 1.6 remains open. We shall discuss at the end of this paper the main difficulty of extending the linear decay property in Theorem 1.5 to the time-decay of the nonlinear system.
7 DISSIPATIVE PROPERTY OF VLASOV-MAXWELL-BOLTZMANN SYSTEM Related work and key points in the proof. As mentioned before, the method of constructing energy functionals or time-frequency functionals to deal with the global existence and time-decay estimates presented in this paper has been also extensively applied in [4, 7, 6, 3, 5]. Specifically, [4] is a starting point of these series of work. In [4], some interactive energy functionals were constructed to consider the dissipation of the macroscopic part of the solution and also the global existence of solutions without any initial layer was proved. Later, the same thing was done in [7] for the one-species Vlasov-Poisson-Boltzmann system, where the additional efforts are made to take care of the coupling effect from the self-consistent potential force through the Poisson equation. In order to investigate the optimal rate of convergence for the one-species Vlasov-Poisson-Boltzmann system, a new method on the basis of the linearized Fourier analysis was developed in [6] to study the time-decay property of the linear solution operator, where the key point is again to construct some proper time-frequency functionals so as to capture the optimal dissipation rate of the system. At the same time, [3] provided another method to study the exponential time-decay for the linear Boltzmann equation with a confining force by using the operator calculations instead of the Fourier analysis. Recently, following a combination of [6] and [17], the optimal large-time behavior of the two-species Vlasov-Maxwell-Boltzmann system was analyzed in [5]. The main finding in [5] is that although the non-homogeneous Maxwell system conserves the energy of the electromagnetic field, the coupling of the Boltzmann equation with the Maxwell system can generate some weak dissipation of the electromagnetic field which is actually of the regularity-loss type. It should be pointed out that even though the form of two-species Vlasov-Maxwell-Boltzmann system looks more complicated than that of the case of one-species, the study of global existence and time-decay rate is much more delicate in the case of one-species because the coupling term in the source of the Maxwell system M 1/ u(t, x, )d = b(t, x) corresponds to the momentum component of the macroscopic part of the solution which is degenerate with respect to the linearized operator L. Essentially, it is this kind of the macroscopic coupling feature that leads to some different dissipation properties between two-species and one-species for the Vlasov-Maxwell-Boltzmann system. For the convenience of readers, let us list a table below to present in a clear way similarity and difference of the dissipative and time-decay properties for three models: Boltzmann equation (BE), Vlasov-Poisson-Boltzmann system (VPB) and Vlasov-Maxwell-Boltzmann system (VMB). In Table 1, 1-s means one-species and -s two-species. Corresponding to different models, E(t, k) stands for some time-frequency functional equivalent with the naturally existing one and D(t, k) denotes the optimal dissipation rate of E(t, k) satisfying d E(t, k) + λd(t, k) 0 dt for all t 0 and k. All estimates are written for the linearized homogeneous equation or system. For more details and proof of other models in the above table, interested readers can refer to [3, 6, 5]. Notice that the 1-s VMB system decays faster than the 1-s VPB system due to the choice of initial data, that is, E 0 L 1 x is assumed in the case of 1-s VMB system, whereas in the case of 1-s VPB system, the
8 8 R.-J. DUAN electric field E 0 = x φ 0 with the potential force φ 0 satisfying the Poisson equation x φ 0 = M 1/ u 0 d may not belong to L 1 x under the assumption of u 0 Z 1 L. We remark that the decay rate (1 + t) 1/4 for the 1-s VPB system can be improved to be (1+t) 3/4 provided that x φ 0 L 1 x is additionally supposed. Finally, it should also be pointed out that the method developed in [3, 6, 5] and this paper could provide a good tool to deal with similar studies for other physical models with the structure involving not only the free transport operator but also the degenerately dissipative operator; see also [0]. E(t, k) D(t, k) = u(t) L BE û L ν 1/ {I P}û L + k 1+ k [â, ˆb, C(1 + t) 3 ĉ] 4 u 0 Z1 L 1-s VPB û L + â k ν 1/ {I P}û L + k 1+ k [â, ˆb, C(1 + t) 1 ĉ] + â 4 u 0 Z1 L 1-s VMB û L ν1/ {I P}û L + [Ê, ˆB] + k 1+ k [â, ˆb, ĉ] + k Ê + k (1+ k ) Ê + k 4 (1+ k ) ˆB 3 C(1 + t) 3 8 ( U 0 Z1 + x U 0 Z ) -s VMB û L ν1/ {I P}û L + [Ê, ˆB] + k 1+ k [â ±, ˆb, ĉ] + k Ê + k (1+ k ) ( Ê + ˆB ) C(1 + t) 3 4 ( U 0 Z1 + xu 0 Z ) Table 1.1 Dissipative and time-decay properties of different models Since the current work is a further development in the study of the Vlasov- Maxwell-Boltzmann system as in [5], we omit the detailed literature review for brevity, and readers can refer to [5] and reference therein. Here, we only mention some of them. The spectral analysis and global existence for the Boltzmann equation with near-equilibrium initial data was given by [19]. For the same topic, thirteen moments method and global existence was found by [13]. The energy method of the Boltzmann equation was developed independently in [9, 10, 11] and [15, 14, 1] by using the different macro-micro decomposition. The almost exponential rate of convergence of the Boltzmann equation on torus for large initial data was obtained in [] under some additional regularity assumption. [18] provided a very simple proof of [] in the framework of small perturbation. The diffusive limit of the two-species Vlasov- Maxwell-Boltzmann system over the torus was studied in [8]. It could be interesting to consider the same issue as in [8] for the one-species Vlasov-Maxwell-Boltzmann system because of its weaker dissipation property. In what follows, let us explain some new technical points in the proof of our main results which are different from previous work. The first one is about the dissipation estimate on the momentum component b(t, x) in the macroscopic part Pu in Theorem 1.1. Recall that in the previous work [7] and [3], the dissipation estimate of b(t, x)
9 DISSIPATIVE PROPERTY OF VLASOV-MAXWELL-BOLTZMANN SYSTEM 9 was based on (.7) 1 and (.6). This fails in the case of one-species Vlasov-Maxwell- Boltzmann system because it is impossible to control a term such as δ ij c( i E j + j E i )dx. ij Instead, the right way is to make estimates on (.9) 4. Therefore, one can get the macroscopic dissipation estimate (3.8) with the dissipation of E multiplied by a small constant on the right-hand side. The second key point is about the dissipation of the electromagnetic field. Different from [5], although there is no cancelation in the case of one-species, one can still design some interactive energy functional E lin, N (U(t)) to capture the weaker dissipation rate α E + α B. 1 The third key point is about the time-decay estimate on the linearized solution operator. In fact, for a general frequency function φ(k) which is of the regularity-loss type as in Lemma 3.1, one can repeatedly apply the Minkowski inequality to make interchanges between frequency and space variables so that the more general L p -L q type time-decay than in Theorem 1.5 can be obtained. The last key point is about the estimate on the nonlinear term Eu in the proof of Lemma 4.1 concerning the a priori estimates of solutions. In particular, it is impossible to bound E(Pu) dxd by using E N (U(t)) 1/ D N (U(t)) up to a constant since both E and Pu do not enter into the dissipation rate D N (U(t)) given in (1.4). Instead, we first take the velocity integration and then use the macroscopic balance laws (.6) so as to obtain an estimate as E(Pu) dxd d dt [ ] b (a + c) dx + C E N (U(t)) 1/ + E N (U(t)) D N (U(t)). We remark that this also has been observed in [7] in the study of the one-species Vlasov-Poisson-Boltzmann system Notations. Throughout this paper, C denotes some positive (generally large) constant and λ denotes some positive (generally small) constant, where both C and λ may take different values in different places. In addition, A B means λa B 1 λa for a generic constant 0 < λ < 1. For any integer m 0, we use Hx, m, Hm x, H m to denote the usual Hilbert spaces H m ( x ), Hm ( x), H m ( ), respectively, and L, L x, L are used for the case when m = 0. When without confusion, we use H m to denote Hx m and use L to denote L x or L x,. For a Banach space X, X denotes the corresponding norm, while always denotes the norm L for simplicity. For r 1, we also define the standard time-space mixed Lebesgue space Z r = L (Lr x) = L ( ; Lr ( x)) with the norm g Zr = ( ( ) ) /r 1/ g(x, ) r dx d, g = g(x, ) Z r. For multi-indices α = [α 1, α, α 3 ] and β = [β 1, β, β 3 ], we denote α β = α1 x 1 α x α3 x 3 β1 1 β β3 3.
10 10 R.-J. DUAN The length of α is α = α 1 + α + α 3 and the length of β is β = β 1 + β + β 3. For simplicity, we also use j to denote xj for each j = 1,, 3. For an integrable function g : R, its Fourier transform is defined by ĝ(k) = Fg(k) = e πix k g(x)dx, x k := 3 x j k j, k, where i = 1 C is the imaginary unit. For two complex vectors a, b C 3, (a b) denotes the dot product of a with the complex conjugate of b over the complex field. x = (1 + x ) 1/ is defined in terms of the Fourier transform. The rest of this paper is arranged as follows. In Section we present some basic property of the linearized collision operator and derive some macroscopic moment equations. In Section 3 we study the linearized non-homogeneous Vlasov-Maxwell- Boltzmann system in order to prove Theorem 1.1, Theorem 1. and Theorem 1.5. Finally, we prove in Section 4 Theorem 1.6 for the global existence of solutions to the nonlinear Cauchy problem.. Moment equations. It is easy to see that Lu and Γ(u, u) are given by Lu = 1 [Q(M, Mu) + Q( ] Mu, M), (.1) M Γ(u, u) = 1 M Q( Mu, Mu). (.) For the linearized collision operator L, one has the following standard facts [1]. L can be split as Lu = ν()u + Ku, where the collision frequency is given by ν() = ( ) ω M( ) dωd. (.3) S Notice that ν() (1 + ) 1/. The null space of L is given by { N = span M 1/, i M 1/ (1 i 3), M 1/}. The linearized collision operator L is non-positive and further L is known to be locally coercive in the sense that there is a constant λ 0 > 0 such that [1]: ulu d λ 0 ν() {I P}u d, (.4) where, for fixed (t, x), P denotes the orthogonal projection from L to N. Given any u(t, x, ), one can write P in (.4) as Pu = {a(t, x) + b(t, x) + c(t, x)( 3)}M 1/. (.5) Since P is a projection, the coefficient functions a(t, x), b(t, x) [b 1 (t, x), b (t, x), b 3 (t, x)] and c(t, x) depend on u(t, x, ) in terms of a = M 1/ ud = M 1/ Pud, R 3 R 3 b i = i M 1/ ud = i M 1/ Pud, 1 i 3, c = 1 ( 3)M 1/ ud = 1 ( 3)M 1/ Pud. 6 R 6 3 j=1
11 DISSIPATIVE PROPERTY OF VLASOV-MAXWELL-BOLTZMANN SYSTEM 11 To derive evolution equations of a, b and c, we start from the local balance laws of the original system (1.1) to obtain t fd + x fd = 0, R 3 R 3 t fd + x fd E fd fd B = 0, R t R fd + x 3 R fd E fd = 0. 3 The above system implies t a + x b = 0, t b + x (a + c) + x M 1/ {I P}ud E(1 + a) b B = 0, t c x b x M 1/ {I P}udx 1 E b = 0. 3 As in [6], define Θ ij (u) = ( i j 1)M 1/ ud, Λ i (u) = 1 ( 5) i M 1/ ud R 10 3 for 1 i, j 3. Applying them to the first equation of (1.3), one has { t [Θ ij ({I P}u) + cδ ij ] + i b j + j b i = Θ ij (l + g), t Λ i ({I P}u) + i c = Λ i (l + g), (.6) (.7) where δ ij means Kronecker delta, and l = x {I P}u + Lu, g = 1 Eu (E + B) u + Γ(u, u). (.8) One can replace t c in (.7) 1 by using (.6) 3 so that t Θ ij ({I P}u) + i b j + j b i 3 δ ij x b 10 3 δ ij x Λ({I P}u) = Θ ij (l + g) 3 δ ije b. In a summary, we obtained the following moment system t a + x b = 0, t b + x (a + c) + x Θ({I P}u) E = Ea + b B, t c x b x Λ({I P}u) = 1 3 E b, t Θ ij ({I P}u) + i b j + j b i 3 δ ij x b 10 3 δ ij x Λ({I P}u) = Θ ij (l + g) 3 δ ije b, 1 i, j 3, t Λ i ({I P}u) + i c = Λ i (l + g), 1 i 3. (.9)
12 1 R.-J. DUAN On the other hand, the Maxwell system in (1.4) is equivalent with t E x B = b, t B + x E = 0, x E = 0, x B = 0. Finally, we should point out that the key analysis of all results in this paper is based on the above moment equations coupled with the Maxwell system. 3. Linear non-homogeneous system. In this section we consider the Cauchy problem (1.5) and (1.6) on the linearized Vlasov-Maxwell-Boltzmann system. For convenience of readers, recall it by t u + x u M 1/ E = Lu + h, t E x B = M 1/ ud, (3.1) t B + x E = 0, E = M 1/ ud, x B = 0, with u(0, x, ) = u 0 (x, ), E(0, x) = E 0 (x), B(0, x) = B 0 (x). (3.) Here, h = h(t, x, ) is a given non-homogenous source term, satisfying Ph = L energy functional and its optimal dissipation rate. In this subsection we shall prove Theorem 1.1. Before that, similar to obtain (.9), one can also derive the following moment equations corresponding to the linear equation (3.1): t a + x b = 0, t b + x (a + c) + x Θ({I P}u) E = 0, t c x b x Λ({I P}u) = 0, t Θ ij ({I P}u) + i b j + j b i 3 δ ij x b 10 3 δ ij x Λ({I P}u) = Θ ij (l + h), t Λ i ({I P}u) + i c = Λ i (l + h), (3.3) where 1 i, j 3, and as in (.8), l still denotes l = x {I P}u + Lu. The Maxwell system also takes the form of t E x B = b, t B + x E = 0, x E = 0, x B = 0. (3.4) Proof of Theorem 1.1: Let N 3. First of all, a usual energy estimate on (3.1) gives 1 d dt ( α u + α [E, B] )+λ ν 1/ α {I P}u C ν 1/ α g. α N α N α N (3.5)
13 DISSIPATIVE PROPERTY OF VLASOV-MAXWELL-BOLTZMANN SYSTEM 13 As in [3, 7] or [5], one can further deduce from (3.3) the dissipation of a, b and c. In fact, let ɛ 1 > 0, ɛ > 0 be arbitrary constants to be chosen later. From (3.3) 5 and (3.3) 3, it follows that d x α c Λ( α {I P}u)dx + λ x α c dt ɛ 1 x α b + C α {I P}u + ν 1/ α g. ɛ 1 α N (3.3) 4 and (3.3) imply d dt 3 ij=1 + λ ( i α b j + j α b i δ ij x α b)θ ij ( α {I P}u)dx 3 x α b ɛ x α [a, c] + x α E + C α {I P}u + ɛ α N α N ν 1/ α g. It holds from (3.3) and (3.3) 1 that d dt x α a α bdx + λ C α N α a x α [b, c] + C x α {I P}u. Define E lin,1 N (U(t)) = + ij=1 + κ 1 x α c Λ( α {I P}u)dx 3 ( i α b j + j α b i 3 δ ij x α b)θ ij ( α {I P}u)dx x α a α bdx, (3.6) for some small constant κ 1 > 0. Therefore, by taking κ 1 > 0 small enough and then letting ɛ 1 > 0 and ɛ > 0 be small enough, the dissipation of a, b and c can be obtained
14 14 R.-J. DUAN by the following inequality d dt E lin,1 N (U(t)) + λ ɛ α N x α [a, b, c] + a x α E + C α {I P}u + ɛ α N ν 1/ α g, where ɛ > 0 is still left to be chosen later on. The key part is to estimate the dissipation of E and B. We claim that d dt 1 ɛ 3 α E α bdx + λ 1 α N + C with ɛ 3 > 0 to be chosen, and d dt 1 α N 1 α x B + C ɛ 3 1 α x B α Edx + λ 1 α N α E x α b (3.7) ( x α [a, c] + x α {I P}u ), (3.8) 1 α N α x E + C α x B 1 α N α b. (3.9) For this time, suppose that (3.8) and (3.9) hold true. Define E lin, N (U(t)) = α E α bdx κ α x B α Edx 1 1 α N (3.10) for some constant κ > 0. Then, by taking κ > 0 and further ɛ 3 > 0 both small enough, it follows from (3.8) and (3.9) that d dt E lin, N (U(t)) + λ 1 C x b + C α E + λ 1 α B ( x α [a, b, c] + x α {I P}u ), (3.11) which is the desired dissipation estimate of the electromagnetic field E and B. Now, define EN lin (U(t)) = ( ) ( α u + α [E, B] )+κ 4 E lin,1 N (U(t)) + κ 3E lin, N (U(t)) (3.1) α N with constants κ 3 > 0 and κ 4 > 0, where E lin,1 lin, N (U(t)), EN (U(t)) are defined in (3.6) and (3.10). In the same way as before, by taking properly small constants κ 3 > 0,
15 DISSIPATIVE PROPERTY OF VLASOV-MAXWELL-BOLTZMANN SYSTEM 15 ɛ > 0 and κ 4 > 0 in turn, (1.9) follows from the linear combination of (3.5), (3.7) and (3.11). Moreover, it is easy to verify that EN lin(u(t)) is the desired L energy functional satisfying (1.7) and DN lin (U(t)) is given by (1.8). Here, one has to check α B C α x B. (3.13) 1 α N In fact, by taking α with α N 1 and using B = 1 x x x B due to x B = 0, it holds α B = α 1 x x x B = α γi i 1 x x x B, for some γ i (1 i 3) with γ i = 1. Since i 1 x j for any 1 i, j 3 is a bounded operator from L p to itself with 1 < p <, α B C α γi x B. Hence, (3.13) follows from taking summation of the above inequality over α N 1. Now, the rest is to prove (3.8) and (3.9). Take α with 1 α N 1. By using (3.3) to replace E and then using (3.4) 1 to replace t E, one can compute α E = α E α Edx (3.14) R 3 = α E α [ t b + x (a + c) + x Θ({I P}u)]dx = d α E α bdx α t E α bdx dt R 3 + α E α [ x (a + c) + x Θ({I P}u)]dx = d α E α bdx + α (b x B) α bdx dt R 3 + α E α [ x (a + c) + x Θ({I P}u)]dx. Then, (3.8) follows from the above identity after taking summation over 1 α N 1 and further applying the Cauchy-Schwarz inequality and integration by parts. In fact, it suffices to consider the second term on the r.h.s. of (3.14). It can be estimated by 1 α (b x B) α bdx = α x B α bdx + 1 α N ɛ 3 1 α N 1 α =N 1 α b α γi x B α+γi bdx α x B + C ɛ 3 x α b, where as before 0 < ɛ 3 1 is small to be chosen, and γ i denotes a multi-index with γ i = 1 for some 1 i 3. To prove (3.9), take α with 1 α N. By using
16 16 R.-J. DUAN (3.4) 1 to replace x B and then using (3.4) to replace t B, one has α x B = α x B α x Bdx = α x B α ( t E + b)dx = d α x B α Edx α x t B α Edx + α x B α bdx dt = d α x B α Edx + α x x E α Edx + α x B α bdx. dt Then, (3.9) follows by applying integration by part to the right-hand second term of the above identity, using the Cauchy-Schwarz inequality and then taking summation over 1 α N. The proof of Theorem 1.1 is complete. 3.. L time-frequency functional and its optimal dissipation rate. In this subsection we shall prove Theorem 1. as well as Corollary 1.3 and Corollary 1.4. For that, we need to consider the solution U = [u, E, B] to the Cauchy problem (3.1)-(3.) in the Fourier space k. By taking the Fourier transform in x from (3.1) 1, (3.3) and (3.4), one has t û + i kû M 1/ Ê = Lû + ĥ, (3.15) t â + ik ˆb = 0, tˆb + ik(â + ĉ) + ik Θ({I P}û) Ê = 0, t ĉ ik ˆb + 5 ik Λ({I P}û) = 0, 3 t Θ ij ({I P}û) + ik iˆbj + ik jˆbi 3 δ ijik ˆb 10 3 δ ijik Λ({I P}û) = Θ ij (ˆl + ĥ), t Λ i ({I P}û) + ik i ĉ = Λ i (ˆl + ĥ), (3.16) and where ˆl is given by t Ê ik ˆB = ˆb, t ˆB + ik Ê = 0, ik Ê = â, ik ˆB = 0, ˆl = ik {I P}û + Lû. (3.17) These equations above are ones to be used through this subsection. Proof of Theorem 1.: It is similar to the proof of Theorem 1.1. The difference is that all calculations are made in the Fourier space. Thus, some details in the following proof will be omitted for simplicity. First of all, as in [3], on one hand, from (3.15) and (3.17), one has 1 t( û L + [Ê, ˆB] ) + λ ν 1/ {I P}û L C ν 1/ ĝ L, (3.18)
17 DISSIPATIVE PROPERTY OF VLASOV-MAXWELL-BOLTZMANN SYSTEM 17 and on the other hand, from (3.16), the following three estimates hold true: t R(ikĉ Λ({I P}û) + λ k ĉ ɛ 1 k ˆb + C (1 + k ) {I P}û L + C ν 1/ ĝ ɛ 1 L (3.19), and t R 3 (ik iˆbj + ik jˆbi 3 δ ijik ˆb Θ ij ({I P}û) + λ k ˆb ij=1 ɛ k [â, ĉ] k + ɛ 1 + k Ê + C (1 + k ) {I P}û L + C ν 1/ ĝ ɛ L, (3.0) t R(ikâ ˆb) + λ(1 + k ) â k ˆb + C k ĉ + C k {I P}û L, (3.1) where constants 0 < ɛ 1, ɛ 1 are to be chosen. Define E lin,1 (Û(t)) = 1 R{(ikĉ Λ({I P}û) 1 + k 3 + (ik iˆbj + ik jˆbi 3 δ ijik ˆb Θ ij ({I P}û) + κ 1 (ikâ ˆb)} (3.) ij=1 for a constant κ 1 > 0. One can take κ 1 and then ɛ 1 both small enough such that the sum of (3.19), (3.0) and κ 1 (3.1) gives lin,1 t E (Û(t)) + λ k 1 + k [â, ˆb, ĉ] + â k ɛ (1 + k ) Ê + C ) ( {I P}û L + ν 1/ ĝ ɛ L. (3.3) For estimates on the dissipation of Ê and ˆB, it is straightforward to deduce from (3.15) and (3.16) the following two identities and t (k Ê k ˆb) + k Ê = k ˆb (ik k ˆB k ˆb) + (k Ê ik (k Θ({I P}û))), (3.4) t (ik ˆB Ê) + k ˆB = k Ê + (ik ˆB ˆb). (3.5) By applying the Cauchy-Schwarz to (3.5) and then multiplying it by k /(1 + k ) 3, one has t k (1 + k ) 3 R(ik ˆB Ê) + λ k k ˆB (1 + k ) 3 k k Ê (1 + k ) 3 + C k ˆb (1 + k ) 3,
18 18 R.-J. DUAN which implies k t (1 + k ) 3 R(ik ˆB Ê) + k ˆB k λ k Ê (1 + k ) 3 (1 + k ) + C k ˆb 1 + k. (3.6) Similarly, after dividing (3.4) by (1 + k ) and then using Cauchy-Schwarz, tr(k Ê k ˆb) k Ê (1 + k ) + λ (1 + k ) (3.7) k ˆb (1 + k ) + k k ˆB k ˆb (1 + k ) + C k 4 Θ({I P}û) (1 + k ) k k ɛ ˆB 3 (1 + k ) 3 + C k ɛ k ˆb + C {I P}û L, where we used the inequality k k ˆB k ˆb (1 + k ) ɛ 3 k k ˆB (1 + k ) 3 + C ɛ 3 k ˆb 1 + k for an arbitrary constant 0 < ɛ 3 1. Then, in terms of (3.6) and (3.7), let us define lin, Ê k ˆb) k R(ik E (Û(t)) = R(k (1 + k ) κ ˆB Ê) (1 + k ) 3, (3.8) where κ > 0 is chosen small enough such that lin, k Ê t E (Û(t)) + λ (1 + k ) + λ k 4 ˆB (1 + k ) 3 C k ˆb 1 + k + C {I P}û L. (3.9) Here, we used k ˆB = k k ˆB due to k ˆB = 0. Now, in terms of (3.18), (3.3) with 0 < ɛ 1 and (3.9), we define E lin (Û(t)) = û L + [Ê, ˆB] ( ) + κ 4 E lin,1 (Û(t)) + κ 3E lin, (Û(t)), (3.30) where E lin,1 (Û(t)), E lin, (Û(t)) are denoted by (3.) and (3.8), and κ 3 > 0, ɛ > 0 and κ 4 > 0 are chosen in turn small enough such that (1.10) and (1.1) hold true and lin D (Û(t)) is given by (1.11). Here, notice that we used k Ê = k Ê + k Ê = â + k Ê. The proof of Theorem 1. is complete. Proof of Corollary 1.3 and Corollary 1.4: First, (1.13) in Corollary 1.3 immediately results from (1.1) by noticing D(Û(t, k)) λ k 4 (1 + k E(Û(t, k)) ) 3 due to the definitions (1.10), (1.11) of E(Û(t, k)), D(Û(t, k)). For Corollary 1.4, it suffices to define Em lin (Û(t, k)) = α =m E( k xu(t, k)) = α =m E((ik) α Û(t, k)). (3.31) Since the system (3.1) satisfied by U = [u, E, B] is linear, it is easy to see from Theorem 1. that Em lin (Û(t, k)) satisfies (1.14) and (1.16) with Dlin m (Û(t, k)) given by (1.15). This hence completes the proof of Corollary 1.4.
19 DISSIPATIVE PROPERTY OF VLASOV-MAXWELL-BOLTZMANN SYSTEM Time-decay estimates. In this subsection we turn to the proof of Theorem 1.5. Theorem 1.5 actually follows from (1.13) in Corollary 1.3. However, we would rather provide a similar much more general result. Recall the definition of A(t) as in (1.18). In what follows, for the Fourier transform Û(t, k) of U = [u, E, B], we set Û(t, k) = û L + [Ê, ˆB] (3.3) for simplicity. Lemma 3.1. Assume that for any initial data U 0, the linear homogeneous solution U(t) = A(t)U 0 obeys the pointwise estimate Û(t, k) Ce φ(k)t Û0(k) (3.33) for all t 0, k, where φ(k) is a strictly positive, continuous and real-valued function over k and satisfies φ(k) { O(1) k σ + as k 0, O(1) k σ as k, (3.34) for two constants σ > σ + > 0. Let m 0 be an integer, 1 p, r q and σ 0. Then, U(t) = A(t)U 0 obeys the time-decay estimate m x U(t) Zq C(1 + t) 3 σ + ( 1 p 1 q ) m σ + U0 Zp + C(1 + t) σ σ m+[σ+3( 1 r 1 q )]+ x U 0 Zr, (3.35) for any t 0, where [ ] + is defined in (1.). Proof. Take a constant R > 0. From the assumptions on φ(k), it is easy to see φ(k) { λ k σ + if k R, λ k σ if k R. Take q and an integer m 0. From Hausdorff-Young inequality, m x U(t) Zq C k m e φ(k)t Û 0 Zq C k m e λ k σ + t Û 0 + C k Zq m e λ k σ t Û 0 := I 1 + I, ( k R) Zq ( k R) (3.36) where 1 q + 1 q = 1. For I 1, by the definition of the norm Zq, I 1 = C k m e λ k σ + tû L 0 + C k m e λ k σ + t [Ê0, ˆB 0 ] (L q ( k R)) Here, note that since 1 q, from the Minkowski inequality, k m e λ k σ + tû L 0 (L q ( k R)) k m e λ k σ + tû 0 L q ( k R;L ) = L q ( k R) k m e λ k σ + t L û 0 L. q ( k R).
20 0 R.-J. DUAN Hence, we arrive at I 1 C k m e λ k σ + t [ û 0 L, Ê0, ˆB 0 ]. L q ( k R) Take 1 p. Further using the Hölder inequality for 1 q = p q p q + 1 p with p given by 1 p + 1 p = 1, I 1 C k m e λ k σ + t p q L p q ( k R) [ û 0 L, Ê0, ˆB 0 ] L p ( k R) Here, the right-hand first term can be estimated in a standard way [1] as k m e λ k σ + t p q C(1 + t) 3 σ ( 1 + p 1 q ) m σ + L p q ( k R) by using change of variable kt 1 σ + k, and the right-hand second term is estimated by Minkowski and Hausdorff-Young inequalities as [ û 0 L, Ê0, ˆB [ 0 ] = û 0 L p L p ( k R;L ( k R) ) + Ê 0, ˆB 0 ] [ û 0 L (L p ( k R)) + Ê 0, ˆB 0 ] L p ( k R) L p ( k R) C( u 0 L (L p x) + [E 0, B 0 ] L p) = C U 0 Zp, where the Minkowski inequality was validly used due to p. Therefore, for I 1, one has I 1 C(1 + t) 3 σ + ( 1 p 1 q ) m σ + U0 Zp. To estimate I, take a constant σ 0 so that I = C k m e λ k σ t Û 0 C sup k Zq σ e λ k σ t k m+σ Û 0. ( k R) k R Zq ( k R) Here, the right-hand first term decays in time as sup k σ e λ k σ t C(1 + t) σ σ. k R We estimate the right-hand second term as follows. Take 1 r with 1 r + 1 r = 1 and take a constant ɛ > 0 small enough. Then, similarly as before, from Minkowski and Hölder inequalities for 1 q = r q + 1 r, one has r q k m+σ Û 0 k Zq m+σ [ û 0 L ( k R), Ê0, ˆB 0 ] L q ( k R) q L q r q r q r q [ û 0 L, Ê0, ˆB 0 ] k 3(1+ɛ) r r k m+σ+3(1+ɛ) r q ( k R) k m+[σ+3( C 1 r 1 q )]+ ɛ [ û 0 L, Ê0, ˆB 0 ] L r ( k R) L r ( k R).
21 DISSIPATIVE PROPERTY OF VLASOV-MAXWELL-BOLTZMANN SYSTEM 1 Here, by Minkowski inequality due to q once again and further by Hausdorff- Young inequality, k m+[σ+3( 1 r 1 q )]+ [ û 0 L, Ê0, ˆB 0 ] k m+[σ+3( 1 r 1 q )]+ L û 0 + (L r ( k R)) Then, it follows that Thus, I is estimated by L r ( k R) k m+[σ+3( 1 r 1 q )]+ [Ê0, ˆB 0 ] k m+σ Û 0 C Zq m+[σ+3( 1 r 1 q )]+ x U 0 Zr. ( k R) I C(1 + t) σ σ m+[σ+3( 1 r 1 q )]+ x U 0 Zr. L r ( k R) C m+[σ+3( 1 r 1 q )]+ x U 0 Zr. Now, (3.35) follows by plugging the estimates of I 1 and I into (3.36). This completes the proof of Lemma 3.1. Proof of Theorem 1.5: To prove (1.0), by letting h = 0 and using Corollary 1.3, E(Û(t, k)) λ k e (1+ k ) 3 t E(Û0(k)) 4 for any t 0 and k, where we have set U = U I for simplicity. Due to (1.10) and (3.3), E(Û(t, k)) Û(t, k) holds so that Û(t, k) λ k Ce (1+ k ) 3 t Û0(k) 4 for any t 0 and k. This shows that corresponding to (3.33) and (3.34) of Lemma 3.1, one has the special situation φ(k) = λ k 4 (1 + k ) 3 with σ + = 4, σ =. Thus, one can apply Lemma 3.1 to obtain (1.0) from (3.35). To prove (1.1), we let U 0 = 0, and also set U = U II for simplicity. Note that (1.13) implies Û(t, k) C t e 0 λ k 4 (1+ k ) 3 (t s) ν 1/ ĥ(s) L ds for any t 0 and k. One can again apply Lemma 3.1 with q = r = so that (1.1) follows. The proof of Theorem 1.5 is complete. 4. Nonlinear system. In this section we are concerned with the global existence of solutions to the Cauchy problem (1.3)-(1.4) of the reformulated nonlinear Vlasov- Maxwell-Boltzmann system. We first devote ourselves to the proof of some uniformin-time a priori estimates on the solution. In what follows, U = [u, E, B] is supposed
22 R.-J. DUAN to be smooth in all arguments and satisfy the system (1.3) over 0 t T for some 0 < T. Lemma 4.1. Under the assumption that sup 0 t T (a+c)(t) L is small enough, there is E N (U(t)) satisfying (1.3) such that d dt E N(U(t)) + λd N (U(t)) C(E N (U(t)) 1/ + E N (U(t)))D N (U(t)) (4.1) for any 0 t T, where D N (U(t)) is defined in (1.4). Proof. First, the zero-order energy estimate implies 1 d dt ( u + [E, B] b (a + c)dx + λ ν 1/ {I P}u C(E N (U(t)) 1/ + E N (U(t)))D N (U(t)). (4.) Here and hereafter, when E N (U(t)) occurs in the right-hand terms of inequalities, it means an equivalent energy functional satisfying (1.3) and its explicit representation will be determined later on. In fact, from the system (1.3), 1 d dt ( u + [E, B] ) + λ ν 1/ {I P}u uγ(u, u)dxd + 1 Eu dxd. (4.3) Here, for the first term on the r.h.s. of (4.3), it is a standard fact as in [9] or [10] that it is bounded by CE N (U(t)) 1/ D N (U(t)). Since all estimates on terms involving the nonlinear term Γ(u, u) in the following can be handled in the similar way, we shall omit the details of their proof for brevity. The right-hand second term of (4.3) can be estimated as in [7, Lemma 4.4]: 1 Eu dxd = 1 E Pu dxd + EPu{I P}udxd + 1 E {I P}u dxd. Here, it is easy to see EPu{I P}udxd + 1 E {I P}u dxd C E H ( x [a, b, c] + ν 1/ {I P}u ) CE N (U(t)) 1/ D N (U(t)). And, as in [7, Lemma 4.4], one can compute 1 E Pu dxd = E b(a + c)dx, where by replacing E by equation (.9), using integration by part in t and then replacing t (a + c) by equations (.9) 1 and (.9) 3, gives E b(a + c)dx = 1 d b (a + c)dx dt + b [ 5 6 x b x Λ({I P}u) 1 3 E b]dx + [ x (a + c) + x Θ({I P}u)] b(a + c)dx Ea b(a + c)dx b B b(a + c)dx.
23 DISSIPATIVE PROPERTY OF VLASOV-MAXWELL-BOLTZMANN SYSTEM 3 Note b B b = 0. Hence, it follows that 1 E Pu dxd 1 d b (a + c)dx dt + C [a, b, c, E] H 1( x [a, b, c] + x {I P}u ) + C E x a x [a, b, c] 1 d b (a + c)dx + C(E N (U(t)) 1/ + E N (U(t)))D N (U(t)). dt Collecting the above estimates and putting them into (4.3) proves (4.). Next, for the estimates on all derivatives including the pure spatial derivatives and space-velocity mixed derivatives, one has and 1 d dt 1 d dt 1 α N N k=1 C k ( α u + α [E, B] ) + λ β =k α + β N CE N (U(t)) 1/ D N (U(t))+C α β {I P}u + λ 1 α N β 1 α + β N ν 1/ α {I P}u CE N (U(t)) 1/ D N (U(t)), (4.4) ν 1/ α β {I P}u α x [a, b, c] +C α N ν 1/ α {I P}u, (4.5) where C k (1 k N) are strictly positive constants. Since the nonlinear term g takes the form as in (.8), the proof of (4.4) and (4.5) is almost the same as in [7] for the case of the Vlasov-Poisson-Boltzmann system and thus details are omitted for brevity. Finally, the key step is to obtain the macroscopic dissipation of a, b, c and E, B for the nonlinear system (1.3). This is similar to the proof of Theorem 1.1 in Subsection 3.1 for the linearized system. Here, the additional efforts should be made to take care of all quadratically nonlinear terms in g defined by (.8). But, these nonlinear estimates once again are almost the same as in [7] for the case of the Vlasov-Poisson-Boltzmann system so we omit details for brevity. Thus, we have the following estimates. Recall the definitions (3.6) and (3.10) of two interactive functionals E lin,1 lin, N (U(t)) and EN (U(t)), where κ 1 > 0, κ > 0 in (3.6) and (3.10) are sufficiently small. It turns out that d ( dt E lin,1 N ) (U(t)) + κ 3E lin, N (U(t)) + λ α a + λ + λ 1 C α N α E + λ α N 1 α N α B α [b, c] α {I P}u + CE N (U(t))D N (U(t)), (4.6) where κ 3 > 0 is small enough. This is the desired estimate on the macroscopic dissipation.
24 4 R.-J. DUAN Now, we are in a position to prove (4.1). Let κ 4 > 0 be taken as in (3.1). Define E N (U(t)) = EN lin (U(t)) b (a + c)dx + κ 5 N C k k=1 β =k α + β N β α {I P}u (4.7), where κ 5 > 0 is a constant to be chosen. Here, notice that by recalling the definition (3.1) of E lin N (U(t)), E N(U(t)) can be rewritten as E N (U(t)) (4.8) = ( α u + α [E, B] ) b (a + c)dx α N +κ 4 x α c Λ( α {I P}u)dx + +κ 4 κ 3 +κ 5 N k=1 1 i,j 3 +κ 1 C k 1 β =k α + β N ( i α b j + j α b i 3 δ ij x α b)θ ij ( α {I P}u)dx x α a α bdx R 3 α E α bdx κ α β {I P}u. 1 α N α x B α Edx R 3 Due to smallness of sup 0 t T (a+c)(t) L by the assumption, from (4.7) and (1.7), it is easy to see which further implies E N (U(t)) E lin N (U(t)) + N k=1 C k β =k α + β N E N (U(t)) u(t) L (HN x ) + [E(t), B(t)] H + N N u(t) H N x, + [E(t), B(t)] H N. k=1 α β {I P}u, C k β =k α + β N α β {I P}u Thus, E N (U(t)) satisfies (1.3) for any κ 5 > 0. Moreover, by taking κ 5 > 0 small enough, the summation of (4.), (4.4), κ 4 (4.6) and then κ 5 (4.6) leads to (4.1) with D N (U(t)) defined in (1.4). This completes the proof of Lemma 4.1. Proof of Theorem 1.6: Let us consider the uniform-in-time a priori estimates of solutions under the smallness assumption that sup 0 t T ( ) u(t) H + [E(t), B(t)] x, N H δ N
25 DISSIPATIVE PROPERTY OF VLASOV-MAXWELL-BOLTZMANN SYSTEM 5 for a sufficiently small constant δ > 0. This smallness assumption also implies that sup 0 t T (a + c)(t) L is small enough since N 4. Thus, it follows from Lemma 4.1 that there are E N (U(t)), D N (U(t)) defined in (4.8) and (1.4) such that (1.3) holds true for E N (U(t)) and for any 0 t T, that is d dt E N(U(t)) + λd N (U(t)) C(δ 1/ + δ)d N (U(t)) E N (U(t)) + λ t 0 D N (U(s))ds E N (U 0 ) for any 0 t T, since δ > 0 is small enough. Now, the rest proof follows from the standard process by combining the above uniform-in-time a priori estimates with the local existence as well as the continuity argument as in [11] or [7] under the assumption that E N (U 0 ) is sufficiently small, and details are omitted for simplicity. The proof of Theorem 1.6 is complete. We conclude this paper with a discussion about the large-time behavior of solutions to the nonlinear system. Although Theorem 1.6 shows the global existence of close-to-equilibrium solutions to the Cauchy problem (1.3)-(1.4) of the nonlinear Vlasov-Maxwell-Boltzmann system, the decay rate of the obtained solution remains open. This issue has been studied in [5] for the case of two-species. However, the approach of [5] by applying the linear decay property together with the Duhamel s principle to the nonlinear system can not be applied to the case of one-species here. Let us explain a little the key difficulty in a formal way. In fact, the linear system in one-species decays as (1 + t) 3/8 which is slower than (1 + t) 3/4 in two-species as pointed out in Table 1. Thus, in one-species case, the quadratic nonlinear source decays as at most (1+t) 3/4, and if the Duhamel s principle was used, the time-integral term generated from the nonhomogeneous source decays as t 0 (1 + t s) 3 8 (1 + s) 3 4 ds C(1 + t) 1 8. Thus, the bootstrap argument breaks down and one can not expect the solution to decay as (1 + t) 3/8 in the nonlinear case. Therefore, the study of the large-time behavior for the one-species Vlasov-Maxwell-Boltzmann system becomes much more difficult than for the two-species case as in [5]. Acknowledgments. This work was supported partially by the Direct Grant 010/011 from CUHK and by the General Research Fund (Project No.: ) from RGC of Hong Kong. The author would like to thank Peter Markowich and Massimo Fornasier for their strong support when he stayed in RICAM, Austrian Academy of Sciences for the Post-doc study during , and thank Professors Tong Yang, Chang-Jiang Zhu and Hui-Jiang Zhao for their continuous encouragement. Some useful discussions from Robert M. Strain are much acknowledged. Thanks also go to the referee for useful suggestions to improve the presentation of the paper. REFERENCES [1] C. Cercignani, R. Illner and M. Pulvirenti, The Mathematical Theory of Dilute Gases, Applied Mathematical Sciences. Springer-Verlag, New York, 1994.
26 6 R.-J. DUAN [] L. Desvillettes and C. Villani, On the trend to global equilibrium for spatially inhomogeneous kinetic systems: The Boltzmann equation, Invent. Math., () (005), pp [3] R.-J. Duan, Hypocoercivity of linear degenerately dissipative kinetic equations, Nonlinearity, 4 (011), pp [4] R.-J. Duan, On the Cauchy problem for the Boltzmann equation in the whole space: Global existence and uniform stability in L (HN x ), Journal of Differential Equations, 44 (008), pp [5] R.-J. Duan and R. M. Strain, Optimal large-time behavior of the Vlasov-Maxwell-Boltzmann system in the whole space, Communications on Pure and Applied Mathematics, 64 (011), pp [6] R.-J. Duan and R. M. Strain, Optimal time decay of the Vlasov-Poisson-Boltzmann system in, Archive for Rational Mechanics and Analysis, 199 (011), pp [7] R.-J. Duan and T. Yang, Stability of the one-species Vlasov-Poisson-Boltzmann system, SIAM Journal on Mathematical Analysis, (6) (010), pp [8] J. Jang, Vlasov-Maxwell-Boltzmann diffusive limit, Arch. Rational Mech. Anal., 194 (009), pp [9] Y. Guo, The Boltzmann equation in the whole space, Indiana Univ. Math. J., 53 (004), pp [10] Y. Guo, The Vlasov-Poisson-Boltzmann system near Maxwellians, Comm. Pure Appl. Math., (9) (00), pp [11] Y. Guo, The Vlasov-Maxwell-Boltzmann system near Maxwellians, Invent. Math., 153 (003), pp [1] S. Kawashima, Systems of a hyperbolic-parabolic composite type, with applications to the equations of magnetohydrodynamics, thesis, Kyoto University (1983). [13] S. Kawashima, The Boltzmann equation and thirteen moments, Japan J. Appl. Math., 7 (1990), pp [14] T. P. Liu, T. Yang and S. H. Yu, Energy method for the Boltzmann equation, Physica D, 188 (3-4) (004), pp [15] T. P. Liu and S. H. Yu, Boltzmann equation: micro-macro decompositions and positivity of shock profiles, Commun. Math. Phys., 46 (1) (004), pp [16] P. A. Markowich, C. A. Ringhofer and C. Schmeiser, Semiconductor Equations, Springer- Verlag, Vienna, x+48 pp. [17] R. M. Strain, The Vlasov-Maxwell-Boltzmann system in the whole space, Commun. Math. Phys., 68 () (006), pp [18] R. M. Strain and Y. Guo, Exponential decay for soft potentials near Maxwellian, Arch. Ration. Mech. Anal., 187 (008), pp [19] S. Ukai, On the existence of global solutions of mixed problem for non-linear Boltzmann equation, Proceedings of the Japan Academy, 50 (1974), pp [0] C. Villani, Hypocoercivity, Memoirs Amer. Math. Soc., 0 (009), iv+141. [1] T. Yang and H. J. Zhao, Global existence of classical solutions to the Vlasov-Poisson- Boltzmann system, Commun. Math. Phys., 68 (006), pp
THE VLASOV-POISSON-BOLTZMANN SYSTEM FOR SOFT POTENTIALS
May 9, 01 11:15 WSPC/INSTRUCTION FILE M3AS-DYZ-s Mathematical Models and Methods in Applied Sciences c World Scientific Publishing Company THE VLASOV-POISSON-BOLTZMANN SYSTEM FOR SOFT POTENTIALS RENJUN
More informationALMOST EXPONENTIAL DECAY NEAR MAXWELLIAN
ALMOST EXPONENTIAL DECAY NEAR MAXWELLIAN ROBERT M STRAIN AND YAN GUO Abstract By direct interpolation of a family of smooth energy estimates for solutions near Maxwellian equilibrium and in a periodic
More informationThe Boltzmann Equation Near Equilibrium States in R n
The Boltzmann Equation Near Equilibrium States in R n Renjun Duan Department of Mathematics, City University of Hong Kong Kowloon, Hong Kong, P.R. China Abstract In this paper, we review some recent results
More informationhal , version 1-22 Nov 2009
Author manuscript, published in "Kinet. Relat. Models 1, 3 8) 355-368" PROPAGATION OF GEVREY REGULARITY FOR SOLUTIONS OF LANDAU EQUATIONS HUA CHEN, WEI-XI LI AND CHAO-JIANG XU Abstract. By using the energy-type
More informationOPTIMAL CONVERGENCE RATES FOR THE COMPRESSIBLE NAVIER-STOKES EQUATIONS WITH POTENTIAL FORCES
OPTIMAL CONVERGENCE RATES FOR THE COMPRESSIBLE NAVIER-STOKES EQUATIONS WITH POTENTIAL FORCES RENJUN DUAN Department of Mathematics, City University of Hong Kong 83 Tat Chee Avenue, Kowloon, Hong Kong,
More informationGLOBAL WELL-POSEDNESS IN SPATIALLY CRITICAL BESOV SPACE FOR THE BOLTZMANN EQUATION
GLOBAL WELL-POSEDNESS IN SPATIALLY CRITICAL BESOV SPACE FOR THE BOLTZMANN EQUATION RENJUN DUAN, SHUANGQIAN LIU, AND JIANG XU Abstract. The unique global strong solution in the Chemin-Lerner type space
More informationarxiv: v1 [math.ap] 28 Apr 2009
ACOUSTIC LIMIT OF THE BOLTZMANN EQUATION: CLASSICAL SOLUTIONS JUHI JANG AND NING JIANG arxiv:0904.4459v [math.ap] 28 Apr 2009 Abstract. We study the acoustic limit from the Boltzmann equation in the framework
More informationHypocoercivity for kinetic equations with linear relaxation terms
Hypocoercivity for kinetic equations with linear relaxation terms Jean Dolbeault dolbeaul@ceremade.dauphine.fr CEREMADE CNRS & Université Paris-Dauphine http://www.ceremade.dauphine.fr/ dolbeaul (A JOINT
More informationUne approche hypocoercive L 2 pour l équation de Vlasov-Fokker-Planck
Une approche hypocoercive L 2 pour l équation de Vlasov-Fokker-Planck Jean Dolbeault dolbeaul@ceremade.dauphine.fr CEREMADE CNRS & Université Paris-Dauphine http://www.ceremade.dauphine.fr/ dolbeaul (EN
More informationAuthor(s) Huang, Feimin; Matsumura, Akitaka; Citation Osaka Journal of Mathematics. 41(1)
Title On the stability of contact Navier-Stokes equations with discont free b Authors Huang, Feimin; Matsumura, Akitaka; Citation Osaka Journal of Mathematics. 4 Issue 4-3 Date Text Version publisher URL
More informationHypocoercivity and Sensitivity Analysis in Kinetic Equations and Uncertainty Quantification October 2 nd 5 th
Hypocoercivity and Sensitivity Analysis in Kinetic Equations and Uncertainty Quantification October 2 nd 5 th Department of Mathematics, University of Wisconsin Madison Venue: van Vleck Hall 911 Monday,
More informationYAN GUO, JUHI JANG, AND NING JIANG
LOCAL HILBERT EXPANSION FOR THE BOLTZMANN EQUATION YAN GUO, JUHI JANG, AND NING JIANG Abstract. We revisit the classical ork of Caflisch [C] for compressible Euler limit of the Boltzmann equation. By using
More informationA LOWER BOUND ON BLOWUP RATES FOR THE 3D INCOMPRESSIBLE EULER EQUATION AND A SINGLE EXPONENTIAL BEALE-KATO-MAJDA ESTIMATE. 1.
A LOWER BOUND ON BLOWUP RATES FOR THE 3D INCOMPRESSIBLE EULER EQUATION AND A SINGLE EXPONENTIAL BEALE-KATO-MAJDA ESTIMATE THOMAS CHEN AND NATAŠA PAVLOVIĆ Abstract. We prove a Beale-Kato-Majda criterion
More informationOptimal L p (1 p ) rates of decay to linear diffusion waves for nonlinear evolution equations with ellipticity and dissipation
Nonlinear Analysis ( ) www.elsevier.com/locate/na Optimal L p (1 p ) rates of decay to linear diffusion waves for nonlinear evolution equations with ellipticity and dissipation Renjun Duan a,saipanlin
More informationGlobal existence and asymptotic behavior of the solutions to the 3D bipolar non-isentropic Euler Poisson equation
Nonlinear Analysis: Modelling and Control, Vol., No. 3, 35 33 ISSN 139-5113 http://dx.doi.org/1.15388/na.15.3.1 Global existence and asymptotic behavior of the solutions to the 3D bipolar non-isentropic
More informationGlobal regularity of a modified Navier-Stokes equation
Global regularity of a modified Navier-Stokes equation Tobias Grafke, Rainer Grauer and Thomas C. Sideris Institut für Theoretische Physik I, Ruhr-Universität Bochum, Germany Department of Mathematics,
More informationGlobal Solutions for a Nonlinear Wave Equation with the p-laplacian Operator
Global Solutions for a Nonlinear Wave Equation with the p-laplacian Operator Hongjun Gao Institute of Applied Physics and Computational Mathematics 188 Beijing, China To Fu Ma Departamento de Matemática
More informationTHE WAVE EQUATION. d = 1: D Alembert s formula We begin with the initial value problem in 1 space dimension { u = utt u xx = 0, in (0, ) R, (2)
THE WAVE EQUATION () The free wave equation takes the form u := ( t x )u = 0, u : R t R d x R In the literature, the operator := t x is called the D Alembertian on R +d. Later we shall also consider the
More informationLARGE TIME BEHAVIOR OF THE RELATIVISTIC VLASOV MAXWELL SYSTEM IN LOW SPACE DIMENSION
Differential Integral Equations Volume 3, Numbers 1- (1), 61 77 LARGE TIME BEHAVIOR OF THE RELATIVISTIC VLASOV MAXWELL SYSTEM IN LOW SPACE DIMENSION Robert Glassey Department of Mathematics, Indiana University
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 informationDecay rate of the compressible Navier-Stokes equations with density-dependent viscosity coefficients
South Asian Journal of Mathematics 2012, Vol. 2 2): 148 153 www.sajm-online.com ISSN 2251-1512 RESEARCH ARTICLE Decay rate of the compressible Navier-Stokes equations with density-dependent viscosity coefficients
More informationEnergy transfer model and large periodic boundary value problem for the quintic NLS
Energy transfer model and large periodic boundary value problem for the quintic NS Hideo Takaoka Department of Mathematics, Kobe University 1 ntroduction This note is based on a talk given at the conference
More informationStability of Mach Configuration
Stability of Mach Configuration Suxing CHEN Fudan University sxchen@public8.sta.net.cn We prove the stability of Mach configuration, which occurs in moving shock reflection by obstacle or shock interaction
More informationINSTITUTE OF MATHEMATICS THE CZECH ACADEMY OF SCIENCES. Note on the fast decay property of steady Navier-Stokes flows in the whole space
INSTITUTE OF MATHEMATICS THE CZECH ACADEMY OF SCIENCES Note on the fast decay property of stea Navier-Stokes flows in the whole space Tomoyuki Nakatsuka Preprint No. 15-017 PRAHA 017 Note on the fast
More informationLOCAL WELL-POSEDNESS FOR AN ERICKSEN-LESLIE S PARABOLIC-HYPERBOLIC COMPRESSIBLE NON-ISOTHERMAL MODEL FOR LIQUID CRYSTALS
Electronic Journal of Differential Equations, Vol. 017 (017), No. 3, pp. 1 8. ISSN: 107-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu LOCAL WELL-POSEDNESS FOR AN ERICKSEN-LESLIE S
More informationin Bounded Domains Ariane Trescases CMLA, ENS Cachan
CMLA, ENS Cachan Joint work with Yan GUO, Chanwoo KIM and Daniela TONON International Conference on Nonlinear Analysis: Boundary Phenomena for Evolutionnary PDE Academia Sinica December 21, 214 Outline
More informationGlobal Weak Solution to the Boltzmann-Enskog equation
Global Weak Solution to the Boltzmann-Enskog equation Seung-Yeal Ha 1 and Se Eun Noh 2 1) Department of Mathematical Science, Seoul National University, Seoul 151-742, KOREA 2) Department of Mathematical
More informationStability of an abstract wave equation with delay and a Kelvin Voigt damping
Stability of an abstract wave equation with delay and a Kelvin Voigt damping University of Monastir/UPSAY/LMV-UVSQ Joint work with Serge Nicaise and Cristina Pignotti Outline 1 Problem The idea Stability
More informationBLOW-UP AND EXTINCTION OF SOLUTIONS TO A FAST DIFFUSION EQUATION WITH HOMOGENEOUS NEUMANN BOUNDARY CONDITIONS
Electronic Journal of Differential Equations, Vol. 016 (016), No. 36, pp. 1 10. ISSN: 107-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu BLOW-UP AND EXTINCTION OF SOLUTIONS TO A FAST
More informationWELL-POSEDNESS FOR HYPERBOLIC PROBLEMS (0.2)
WELL-POSEDNESS FOR HYPERBOLIC PROBLEMS We will use the familiar Hilbert spaces H = L 2 (Ω) and V = H 1 (Ω). We consider the Cauchy problem (.1) c u = ( 2 t c )u = f L 2 ((, T ) Ω) on [, T ] Ω u() = u H
More informationVector Spaces. Vector space, ν, over the field of complex numbers, C, is a set of elements a, b,..., satisfying the following axioms.
Vector Spaces Vector space, ν, over the field of complex numbers, C, is a set of elements a, b,..., satisfying the following axioms. For each two vectors a, b ν there exists a summation procedure: a +
More informationThe 2D Magnetohydrodynamic Equations with Partial Dissipation. Oklahoma State University
The 2D Magnetohydrodynamic Equations with Partial Dissipation Jiahong Wu Oklahoma State University IPAM Workshop Mathematical Analysis of Turbulence IPAM, UCLA, September 29-October 3, 2014 1 / 112 Outline
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 informationDecay in Time of Incompressible Flows
J. math. fluid mech. 5 (23) 231 244 1422-6928/3/3231-14 c 23 Birkhäuser Verlag, Basel DOI 1.17/s21-3-79-1 Journal of Mathematical Fluid Mechanics Decay in Time of Incompressible Flows Heinz-Otto Kreiss,
More informationSHARP BOUNDARY TRACE INEQUALITIES. 1. Introduction
SHARP BOUNDARY TRACE INEQUALITIES GILES AUCHMUTY Abstract. This paper describes sharp inequalities for the trace of Sobolev functions on the boundary of a bounded region R N. The inequalities bound (semi-)norms
More informationarxiv:math/ v1 [math.ap] 28 Oct 2005
arxiv:math/050643v [math.ap] 28 Oct 2005 A remark on asymptotic completeness for the critical nonlinear Klein-Gordon equation Hans Lindblad and Avy Soffer University of California at San Diego and Rutgers
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 informationGlobal existence of smooth solutions to two-dimensional compressible isentropic Euler equations for Chaplygin gases
Global existence of smooth solutions to two-dimensional compressible isentropic Euler equations for Chaplygin gases De-Xing Kong a Yu-Zhu Wang b a Center of Mathematical Sciences, Zhejiang University Hangzhou
More informationNULL CONDITION FOR SEMILINEAR WAVE EQUATION WITH VARIABLE COEFFICIENTS. Fabio Catalano
Serdica Math J 25 (999), 32-34 NULL CONDITION FOR SEMILINEAR WAVE EQUATION WITH VARIABLE COEFFICIENTS Fabio Catalano Communicated by V Petkov Abstract In this work we analyse the nonlinear Cauchy problem
More informationPresenter: Noriyoshi Fukaya
Y. Martel, F. Merle, and T.-P. Tsai, Stability and Asymptotic Stability in the Energy Space of the Sum of N Solitons for Subcritical gkdv Equations, Comm. Math. Phys. 31 (00), 347-373. Presenter: Noriyoshi
More informationBLOW-UP OF SOLUTIONS FOR VISCOELASTIC EQUATIONS OF KIRCHHOFF TYPE WITH ARBITRARY POSITIVE INITIAL ENERGY
Electronic Journal of Differential Equations, Vol. 6 6, No. 33, pp. 8. ISSN: 7-669. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu BLOW-UP OF SOLUTIONS FOR VISCOELASTIC EQUATIONS OF KIRCHHOFF
More informationSemigroup factorization and relaxation rates of kinetic equations
Semigroup factorization and relaxation rates of kinetic equations Clément Mouhot, University of Cambridge Analysis and Partial Differential Equations seminar University of Sussex 24th of february 2014
More informationHydrodynamic Limits for the Boltzmann Equation
Hydrodynamic Limits for the Boltzmann Equation François Golse Université Paris 7 & Laboratoire J.-L. Lions golse@math.jussieu.fr Academia Sinica, Taipei, December 2004 LECTURE 2 FORMAL INCOMPRESSIBLE HYDRODYNAMIC
More informationLECTURE 1: SOURCES OF ERRORS MATHEMATICAL TOOLS A PRIORI ERROR ESTIMATES. Sergey Korotov,
LECTURE 1: SOURCES OF ERRORS MATHEMATICAL TOOLS A PRIORI ERROR ESTIMATES Sergey Korotov, Institute of Mathematics Helsinki University of Technology, Finland Academy of Finland 1 Main Problem in Mathematical
More informationNONLINEAR FREDHOLM ALTERNATIVE FOR THE p-laplacian UNDER NONHOMOGENEOUS NEUMANN BOUNDARY CONDITION
Electronic Journal of Differential Equations, Vol. 2016 (2016), No. 210, pp. 1 7. ISSN: 1072-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu NONLINEAR FREDHOLM ALTERNATIVE FOR THE p-laplacian
More informationNonlinear Analysis 71 (2009) Contents lists available at ScienceDirect. Nonlinear Analysis. journal homepage:
Nonlinear Analysis 71 2009 2744 2752 Contents lists available at ScienceDirect Nonlinear Analysis journal homepage: www.elsevier.com/locate/na A nonlinear inequality and applications N.S. Hoang A.G. Ramm
More informationThe existence of Burnett coefficients in the periodic Lorentz gas
The existence of Burnett coefficients in the periodic Lorentz gas N. I. Chernov and C. P. Dettmann September 14, 2006 Abstract The linear super-burnett coefficient gives corrections to the diffusion equation
More informationExponential tail behavior for solutions to the homogeneous Boltzmann equation
Exponential tail behavior for solutions to the homogeneous Boltzmann equation Maja Tasković The University of Texas at Austin Young Researchers Workshop: Kinetic theory with applications in physical sciences
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 informationALEKSANDROV-TYPE ESTIMATES FOR A PARABOLIC MONGE-AMPÈRE EQUATION
Electronic Journal of Differential Equations, Vol. 2005(2005), No. 11, pp. 1 8. ISSN: 1072-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu ftp ejde.math.txstate.edu (login: ftp) ALEKSANDROV-TYPE
More informationOn the Three-Phase-Lag Heat Equation with Spatial Dependent Lags
Nonlinear Analysis and Differential Equations, Vol. 5, 07, no., 53-66 HIKARI Ltd, www.m-hikari.com https://doi.org/0.988/nade.07.694 On the Three-Phase-Lag Heat Equation with Spatial Dependent Lags Yang
More informationEXISTENCE OF GLOBAL WEAK SOLUTIONS TO A HYBRID VLASOV-MHD MODEL FOR MAGNETIZED PLASMAS
EXISTENCE OF GLOBAL WEAK SOLUTIONS TO A HYBRID VLASOV-MHD MODEL FOR MAGNETIZED PLASMAS BIN CHENG, ENDRE SÜLI, AND CESARE TRONCI Abstract. We prove the global-in-time existence of large-data finite-energy
More informationSeparation of Variables in Linear PDE: One-Dimensional Problems
Separation of Variables in Linear PDE: One-Dimensional Problems Now we apply the theory of Hilbert spaces to linear differential equations with partial derivatives (PDE). We start with a particular example,
More informationON WEAKLY NONLINEAR BACKWARD PARABOLIC PROBLEM
ON WEAKLY NONLINEAR BACKWARD PARABOLIC PROBLEM OLEG ZUBELEVICH DEPARTMENT OF MATHEMATICS THE BUDGET AND TREASURY ACADEMY OF THE MINISTRY OF FINANCE OF THE RUSSIAN FEDERATION 7, ZLATOUSTINSKY MALIY PER.,
More informationMildly degenerate Kirchhoff equations with weak dissipation: global existence and time decay
arxiv:93.273v [math.ap] 6 Mar 29 Mildly degenerate Kirchhoff equations with weak dissipation: global existence and time decay Marina Ghisi Università degli Studi di Pisa Dipartimento di Matematica Leonida
More informationGLOBAL EXISTENCE AND ENERGY DECAY OF SOLUTIONS TO A PETROVSKY EQUATION WITH GENERAL NONLINEAR DISSIPATION AND SOURCE TERM
Georgian Mathematical Journal Volume 3 (26), Number 3, 397 4 GLOBAL EXITENCE AND ENERGY DECAY OF OLUTION TO A PETROVKY EQUATION WITH GENERAL NONLINEAR DIIPATION AND OURCE TERM NOUR-EDDINE AMROUN AND ABBE
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 informationON THE DYNAMICAL SYSTEMS METHOD FOR SOLVING NONLINEAR EQUATIONS WITH MONOTONE OPERATORS
PROCEEDINGS OF THE AMERICAN MATHEMATICAL SOCIETY Volume, Number, Pages S -9939(XX- ON THE DYNAMICAL SYSTEMS METHOD FOR SOLVING NONLINEAR EQUATIONS WITH MONOTONE OPERATORS N. S. HOANG AND A. G. RAMM (Communicated
More information2 Formal derivation of the Shockley-Read-Hall model
We consider a semiconductor crystal represented by the bounded domain R 3 (all our results are easily extended to the one and two- dimensional situations) with a constant (in space) number density of traps
More informationStabilization for the Wave Equation with Variable Coefficients and Balakrishnan-Taylor Damping. Tae Gab Ha
TAIWANEE JOURNAL OF MATHEMATIC Vol. xx, No. x, pp., xx 0xx DOI: 0.650/tjm/788 This paper is available online at http://journal.tms.org.tw tabilization for the Wave Equation with Variable Coefficients and
More informationAW -Convergence and Well-Posedness of Non Convex Functions
Journal of Convex Analysis Volume 10 (2003), No. 2, 351 364 AW -Convergence Well-Posedness of Non Convex Functions Silvia Villa DIMA, Università di Genova, Via Dodecaneso 35, 16146 Genova, Italy villa@dima.unige.it
More informationExistence of global solutions to the Cauchy problem for the inelastic Boltzmann equation with near-vacuum data
Existence of global solutions to the Cauchy problem for the inelastic Boltzmann equation with near-vacuum data Ricardo J. Alonso July 12, 28 Abstract The Cauchy problem for the inelastic Boltzmann equation
More informationA COUNTEREXAMPLE TO AN ENDPOINT BILINEAR STRICHARTZ INEQUALITY TERENCE TAO. t L x (R R2 ) f L 2 x (R2 )
Electronic Journal of Differential Equations, Vol. 2006(2006), No. 5, pp. 6. ISSN: 072-669. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu ftp ejde.math.txstate.edu (login: ftp) A COUNTEREXAMPLE
More informationA REMARK ON AN EQUATION OF WAVE MAPS TYPE WITH VARIABLE COEFFICIENTS
A REMARK ON AN EQUATION OF WAVE MAPS TYPE WITH VARIABLE COEFFICIENTS DAN-ANDREI GEBA Abstract. We obtain a sharp local well-posedness result for an equation of wave maps type with variable coefficients.
More informationDecay rates for partially dissipative hyperbolic systems
Outline Decay rates for partially dissipative hyperbolic systems Basque Center for Applied Mathematics Bilbao, Basque Country, Spain zuazua@bcamath.org http://www.bcamath.org/zuazua/ Numerical Methods
More informationRegularity of the density for the stochastic heat equation
Regularity of the density for the stochastic heat equation Carl Mueller 1 Department of Mathematics University of Rochester Rochester, NY 15627 USA email: cmlr@math.rochester.edu David Nualart 2 Department
More informationScientiae Mathematicae Japonicae Online, Vol. 5, (2001), Ryo Ikehata Λ and Tokio Matsuyama y
Scientiae Mathematicae Japonicae Online, Vol. 5, (2), 7 26 7 L 2 -BEHAVIOUR OF SOLUTIONS TO THE LINEAR HEAT AND WAVE EQUATIONS IN EXTERIOR DOMAINS Ryo Ikehata Λ and Tokio Matsuyama y Received November
More informationClassical inequalities for the Boltzmann collision operator.
Classical inequalities for the Boltzmann collision operator. Ricardo Alonso The University of Texas at Austin IPAM, April 2009 In collaboration with E. Carneiro and I. Gamba Outline Boltzmann equation:
More informationMP463 QUANTUM MECHANICS
MP463 QUANTUM MECHANICS Introduction Quantum theory of angular momentum Quantum theory of a particle in a central potential - Hydrogen atom - Three-dimensional isotropic harmonic oscillator (a model of
More informationLecture No 1 Introduction to Diffusion equations The heat equat
Lecture No 1 Introduction to Diffusion equations The heat equation Columbia University IAS summer program June, 2009 Outline of the lectures We will discuss some basic models of diffusion equations and
More informationChapter 8 Integral Operators
Chapter 8 Integral Operators In our development of metrics, norms, inner products, and operator theory in Chapters 1 7 we only tangentially considered topics that involved the use of Lebesgue measure,
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 informationThe propagation of chaos for a rarefied gas of hard spheres
The propagation of chaos for a rarefied gas of hard spheres Ryan Denlinger 1 1 University of Texas at Austin 35th Annual Western States Mathematical Physics Meeting Caltech February 13, 2017 Ryan Denlinger
More informationSYMMETRY RESULTS FOR PERTURBED PROBLEMS AND RELATED QUESTIONS. Massimo Grosi Filomena Pacella S. L. Yadava. 1. Introduction
Topological Methods in Nonlinear Analysis Journal of the Juliusz Schauder Center Volume 21, 2003, 211 226 SYMMETRY RESULTS FOR PERTURBED PROBLEMS AND RELATED QUESTIONS Massimo Grosi Filomena Pacella S.
More informationON THE STRONG SOLUTIONS OF THE INHOMOGENEOUS INCOMPRESSIBLE NAVIER-STOKES EQUATIONS IN A THIN DOMAIN
ON THE STRONG SOLUTIONS OF THE INHOMOGENEOUS INCOMPRESSIBLE NAVIER-STOKES EQUATIONS IN A THIN DOMAIN XIAN LIAO Abstract. In this work we will show the global existence of the strong solutions of the inhomogeneous
More informationHilbert Spaces. Hilbert space is a vector space with some extra structure. We start with formal (axiomatic) definition of a vector space.
Hilbert Spaces Hilbert space is a vector space with some extra structure. We start with formal (axiomatic) definition of a vector space. Vector Space. Vector space, ν, over the field of complex numbers,
More informationSTABILITY OF THE RELATIVISTIC MAXWELLIAN IN A COLLISIONAL PLASMA
STABILITY OF THE RELATIVISTIC MAXWELLIAN IN A COLLISIONAL PLASMA ROBERT M. STRAIN AND YAN GUO Abstract. The relativistic Landau-Maxwell system is among the most fundamental and complete models for describing
More informationTong Sun Department of Mathematics and Statistics Bowling Green State University, Bowling Green, OH
Consistency & Numerical Smoothing Error Estimation An Alternative of the Lax-Richtmyer Theorem Tong Sun Department of Mathematics and Statistics Bowling Green State University, Bowling Green, OH 43403
More informationSome Applications of an Energy Method in Collisional Kinetic Theory
Some Applications of an Energy Method in Collisional Kinetic Theory by Robert Mills Strain III B. A., New York University, Sc. M., Brown University, A Dissertation submitted in partial fulfillment of the
More informationAnalysis in weighted spaces : preliminary version
Analysis in weighted spaces : preliminary version Frank Pacard To cite this version: Frank Pacard. Analysis in weighted spaces : preliminary version. 3rd cycle. Téhéran (Iran, 2006, pp.75.
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 informationDissipative quasi-geostrophic equations with L p data
Electronic Journal of Differential Equations, Vol. (), No. 56, pp. 3. ISSN: 7-669. URL: http://ejde.math.swt.edu or http://ejde.math.unt.edu ftp ejde.math.swt.edu (login: ftp) Dissipative quasi-geostrophic
More informationSmoluchowski Navier-Stokes Systems
Smoluchowski Navier-Stokes Systems Peter Constantin Mathematics, U. of Chicago CSCAMM, April 18, 2007 Outline: 1. Navier-Stokes 2. Onsager and Smoluchowski 3. Coupled System Fluid: Navier Stokes Equation
More informationMathematical modelling of collective behavior
Mathematical modelling of collective behavior Young-Pil Choi Fakultät für Mathematik Technische Universität München This talk is based on joint works with José A. Carrillo, Maxime Hauray, and Samir Salem
More information16. Local theory of regular singular points and applications
16. Local theory of regular singular points and applications 265 16. Local theory of regular singular points and applications In this section we consider linear systems defined by the germs of meromorphic
More informationarxiv: v4 [math.ap] 14 Apr 2012
OPTIMAL TIME DECAY OF THE NON CUT-OFF BOLTZMANN EQUATION IN THE WHOLE SPACE arxiv:10115561v4 [mathap] 14 Apr 01 ROBERT M STRAIN Abstract In this paper we study the large-time behavior of perturbative classical
More informationBIHARMONIC WAVE MAPS INTO SPHERES
BIHARMONIC WAVE MAPS INTO SPHERES SEBASTIAN HERR, TOBIAS LAMM, AND ROLAND SCHNAUBELT Abstract. A global weak solution of the biharmonic wave map equation in the energy space for spherical targets is constructed.
More informationEntropy Methods for Reaction-Diffusion Equations with Degenerate Diffusion Arising in Reversible Chemistry
1 Entropy Methods for Reaction-Diffusion Equations with Degenerate Diffusion Arising in Reversible Chemistry L. Desvillettes CMLA, ENS Cachan, IUF & CNRS, PRES UniverSud, 61, avenue du Président Wilson,
More informationON LIQUID CRYSTAL FLOWS WITH FREE-SLIP BOUNDARY CONDITIONS. Chun Liu and Jie Shen
DISCRETE AND CONTINUOUS Website: http://aimsciences.org DYNAMICAL SYSTEMS Volume 7, Number2, April2001 pp. 307 318 ON LIQUID CRYSTAL FLOWS WITH FREE-SLIP BOUNDARY CONDITIONS Chun Liu and Jie Shen Department
More informationABOUT SOME TYPES OF BOUNDARY VALUE PROBLEMS WITH INTERFACES
13 Kragujevac J. Math. 3 27) 13 26. ABOUT SOME TYPES OF BOUNDARY VALUE PROBLEMS WITH INTERFACES Boško S. Jovanović University of Belgrade, Faculty of Mathematics, Studentski trg 16, 11 Belgrade, Serbia
More informationWaves in Flows. Global Existence of Solutions with non-decaying initial data 2d(3d)-Navier-Stokes ibvp in half-plane(space)
Czech Academy of Sciences Czech Technical University in Prague University of Pittsburgh Nečas Center for Mathematical Modelling Waves in Flows Global Existence of Solutions with non-decaying initial data
More informationS. Lototsky and B.L. Rozovskii Center for Applied Mathematical Sciences University of Southern California, Los Angeles, CA
RECURSIVE MULTIPLE WIENER INTEGRAL EXPANSION FOR NONLINEAR FILTERING OF DIFFUSION PROCESSES Published in: J. A. Goldstein, N. E. Gretsky, and J. J. Uhl (editors), Stochastic Processes and Functional Analysis,
More informationSome asymptotic properties of solutions for Burgers equation in L p (R)
ARMA manuscript No. (will be inserted by the editor) Some asymptotic properties of solutions for Burgers equation in L p (R) PAULO R. ZINGANO Abstract We discuss time asymptotic properties of solutions
More informationSufficient conditions for functions to form Riesz bases in L 2 and applications to nonlinear boundary-value problems
Electronic Journal of Differential Equations, Vol. 200(200), No. 74, pp. 0. ISSN: 072-669. URL: http://ejde.math.swt.edu or http://ejde.math.unt.edu ftp ejde.math.swt.edu (login: ftp) Sufficient conditions
More informationResearch Article On a Quasi-Neutral Approximation to the Incompressible Euler Equations
Applied Mathematics Volume 2012, Article ID 957185, 8 pages doi:10.1155/2012/957185 Research Article On a Quasi-Neutral Approximation to the Incompressible Euler Equations Jianwei Yang and Zhitao Zhuang
More informationSpectral analysis for rank one perturbations of diagonal operators in non-archimedean Hilbert space
Comment.Math.Univ.Carolin. 50,32009 385 400 385 Spectral analysis for rank one perturbations of diagonal operators in non-archimedean Hilbert space Toka Diagana, George D. McNeal Abstract. The paper is
More informationNew Identities for Weak KAM Theory
New Identities for Weak KAM Theory Lawrence C. Evans Department of Mathematics University of California, Berkeley Abstract This paper records for the Hamiltonian H = p + W (x) some old and new identities
More informationPositive solutions of BVPs for some second-order four-point difference systems
Positive solutions of BVPs for some second-order four-point difference systems Yitao Yang Tianjin University of Technology Department of Applied Mathematics Hongqi Nanlu Extension, Tianjin China yitaoyangqf@63.com
More informationAbstract. 1. Introduction
Journal of Computational Mathematics Vol.28, No.2, 2010, 273 288. http://www.global-sci.org/jcm doi:10.4208/jcm.2009.10-m2870 UNIFORM SUPERCONVERGENCE OF GALERKIN METHODS FOR SINGULARLY PERTURBED PROBLEMS
More informationON THE CONSTRUCTION OF DUALLY FLAT FINSLER METRICS
Huang, L., Liu, H. and Mo, X. Osaka J. Math. 52 (2015), 377 391 ON THE CONSTRUCTION OF DUALLY FLAT FINSLER METRICS LIBING HUANG, HUAIFU LIU and XIAOHUAN MO (Received April 15, 2013, revised November 14,
More information