Small BGK waves and nonlinear Landau damping (higher dimensions)

Size: px
Start display at page:

Download "Small BGK waves and nonlinear Landau damping (higher dimensions)"

Transcription

1 Small BGK waves and nonlinear Landau damping higher dimensions Zhiwu Lin and Chongchun Zeng School of Mathematics Georgia Institute of Technology Atlanta, GA, USA Abstract Consider Vlasov-Poisson system with a fied ion background and periodic boundary conditions on the space variables, in dimension d =,. First, we show that for general homogeneous equilibrium and any periodic bo, within any small neighborhood in the Sobolev space W,v s,p p >, s < + of the steady distribution function, there eist nontrivial travelling wave solutions BGK waves with arbitrary trav- p eling speed. This implies that nonlinear Landau damping is not true in W s s,p < + space for any homogeneous equilibria and in any period p bo. The BGK waves constructed are one dimensional, that is, depending only on one space variable. Higher dimensional BGK waves are shown to not eist. Second, we prove that there eist no nontrivial invariant structures in some neighborhood of stable homogeneous equilibria, in the + v b -weighted H s,v b > d, s > 4 space. Since arbitrarilly small BGK waves can also be constructed near any homogeneous equilibria in such weighted H,v s s < norm, this shows that s = is the critical regularity for the eistence of nontrivial invariant structures near stable homogeneous equilibria. These generalize our previous results in the one dimensional case. Introduction Consider a collisionless electron plasma with a fied homogeneous neutralizing ion background. The Vlasov-Poisson system in d dimension is t f + v f E v f =, E = φ, φ = f dv +, d a b where f t,, v is the distribution function, E, t is the electrical field and φ, t is the electrical potential. We consider the Vlasov-Poisson system

2 in a periodic bo, with periods T i in i. In 946, Landau [6], looking for analytical solutions of the linearized Vlasov-Poisson system around Mawellian e v,, pointed out that the electric field is subject to time decay even in the absence of collisions. The effect of this Landau damping, as it is subsequently called, plays a fundamental role in the study of plasma physics. However, Landau s treatment is in the linear regime; that is, only for infinitesimally small initial perturbations. ecently, nonlinear Landau damping was shown [] for analytical perturbations of stable equilibria with linear eponential decay. For general perturbations in Sobolev spaces, the proof of nonlinear damping remains open. We refer to [9] [] for more discussions and references on this topic. In [9], the following results were obtained for D Vlasov-Poisson system. First, we show that for general homogeneous equilibria, within any small neighborhood in the Sobolev space W p s,p >, s < + p of the steady distribution function, there eist nontrivial travelling wave solutions BGK waves with arbitrary minimal period and traveling speed. This implies that nonlinear Landau damping is not true in W s s,p < + p space for any homogeneous equilibria and any spatial period. Second, it is shown that for homogeneous equilibria satisfying Penrose s linear stability condition, there eist no nontrivial travelling BGK waves and unstable homogeneous states in some W p s,p >, s > + p neighborhood. Furthermore, we prove that there eist no nontrivial invariant structures in the H s s > neighborhood of stable homogeneous states. In particular, these results suggest the contrasting long time dynamics in the H s s > and H s s < neighborhoods of a stable homogeneous state. In this paper, we generalize above results to higher dimensions d =,. Denote the fractional order Sobolev spaces by W s,p d or W s,p,v, T d with p >, s. These spaces are the comple interpolation of of L p space and Sobolev spaces W m,p m positive integer. Our first result is to construct D BGK waves in W s,p,v s < + p spaces. Theorem. Assume the homogeneous distribution function f v W s,p d d, p >, s [, + p satisfies f v, f v dv =, v f v dv < +. Fi T > and c. Then for any ε >, there eist travelling BGK wave solutions of the form f = f ε ct, v, E = Eε ct e to??, such that f ε, v, E ε has minimal period T in, f ε, v, E ε is not identically zero, and f ε f L,v + T d v f ε f d dv + f ε f W s,p,v < ε.

3 In Proposition., we show that there eist no D and D BGK waves. Therefore, the form of D BGK waves in Theorem. is in some sense necessary. For any b > d 4, we denote Hs,b v d to be the + v b weighted H s space, that is, Hv s,b d { = f + v } b s f <. L d and f H s,b v = + v b s f L d. Let T d be a periodic bo with periods T i in i i =,, d, and { π Z d = j,, π } j d j,, j d are integers. T T d We define the space H s Hv sv,b T d d by h = e i k h k v H s Hv sv,b if h H s Hv k Z d = h Hv + k Z d k s h k Hv <. Given a homogeneous profile f v H s,b d d, s >, b > d 4, we have the following generalization of Penrose s stability condition in higher dimensions k ma v i S k/ k f, k/ k α α v i dα >, for any k Z d, 4 where f, k/ k α is the projection of f v to k direction, as defined by and S k/ k is the set of all critical points of f, k/ k α. By Corollary., one only need to check the condition 4 for finitely many k Z d with k C d, s, b f H s,b d. In particular, for a single humped profile f v = µ v with µ e <, the stability condition 4 is satisfied for any period set T,, T d. The following Theorem ecludes any nontrivial invariant structures steady, time periodic, quasi-periodic etc. near stable homogeneous equilibria satisfying 4, in the H s Hv sv,b spaces of high v regularity.

4 Theorem. Consider a homogeneous profile f v H s,b d d, s >, b > d, 5 4 with f v and f v dv =. Let T d be a periodic bo with periods T i in i i =,, d. Assume that f v satisfies the Penrose stability condition 4 for T,, T d. Let f, v, t, E, v, t be the solution of?? in T d. For any s, s v satisfying there eists ε >, such that if s, s > d, and < s v s, 6 f t f H s then E t for all t. Hv < ε, for all t, 7 In the above Theorem, the assumption s > d, s is to make H s an algebra which is needed in the proof of Lemma.. The use of weighted Sobolev space Hv s,b in Theorem. is rather natural in higher dimensions. Indeed, even to state the Penrose s stability condition 4, we need to assume that the homogeneous equilibrium f v H s,b d with s, b satisfying 5. This is because that linear instability stability of homogeneous equilibria of Vlasov- Poisson is longitudinal along the wave direction of perturbation. The weighted Sobolev space 5 is needed to ensure that the projected steady distribution function in any wave direction is in H s s > which is necessary to get the D Penrose stability criterion. Moreover, in Theorem. we also construct D BGK waves arbitrarily near any homogeneous equilibrium in H s Hv sv,b d, b > d 4 for any s > and s v <. Combined with Theorem., this shows that for weighted Sobolev spaces H s Hv sv,b, the critical v regularity for the eistence of nontrivial invariant structures near a stable homogeneous equilibrium is s v =. This generalizes our D results in [9] to higher dimensions. We note that the critical regularity s v = does not depend on the dimension. This illustrates again the longitudinal D nature of Landau damping, which is obvious in the linear regime. We briefly mention some differences of the long time behaviors of Vlasov- Poisson in D and higher dimensions. For the D case, numerical simulations e.g. [4] [] indicated that for certain small initial data near a stable homogeneous state including Mawellian, there is no decay of electric fields and the asymptotic state is a BGK wave or the superposition of BGK waves. Moreover, BGK waves also appear as the asymptotic states for the saturation of an unstable homogeneous state [] in D. These suggest that small BGK waves play an important role in understanding the long time behaviors of D Vlasov-Poisson system. However, for D and D Vlasov-Poisson system, numerical simulations [] [] suggested that when starting near a homogeneous state both stable 4

5 and unstable, the electric fields decay eventually. Our Theorems. and. on eistence of D BGK waves show that such decay of electric field is not true for general initial data near homogeneous states. But the numerical simulations seem to suggest that these D BGK waves do not appear in the long time dynamics in D and D cases. To eplain these phenomena, it will be important to understand the transversal instability of D BGK waves. This paper is organized as follows. In Section, we prove the eistence of D BGK waves in W s s,p < + p neighborhoods of homogeneous states. In Section, we use the linear decay estimate to show that all invariant structures near stable homogeneous equilibria in H s Hv sv,b spaces satisfying 6 are trivial. Throughout this paper, we use C to denote a generic constant in the estimates and the dependence of C is indicated only when it matters in the proof. Eistence of BGK waves in W s,p s < + p In this Section, we construct nontrivial steady states BGK waves near any homogeneous state in the space W,v s,p s < + p. We consider d = only, since the proof is almost the same for d =. The BGK waves we construct are one-dimensional, that is, the steady distribution f = f, v, v and the electric field E = E e. We will show that such a restriction is necessary by ecluding D and D BGK waves. Proof of Theorem.. We adapt the line of proof in [9] to construct BGK wave solutions for D Vlasov-Poisson equations. First, we modify f v to a smooth function f v with some additional properties. Let η v v be the standard mollifier function. For > define f v = η v f v, where η v = η and when is small enough v. Then by the properties of mollifiers, we have f C, f v, f v dv =, f f L + v f f dv + f f W s,p ε 6. Modifying f v near infinity by cut-off, we can assume in addition that f v H,b defined in. In the second step, let σ = σ be the D cut-off function. Let > be a small number, and define v f v, v + f v, v v f, v, v = f v, v σ + σ f v, v f v, v v = f v, v σ. 5

6 Then, f, C, f, v >, f, v dv = f v dv =, and f, v, v is even in v when v [, ]. We show that: when is small enough f, f L + v f, f dv+ f, f W s,p ε 6. 8 Fist, a minor modification of the proof of Lemma. in [9] yields that: when, f f, L + v f f, dv + f f, W,p. It remains to show that f f, W s,p, when. 9 We have v f v, v v f v, v v f f, = σ v, and v f v, v + v f v, v v f f, = σ + σ v v f v, v f v, v. v v By a scaling argument as in the proof of Lemma. of [9], f v, v f v, v v C f W s,p. W s,p So 9 follows from Lemma. below. Thus for any fied ε >, by choosing, small enough, we get f, f L + v f, f dv + f, f W s,p ε. We set f v, v = f, v, v, then f v >, f v C H,b, f v dv =, f v is even for v in [, ] and within ε distance of f v in the norm of. Below, we denote a = /. 6

7 Fi the period T >, we only consider the travel speed c = since the construction for any c follows by the Galilean transform as in [9]. Our strategy is to construct BGK wave solutions of the form f ε, v, v, E ε e by bifurcation at a modified homogeneous profile near f v, v. Denote σ = σ to be the cut-off function such that σ C, σ ; σ = when ; σ = when. Similar to Lemma. in [9], there eists g, C, g = when 4a, such that v f v, v σ = g v a, v. Define g +,, g, C by f ±, σ a g ±, = Then + g, if > a g, if 4a < a. if 4a { g+ v f v, v =, v g v, v if v > if v. Since v f, v =, f C H,b, we have v f v dv v <. Indeed, let f v = f v, v dv, then since f =, by Corollary., v f v dv v = f v dv v C f H,b. We consider three cases. Case : F v = ep v f v v dv < v v π T. Let + ep v + v = G v, and F v = e v, where v is a large positive constant such that F v dv >. v Let γ, > be two small parameters to be fied, define [ f γ, v, v = f v, v + γ F + C γ 7 v γ ] F v,

8 where C = F v F v dv >. The rest of the proof is similar to the proof of Proposition. in [9]. We sketch it below. There eists < < such that for γ > small enough < v f γ, v, v v dv < π < T v f γ, v, v v dv, when < γ < γ. Let β be a T periodic function and denote e = v β. We look for D BGK wave solution f = f β γ,, v, E = E e near f γ,,, where f β γ, +C, v = γ [g + e, v + γ G +C γ [g e, v + γ G e γ e γ ] F v ] F v if v > if v and E = β. The steady Vlasov-Poisson equation is reduced to the ODE β = f β γ,, v dv 4 = + C γ { g + e, v dv + g e, v dv v > v γ + G e γ F v dv} := h γ, β. Since and h γ, = f γ, v dv = d = h γ, + C γ { v > g + v, v dv + + γ γ G v v γ g v, v dv F v dv } v f γ, v, v = dv <, when < γ < γ, < <, v so β = is a center for the ODE 4 and there eist bifurcation of periodic solutions. More precisely, for any fied γ, γ, there eists r > independent 8

9 of,, such that for each < r < r, there eists a T γ, ; r periodic solution β γ,;r to the ODE 4 with β γ,;r H,T γ,;r = r. Moreover, π v f γ, v, v dv, when r. T γ, ; r v To get a solution with the given period T, we adjust [, ] by using the inequality and the fact that T γ, ; r is continuous to. So for each γ, r > small enough, there eists T γ, r,, such that T γ, T ; r = T. Define f γ,r, v = f β γ, T, v, β γ,r = β γ,t ;r and let E γ,r = β γ,r e. Then f γ,r, v, E γ,r is a nontrivial steady solution to?? with period T. For any fied γ >, let γ = lim r T γ, r [, ]. By the dominant convergence theorem, it is easy to show that fγ,r, v f γ,γ v T L + v fγ,r, v f γ,γ v d dv,v + fγ,r, v f γ,γ v W,p,,v when r = β γ,r H,T. So for any γ > and ε >, there eists r = r γ, ε > such that fγ,r, v f γ,γ v L,v + Since f v f γ,γ v = T v fγ,r, v f γ,γ v d dv + fγ,r, v f γ,γ v W,p < ε,v. [ + C γ C γ f v γ F and γ [, ], by using Lemma., for s < + p, f v f γ,γ v W, when γ. s,p ] v F v. γ It is also easy to show that f v f γ,γ v L + v f v f γ,γ v dv, when γ. Thus we can choose γ > small enough such that T f v f γ,γ v L + T v f v f γ,γ v dv 9

10 + f v f γ,γ v < ε W s,p. So the nontrivial steady solution f γ,r, v, E γ,r is within ε distance of the homogeneous state f v, in the norm of. Case : v f v. π v dv < T Choose F v = ep v and F v is the same as before. Define f γ, v as in Case see. Then there eists < < such that v f γ, < v, v π v f γ, dv < < v, v dv. v T v The rest of the proof is the same as in Case. v f v π T. For >, define f v, v = f Case : v dv = For any ε >, there eist < ε < < ε such that < and when ε, ε, v f v v dv < T f v f v L + T π < T v f v dv, v v f v f v dv + f v f v W s,p < ε. We construct steady BGK waves near f v,, which are of the form { f β, v = g + g e e, v, v v, v. if v > if v, e = v β, 5 and E = β e. The eistence of BGK waves is then reduced to solve the ODE β = f β, v dv := h β. 6 As in Case, for any ε, ε, r ε > independent of such that for each < r < r, there eists a T ; r periodic solution β ;r to the ODE 6, satisfying β ;r H,T ;r = r and π v f v dv, when r. T ; r v For r small enough, again there eists T r, ε ε, ε such that T T ; r = T. Define f r,ε, v = f β T, v and E r,ε = β T ;r e.

11 Then f r,ε, v, E r,ε is a nontrivial steady solution to?? with period T. As in Cases and, by choosing r small enough, f r,ε, v is within ε distance of the homogeneous state f v, in the norm of. This finishes the proof of the Theorem.. Lemma. i Given f W p,p L,and Then for >, f g W s,p p >, s < p v ii Given f, g W s,p f v. g v, v, when. 7 W s,p p >, s < p. Then for >, g v, when. W s,p Proof. Proof of i: First we consider g C for W s,p norm see [5], we have v f g v, v W s,p v C f g v, v v + f L p v W s,p v. By Fubini Theorem g v, v W s,p v. L p v By the estimates in the proof of Theorem. of [5], for any p >, s < p, when h W s,p, h W p,p L, we have h h W s,p C h W s,p h W p,p + h L. So f v C f C f C f g v, v v v W s,p v L p v + g L v W s,p g p,p Wv W s,p g Wv,p L p v W g s,p W,p, v L p v

12 when.since f v W s,p under the assumption s < p see [9] for a proof. By the trace Theorem, we also have v v Lp f g v, v W s,p v C f g L p v W,p, when. This proves 7 for g C. When g W s,p, 7 can be proved by using C functions as approimations. Proof of ii: By Fubini Theorem for W s,p norm, v f g v W s,p f v W v L C g v s,p L p v v + g v W s,p v f p v, when. By the similar proof of Theorem., we can get the following. Theorem. Assume the homogeneous distribution function satisfies f v H sv,b d d, b > d 4, s v [, f v, f v dv =, v f v dv < +. Fi T > and c. Then for any ε >, s, there eist travelling wave solutions of the form f = f ε ct, v, E = E ε ct e to??, such that f ε, v, E ε has minimal period T in, f ε, v, E ε is not identically zero, and T f ε f + L v f ε, v f v d dv+ f ε f,v H s H < ε. d v 8 Proof. The construction of BGK waves follows the same line of the proof of Theorem.. First, we modify f v to a smooth profile f v. Then by adding proper perturbations in a scaling form to f v, we get the modified profile f γ, v. The BGK waves f ε, v, E ε are obtained by bifurcation near f γ, v,. To show the estimate 8, we need to control three deviations in the norm of 8: i f ε, v f γ, v ; ii f γ, v f v ; and iii f v f v. For the estimate of i, we choose integers s s, s v s v, and b b and it is easy to show that f ε, v f γ, v H s Hv C f ε, v f γ, v H s sv, b Hv and the right hand side can be made arbitrarily small by using the dominant convergence Theorem. For estimates of ii and iii, we use the following analogue of Lemma..

13 Lemma. i Given f v H L, g v, v H s,b s <, b > 4. For >, f v g v, v, when. H s,b ii Given f, g H s,b s <, b > 4. For >, v f g v, when. H s,b Proof. First, we show that for any function h H s,b d d, s, b > 4, the norm h H s,b d defined by is equivalent to both + v b H f 9 s d and + v b + + v d b f H s d. We only need to prove the equivalence of the norms and 9 for s = and s =, since then for < s < it follows from interpolation. For s =, it is trivial. For s =, by choosing a > small enough, we have + a v b f + a v b L f d = f + a v b + f + a v b L + a v b L f. Thus + a v b L f + a v b L f + a v b L f. The equivalence of and can be proved in the same way. Proof of i: By Lemma. i, v f g v, v C H + v b v f g v, v s,b, Hs when. Since f v H, L and + v b g v, v C g H s,b <. H s

14 Proof of ii: By using the equivalent norm and Lemma. ii, v f g v H s,b C + v b + v b v f g v Hs v C{ f g v + b v b v f g v H s Hs v + f v b g v } H s, when. In the following, we show that there eist no truly D or D BGK solutions. Therefore, the D BGK form of solutions constructed in Theorem. is in some sense necessary. Proposition. i d = Assume µ C L +, µ. If f, v = µ v β, E = β is a solution of the Vlasov-Poisson system, then E. ii d = Assume µ C L +, µ. If f,, v, v, v = µ v + v β,, v, E, = β, β, is a solution of the Vlasov-Poisson system, then E. iii d = Assume µ C, µ, µ r r L +. If f, v = µ v β, E = β is a solution of the Vlasov-Poisson system, then E. Proof. We only prove i since the proof of ii and iii is similar. The electric potential β satisfies β = v β dv + = g β. µ By the assumptions on µ, we have g β C and g β = v β dv = π µ µ s β ds = πµ β, 4

15 Taking derivative of and integrating with β, we have β T d = g β β d. T So T β d = and β is a constant C. By the periodic assumption of β, C = and thus β. Similarly, β. emark. In D and D, the function g β defined in always satisfies g β and thus the elliptic problem only has trivial solutions. For D, the function g β can change signs and thus the eistence of D BGK waves is possible. We note that Proposition. does not eclude steady travelling wave solutions not of BGK types. It would be interesting to construct or eclude nontrivial steady solutions not of BGK types. Non-eistence of Invariant structures In this section, we prove that there eist no nontrivial invariant structures near stable equilibria in spaces H s Hv sv,b sv >. First, we derive the optimal linear decay estimate in a space-time norm under the Penrose stability condition 4. By assuming the uniform in time bound close to the stable homogeneous state above the regularity threshold s v >, it is shown that the solutions of nonlinear equation is dominated by the linearized equation and a decay estimate in space-time is obtained for the nonlinear solution. By the time translation property of Vlasov-Poisson system, such a space-time decay implies that the invariant solution must be trivial, that is, homogeneous. Lemma. Given f v H s,b d d, s, b > d 4. For any unit vector e d, let v = α e + w where v d and w e. Define f e α = f α e + w dw. d Then for some constant C independent of e. f e α H s C f H s,b d Proof. To simplify notations, we only consider e =,,,. Then α = v and f e v = f v dv dv d. d 5

16 Let ξ = ξ, ξ d be the Fourier variable. Then f e v H s = + ξ s ˆf e ξ L = + ξ s ˆf ξ,,, L C + ξ s H ˆf ξ b d = C + v b s f = C f H s,b d. L d Here, the first inequality above is due to the trace theorem and that b > d. Corollary. Given f v H s,b d d, s >, b > d 4. For any unit vector e d, we have i If α is a critical point of f e α, then f e α α α dα C d, s, b f H s,b d. ii For any α, P f e α α α dα C d, s, b f H s,b d, where P is the principal value integral. Proof. i follows from Lemma. and the following Hardy inequality see Lemma. of [9]: If u v W s,p p >, s > p, and u =, then u v v dv C u W s,p, for some constant C. Proof of ii: Since P f e α α α dα = d dα f e α + α + f e α + α α by Hardy inequality P f e α α α dα f C e α + α Hs + f e α + α Hs C f e α H s C f H s,b d dα, by Lemma.. 6

17 The net lemma generalizes the one dimensional result in [9] about linear decay estimates in a space-time norm. The linearized Vlasov-Poisson system at an homogeneous state f, E = f v, is t f + v f E v f =, E = φ, φ = f dv, d 4a 4b Lemma. Assume f v H s,b d d, s >, b > d 4 satisfies the Penrose stability condition 4 for period tuple T,, T d. Let f, v, t, E, t be a solution of the linearized Vlasov-Poisson system 4a-4b with period tuple T,, T d. If f, v, H s Hv sv,b with s v s, then t sv E, t C f, v, H s, 5 L t H +s+sv Hv Proof. First, we reduce the linearized problem to the one dimensional case. Since the homogeneous component of f, v, t remains steady for the linearized equation and therefore has no effect on E, t, we assume that f has no homogeneous component. Let f, v, t = e i k f k v, t k Z d and the electric potential Then E, t = φ = φ, t = k Z d e i k φ k t. k Z d i kφ k t e i k = where E k t = i kφ k t. Denote e = k/ k, then k Z d E k t e i k, E k t = i kφ k t = Ẽ k t e where Ẽ k t = i k φ k t. Let v = α e + w where α, w e, and f k α, t = f k, e α, t = f k α e + w, t dw d 7

18 The linearized Vlasov equation implies that = t f k + v i k f k E k v f = t f k + iα k f k Ẽk α f An integration of the w variable on above equation yields t f k α, t + iα k f k α, t Ẽkf, e α =. 6 The Poisson equation implies k φ k t = f k v, t dv, d and thus i k Ẽ k t = f k α, t dα. 7 Equations 6 and 7 imply that f k α, t, Ẽk t e i k solves the linearized D Vlasov-Poisson equations at the homogeneous profile f, e α. Thus by the D representation formula in [9] and the Penrose stability condition 4, we have k G k y + i Ẽ k t = π k e i k yt dy. F e y + i Here, and G k y + i = P F e y + i = P f k α, α y dα + iπ f k y, f, e α α y dα + iπf, e y. By the Penrose stability condition 4 and P f, e α α y dα C d, s, b f H s,b d Corollary., there eists c > independent of k, such that k F e y + i c k. Then by the same proof of Proposition 4. in [9], t svẽ k t C k sv f k α, L Hv sv C k sv f k v,. H s,b d 8

19 So t sv E, t = C k Z d k Z d L t H +s+sv k +sv+s t s v Ẽ k t C f, v, H s L k s f k v, H s,b d Hv. Lemma. Assume f v H s,b d d, s >, b > d 4 and the Penrose stability condition 4 is satisfied for period tuple T,, T d. Let f, v, t, E, t be a solution of the Vlasov-Poisson system a-b with period tuple T,, T d. For any s, s v satisfying 6, there eists ε >, such that if then f t f H s Hv + t s v E, t L < ε, for all t, {t } H +s Cε. 8 Proof. Denote L to be the linearized operator corresponding to the linearized Vlasov-Poisson equation at f v,, and E is the mapping from f, v to E by the Poisson equation 4b.It follows from Lemma. that: For any s v s, if h, v H s,then + t sv E e tl h L t H C h, v +s H s Denote f t = f t f, then Thus f t = e tl f + and correspondingly t Hv sv,b t f = L f + E v f. Hv. 9 e t ul E v f u du = f lin t + f non t, E t = E f lin t + E f non t = E lin t + E non t. By the linear estimate 9, + t sv Elin, t L {t } H +s = + t sv E e tl f L t H C f H s Hv, +s 9

20 and + t sv Enon, t = C = C C C t L {t } H +s + t sv Enon, t t + t sv E H +s dt [e t ul E v H f u] du dt +s t + t sv + t u sv + u sv du t + u sv + t u sv E [ e t ul E v f u ] + u sv + t u sv E [ e t ul E v f u + u sv u + u sv E v f u + u sv E u Cε + t sv E, t + t u sv E [ e t ul E v f u H +s L {t } H +s H s In the above estimate, we use the fact that t du sv,b Hv f u H s. Hv + t u sv + u sv du C + t sv du H +s ] dudt H +s because of the assumption that s v >. The assumption 6 ensures that the following inequality is true E v f u H s Thus + t sv E L, t + t sv Elin, t C f H s Hv sv,b Hv {t } H +s L {t } H +s C E u H +s f u H s + + t sv Enon, t + Cε + t s v E L, t {t } H. +s By taking ε = C, we get the estimate 8. Hv. ] L {t } H +s H +s dudt dtdu

21 Theorem. follows from Lemma. and the time translation symmetry of the Vlasov-Poisson equation. Since the arguments are eactly the same as in the D case [9], we skip the details. As a corollary of Theorem., we get the following nonlinear instability result. Corollary. Assume f v H s,b d d, s >, b > d 4 and the Penrose stability condition 4 is satisfied for the period tuple T,, T d. For any s, s v satisfying 6, there eists ε > such that for any solution f, v, t, E, t of the Vlasov-Poisson system?? with period tuple T,, T d and E, not identically zero, the following is true: f T f H s Hv ε, for some T. We can also study the positive negative invariant structures near f v,, which are solutions f t, E t of nonlinear Vlasov-Poisson equation satisfying the conditions 7 for all t t. The net theorem shows that the electric field of these semi-invariant structures must decay when t + t. Theorem. Given a homogeneous profile f v H s,b d d, s >, b > d. 4 Assume that f v satisfies the Penrose stability condition 4 for T,, T d. Let f, v, t, E, v, t be a solution of?? in T d. For any s, s v satisfying 6, there eists ε >, such that if f t f H s Hv < ε, for all t or t, with f L,v <, v f,, v dvd <, T d d then E t, when t + or t. L Proof. By energy conservation, v f, v, t dvd + E, t T d d L = v f, v, dvd + E, L < C. T d d

22 Let j = vf dv. When d =,we have j t = vf t dv v dv f t L,v + v fdv v A A v A C f L,v A + v fdv C f 4L,v v 4 fdv, A by choosing A = v fdv/ f L,v 4. Thus j, t L 4 C v f dvd C. Since d E dt, t = L T d j, t E, t d j, t L 4 E, t L 4 C E, t H, and by Lemma. E, t H dt Cε, + t s dt + t s E, t dt H thus lim t E L, t eists and equals zero. When d =, the proof is very similar. The estimates become j, t 5 L 4 C, and d E dt, t j, t 5 L L 4 E, t L 5 C E, t H +s. The rest is the same. Acknowledgement This work is supported partly by the NSF grants DMS-9875 Lin and DMS-89 Zeng.

23 eferences [] Armstrong, T., Montgomery, D., Asymptotic state of the two-stream instability, J.Plasma. Physics,, part 4, [] Bernstein, I., Greene, J., Kruskal, M., Eact nonlinear plasma oscillations. Phys. ev. 8,, [] Danielson, J.., Anderegg, F. and Driscoll, C. F. Measurement of Landau Damping and the Evolution to a BGK Equilibrium, Physical eview Letters 9, [4] Demeio, L. and Zweifel, P. F. Numerical simulations of perturbed Vlasov equilibria, Phys. Fluids B, [5] Isichenko, M. B., Nonlinear Landau Damping in Collisionless Plasma and Inviscid Fluid, Phys. ev. Lett. 78, [6] Landau, L. On the vibration of the electronic plasma. J. Phys. USS, [7] Lin, Zhiwu, Instability of periodic BGK waves, Math. es. Letts., 8, [8] Lin, Zhiwu and Zeng, Chongchun, Inviscid dynamical structures near Couette flow, Arch. ational Mech. Anal.,, [9] Lin, Zhiwu and Zeng, Chongchun, Small BGK waves and nonlinear Landau damping, Commun. Math. Phys, 6, 9-,. [] Manfredi, Giovanni, Long-Time Behavior of Nonlinear Landau Damping, Physical eview Letters []. Morse, C. Nielson. One-, Two-, and Three-Dimensional Numerical Simulation of Two-Beam Plasmas, Physical eview Letters, : [] Mouhot, C., and Villani, C., On Landau damping, Acta Mathematica, to appear. [] M. M. Oppenheim, G. Vetoulis, D. L. Newman, M. V. Goldman, Evolution of electron phase-space holes in D, Geophysical esearch Letters 8, [4] Penrose, O., Electrostatic instability of a non-mawellian plasma. Phys. Fluids, [5] Strichartz, obert S, Multipliers on fractional Sobolev spaces, J. Math. Mech. 6,

Nonlinear instability of periodic BGK waves for Vlasov-Poisson system

Nonlinear instability of periodic BGK waves for Vlasov-Poisson system Nonlinear instability of periodic BGK waves for Vlasov-Poisson system Zhiwu Lin Courant Institute Abstract We investigate the nonlinear instability of periodic Bernstein-Greene-Kruskal(BGK waves. Starting

More information

Bernstein-Greene-Kruskal (BGK) Modes in a Three Dimensional Unmagnetized Plasma

Bernstein-Greene-Kruskal (BGK) Modes in a Three Dimensional Unmagnetized Plasma Bernstein-Greene-Kruskal (BGK) Modes in a Three Dimensional Unmagnetized Plasma C. S. Ng & A. Bhattacharjee Space Science Center University of New Hampshire Space Science Center Seminar, March 9, 2005

More information

Piecewise Smooth Solutions to the Burgers-Hilbert Equation

Piecewise 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 information

LANDAU DAMPING IN SOBOLEV SPACES FOR THE VLASOV-HMF MODEL

LANDAU DAMPING IN SOBOLEV SPACES FOR THE VLASOV-HMF MODEL LANDAU DAMPING IN SOBOLEV SPACES FOR THE VLASOV-HMF MODEL ERWAN FAOU AND FRÉDÉRIC ROUSSET Abstract. We consider the Vlasov-HMF (Hamiltonian Mean-Field) model. We consider solutions starting in a small

More information

Small BGK waves and nonlinear Landau damping

Small BGK waves and nonlinear Landau damping Small BGK waes and nonlinear Landau damping hiwu Lin and Chongchun eng School of Mathematics Georgia Institute of Technology Atlanta, GA 333, USA Abstract Consider D Vlaso-poisson system with a ed ion

More information

hal , version 1-22 Nov 2009

hal , 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 information

LARGE TIME BEHAVIOR OF THE RELATIVISTIC VLASOV MAXWELL SYSTEM IN LOW SPACE DIMENSION

LARGE 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 information

Stability for a class of nonlinear pseudo-differential equations

Stability for a class of nonlinear pseudo-differential equations Stability for a class of nonlinear pseudo-differential equations Michael Frankel Department of Mathematics, Indiana University - Purdue University Indianapolis Indianapolis, IN 46202-3216, USA Victor Roytburd

More information

A COUNTEREXAMPLE TO AN ENDPOINT BILINEAR STRICHARTZ INEQUALITY TERENCE TAO. t L x (R R2 ) f L 2 x (R2 )

A 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 information

Accurate representation of velocity space using truncated Hermite expansions.

Accurate representation of velocity space using truncated Hermite expansions. Accurate representation of velocity space using truncated Hermite expansions. Joseph Parker Oxford Centre for Collaborative Applied Mathematics Mathematical Institute, University of Oxford Wolfgang Pauli

More information

Anomalous transport of particles in Plasma physics

Anomalous 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 information

Nonlinear Modulational Instability of Dispersive PDE Models

Nonlinear Modulational Instability of Dispersive PDE Models Nonlinear Modulational Instability of Dispersive PDE Models Jiayin Jin, Shasha Liao, and Zhiwu Lin Georgia Tech ICERM workshop on water waves, 4/28/2017 Jiayin Jin, Shasha Liao, and Zhiwu Lin Georgia Tech

More information

u( x) = g( y) ds y ( 1 ) U solves u = 0 in U; u = 0 on U. ( 3)

u( x) = g( y) ds y ( 1 ) U solves u = 0 in U; u = 0 on U. ( 3) M ath 5 2 7 Fall 2 0 0 9 L ecture 4 ( S ep. 6, 2 0 0 9 ) Properties and Estimates of Laplace s and Poisson s Equations In our last lecture we derived the formulas for the solutions of Poisson s equation

More information

Traveling waves of a kinetic transport model for the KPP-Fisher equation

Traveling waves of a kinetic transport model for the KPP-Fisher equation Traveling waves of a kinetic transport model for the KPP-Fisher equation Christian Schmeiser Universität Wien and RICAM homepage.univie.ac.at/christian.schmeiser/ Joint work with C. Cuesta (Bilbao), S.

More information

Time-asymptotic wave propagation in collisionless plasmas

Time-asymptotic wave propagation in collisionless plasmas PHYSICAL EVIEW E 68, 646 3 Time-asymptotic wave propagation in collisionless plasmas Carlo Lancellotti Department of Mathematics, City University of New York CSI, Staten Island, New Yor34, USA J. J. Dorning

More information

Waves in plasma. Denis Gialis

Waves in plasma. Denis Gialis Waves in plasma Denis Gialis This is a short introduction on waves in a non-relativistic plasma. We will consider a plasma of electrons and protons which is fully ionized, nonrelativistic and homogeneous.

More information

We denote the space of distributions on Ω by D ( Ω) 2.

We denote the space of distributions on Ω by D ( Ω) 2. Sep. 1 0, 008 Distributions Distributions are generalized functions. Some familiarity with the theory of distributions helps understanding of various function spaces which play important roles in the study

More information

Downloaded 08/07/18 to Redistribution subject to SIAM license or copyright; see

Downloaded 08/07/18 to Redistribution subject to SIAM license or copyright; see SIAM J. MATH. ANAL. Vol. 5, No. 4, pp. 43 4326 c 28 Society for Industrial and Applied Mathematics Downloaded 8/7/8 to 3.25.254.96. Redistribution subject to SIAM license or copyright; see http://www.siam.org/journals/ojsa.php

More information

ROTH S THEOREM THE FOURIER ANALYTIC APPROACH NEIL LYALL

ROTH S THEOREM THE FOURIER ANALYTIC APPROACH NEIL LYALL ROTH S THEOREM THE FOURIER AALYTIC APPROACH EIL LYALL Roth s Theorem. Let δ > 0 and ep ep(cδ ) for some absolute constant C. Then any A {0,,..., } of size A = δ necessarily contains a (non-trivial) arithmetic

More information

Сollisionless damping of electron waves in non-maxwellian plasma 1

Сollisionless damping of electron waves in non-maxwellian plasma 1 http:/arxi.org/physics/78.748 Сollisionless damping of electron waes in non-mawellian plasma V. N. Soshnio Plasma Physics Dept., All-Russian Institute of Scientific and Technical Information of the Russian

More information

Theory of PDE Homework 2

Theory of PDE Homework 2 Theory of PDE Homework 2 Adrienne Sands April 18, 2017 In the following exercises we assume the coefficients of the various PDE are smooth and satisfy the uniform ellipticity condition. R n is always an

More information

The 2D Magnetohydrodynamic Equations with Partial Dissipation. Oklahoma State University

The 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 information

A Quasi-Linear Parabolic Partial Differential Equation with Accretive Property

A Quasi-Linear Parabolic Partial Differential Equation with Accretive Property ONLINE ISSN 8-749 : Volume 3, Issue, 433-438 A Quasi-Linear Parabolic Partial Differential Equation with Accretive Property Aminu U. Bawa *, Micheal O. Egwurube and Murtala M. Ahmad 3 Department of Computer

More information

On the local existence for an active scalar equation in critical regularity setting

On the local existence for an active scalar equation in critical regularity setting On the local existence for an active scalar equation in critical regularity setting Walter Rusin Department of Mathematics, Oklahoma State University, Stillwater, OK 7478 Fei Wang Department of Mathematics,

More information

Dissipative quasi-geostrophic equations with L p data

Dissipative 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 information

ALMOST EXPONENTIAL DECAY NEAR MAXWELLIAN

ALMOST 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 information

Properties of the Scattering Transform on the Real Line

Properties of the Scattering Transform on the Real Line Journal of Mathematical Analysis and Applications 58, 3 43 (001 doi:10.1006/jmaa.000.7375, available online at http://www.idealibrary.com on Properties of the Scattering Transform on the Real Line Michael

More information

This article appeared in a journal published by Elsevier. The attached copy is furnished to the author for internal non-commercial research and

This article appeared in a journal published by Elsevier. The attached copy is furnished to the author for internal non-commercial research and This article appeared in a journal published by Elsevier. The attached copy is furnished to the author for internal non-commercial research and education use, including for instruction at the authors institution

More information

IMA Preprint Series # 2260

IMA Preprint Series # 2260 LONG-TIME BEHAVIOR OF HYDRODYNAMIC SYSTEMS MODELING THE NEMATIC LIUID CRYSTALS By Hao Wu Xiang Xu and Chun Liu IMA Preprint Series # 2260 ( June 2009 ) INSTITUTE FOR MATHEMATICS AND ITS APPLICATIONS UNIVERSITY

More information

Chapter 1. Introduction to Nonlinear Space Plasma Physics

Chapter 1. Introduction to Nonlinear Space Plasma Physics Chapter 1. Introduction to Nonlinear Space Plasma Physics The goal of this course, Nonlinear Space Plasma Physics, is to explore the formation, evolution, propagation, and characteristics of the large

More information

arxiv: v1 [math.ap] 28 Mar 2014

arxiv: v1 [math.ap] 28 Mar 2014 GROUNDSTATES OF NONLINEAR CHOQUARD EQUATIONS: HARDY-LITTLEWOOD-SOBOLEV CRITICAL EXPONENT VITALY MOROZ AND JEAN VAN SCHAFTINGEN arxiv:1403.7414v1 [math.ap] 28 Mar 2014 Abstract. We consider nonlinear Choquard

More information

OPTIMAL CONVERGENCE RATES FOR THE COMPRESSIBLE NAVIER-STOKES EQUATIONS WITH POTENTIAL FORCES

OPTIMAL 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 information

Oscillation Criteria for Certain nth Order Differential Equations with Deviating Arguments

Oscillation Criteria for Certain nth Order Differential Equations with Deviating Arguments Journal of Mathematical Analysis Applications 6, 601 6 001) doi:10.1006/jmaa.001.7571, available online at http://www.idealibrary.com on Oscillation Criteria for Certain nth Order Differential Equations

More information

Wave breaking in the short-pulse equation

Wave breaking in the short-pulse equation Dynamics of PDE Vol.6 No.4 29-3 29 Wave breaking in the short-pulse equation Yue Liu Dmitry Pelinovsky and Anton Sakovich Communicated by Y. Charles Li received May 27 29. Abstract. Sufficient conditions

More information

DIFFERENTIAL EQUATIONS WITH A DIFFERENCE QUOTIENT

DIFFERENTIAL EQUATIONS WITH A DIFFERENCE QUOTIENT Electronic Journal of Differential Equations, Vol. 217 (217, No. 227, pp. 1 42. ISSN: 172-6691. URL: http://ejde.math.tstate.edu or http://ejde.math.unt.edu DIFFERENTIAL EQUATIONS WITH A DIFFERENCE QUOTIENT

More information

Convergence Rate of Nonlinear Switched Systems

Convergence Rate of Nonlinear Switched Systems Convergence Rate of Nonlinear Switched Systems Philippe JOUAN and Saïd NACIRI arxiv:1511.01737v1 [math.oc] 5 Nov 2015 January 23, 2018 Abstract This paper is concerned with the convergence rate of the

More information

Existence of Secondary Bifurcations or Isolas for PDEs

Existence of Secondary Bifurcations or Isolas for PDEs Existence of Secondary Bifurcations or Isolas for PDEs Marcio Gameiro Jean-Philippe Lessard Abstract In this paper, we introduce a method to conclude about the existence of secondary bifurcations or isolas

More information

= 1 2 x (x 1) + 1 {x} (1 {x}). [t] dt = 1 x (x 1) + O (1), [t] dt = 1 2 x2 + O (x), (where the error is not now zero when x is an integer.

= 1 2 x (x 1) + 1 {x} (1 {x}). [t] dt = 1 x (x 1) + O (1), [t] dt = 1 2 x2 + O (x), (where the error is not now zero when x is an integer. Problem Sheet,. i) Draw the graphs for [] and {}. ii) Show that for α R, α+ α [t] dt = α and α+ α {t} dt =. Hint Split these integrals at the integer which must lie in any interval of length, such as [α,

More information

Kinetic theory of gases

Kinetic theory of gases Kinetic theory of gases Toan T. Nguyen Penn State University http://toannguyen.org http://blog.toannguyen.org Graduate Student seminar, PSU Jan 19th, 2017 Fall 2017, I teach a graduate topics course: same

More information

Lectures on basic plasma physics: Kinetic approach

Lectures on basic plasma physics: Kinetic approach Lectures on basic plasma physics: Kinetic approach Department of applied physics, Aalto University April 30, 2014 Motivation Layout 1 Motivation 2 Boltzmann equation (a nasty bastard) 3 Vlasov equation

More information

Partial Differential Equations, 2nd Edition, L.C.Evans Chapter 5 Sobolev Spaces

Partial Differential Equations, 2nd Edition, L.C.Evans Chapter 5 Sobolev Spaces Partial Differential Equations, nd Edition, L.C.Evans Chapter 5 Sobolev Spaces Shih-Hsin Chen, Yung-Hsiang Huang 7.8.3 Abstract In these exercises always denote an open set of with smooth boundary. As

More information

arxiv: v1 [nlin.ps] 18 Sep 2008

arxiv: v1 [nlin.ps] 18 Sep 2008 Asymptotic two-soliton solutions solutions in the Fermi-Pasta-Ulam model arxiv:0809.3231v1 [nlin.ps] 18 Sep 2008 Aaron Hoffman and C.E. Wayne Boston University Department of Mathematics and Statistics

More information

Taylor Series and Asymptotic Expansions

Taylor Series and Asymptotic Expansions Taylor Series and Asymptotic Epansions The importance of power series as a convenient representation, as an approimation tool, as a tool for solving differential equations and so on, is pretty obvious.

More information

HARDY INEQUALITIES WITH BOUNDARY TERMS. x 2 dx u 2 dx. (1.2) u 2 = u 2 dx.

HARDY INEQUALITIES WITH BOUNDARY TERMS. x 2 dx u 2 dx. (1.2) u 2 = u 2 dx. Electronic Journal of Differential Equations, Vol. 003(003), No. 3, pp. 1 8. ISSN: 107-6691. UL: http://ejde.math.swt.edu or http://ejde.math.unt.edu ftp ejde.math.swt.edu (login: ftp) HADY INEQUALITIES

More information

Author(s) Huang, Feimin; Matsumura, Akitaka; Citation Osaka Journal of Mathematics. 41(1)

Author(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 information

Vlasov simulations of electron holes driven by particle distributions from PIC reconnection simulations with a guide field

Vlasov simulations of electron holes driven by particle distributions from PIC reconnection simulations with a guide field GEOPHYSICAL RESEARCH LETTERS, VOL. 35, L22109, doi:10.1029/2008gl035608, 2008 Vlasov simulations of electron holes driven by particle distributions from PIC reconnection simulations with a guide field

More information

Classical solutions for the quasi-stationary Stefan problem with surface tension

Classical solutions for the quasi-stationary Stefan problem with surface tension Classical solutions for the quasi-stationary Stefan problem with surface tension Joachim Escher, Gieri Simonett We show that the quasi-stationary two-phase Stefan problem with surface tension has a unique

More information

Decay in Time of Incompressible Flows

Decay 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 information

NONLOCAL DIFFUSION EQUATIONS

NONLOCAL DIFFUSION EQUATIONS NONLOCAL DIFFUSION EQUATIONS JULIO D. ROSSI (ALICANTE, SPAIN AND BUENOS AIRES, ARGENTINA) jrossi@dm.uba.ar http://mate.dm.uba.ar/ jrossi 2011 Non-local diffusion. The function J. Let J : R N R, nonnegative,

More information

analysis for transport equations and applications

analysis for transport equations and applications Multi-scale analysis for transport equations and applications Mihaï BOSTAN, Aurélie FINOT University of Aix-Marseille, FRANCE mihai.bostan@univ-amu.fr Numerical methods for kinetic equations Strasbourg

More information

Long time behavior of solutions of Vlasov-like Equations

Long time behavior of solutions of Vlasov-like Equations Long time behavior of solutions of Vlasov-like Equations Emanuele Caglioti "Sapienza" Università di Roma (Mathematics) Journée Claude Bardos - Laboratoire Jacques Louis Lions - 2011 Equations Vlasov Poisson

More information

New approach to study the van der Pol equation for large damping

New approach to study the van der Pol equation for large damping Electronic Journal of Qualitative Theor of Differential Equations 2018, No. 8, 1 10; https://doi.org/10.1422/ejqtde.2018.1.8 www.math.u-szeged.hu/ejqtde/ New approach to stud the van der Pol equation for

More information

REMARKS ON THE VANISHING OBSTACLE LIMIT FOR A 3D VISCOUS INCOMPRESSIBLE FLUID

REMARKS ON THE VANISHING OBSTACLE LIMIT FOR A 3D VISCOUS INCOMPRESSIBLE FLUID REMARKS ON THE VANISHING OBSTACLE LIMIT FOR A 3D VISCOUS INCOMPRESSIBLE FLUID DRAGOŞ IFTIMIE AND JAMES P. KELLIHER Abstract. In [Math. Ann. 336 (2006), 449-489] the authors consider the two dimensional

More information

REGULARITY OF THE MINIMIZER FOR THE D-WAVE GINZBURG-LANDAU ENERGY

REGULARITY OF THE MINIMIZER FOR THE D-WAVE GINZBURG-LANDAU ENERGY METHODS AND APPLICATIONS OF ANALYSIS. c 2003 International Press Vol. 0, No., pp. 08 096, March 2003 005 REGULARITY OF THE MINIMIZER FOR THE D-WAVE GINZBURG-LANDAU ENERGY TAI-CHIA LIN AND LIHE WANG Abstract.

More information

TADAHIRO OH 0, 3 8 (T R), (1.5) The result in [2] is in fact stated for time-periodic functions: 0, 1 3 (T 2 ). (1.4)

TADAHIRO OH 0, 3 8 (T R), (1.5) The result in [2] is in fact stated for time-periodic functions: 0, 1 3 (T 2 ). (1.4) PERIODIC L 4 -STRICHARTZ ESTIMATE FOR KDV TADAHIRO OH 1. Introduction In [], Bourgain proved global well-posedness of the periodic KdV in L T): u t + u xxx + uu x 0, x, t) T R. 1.1) The key ingredient

More information

ASYMPTOTIC STRUCTURE FOR SOLUTIONS OF THE NAVIER STOKES EQUATIONS. Tian Ma. Shouhong Wang

ASYMPTOTIC STRUCTURE FOR SOLUTIONS OF THE NAVIER STOKES EQUATIONS. Tian Ma. Shouhong Wang DISCRETE AND CONTINUOUS Website: http://aimsciences.org DYNAMICAL SYSTEMS Volume 11, Number 1, July 004 pp. 189 04 ASYMPTOTIC STRUCTURE FOR SOLUTIONS OF THE NAVIER STOKES EQUATIONS Tian Ma Department of

More information

1 Lyapunov theory of stability

1 Lyapunov theory of stability M.Kawski, APM 581 Diff Equns Intro to Lyapunov theory. November 15, 29 1 1 Lyapunov theory of stability Introduction. Lyapunov s second (or direct) method provides tools for studying (asymptotic) stability

More information

The 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 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 information

Microlocal Methods in X-ray Tomography

Microlocal Methods in X-ray Tomography Microlocal Methods in X-ray Tomography Plamen Stefanov Purdue University Lecture I: Euclidean X-ray tomography Mini Course, Fields Institute, 2012 Plamen Stefanov (Purdue University ) Microlocal Methods

More information

SYLLABUS FOR ENTRANCE EXAMINATION NANYANG TECHNOLOGICAL UNIVERSITY FOR INTERNATIONAL STUDENTS A-LEVEL MATHEMATICS

SYLLABUS FOR ENTRANCE EXAMINATION NANYANG TECHNOLOGICAL UNIVERSITY FOR INTERNATIONAL STUDENTS A-LEVEL MATHEMATICS SYLLABUS FOR ENTRANCE EXAMINATION NANYANG TECHNOLOGICAL UNIVERSITY FOR INTERNATIONAL STUDENTS A-LEVEL MATHEMATICS STRUCTURE OF EXAMINATION PAPER. There will be one -hour paper consisting of 4 questions..

More information

Decay rates for partially dissipative hyperbolic systems

Decay 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 information

Landau Damping Simulation Models

Landau Damping Simulation Models Landau Damping Simulation Models Hua-sheng XIE (u) huashengxie@gmail.com) Department of Physics, Institute for Fusion Theory and Simulation, Zhejiang University, Hangzhou 310027, P.R.China Oct. 9, 2013

More information

2. Higher-order Linear ODE s

2. Higher-order Linear ODE s 2. Higher-order Linear ODE s 2A. Second-order Linear ODE s: General Properties 2A-1. On the right below is an abbreviated form of the ODE on the left: (*) y + p()y + q()y = r() Ly = r() ; where L is the

More information

Nonlinear eigenvalue problems for higher order model equations

Nonlinear eigenvalue problems for higher order model equations Nonlinear eigenvalue problems for higher order model equations L.A. Peletier Mathematical Institute, Leiden University PB 951, 3 RA Leiden, The Netherlands e-mail: peletier@math.leidenuniv.nl 1 Introduction

More information

Integration of Vlasov-type equations

Integration of Vlasov-type equations Alexander Ostermann University of Innsbruck, Austria Joint work with Lukas Einkemmer Verona, April/May 2017 Plasma the fourth state of matter 99% of the visible matter in the universe is made of plasma

More information

Lecture 2: A Strange Term Coming From Nowhere

Lecture 2: A Strange Term Coming From Nowhere Lecture 2: A Strange Term Coming From Nowhere Christophe Prange February 9, 2016 In this lecture, we consider the Poisson equation with homogeneous Dirichlet boundary conditions { u ε = f, x ε, u ε = 0,

More information

arxiv:math/ v1 [math.ap] 28 Oct 2005

arxiv: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 information

On semilinear elliptic equations with nonlocal nonlinearity

On semilinear elliptic equations with nonlocal nonlinearity On semilinear elliptic equations with nonlocal nonlinearity Shinji Kawano Department of Mathematics Hokkaido University Sapporo 060-0810, Japan Abstract We consider the problem 8 < A A + A p ka A 2 dx

More information

Rate of convergence for certain families of summation integral type operators

Rate of convergence for certain families of summation integral type operators J Math Anal Appl 296 24 68 618 wwwelseviercom/locate/jmaa Rate of convergence for certain families of summation integral type operators Vijay Gupta a,,mkgupta b a School of Applied Sciences, Netaji Subhas

More information

Interfacial waves in steady and oscillatory, two-layer Couette flows

Interfacial waves in steady and oscillatory, two-layer Couette flows Interfacial waves in steady and oscillatory, two-layer Couette flows M. J. McCready Department of Chemical Engineering University of Notre Dame Notre Dame, IN 46556 Page 1 Acknowledgments Students: M.

More information

Instability, index theorem, and exponential trichotomy for Linear Hamiltonian PDEs

Instability, index theorem, and exponential trichotomy for Linear Hamiltonian PDEs Instability, index theorem, and exponential trichotomy for Linear Hamiltonian PDEs Zhiwu Lin and Chongchun Zeng School of Mathematics Georgia Institute of Technology Atlanta, GA 30332, USA Abstract Consider

More information

Asymptotic behavior of infinity harmonic functions near an isolated singularity

Asymptotic behavior of infinity harmonic functions near an isolated singularity Asymptotic behavior of infinity harmonic functions near an isolated singularity Ovidiu Savin, Changyou Wang, Yifeng Yu Abstract In this paper, we prove if n 2 x 0 is an isolated singularity of a nonegative

More information

arxiv: v1 [math.ap] 28 Apr 2009

arxiv: 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 information

Dynamical Systems. 1.0 Ordinary Differential Equations. 2.0 Dynamical Systems

Dynamical Systems. 1.0 Ordinary Differential Equations. 2.0 Dynamical Systems . Ordinary Differential Equations. An ordinary differential equation (ODE, for short) is an equation involves a dependent variable and its derivatives with respect to an independent variable. When we speak

More information

Conservation law equations : problem set

Conservation law equations : problem set Conservation law equations : problem set Luis Silvestre For Isaac Neal and Elia Portnoy in the 2018 summer bootcamp 1 Method of characteristics For the problems in this section, assume that the solutions

More information

Week 6 Notes, Math 865, Tanveer

Week 6 Notes, Math 865, Tanveer Week 6 Notes, Math 865, Tanveer. Energy Methods for Euler and Navier-Stokes Equation We will consider this week basic energy estimates. These are estimates on the L 2 spatial norms of the solution u(x,

More information

arxiv: v1 [math.ap] 28 May 2009

arxiv: v1 [math.ap] 28 May 2009 Wave breaking in the short-pulse equation arxiv:95.4668v [math.ap] 8 May 9 Yue Liu, Dmitry Pelinovsky, and Anton Sakovich Department of Mathematics, University of Teas at Arlington, Arlington, TX, 769,

More information

Math 256: Applied Differential Equations: Final Review

Math 256: Applied Differential Equations: Final Review Math 256: Applied Differential Equations: Final Review Chapter 1: Introduction, Sec 1.1, 1.2, 1.3 (a) Differential Equation, Mathematical Model (b) Direction (Slope) Field, Equilibrium Solution (c) Rate

More information

Non-linear Stability of Steady Geophysical Flows

Non-linear Stability of Steady Geophysical Flows Non-linear Stability of Steady Geophysical Flows Di Qi, and Andrew J. Majda Courant Institute of Mathematical Sciences Fall 2016 Advanced Topics in Applied Math Di Qi, and Andrew J. Majda (CIMS) Non-linear

More information

RENORMALIZATION OF DYSON S VECTOR-VALUED HIERARCHICAL MODEL AT LOW TEMPERATURES

RENORMALIZATION OF DYSON S VECTOR-VALUED HIERARCHICAL MODEL AT LOW TEMPERATURES RENORMALIZATION OF DYSON S VECTOR-VALUED HIERARCHICAL MODEL AT LOW TEMPERATURES P. M. Bleher (1) and P. Major (2) (1) Keldysh Institute of Applied Mathematics of the Soviet Academy of Sciences Moscow (2)

More information

You may not start to read the questions printed on the subsequent pages until instructed to do so by the Invigilator.

You may not start to read the questions printed on the subsequent pages until instructed to do so by the Invigilator. MATHEMATICAL TRIPOS Part III Thursday 1 June 2006 1.30 to 4.30 PAPER 76 NONLINEAR CONTINUUM MECHANICS Attempt FOUR questions. There are SIX questions in total. The questions carry equal weight. STATIONERY

More information

AVERAGE RECIPROCALS OF THE ORDER OF a MODULO n

AVERAGE RECIPROCALS OF THE ORDER OF a MODULO n AVERAGE RECIPROCALS OF THE ORDER OF a MODULO n KIM, SUNGJIN Abstract Let a > be an integer Denote by l an the multiplicative order of a modulo integers n We prove that l = an Oa ep 2 + o log log, n,n,a=

More information

Study of Laser Plasma Interactions Using an Eulerian Vlasov Code

Study of Laser Plasma Interactions Using an Eulerian Vlasov Code PSFC/JA-04-6 Study of Laser Plasma Interactions Using an Eulerian Vlasov Code D. J. Strozzi, M. M. Shoucri*, and A. Bers March 2004 Plasma Science and Fusion Center Massachusetts Institute of Technology

More information

Identification Methods for Structural Systems. Prof. Dr. Eleni Chatzi System Stability - 26 March, 2014

Identification Methods for Structural Systems. Prof. Dr. Eleni Chatzi System Stability - 26 March, 2014 Prof. Dr. Eleni Chatzi System Stability - 26 March, 24 Fundamentals Overview System Stability Assume given a dynamic system with input u(t) and output x(t). The stability property of a dynamic system can

More information

Gaussian beams for high frequency waves: some recent developments

Gaussian beams for high frequency waves: some recent developments Gaussian beams for high frequency waves: some recent developments Jianliang Qian Department of Mathematics, Michigan State University, Michigan International Symposium on Geophysical Imaging with Localized

More information

1 Introduction and Notations he simplest model to describe a hot dilute electrons moving in a fixed ion background is the following one dimensional Vl

1 Introduction and Notations he simplest model to describe a hot dilute electrons moving in a fixed ion background is the following one dimensional Vl Numerical Study on Landau Damping ie Zhou 1 School of Mathematical Sciences Peking University, Beijing 100871, China Yan Guo 2 Lefschetz Center for Dynamical Systems Division of Applied Mathematics Brown

More information

A global solution curve for a class of free boundary value problems arising in plasma physics

A global solution curve for a class of free boundary value problems arising in plasma physics A global solution curve for a class of free boundary value problems arising in plasma physics Philip Korman epartment of Mathematical Sciences University of Cincinnati Cincinnati Ohio 4522-0025 Abstract

More information

An introduction to Birkhoff normal form

An introduction to Birkhoff normal form An introduction to Birkhoff normal form Dario Bambusi Dipartimento di Matematica, Universitá di Milano via Saldini 50, 0133 Milano (Italy) 19.11.14 1 Introduction The aim of this note is to present an

More information

Laplace s Equation. Chapter Mean Value Formulas

Laplace s Equation. Chapter Mean Value Formulas Chapter 1 Laplace s Equation Let be an open set in R n. A function u C 2 () is called harmonic in if it satisfies Laplace s equation n (1.1) u := D ii u = 0 in. i=1 A function u C 2 () is called subharmonic

More information

Excitation and Decay of Electron Acoustic Waves

Excitation and Decay of Electron Acoustic Waves Excitation and Decay of Electron Acoustic Waves Francesco Valentini, Thomas M. O Neil and Daniel H. E. Dubin Dipt. di Fisica and INFM, Univ. della Calabria, 87036 Rende, Italy Department of Physics, University

More information

Summer College on Plasma Physics. 30 July - 24 August, The particle-in-cell simulation method: Concept and limitations

Summer College on Plasma Physics. 30 July - 24 August, The particle-in-cell simulation method: Concept and limitations 1856-30 2007 Summer College on Plasma Physics 30 July - 24 August, 2007 The particle-in-cell M. E. Dieckmann Institut fuer Theoretische Physik IV, Ruhr-Universitaet, Bochum, Germany The particle-in-cell

More information

Courant Instructor, Courant Institute of New York University,

Courant Instructor, Courant Institute of New York University, Zhiwu Lin Contact Information School of Mathematics Offi ce: (404)-8949232 Georgia Institute of Technology Atlanta, GA 30332, USA E-mail: zlin@math.gatech.edu Research Interests Fluid dynamics, kinetic

More information

High Order Semi-Lagrangian WENO scheme for Vlasov Equations

High Order Semi-Lagrangian WENO scheme for Vlasov Equations High Order WENO scheme for Equations Department of Mathematical and Computer Science Colorado School of Mines joint work w/ Andrew Christlieb Supported by AFOSR. Computational Mathematics Seminar, UC Boulder

More information

Chapter One. The Calderón-Zygmund Theory I: Ellipticity

Chapter One. The Calderón-Zygmund Theory I: Ellipticity Chapter One The Calderón-Zygmund Theory I: Ellipticity Our story begins with a classical situation: convolution with homogeneous, Calderón- Zygmund ( kernels on R n. Let S n 1 R n denote the unit sphere

More information

Chapter 9. Derivatives. Josef Leydold Mathematical Methods WS 2018/19 9 Derivatives 1 / 51. f x. (x 0, f (x 0 ))

Chapter 9. Derivatives. Josef Leydold Mathematical Methods WS 2018/19 9 Derivatives 1 / 51. f x. (x 0, f (x 0 )) Chapter 9 Derivatives Josef Leydold Mathematical Methods WS 208/9 9 Derivatives / 5 Difference Quotient Let f : R R be some function. The the ratio f = f ( 0 + ) f ( 0 ) = f ( 0) 0 is called difference

More information

On Pressure Stabilization Method and Projection Method for Unsteady Navier-Stokes Equations 1

On Pressure Stabilization Method and Projection Method for Unsteady Navier-Stokes Equations 1 On Pressure Stabilization Method and Projection Method for Unsteady Navier-Stokes Equations 1 Jie Shen Department of Mathematics, Penn State University University Park, PA 1682 Abstract. We present some

More information

ATTRACTORS FOR SEMILINEAR PARABOLIC PROBLEMS WITH DIRICHLET BOUNDARY CONDITIONS IN VARYING DOMAINS. Emerson A. M. de Abreu Alexandre N.

ATTRACTORS FOR SEMILINEAR PARABOLIC PROBLEMS WITH DIRICHLET BOUNDARY CONDITIONS IN VARYING DOMAINS. Emerson A. M. de Abreu Alexandre N. ATTRACTORS FOR SEMILINEAR PARABOLIC PROBLEMS WITH DIRICHLET BOUNDARY CONDITIONS IN VARYING DOMAINS Emerson A. M. de Abreu Alexandre N. Carvalho Abstract Under fairly general conditions one can prove that

More information

ME 391 Mechanical Engineering Analysis

ME 391 Mechanical Engineering Analysis Solve the following differential equations. ME 9 Mechanical Engineering Analysis Eam # Practice Problems Solutions. e u 5sin() with u(=)= We may rearrange this equation to e u 5sin() and note that it is

More information

Introduction. Chapter Plasma: definitions

Introduction. Chapter Plasma: definitions Chapter 1 Introduction 1.1 Plasma: definitions A plasma is a quasi-neutral gas of charged and neutral particles which exhibits collective behaviour. An equivalent, alternative definition: A plasma is a

More information

Math Ordinary Differential Equations

Math Ordinary Differential Equations Math 411 - Ordinary Differential Equations Review Notes - 1 1 - Basic Theory A first order ordinary differential equation has the form x = f(t, x) (11) Here x = dx/dt Given an initial data x(t 0 ) = x

More information