ON THE VANISHING VISCOSITY LIMIT IN A DISK
|
|
- Silvester Houston
- 5 years ago
- Views:
Transcription
1 ON THE VANISHING VISCOSITY LIMIT IN A DISK JAMES P. KELLIHER Abstract. Let u be a solution to the Navier-Stokes equations in the unit disk with no-slip boundary conditions and viscosity >, and let u be a smooth solution to the Euler equations. We say that the vanishing viscosity limit holds on [, T ] if u converges to u in L ([, T ]; L 2 ). We show that a necessary and sufficient condition for the vanishing viscosity limit to hold is the vanishing with the viscosity of the time-space average of the energy of u in a boundary layer of width proportional to due to the modes (eigenfunctions of the Stokes operator) whose frequencies in the radial or the tangential direction lie between L() and M(). Here, L() must be of order less than / and M() must be of order greater than /. Active links added 2 Jan 2; sign error in (4.3) corrected 8 Sept 24. Introduction In the presence of a boundary, the question of whether solutions of the Navier-Stokes equations with no-slip boundary conditions converge to a solution of the Euler equations as the viscosity vanishes the so-called vanishing viscosity limit is very difficult. The convergence of most interest is of the velocities, uniformly over finite time and L 2 in space. Except in the very special case of radially symmetric initial vorticity in a disk, where convergence is known to hold (see Theorem 6.), the question of convergence or the lack thereof is unresolved for nonzero initial velocity in a bounded domain. (For a half-space with analytic initial data, the vanishing viscosity limit is shown to hold in [4].) Tosio Kato in [6] gave necessary and sufficient conditions on the velocity u of the Navier-Stokes equations for the vanishing viscosity limit to hold. The most interesting of these is that T u(t) 2 L 2 (Γ c) dt as, where Γ c is the boundary strip of width c with c > fixed but arbitrary. Making only a small change to Kato s proof, it is possible to replace u with the vorticity ω = ω(u) = u 2 2 u, giving Equation (2.3) (see [7]). (The necessity of Equation (2.3) is immediate from Kato s condition, but 2 Mathematics Subject Classification. Primary 76D5, 76B99, 76D99. Key words and phrases. Vanishing viscosity, Navier-Stokes equations, Euler equations.
2 2 JAMES P. KELLIHER because we do not have a boundary condition on the inner boundary of Γ c the sufficiency of the condition requires proof.) Other necessary and sufficient conditions were established by Teman and Wang in [6] and [7]. These are the conditions in Equation (2.5) and Equation (2.6) of Theorem 2.3, and involve only the derivatives in the directions tangential to the boundary of either the tangential or normal components of the velocity, though for a slightly larger boundary layer. Finally, a condition that requires that the average energy density in the boundary layer of the same width as Kato s vanish with viscosity, Equation (2.7), is proven in [7]. All these conditions (which apply to a bounded domain in dimensions 2 and higher) are summarized in Theorem 2.3. We consider the issue of vanishing viscosity in the (unit) disk and look for weaker necessary and sufficient conditions for the limit to hold. The reason for working in the disk is that the simple geometry allows us to make quite explicit calculations using the eigenfunctions of the Stokes operator, which are composed of Bessel functions of the first kind. In a sense, this connects the energy method with the geometry. What we find is that we need only consider certain ranges of frequencies (or equivalently, length scales) in the various conditions: this is Theorem 2.4. Although Theorem 2.4 is specific to the disk, there is no hydrodynamical reason to expect the disk to be special as regards the vanishing viscosity limit, so one would expect a version of the theorem to apply to all sufficiently smooth bounded domains in R 2, and probably in higher dimensions as well. We discuss this issue more fully in Remark (2.). In [3], Cheng and Wang obtain a result regarding vanishing viscosity in two dimensions analogous to Equation (2.5) and Equation (2.6). Their result applies to an approximating sequence to a solution of the Navier- Stokes equations as the viscosity vanishes, whereas our result applies to the necessary and sufficient condition for the vanishing viscosity limit to hold. While for the other conditions in Theorem 2.4 we use very different techniques than those in [3], our proof of the necessity and sufficiency of Equation (2.5) and Equation (2.6) uses the key inequality in their paper. Section 7 contains a brief comparison between the two results. In [3], the authors consider the Stokes problem (linearized Navier-Stokes equations) external to a disk with time-varying Dirichlet boundary conditions, showing that the vanishing viscosity limit holds. In fact, they do much more than this, giving an explicit construction of the solution to the Stokes problem and showing that it can be decomposed into the sum of the solution to the linearized Euler equations, the solution to the associated Prantdl equations, and a small correction term. The symmetry of the geometry allows the authors of [3] to construct the solutions in an explicit form (involving Bessel functions of the first and second kind). The nonlinear term in the Navier-Stokes equations makes an explicit solution impossible for us; however, we can expand the solution in terms of eigenfunctions of the Stokes operator for which we have an explicit form (in terms of Bessel
3 ON VANISHING VISCOSITY IN A DISK 3 functions of the first kind) which we can use to obtain finer estimates on the behavior of the Navier-Stokes equations in the boundary layer than would be possible for a general domain. A word on notation: We use C to represent an unspecified constant that always has the same value on both sides of an equality but may have a different value on each side of an inequality. 2. Definitions and Kato-type conditions We now give definitions of the Euler and Navier-Stokes equations, and state the results from [6], [7], [6], and [7] that we will need. In Section 4 we will specialize to the unit disk, but for now we assume only that Ω is a bounded domain in R 2 with C 2 -boundary Γ, and we let n be the outward normal vector to Γ. A classical solution (u, p) to the Euler equations satisfies, for fixed T >, { t u + u u + p = f and div u = on [, T ] Ω, (E) u n = on [, T ] Γ, and u = u on {} Ω, where div u =. These equations describe the motion of an incompressible fluid of constant density and zero viscosity. We assume that u is in C k+ɛ (Ω), ɛ >, and that f is in C k ([, t] Ω) for all t >, where k = or 2. Then as shown in [8] (Theorem and the remarks on p ), there exists a unique solution u in C loc ([, ); Ck+ɛ (Ω)). The Navier-Stokes equations describe the motion of an incompressible fluid of constant density and positive viscosity. A classical solution to the Navier-Stokes equations can be defined in analogy with (E) by { t u + u u + p = u + f and div u = on [, T ] Ω, u = on [, T ] Γ, and u = u on {} Ω. We will work, however, with weak solutions to the Navier-Stokes equations. Definition 2. (Weak Navier-Stokes Solutions). Given T >, viscosity >, and initial velocity u in H, u in L 2 ([, T ]; V ) with t u in L 2 ([, T ]; V ) is a weak solution to the Navier-Stokes equations if u() = u and (NS) t u v + (u u) v + u v = fv Ω for all v in V. (The spaces H and V are defined in Section 3.) Ω Definition 2.2. We say that the vanishing viscosity limit holds if u u in L ([, T ]; L 2 (Ω)) as. (2.) Theorem 2.3 applies to a bounded domain with C 2 -boundary in R d, d 2. The conditions in Equation (2.2) and Equation (2.4) are due to Kato ([6]), the conditions in Equation (2.3) and Equation (2.7) appear in [7], and the conditions in Equation (2.5) and Equation (2.6) are due to Temam and Wang ([6], [7]). Ω Ω
4 4 JAMES P. KELLIHER Theorem 2.3. Let T > and assume that u is in H and that u is in C k+ɛ (Ω), ɛ > with k = or 2. In addition, assume that (a) u u in L 2 (Ω) as, (b) f is in L ([, T ]; L 2 (Ω)), (c) f f L ([,T ];L 2 (Ω)) as. Let δ : [, ) [, ) be such that δ() converges to while δ()/ diverges to as. Then the the vanishing viscosity limit (Definition 2.2) holds if and only if any of the following conditions holds: T T T T T ω(s) 2 L 2 (Ω) ds as, (2.2) ω(s) 2 L 2 (Γ c) ds as, (2.3) u(s) 2 L 2 (Γ c) ds as, (2.4) τ uτ (s) 2 L 2 (Γ δ() ) ds as, (2.5) τ u n (s) 2 L 2 (Γ δ() ) ds as. (2.6) Here τ represents the derivatives in the boundary layer in the directions tangential to the boundary, uτ is the projection of u in the direction tangential to the boundary, and u n is the projection of u in the direction normal to the boundary. When k = 2, these conditions are also equivalent to T u(s) 2 L 2 (Γ c) ds as. (2.7) The quantity in Equation (2.7) is proportional to the space-time average of the energy in the boundary layer. We show (see Remark (5.3)) that in Equation (2.2), Equation (2.4), and Equation (2.7), contributions from the high frequency modes can be ignored. This result applies to an arbitrary bounded domain in R d, d 2, with a C 2 -boundary. Our main result is Theorem 2.4, which is an improvement of Theorem 2.3 in the special case of the unit disk. In what follows we decompose the solution u in the form u(t, x) = g mj (t)u mj (x), m= j=
5 ON VANISHING VISCOSITY IN A DISK 5 where (u mj ) are the eigenfunctions of the Stokes operator described in Section 4, and let N N u N (t, x) = g mj (t)u mj (x) (2.8) and ũ N (t, x) = m= N j= N m= N j= g mj (t)u mj (x) (2.9) with vorticities ω N (t, x) = ω(u N (t, x)) and ω N (t, x) = ω(ũ N (t, x)). As we will see in Section 4, the frequency of u mk in the tangential direction is m and the radial frequency of u mk is, in effect, k. Thus, u N includes the contributions from all modes with both frequencies less than N, while ũ N includes the contributions from all modes with tangential frequency less than N. Theorem 2.4. Assume that Ω is the unit disk and make the same assumptions on the initial data, forcing, and the function δ as in Theorem 2.3. Let L and M be any functions mapping (, ) to Z + with L(), M() as. (2.) Then the the vanishing viscosity limit (Definition 2.2) holds if and only if any of the following conditions holds: T T T T T T ω(s) M() ω L() (s) 2 L 2 (Ω) ds as, (2.) ω(s) ω L() (s) 2 L 2 (Ω) ds as, (2.2) ω(s) ω L() (s) 2 L 2 (Γ c) ds as, (2.3) u M() (s) u L() (s) 2 L 2 (Γ c) ds as, (2.4) τ uτ (s) τ ũ L(δ) τ (s) 2 L 2 (Γ δ() ) ds as, (2.5) τ u n (s) τ ũ L(δ) n (s) 2 L 2 (Γ δ() ) ds as. (2.6) When k = 2, these conditions are also equivalent to T u M() (s) u L() (s) 2 L 2 (Γ c) ds as. (2.7)
6 6 JAMES P. KELLIHER Observe, for instance, that u M() u L() in Equation (2.7) represents the contribution from all modes whose frequencies in the radial or the tangential direction lie between L() and M(). Remark 2.. By Lemma A.3 and Equation (4.7), u N is essentially the contributions of all the modes with eigenvalues less than CN 2. In fact, suppose that we replace the definition of u N in Equation (2.8) with u N (t, x) = g j (t)u j (x), (2.8) {j:λ j <N 2 } the single subscripts in Equation (2.8) referring to the eigenfunctions and eigenvalues of the Stokes operator on a general domain in R d, d 2, defined in Section 3. It follows easily from Theorem 2.4 that the conditions in Equation (2.), Equation (2.3), Equation (2.4), and Equation (2.7) continue to be equivalent to the vanishing viscosity limit. It is in this form that we would expect Theorem 2.4 to generalize to fairly arbitrary smooth domains in R 2 and with N 2 in Equation (2.8) replaced by N raised to some other power to domains in R d, d 3. The obstacle to establishing this generalization is the difficulty of obtaining the equivalents of Lemma A.8 and Lemma A.9 along with an approximate form of Lemma A. for high frequencies. 3. The Stokes operator in a bounded domain Before specializing to the case of a disk, we discuss first some general properties related to the Stokes operator. We define the function spaces H and V as follows (see Section I..4 of [5] for more details). First let V = { u (D(Ω)) 2 : div u = } be the space of vector-valued divergence-free distributions on Ω. We let H be the closure of V in L 2 (Ω) and V be the closure of V in H (Ω). Alternate characterizations of H and V are H = { u (L 2 (Ω)) 2 : div u = in Ω, u n = on Γ }, V = { u (L 2 (Ω)) 2 : div u = in Ω, u = on Γ }, the boundary conditions applying in terms of a trace. By, we mean the inner product in L 2 (Ω): f, g = Ω fg. (It will be convenient to use complex-valued eigenfunctions, so the complex conjugate is required in this definition. Our velocity fields and vorticities, however, are real, so conjugation will not always appear in our calculations.) Then u, v H = u, v and u, v V = u, v. Although V is dense in H it is not dense in H H (Ω) (with the H -norm). This is because if v is a vector in H H (Ω) that does not vanish on the boundary it cannot be approximated in the H /2 (Γ)-norm by a sequence of vectors in V. By the continuity of the trace operator from H (Ω) to H /2 (Γ),
7 ON VANISHING VISCOSITY IN A DISK 7 then, v cannot be approximated in the H (Ω)-norm by any sequence of vectors in V. Informally, this is because each element of V is zero on Γ and so the limit of a sequence of elements in V cannot become nonzero on the boundary without the gradient near the boundary becoming indefinitely large. More formally, we have Lemma 3.. Lemma 3.. The space V is not dense in H H (Ω). Proof. Let u be any element of V. Then its vorticity ω is in L 2 (Ω) L (Ω) and must satisfy ω = ψ = div ψ = ψ n = u n =, (3.) Ω Ω Ω Γ Γ where ψ is the stream function: u = ψ = ( 2 ψ, ψ) and ω = ψ. Because only u n = for u in H, the same cannot be said for u in H H (Ω): Let u in H H (Ω) have vorticity ω with nonzero total mass. Then for any sequence {v j } in V, ω ω(v j ) L 2 (Ω) C ω ω(v j) L (Ω) C (ω ω(v j )) Ω = C ω >, so V cannot be dense in H H (Ω). Ω We now briefly describe the properties we will need of the Stokes operator A on Ω, referring the reader, for instance, to Section I.2 of [5] for more details. One way to define A is that given u in V H 2 (Ω), Au in H satisfies Au = u + p for some harmonic scalar field p. We have D(A) = V H 2 (Ω) with A mapping D(A) onto H, and there exists a set of eigenfunctions {u j } for A, complete in H and in V, with corresponding eigenvalues {λ j }, < λ λ 2, and each u j is in H 2 (Ω) since we are assuming that Γ is C 2. (When we specialize to the disk, the eigenfunctions will be in C (Ω).) An eigenfunction u j of A satisfies Au j = λ j u j or, equivalently, { uj + λ j u j = p j, p j =, div u j = on Ω, (3.2) u j = on Γ. The eigenfunctions are orthogonal in both H and V. The usual convention is to make the eigenvectors orthonormal in H, but we will find it more convenient to normalize them to be orthonormal in V so that u j 2 L 2 (Ω) = ω j 2 L 2 (Ω) = and u j 2 L 2 (Ω) = u j, u j = λ j u j, Au j = λ j u j, u j = λ j. (3.3) Moreover, we have Lemma 3.2.
8 8 JAMES P. KELLIHER Lemma 3.2. If u is in V with ω = ω(u) then u = ω, ω j u j, (3.4) j= with the sum converging in both V and H. Proof. Let u be in V and let u n = n j= ( u, u j H / u j, u j H )u j. Then u n converges in H to u because {u j } is complete in H. But, u, u j H u j, u j H = λ j u, u j λ j u j, u j = u, Au j u j, Au j = u, u j u j, u j = u, u j = ω, ω j, so the expansion of u in V in terms of the eigenfunctions of A is the same as the expansion of u in H (and the coefficients are as given in Equation (3.4)), meaning that u n converges in V to u as well. In the proof of Lemma 3.2 we used the identity u, v = ω(u), ω(v) for all u, v in V, which follows by integrating by parts. Were we to use the definition of ω(u) as the antisymmetric matrix ( u ( u) T )/2, which is usual in higher dimensions, this would have introduced a factor of 2 into Equation (3.4). Corollary 3.3. If u is in V then u = ω, ω j u j and ω = j= with the sums converging in L 2 (Ω). ω, ω j ω j, Since the solution u to (NS) lies in V for all positive time, we can write ω(t) = g j (t)ω j, u(t) = g j (t)u j, j= j= j= ω(t) 2 L 2 (Ω) = g j (t) 2, u(t) 2 L 2 (Ω) = j= j= g j (t) 2 λ j, (3.5) where g j are functions of time. The expansion of u will converge for all t and that of ω for t > and also for t = if and only if the initial velocity is in V ; in general, we only assume that it is H. Because u(t) u in L 2 (Ω) as t, each g j (t) is continuous at t =, though this does not mean that ω(t) is continuous in L 2 (Ω) at t =. Also, note that g j (t) is complex-valued since the eigenvectors are complex-valued, but u(t) and ω(t) are real-valued.
9 ON VANISHING VISCOSITY IN A DISK 9 4. Eigenfunctions of the Stokes operator in the unit disk We now fix Ω to be the unit disk in R 2 centered at the origin. In [], a complete set of eigenfunctions for the annulus is derived in terms of Bessel functions of the first and second kind, J n and Y n. By ignoring the terms involving Y n and modifying somewhat the calculation of the eigenvalues, one can easily obtain the eigenfunctions for a disk. We will, however, derive the vorticity of the eigenfunctions directly, as this is quite easy. In determining the eigenvalues and the velocity of the eigenfunctions, which is more difficult, we will rely on the results in []. Taking the curl of Equation (3.2) (with u = u j ), we see that the vorticity ω = ω(u) satisfies { ω + λω = on Ω, (4.) u = on Γ. That is, ω is an eigenfunction of the negative Laplacian, but with boundary conditions on the velocity u. Ignoring for the moment the issue of boundary conditions, we use separation of variables to look for a complete set of solutions to ω + λω = on Ω. Writing ω(r, θ) = f n (r)e inθ in polar coordinates, n =,, 2,..., and using gives ( d 2 f n dr 2 = 2 r 2 + r r + 2 r 2 θ 2 + r df n dr + (λ n2 r 2 ) f n ) e inθ =. (4.2) Since J n, the Bessel function of the first kind of order n, is a solution of d 2 J n (s) ds 2 + ) dj n (s) + ( n2 s ds s 2 J n (s) =, (4.3) making the change of variables s = λ /2 r, we see that Equation (4.2) holds with f n (r) = J n (λ /2 r). Thus, the eigenfunctions have vorticity of the form J n (λ /2 r)e inθ and it remains to determine the eigenvalues that satisfy u = on the boundary. The easiest way to do this is to use the expressions in []. For n =, we drop the term involving Y in Equation (39) p. 46 of [], giving u k (r, θ) = λ /2 k J (λ /2 k r)ê θ, where λ /2 k, k =, 2,..., are the eigenvalues described below. For n, dropping the terms involving the Bessel functions of the second kind from
10 JAMES P. KELLIHER the last equation on p. (46) of [], we have ( in u nk (r, θ) = λ nk r J n(λ /2 nk r) + D nkin 2 ) λ 2 r n e inθ ê r nk ( ( ) + J n+ (λ /2 nk r) J n (λ /2 nk r) 2λ /2 nk ) D nkn 2 λ 2 r n e inθ ê θ, nk the eigenvalues λ nk, k =, 2,..., being described below and the D nk being undetermined constants. In both cases we scaled the eigenfunctions differently than in []. A direct calculation shows that ω nk (r, θ) def = ω(u nk )(r, θ) = C nk J n (λ /2 nk r)einθ, where C nk is a normalization constant. A direct calculation also shows that div u nk =. For n =, λ /2 k = j k, where j nk is the k-th positive root of J n+ (x) =, (4.4) as this gives u k (, θ) =. Setting u nk (, θ) = we obtain two equations in the two unknowns λ nk and D nk. We eliminate D nk from the two equations to obtain a single equation for λ nk. Then using the identity in Equation (A.2) we obtain the equation Thus, λ /2 nk λ /2 nk J n(λ /2 nk ) nj n(λ /2 nk ) =. is the k-th positive root of xj n(x) nj n (x) = xj n+ (x) =, (4.5) k =, 2,..., where we used Equation (A.4). That is, λ /2 nk = j n+,k. It follows then that D nk = λ nkj n (λ /2 nk ) = j2 n+,k J n(j n+,k ). (4.6) n n Since we are normalizing the eigenfunctions so that ω mj, ω nk = δ mn δ jk, we must choose C nk so that C 2 nk = J n(j n+,k r)e inθ 2 L 2 (Ω) = 2π rj n (j n+,k r) 2 dr = 2π r2 [ Jn (j n+,k r) 2 J n (j n+,k r)j n+ (j n+,k r) ] 2 = πj n (j n+,k ) 2. Here we used Equation (A.9). Because J n (x) = ( ) n J n (x) while e inθ is linearly independent of e inθ for n =, 2,..., we must also include the eigenfunctions for negative integers n. This gives simple eigenvalues when n = and double eigenvalues corresponding to the pairs (n, k) and ( n, k) for n.
11 ON VANISHING VISCOSITY IN A DISK To summarize, the vorticity of the eigenfunctions is given by with eigenvalue and where ω nk (r, θ) = C nk J n (λ /2 nk r)einθ, C nk = λ nk = j 2 n +,k, (4.7) π /2 J n (j n +,k ), for n Z and k N := {, 2,...}. With our choice of normalization of the eigenfunctions (Equation (3.3)), the velocity becomes where α = j n+,k. u nk (r, θ) = J n(αr) J n (α)r n π /2 α 2 J n (α) r ineinθ ê r + α (J n+(αr) J n (αr)) + 2nJ n (α)r n 2π /2 α 2 e inθ ê θ, J n (α) 5. Proof of Theorem 2.4 From the fundamental energy equality for (NS) we have for all t in [, T ], u 2 L 2 (Ω) = ω 2 L 2 (Ω) 2 u 2 H + 4 f 2 L ([,T ];L 2 (Ω)). (4.8) It follows from Equation (3.5) and assumptions (a) and (b) of Theorem 2.3 that for all sufficiently small >, ω 2 L 2 (Ω) = g mj (s) 2 ds C. (5.) m= j= Remark 5.. In Equation (5.) and most of what follows we use the doubly indexed expressions for the eigenfunctions and eigenvalues of Section 4, though we will occasionally find it convenient to refer to the singly indexed version of Section 3 instead. Theorem 5.. With the assumptions of Theorem 2.4, Proof. Using Lemma A., ω L() 2 L 2 (Γ c) = L() lim ω L() 2 L 2 (Γ c) =. (5.2) L() m= L() j= n= L() k= L() L() g mj (s)g nk (s) ds ω mj, ω nk L 2 (Γ c)
12 2 JAMES P. KELLIHER = = L() L() L() g nj (s)g nk (s) ds ω nj, ω nk L 2 (Γ c) n= L() j= k= L() L() L() g nj (s) g nk (s) ds ω nj L 2 (Γ ω c) nk L 2 (Γ c) n= L() j= k= L() n= L() L() n= L() j= L() g nj (s) ω nj L 2 (Γ c) j= 2 ds L() L() g nj (s) 2 ω nj 2 L 2 (Γ ds, c) where we used the Cauchy-Schwarz inequality in the last step. By Lemma A.3, j= /L() < C/(L() + 2) C/j L()+, = Cλ /2 L(),. Since L() as, for all sufficiently small we have c < C/L() 2πλ /2 L(), 2πλ /2 j, for all j L(). So by Lemma A.8, L() ω nj 2 L 2 (Γ c) CL(). (5.3) j= Then using Equation (5.), ω L() 2 L 2 (Γ CL() c) L() n= L() j= L() g nj (s) 2 ds CL(), which vanishes with by the condition in Equation (2.), and Equation (5.2) therefore holds. Remark 5.2. We could try to improve Theorem 5. by using ω(ũ N ) of Equation (2.9) in place of ω N, thereby incorporating all of the frequencies in the radial direction for a given angular frequency. Unfortunately, the best bound that one can achieve on ω nj 2 L 2 (Γ δ ) for j > n is the extension of Lemma A.8 described in Remark (A.), and this is very much insufficient to bound the terms with j > n. Another possible approach is to try to incorporate the destructive interference that occurs in the inner product of two eigenfunctions in the boundary layer that the use of Hölder s inequality in our proof of Theorem 5. ignored. The best bound one can hope to obtain is that ω nj, ω nk L 2 (Γ δ ) C k j
13 ON VANISHING VISCOSITY IN A DISK 3 for all δ in [, ] and without restriction on n, j, or k except that j k. We could then follow the obvious approach of decomposing the equivalent of the first sum in the proof of Theorem 5. into four pieces: a diagonal term where m = n and j = k and three terms containing low frequencies in j and k, low frequencies in j and high frequencies in k, and high frequencies in both j and k. If we do this, however, we will find that the factor of / k j is just insufficient to obtain convergence. Corollary 5.2. The conditions in Equation (2.) and Equation (2.3) of Theorem 2.4 are equivalent. Proof. That Equation (2.) implies Equation (2.3) follows directly from Theorem 2.3. So assume that Equation (2.3) holds. Because A + B 2 2 A B 2 for any norm, ω 2 L 2 (Γ 2 c) ω L() 2 L 2 (Γ + 2 c) ω ω L() 2 L 2 (Γ. c) This vanishes with by Theorem 5. and Equation (2.3), showing that Equation (2.3) holds and hence by Theorem 2.3 that Equation (2.) holds. Theorem 5.3. With the assumptions of Theorem 2.4, and lim lim u(s) u M() (s) 2 L 2 (Γ c) ds = (5.4) u L() (s) 2 L 2 (Γ c) ds =. (5.5) Proof. We can write u(t) u M() (t) = A(t) + B(t), where A(t) = and Now, M() m= M() j=m()+ g mj (t)u mj (x), B(t) = m >M() j= g mj (t)u mj (x) u(t) u M() (t) 2 L 2 (Γ c) 2 A(t) 2 L 2 (Γ c) + 2 B(t) 2 L 2 (Γ c). A(t) 2 L 2 (Γ c) A(t) 2 L 2 (Ω) = = M() m= M() j=m()+ M() m= M() j=m()+ g mj (t) 2 λ mj λ M() ω(t) ω M() (t) 2 L 2 (Ω), g mj (t) 2 u mj 2 L 2 (Ω) M() λ M() m= M() j=m()+ g mj (t) 2
14 4 JAMES P. KELLIHER where we used Equation (3.3). Similarly, B(t) 2 L 2 (Γ c) B(t) 2 L 2 (Ω) λ M() ω(t) ω M() (t) 2 L 2 (Ω). By Equation (4.7) and Lemma A.3, λ M() and λ M() are both bounded below (and above) by CM() 2, so u(t) u M() (t) 2 L 2 (Γ C c) M() 2 ω(t) ωm() (t) 2 L 2 (Ω). Then, = u(s) u M() (s) 2 L 2 (Γ c) ds C M() 2 C M() 2 C 2 M() 2 ω(s) ω M() (s) 2 L 2 (Ω) ds ω(s) 2 L 2 (Ω) ds ω(s) 2 L 2 (Ω) ds C 2 M() 2, where in the last inequality we used Equation (5.). This vanishes with by the assumption on M in Equation (2.) giving Equation (5.4). Arguing as in the proof of Theorem 5., u L (s) 2 L 2 (Γ c) ds CL()3 L() L() n= L() j= L() L() n= L() j= g nj (s) 2 ds CL() L() g nj (s) 2 u nj 2 L 2 (Γ ds c) for all sufficiently small. In the second inequality we used Lemma A.9 and in the last inequality we used Equation (5.). This integral also vanishes with by the assumption on L in Equation (2.) giving Equation (5.5). Corollary 5.4. The conditions in Equation (2.), Equation (2.), Equation (2.4), and Equation (2.7) of Theorem 2.4 are equivalent. j= Proof. For sufficiently large, L() M(), and we have u(s) 2 L 2 (Γ c) 3 um() (s) u L() (s) 2 L 2 (Γ c) + 3 u L() (s) 2 L 2 (Γ c) + 3 u(s) um() (s) 2 L 2 (Γ c). It follows from Theorem 5.3 that t lim sup u(s) 2 L 2 (Γ c) 3 lim sup u M() u L() 2 L 2 (Γ c). In particular, the first limsup is zero if and only if the second limsup is zero (the reverse inequality without the factor of 3 being trivial). Then
15 ON VANISHING VISCOSITY IN A DISK 5 Equation (2.7) of Theorem 2.3 shows that Equation (2.7) holds if and only if Equation (2.) holds. The sufficiency of Equation (2.) and Equation (2.4) for Equation (2.) to hold then follows from Poincaré s inequality in the form u M() (s) u L() (s) 2 L 2 (Γ c) C2 u M() (s) u L() (s) 2 L 2 (Γ c) C 2 u M() (s) u L() (s) 2 L 2 (Ω) = C 2 ω M() (s) ω L() (s) 2 L 2 (Ω). The necessity of Equation (2.) and Equation (2.4) follow immediately from Theorem 2.3. Remark 5.3. If we replace the definition of u N in Equation (2.8) with that in Equation (2.8), then it is clear that Equation (5.4) continues to hold in any bounded domain in R 2 with a C 2 -boundary. It follows as in Corollary 5.4 that the vanishing viscosity limit of Definition 2.2 holds if and only if the condition in Equation (2.), Equation (2.4), or (when k = 2) Equation (2.7) holds with the term involving u L() in each of these conditions removed. A similar result would hold in any dimension for an arbitrary bounded domain with a C 2 -boundary. Theorem 5.5. With the assumptions of Theorem 2.4, and T T τ ũ L(δ) τ (s) 2 L 2 (Γ δ() ) ds as (5.6) τ ũ L(δ) n (s) 2 L 2 (Γ δ() ) ds as. (5.7) Proof. In the unit disk, uτ = u θ and τ = σ, where σ is arc length along the circle of radius r, in which r is held constant. Thus, τ uτ = uθ σ = u θ r θ and for any positive integer N it follows from Poincaré s inequality that τ ũ N τ (s) 2 L 2 (Γ δ ) = r θ (ũn (s)) θ 2 L 2 (Γ δ ) ( δ) 2 θ (ũn (s)) θ 2 L 2 (Γ δ ) Cδ2 2 2 ( δ) 2 r θ (ũn (s)) θ Cδ2 2 2 ( δ) 2 r θ (ũn (s)) θ. But, 2 r θ (ũn (s)) θ = = i N L 2 (Γ δ ) m= N j= N m= N m g mj (s) 2 r θ u mj(r, θ) g mj (s) r u mj(r, θ), j= L 2 (Ω)
16 6 JAMES P. KELLIHER the last equality following from the simple dependence of u mj on θ in Equation (4.8). Thus, r θ (ũn (s)) θ N L 2 (Ω) i m g mj (s) u mj (r, θ) = N m= N m= N j= N m 2 g mj (s) 2 N 2 j= m= N j= L 2 (Ω) g mj (s) 2 N 2 u 2 L 2 (Ω), where we used the orthonormality of the eigenfunctions in V. Combining these two inequalities gives Then using Equation (5.), T τ ũ N τ (s) 2 L 2 (Γ δ ) CN 2 δ 2 τ ũ L(δ) τ (s) 2 L 2 (Γ δ() ) CL(δ)2 δ 2 ( δ) 2 ( δ) 2 u 2 L 2 (Ω). T u 2 L 2 (Ω) CL(δ)2 δ 2 ( δ) 2. This vanishes with δ by the assumption on L in Equation (2.) and hence vanishes with since δ vanishes with, giving Equation (5.6). The proof of Equation (5.7) is entirely analogous. The technique used in the proof of Theorem 5.5 comes from the key inequality following Equation (3.2) in [3]. Corollary 5.6. The conditions in Equation (2.), Equation (2.2), Equation (2.5), and Equation (2.6) of Theorem 2.4 are equivalent. Proof. This corollary can be proved much along the lines of the proofs of Corollary 5.2 and Corollary 5.4. (It is here that we use the assumption that δ()/ diverges to as, which is needed in applying Theorem 2.3.) Together, Corollary 5.2, Corollary 5.4, and Corollary 5.6 establish Theorem Radially symmetric initial vorticity It follows from Lemma 3. that if an initial velocity, no matter how smooth, lies in H but not in V and has a vorticity whose total mass is nonzero, then the velocity of the corresponding solution to (N S) will be discontinuous in H (Ω) at time zero. In the same way, the vanishing viscosity limit of the vorticity cannot hold in L ([, T ]; L 2 (Ω)) for such an initial velocity, since the total mass of the vorticity for the solution to (E) is conserved over time. This means that we might expect a different character to the vanishing viscosity limit when the initial vorticity has zero total mass versus when it has nonzero total mass.
17 ON VANISHING VISCOSITY IN A DISK 7 Indeed, this is what happens in the special case of radially symmetric initial vorticity where, when the initial velocity is in V (which is equivalent for radially symmetric vorticity to the total mass of the vorticity being zero) we obtain convergence in L ([, T ]; L 2 (Ω)) of both the velocity and the vorticity (see []), whereas for initial velocity in H we obtain only the convergence of the velocity in this space, as in Theorem 6.. Such convergence follows immediately from the conditions in Equation (2.2), Equation (2.5), or Equation (2.6) of Theorem 2.4. The convergence also follows from the sufficiency of the conditions in Equation (2.5) and Equation (2.6) as established in [6], since both conditions are satisfied (the gradients in the tangential direction being zero) as pointed out in [7]. When the forcing is zero, however, there is a simple proof that uses only Kato s original conditions. Theorem 6.. Assume that u and u are as in Theorem 2.3 with (for simplicity) u = u, that f = f =, and that ω = ω(u ) is radially symmetric. Then the vanishing viscosity limit of Equation (2.) holds. Proof. Because ω is radially symmetric, ω remains radially symmetric for all time, so ω(u u) = u ω =. Then because Ω is simply connected, u u = q for some scalar field q, and the nonlinear term in (NS) disappears. Thus, (NS) reduces to u () = u and t u v + u v = (6.) Ω for all v in V. This is the heat equation in weak form, which is invariant under the transformation (, t, x) (, t, x). That is, if u is a solution to Equation (6.) with =, then u (t, x) = u (t, x) is a solution to Equation (6.) because u () = u () = u and Ω t u (t, x) v(x) dx + u (t, x) v(x) dx Ω [ ] = ( t u )(t, x) v(x) dx + u (t, x) v(x) dx =. Ω Ω It follows that t ω(s) 2 L 2 (Ω) ds = ω (s) 2 L 2 (Ω) ds = ω (τ) 2 L 2 (Ω) dτ. This vanishes as by the continuity of the integral, because u is in L 2 ([, T ]; V ). The limit in Equation (2.) then follows from the condition in Equation (2.2) of Theorem 2.3. The proof of of Theorem 6. does not yield a bound on the rate of convergence in Equation (2.). Also, without assuming that the initial vorticity is radially symmetric, the argument in the proof of Theorem 6. can be applied to solutions to the Stokes problem (the linearized Navier-Stokes equations) to show that they converge in the vanishing viscosity limit to a solution to the linearized Euler equations (which is just the steady state solution Ω
18 8 JAMES P. KELLIHER u = u ). This would be more interesting, though, if time-varying Dirichlet boundary conditions, for instance, could be incorporated, as in [3]. In the simpler case of u also lying in V, the solution to (E), which is steady state, is zero on the boundary. This eliminates the troublesome boundary term that appears in the direct energy argument bounding u(t) u(t) L 2 (Ω), giving an extremely simple proof of Equation (2.). We give, however, a longer proof of convergence using Theorem 2.3, because it suggests how we might treat more general initial velocities. For convenience, we assume zero forcing and more regularity on the initial velocity than is strictly necessary. Theorem 6.2. Assume that ω is radially symmetric, u = u is in H 3 (Ω) V, and there is no forcing. Then the vanishing viscosity limit of Equation (2.) holds. Proof. Because u is in V H 3 (Ω), u is in L loc ([, ); V H3 (Ω)) (see, for instance, Theorem III.3.6, Remark III.3.7, and Theorem III.3. of [5]). As observed in the proof of Theorem 6., u u = q for some scalar field q. Then t u + (p + q) = u and taking the divergence of both sides we conclude that (p + q) =. Because of the radial symmetry, however, p + q is constant on Γ and hence is constant on Ω. Thus, (p + q) = and t u = u. But u = on Γ so t u = on Γ, and it follows that u = t u is in V. Because u is in V at time zero, we can use the expansion for ω in Equation (3.5) for all time, including for time zero. Because ω remains radially symmetric over time, the expansion reduces to ω(t, r, θ) = g k (t)ω k (r) = g k (t)c k J (j k r). (6.2) k= Since u is in V, it follows from Corollary 3.3 that ω has an expansion like that of Equation (6.2): ω(t, r, θ) = h k (t)c k J (j k r). (6.3) k= (The analogous expansion of ω including all the eigenvectors {ω nk } fails to converge at t = for non-radially symmetric solutions because u is not, in general, in V.) Since J (j k r) = λ k J (j k r), it follows that h k (t) = λ k (t)g k (t) and then by Equation (6.) in strong vorticity form, k= t ω = ω, (6.4) that g k (t) = λ kg k (t). Thus, g k (t) = g k ()e λkt and ω(t, r, θ) = C k g k ()e λkt J (j k r). (6.5) k=
19 But then ON VANISHING VISCOSITY IN A DISK 9 ω(t) 2 L 2 (Ω) = Ck 2 g k() 2 e 2λkt J (j k r) 2 L 2 (Ω) k= g k () 2 C k J (j k r) 2 L 2 (Ω) = ω 2 L 2 (Ω) k= (6.6) and convergence in the vanishing viscosity limit follows from the condition in Equation (2.2) of Theorem 2.3. (Examining the argument in [6] shows that the convergence rate in Equation (2.) is bounded by C(t) /2.) We could have concluded the proof another, more indirect way, as follows. From Equation (6.5) we see that ω(t) is continuous in the L 2 (Ω)-norm at time zero. Thus, we can multiply Equation (6.4) by ω and integrate over time to give d 2 dt ω(t) 2 L 2 (Ω) + ω(t) 2 L 2 (Ω) = ( ω n)ω. (6.7) But ω = ( u) = on Γ because u is in V as we showed, and we conclude by integrating over time that ω is in L ([, ); L 2 (Ω)), and Equation (2.) follows from the condition in Equation (2.2) of Theorem 2.3. This argument shows that whenever the vorticity is continuous at time zero in the L 2 (Ω)-norm and u(t) lies in V for all t >, Equation (2.) holds. The latter condition, however, is very special. In fact, for u with regularity as in Theorem 6.2, t u will always be in V, so u = ( t u + u u + p)/ is in V if and only if u u + p is in V. This, in turn, can hold only if p = u u on Γ. This suggests that for initial velocities for which ω(t) is continuous in the L 2 (Ω)-norm at time zero we might attempt to make an argument using Equation (6.7), though now without the right-hand side vanishing. This would require either control on ω, as we had above, or on ω itself. (For solutions to (NS) with the boundary condition ω = u n = rather than u = this term would vanish and one can obtain the vanishing viscosity limit easily though not using Theorem 2.3, which one would need to show applies to such solutions as done, for instance, in [] and [2].) Another line of attack is also suggested by the proof of Theorem 6.2 to gain control on the coefficients g(t) either by assumptions on the initial velocity or on the solution. Unless bounds as remarkably strong as those obtained in the proof of Theorem 6.2 are achieved, though, something more sophisticated must be employed to obtain Equation (2.). 7. Interpretation in terms of length scales In [3], Cheng and Wang consider the vanishing viscosity limit in the setting of a two-dimensional rectangular channel R, periodic in the x direction with period L and with characteristic boundary conditions (which include no-slip Γ
20 2 JAMES P. KELLIHER boundary conditions as a special case). They decompose any vector u on R of sufficient regularity as u = j= e2πijx/l u j and define the projection P k u = k j= e2πijx/l u j onto the space spanned by the first k modes. This in effect allows one to isolate successively finer-scale spatial variations in the direction tangential to the boundary. They then construct an approximation sequence {v L } to u by letting v L be the solution to the equation that results after projecting each term in (NS) using P N. (We have changed their notation somewhat.) Their v L is the approximate-solution analog of the exact solution truncation represented by ũ L in Equation (2.9). The main result in [3] is that v L() converges to u in L ([, T ]; L 2 (Ω)) as. The requirement on L() is the same as our condition on L in Equation (2.) (with the additional condition that L() as as one would expect), so convergence of v L() to u occurs when only tangential length scales of order larger than are included in the approximations. (All length scales in the normal direction, however, are included. See Remark (5.2) concerning this issue in regards to the vorticity.) The result in [3] makes an important observation about the difficulty of determining numerically whether or not the vanishing viscosity limit holds. Our method of decomposing the solution using the eigenfunctions of the Stokes operator, on the other hand, says little about computation, since approximating this decomposition numerically is probably as least as hard as approximating the solution itself. Nonetheless, it more directly characterizes the properties of the solution itself at different length scales. The analog to the result in [3] is Theorem 5.5, which shows that Temam and Wang s conditions in Equation (2.5) and Equation (2.6), when applied only to the modes with tangential wavelengths of CL() or higher, holds as long as the condition on N in Equation (2.) hold. This does not, however, imply that u L() converges to u in the vanishing viscosity limit, only that if the vanishing viscosity limit fails to hold, the failure originates in the behavior of the tangential component of the gradient projected into the space spanned by the modes with tangential frequencies of order L() or higher; that is, at length scales of order or lower. The other conditions in Theorem 2.4 give alternative ways to measure the behavior of the solution at different length scales or frequencies. They show that we cannot simply say that if the vanishing viscosity fails to hold then the failure lies in the behavior of the solution at any particular range of length scales, but rather that the pertinent range of length scales varies with the measure of behavior. Whether any of these conditions brings us any closer to proving that the vanishing viscosity limit holds in general for smooth initial data in a bounded domain or to proving that it fails to hold in at least one instance remains completely unclear.
21 ON VANISHING VISCOSITY IN A DISK 2 Appendix A. Bounds on the Eigenfunctions In Lemma A. we state the basic identities involving the Bessel functions that we use. We then give a series of lemmas that lead to the bounds on the velocity and vorticity of the eigenfunctions in the boundary layer that we used in the proof of Theorem 2.4. It is perhaps important to note that in the proofs that follow we avoid the use of asymptotic formulas for the Bessel functions, even when such formulas might appear to be useful. This is because we need to deal with the relative values of Bessel functions of different orders near a zero of one of the Bessel functions, and it is precisely in these situations that the errors in the asymptotic formulas dominate. Also, most of the following lemmas apply without change to their proofs with n being any nonnegative real value. Lemma A.. For all nonnegative real numbers n and x, 2nJ n (x) xj n (x) = xj n+ (x), (A.) 2J n(x) = J n (x) J n+ (x), J n (x) = n x J n(x) + J n(x), (A.2) (A.3) (β 2 α 2 ) J n+ (x) = n x J n(x) J n(x), x n J n (αx) α = J n (αx)x n = α x n J n (αx) dx, J n+ (αx)x n dx, (A.4) (A.5) (A.6) xj n (αx)j n (βx) dx = x [ αj n(αx)j n (βx) βj n(βx)j n (αx) ], xj n (ax) 2 dx = 2 xj n (ax) 2 dx = x2 2 (A.7) [ ) ] x 2 J n(ax) 2 + (x 2 n2 J n (ax) 2, (A.8) a 2 [ Jn (ax) 2 J n (ax)j n+ (ax) ]. (A.9) Proof. These are standard identities for Bessel functions. For instance, see Equations (6.28), (6.29), (6.3), (6.3), (6.38), (6.39), (6.5), (6.52), and (6.53) of [2].
22 22 JAMES P. KELLIHER Lemma A.2. For all nonnegative integers n and all positive integers k, < j n+,k j nk < π 2, where j nk is defined in Equation (4.4). Proof. Let j k be the k-th positive zero of J, where we now allow to be a real number in the interval [, ). It is shown in [4] and [5] that for all k, j k is strictly concave as a function of and that dj k /d > (see also [9]). Thus, the function n j n+,k j nk is strictly decreasing as a function of n. But by Equation (2.9) of [5], j n+,k j nk as n, so j n+,k j nk >. The positive zeros of J lie in the intervals (mπ π, mπ π), m =,,..., and the positive zeros of J lie in the intervals (m π + 8 π, m π + 4 π), m =, 2,.... That the zeros lie in only these intervals is shown in Section 5.32 p. 489 and Section 5.34 p. 49 of [8] using an approach of Schafheitlin s. That each of these intervals contains at least one zero is shown on p. 4 of [2]. But j k j k > as we showed above so each interval contains precisely one zero. Because the zeros of J and J are interleaved (see p. 6 of [2], for instance) we can then conclude that j k j k < π 2. But as we observed above, the function n j n+,k j nk is strictly decreasing as a function of n, so j n+,k j nk < π 2 holds for all n. Lemma A.3. For all n =,,... and k =, 2,..., n + k < j nk < π(n/2 + k) π(n + k). Proof. By Lemma A.2, for all n and j, n j nk = j k + (j mk j m,k ) j k + n > n + k, m= because j k > k (which follows directly from Equation (4.3); see p of [8], for instance). By an observation in the proof of Lemma A.2 it follows that j k < πk, and a similar argument using the inequality j n+,k j nk < π 2 from Lemma A.2 gives the upper bound on j nk. Lemma A.4. Let α = j n+,k and β = j nk. For n =,, 2,... and k =, 2,..., J n (αx) J n (α) if β α < x <. Proof. Let g(x) = J n (αx)/ J n (α). From Equation (A.4), J n(α) = (n/α)j n (α), so J n(α) has the same sign as J n (α). From this we conclude that g is increasing in a left-neighborhood N of. Between each zero of J n there is exactly one zero of J n+ (see p. 6 of [2], for instance). Between each zero of J n there is also exactly one zero of J n, because the maximum values of J n are all positive and the minimum values are all negative (see, for instance, p. 7 of [2]) and J n has no repeated
23 ON VANISHING VISCOSITY IN A DISK 23 positive roots (this follows from the defining equation Equation (4.3)). Thus, the neighborhood N includes all x such that β < αx < α. Since g() = it follows that g(x) for all such x. Lemma A.5. Let α = j n+,k and β = j nk. There exists a constant C such that for all n =,,... and k =, 2... n, J n+ (αx) J n (α) Cn( x) if β α < x <. Proof. Since J n+ (α) =, Equation (A.5) with n + in place of n gives J n+ (αx) = α x n+ x t n+ J n (αt) dt. As long as β < αx < α, J n (αt) does not change sign on the interval (x, ] and has its maximum value on this interval at, as observed in the proof of Lemma A.4. Thus, J n+ (αx) α x n+ J n(α) But by Lemma A.2 and Lemma A.3, so β α = α β α (β/α) (n+) x t n+ dt α J n(α) ( x). (β/α) n+ π/2 n + 2 = β α π 2n + 4 ( π ) (n+) e π/2, 2n + 4 the last inequality following from elementary calculus. We conclude that which completes the proof. J n+ (αx) Cn J n (α) ( x), (A.) Lemma A.6. Let α = j n+,k and β = j nk. There exists a constant C such that for all n =,,... and k =, 2..., n, J n (αx) J n (α) C if β α < x <. Proof. Because the positive zeros of J n are interlaced with those of J n, J n does not change sign on the interval [β, α]. From Equation (A.4) with n in place of n, J n (β) = ((n )/β)j n (β), so J n (β) has the same sign as J n (β), and we conclude that J n reaches it maximum value on the interval [β, α] at β. Therefore, for β/α < x <, J n (αx) J n (α) J n (β) J n (α). But, by Equation (A.), J n+ (β) = 2(n/β)J n (β) J n (β) = J n (β), so J n (αx) J n (α) J n+ (β) J n (α) Cn( β/α) C,
24 24 JAMES P. KELLIHER where we used Lemma A.5 and Lemma A.3. Lemma A.7. For all j in N, ω j 2 L 2 (Γ δ ) 2δ when δ λ /2 j. Proof. Let ω j = ω nk and α = j n+,k = λ /2 j, where without loss of generality we assume that n. Then ω j 2 L 2 (Γ δ ) = 2πC2 nk rj n (αr) 2 dr = 2 r J n(αr) 2 δ δ J n (α) 2 dr. In the integrals above, with β = j nk, β/α = (β α)/α /α = λ /2 j δ < r <, where we used Lemma A.2, and the lemma follows from Lemma A.4. Employing Lemma A.7, we can extend its range of applicability, though with a higher bound on the width of the boundary layer. Lemma A.8. For all n in Z, k =,..., n, and all δ < 2πλ /2 n, ω nk 2 L 2 (Γ δ ) 2δ. Proof. Without loss of generality we assume that n. It follows from Lemma A.3 that j n,n /j n, π(n + n)/(n + ) 2π; the lemma follows from this inequality and Lemma A.7. Remark A.. It is possible to extend Lemma A.8 to include all values of k. The idea of the proof is that for k > n, ω nk passes through k complete half-periods (annuli in the unit disk lying between successive nonnegative zeroes of J n (j n+,k r)) and ends with a partial period. Since J n (x) decays like x /2 and the spacing between consecutive zeros of J n approaches a constant, the L 2 -norms of ω nk on each of those half-periods converges to a constant, and since the L 2 -norm of ω nk on the entire unit disk is, the square of the L 2 -norm of ω nk on the last half-period is less than C/k (with C near ). But the last half-period has a width that is greater than C/k. Extending this argument to m periods, what we have shown is that ω nk 2 L 2 (Γ Cm/k ) m/k. With the assumed bound on δ, we choose m so that m/k is of the same order as δ, and the proof is essentially complete. Lemma A.9. There exist positive constants C and C 2 with C 2 < such that for all n in Z and all k =,..., n, when δ < C 2 λ /2 n. u nk 2 L 2 (Γ δ ) C δ 3
25 ON VANISHING VISCOSITY IN A DISK 25 Proof. In the proof that follows, we will often use Lemma A.3 without explicit mention. Also, without loss of generality we assume that n. Let α = j n+,k. We bound first the radial component of u nk. We have, where g n (r) = J n(αr)r n J n (α) and J n (αr) J n (α)r n J n (α)r = α J n (α) r = r n g n (r), B nk (x) = J n+(αx). J n (α) J n+ (αx) B nk (x) x n dx = α r x n dx, To verify the second equality we use the identity in Equation (A.6), from which it follows that the second expression for g n is J n (αr)r n + C J n (α) for some constant C. But all three expressions for g n are zero at r =, so we have the correct limits of integration in the second expression. It follows from Lemma A.5 and our third expression for g n that g n (r) Cnα for all r /α. From Equation (4.8), r x n ( x) dx Cn2 ( r)2 rn u r g(r)rn nk (r, θ) = α 2 π /2 ineinθ ê r, so when δ /α we have u r nk 2 L 2 (Γ δ ) = 2π δ In the last inequality we used so ( δ) 2n r u r nk 2 dr Cn6 α 4 Cn 2 ( δ) 2n δ 5 Cδ 3. δ δ < C 2 λ /2 n = C 2 j n+,k C 2 n + ( r) 4 r 2n ( C ) 2n ( 2 = (G(C 2, n + )) 2 C ) 3 2 n + n + ( C 2 ) 2 ( C 2 ) 3 = ( C 2 ) = C >, where G is the function of Lemma A.. dr
26 26 JAMES P. KELLIHER For the angular component of u nk, we write αj n+ (αr) αj n (αr) + 2nJ n (α)r n = [2nJ n (αr) αrj n (αr)] + αrj n (αr) 2nJ n (αr) + αj n+ (αr) αj n (αr) + 2nJ n (α)r n = α(r + )J n+ (αr) + 2n [ J n (α)r n J n (αr) ] + αj n (αr)(r ). From Equation (4.8) we then have u θ nk 2 C α2 (r + ) 2 J n+ (αr) 2 α 4 J n (α) 2 + C n2 α 4 + C α2 J n (αr) 2 ( r) 2 α 4 J n (α) 2 ( Jn (α)r n J n (αr) J n (α) ) 2 (A.) C( r) 2 + C n 2 ( r n 2 g n (r) ) 2, where we applied both Lemma A.5 and Lemma A.6. The first term in Equation (A.) contributes no more than C δ r( r) 2 dr C δ ( r) 2 dr Cδ 3, and the same is true of the second term in Equation (A.) arguing as for, and this completes the proof. u r nk Lemma A.. When m n, u mj, u nk L 2 (Γ c) = ω mj, ω nk L 2 (Γ c) =. Proof. We have, 2π ω mj, ω nk L 2 (Γ = c) rf(r) e i(m n)θ dθ dr, c where f(r) is a product of two Bessel functions. When m n, the inner integral is zero. A similar argument gives u mj, u nk L 2 (Γ =. c) Lemma A.. Let α be in (, ) and define G α : [, ) [, ) by ( G α (x) = α ) x. x Then for all x >, Proof. The proof is elementary. α G α (x) < e α.
27 ON VANISHING VISCOSITY IN A DISK 27 References [] Jerry L. Bona and Jiahong Wu. The zero-viscosity limit of the 2D Navier-Stokes equations. Stud. Appl. Math., 9(4): , [2] Frank Bowman. Introduction to Bessel functions. Dover Publications Inc., New York, , 22 [3] Wenfang Cheng and Xiaoming Wang. A discrete Kato type theorem on inviscid limit of Navier-Stokes flow. Journal of Mathematical Physics, 48(6), 27. 2, 6, 9, 2 [4] Á. Elbert. Concavity of the zeros of Bessel functions. Studia Sci. Math. Hungar., 2(-2):8 88 (98), [5] Árpád Elbert and Andrea Laforgia. On the square of the zeros of Bessel functions. SIAM J. Math. Anal., 5():26 22, [6] Tosio Kato. Remarks on zero viscosity limit for nonstationary Navier-Stokes flows with boundary. In Seminar on nonlinear partial differential equations (Berkeley, Calif., 983), volume 2 of Math. Sci. Res. Inst. Publ., pages Springer, New York, 984., 3, 9 [7] James P. Kelliher. On Kato s conditions for vanishing viscosity. To appear in Indiana University Mathematics Journal., 2, 3 [8] Herbert Koch. Transport and instability for perfect fluids. Math. Ann., 323(3):49 523, [9] Andrea Laforgia and Martin E. Muldoon. Monotonicity and concavity properties of zeros of Bessel functions. J. Math. Anal. Appl., 98(2):47 477, [] D.-S. Lee and B. Rummler. The eigenfunctions of the Stokes operator in special domains. III. ZAMM Z. Angew. Math. Mech., 82(6):399 47, 22. 9, [] J.-L. Lions. Quelques méthodes de résolution des problèmes aux limites non linéaires. Dunod, [2] Pierre-Louis Lions. Mathematical topics in fluid mechanics. Vol., volume 3 of Oxford Lecture Series in Mathematics and its Applications. The Clarendon Press Oxford University Press, New York, [3] Maria Carmela Lombardo, Russel E. Caflisch, and Marco Sammartino. Asymptotic analysis of the linearized Navier-Stokes equation on an exterior circular domain: explicit solution and the zero viscosity limit. Comm. Partial Differential Equations, 26(-2): , 2. 2, 8 [4] Marco Sammartino and Russel E. Caflisch. Zero viscosity limit for analytic solutions of the Navier-Stokes equation on a half-space. II. Construction of the Navier-Stokes solution. Comm. Math. Phys., 92(2):463 49, 998. [5] Roger Temam. Navier-Stokes equations. AMS Chelsea Publishing, Providence, RI, 2. Theory and numerical analysis, Reprint of the 984 edition. 6, 7, 8 [6] Roger Temam and Xiaoming Wang. On the behavior of the solutions of the Navier- Stokes equations at vanishing viscosity. Ann. Scuola Norm. Sup. Pisa Cl. Sci. (4), 25(3-4): (998), 997. Dedicated to Ennio De Giorgi. 2, 3, 7 [7] Xiaoming Wang. A Kato type theorem on zero viscosity limit of Navier-Stokes flows. Indiana Univ. Math. J., 5(Special Issue):223 24, 2. Dedicated to Professors Ciprian Foias and Roger Temam (Bloomington, IN, 2). 2, 3, 7 [8] G. N. Watson. A treatise on the theory of Bessel functions. Cambridge Mathematical Library. Cambridge University Press, Cambridge, 995. Reprint of the second (944) edition. 22 Department of Mathematics, Brown University, Box 97, Providence, RI 292 Current address: Department of Mathematics, Brown University, Box 97, Providence, RI 292 address: kelliher@math.brown.edu
CONNECTIONS BETWEEN A CONJECTURE OF SCHIFFER S AND INCOMPRESSIBLE FLUID MECHANICS
CONNECTIONS BETWEEN A CONJECTURE OF SCHIFFER S AND INCOMPRESSIBLE FLUID MECHANICS JAMES P. KELLIHER Abstract. We demonstrate connections that exists between a conjecture of Schiffer s (which is equivalent
More informationRecent progress on the vanishing viscosity limit in the presence of a boundary
Recent progress on the vanishing viscosity limit in the presence of a boundary JIM KELLIHER 1 1 University of California Riverside The 7th Pacific RIM Conference on Mathematics Seoul, South Korea 1 July
More informationVANISHING VISCOSITY IN THE PLANE FOR VORTICITY IN BORDERLINE SPACES OF BESOV TYPE
VANISHING VISCOSITY IN THE PLANE FOR VORTICITY IN BORDERLINE SPACES OF BESOV TYPE ELAINE COZZI AND JAMES P. KELLIHER Abstract. The existence and uniqueness of solutions to the Euler equations for initial
More informationASYMPTOTIC 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 informationLACK OF HÖLDER REGULARITY OF THE FLOW FOR 2D EULER EQUATIONS WITH UNBOUNDED VORTICITY. 1. Introduction
LACK OF HÖLDER REGULARITY OF THE FLOW FOR 2D EULER EQUATIONS WITH UNBOUNDED VORTICITY JAMES P. KELLIHER Abstract. We construct a class of examples of initial vorticities for which the solution to the Euler
More informationON A FAILURE TO EXTEND YUDOVICH S UNIQUENESS THEOREM FOR 2D EULER EQUATIONS
ON A FAILURE TO EXTEND YUDOVICH S UNIQUENESS THEOREM FOR 2D EULER EQUATIONS JAMES P. KELLIHER Abstract. In 995, Yudovich extended his own 963 uniqueness result for solutions to the 2D Euler equations with
More informationTHE INVISCID LIMIT FOR TWO-DIMENSIONAL INCOMPRESSIBLE FLUIDS WITH UNBOUNDED VORTICITY. James P. Kelliher
Mathematical Research Letters 11, 519 528 (24) THE INVISCID LIMIT FOR TWO-DIMENSIONAL INCOMPRESSIBLE FLUIDS WITH UNBOUNDED VORTICITY James P. Kelliher Abstract. In [C2], Chemin shows that solutions of
More informationOn the zero viscosity limit in incompressible fluids
On the zero viscosity limit in incompressible fluids A L Mazzucato Department of Mathematics, Penn State University, University Park, PA 1682, U.S.A. E-mail: alm24@psu.edu Abstract. We discuss the convergence
More informationA RECURRENCE THEOREM ON THE SOLUTIONS TO THE 2D EULER EQUATION
ASIAN J. MATH. c 2009 International Press Vol. 13, No. 1, pp. 001 006, March 2009 001 A RECURRENCE THEOREM ON THE SOLUTIONS TO THE 2D EULER EQUATION Y. CHARLES LI Abstract. In this article, I will prove
More informationWeek 6 Notes, Math 865, Tanveer
Week 6 Notes, Math 865, Tanveer. Energy Methods for Euler and Navier-Stokes Equation We will consider this week basic energy estimates. These are estimates on the L 2 spatial norms of the solution u(x,
More informationRemarks on the blow-up criterion of the 3D Euler equations
Remarks on the blow-up criterion of the 3D Euler equations Dongho Chae Department of Mathematics Sungkyunkwan University Suwon 44-746, Korea e-mail : chae@skku.edu Abstract In this note we prove that the
More informationREMARKS 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 informationBOUNDARY LAYER ASSOCIATED WITH A CLASS OF 3D NONLINEAR PLANE PARALLEL CHANNEL FLOWS DONGJUAN NIU 2 3 XIAOMING WANG 4
BOUNDARY LAYER ASSOCIATED WITH A CLASS OF 3D NONLINEAR PLANE PARALLEL CHANNEL FLOWS ANNA L MAZZUCATO 1 Department of Mathematics, Penn State University McAllister Building University Park, PA 16802, U.S.A.
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 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 informationThe Navier-Stokes Equations with Time Delay. Werner Varnhorn. Faculty of Mathematics University of Kassel, Germany
The Navier-Stokes Equations with Time Delay Werner Varnhorn Faculty of Mathematics University of Kassel, Germany AMS: 35 (A 35, D 5, K 55, Q 1), 65 M 1, 76 D 5 Abstract In the present paper we use a time
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 informationREMARKS ON THE MEMBRANE AND BUCKLING EIGENVALUES FOR PLANAR DOMAINS. Leonid Friedlander
REMARKS ON THE MEMBRANE AND BUCKLING EIGENVALUES FOR PLANAR DOMAINS Leonid Friedlander Abstract. I present a counter-example to the conjecture that the first eigenvalue of the clamped buckling problem
More informationMATH 220: INNER PRODUCT SPACES, SYMMETRIC OPERATORS, ORTHOGONALITY
MATH 22: INNER PRODUCT SPACES, SYMMETRIC OPERATORS, ORTHOGONALITY When discussing separation of variables, we noted that at the last step we need to express the inhomogeneous initial or boundary data as
More informationarxiv: v1 [math.ap] 21 Dec 2016
arxiv:1612.07051v1 [math.ap] 21 Dec 2016 On the extension to slip boundary conditions of a Bae and Choe regularity criterion for the Navier-Stokes equations. The half-space case. H. Beirão da Veiga, Department
More informationVANISHING VISCOSITY LIMIT FOR THE 3D NONHOMOGENEOUS INCOMPRESSIBLE NAVIER-STOKES EQUATION WITH SPECIAL SLIP BOUNDARY CONDITION
Electronic Journal of Differential Equations, Vol. 2017 (2017), No. 169, pp. 1 13. ISSN: 1072-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu VANISHING VISCOSITY LIMIT FOR THE 3D NONHOMOGENEOUS
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 informationWeek 2 Notes, Math 865, Tanveer
Week 2 Notes, Math 865, Tanveer 1. Incompressible constant density equations in different forms Recall we derived the Navier-Stokes equation for incompressible constant density, i.e. homogeneous flows:
More information1 Math 241A-B Homework Problem List for F2015 and W2016
1 Math 241A-B Homework Problem List for F2015 W2016 1.1 Homework 1. Due Wednesday, October 7, 2015 Notation 1.1 Let U be any set, g be a positive function on U, Y be a normed space. For any f : U Y let
More informationVORTICITY LAYERS OF THE 2D NAVIER-STOKES EQUATIONS WITH A SLIP TYPE BOUNDARY CONDITION
VORTICITY LAYERS OF THE 2D NAVIER-STOKES EQUATIONS WITH A SLIP TYPE BOUNDARY CONDITION GUNG-MIN GIE 1 AND CHANG-YEOL JUNG 2 Abstract. We study the asymptotic behavior, at small viscosity ε, of the Navier-
More informationOn the well-posedness of the Prandtl boundary layer equation
On the well-posedness of the Prandtl boundary layer equation Vlad Vicol Department of Mathematics, The University of Chicago Incompressible Fluids, Turbulence and Mixing In honor of Peter Constantin s
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 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 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 informationOn convergence criteria for incompressible Navier-Stokes equations with Navier boundary conditions and physical slip rates
On convergence criteria for incompressible Navier-Stokes equations with Navier boundary conditions and physical slip rates Yasunori Maekawa Department of Mathematics, Graduate School of Science, Kyoto
More informationON THE UNIQUENESS IN THE 3D NAVIER-STOKES EQUATIONS
ON THE UNIQUENESS IN THE 3D NAVIER-STOKES EQUATIONS Abdelhafid Younsi To cite this version: Abdelhafid Younsi. ON THE UNIQUENESS IN THE 3D NAVIER-STOKES EQUATIONS. 4 pages. 212. HAL Id:
More informationAppendix A Functional Analysis
Appendix A Functional Analysis A.1 Metric Spaces, Banach Spaces, and Hilbert Spaces Definition A.1. Metric space. Let X be a set. A map d : X X R is called metric on X if for all x,y,z X it is i) d(x,y)
More informationDYNAMIC BIFURCATION THEORY OF RAYLEIGH-BÉNARD CONVECTION WITH INFINITE PRANDTL NUMBER
DYNAMIC BIFURCATION THEORY OF RAYLEIGH-BÉNARD CONVECTION WITH INFINITE PRANDTL NUMBER JUNGHO PARK Abstract. We study in this paper the bifurcation and stability of the solutions of the Rayleigh-Bénard
More informationDIRECTION OF VORTICITY AND A REFINED BLOW-UP CRITERION FOR THE NAVIER-STOKES EQUATIONS WITH FRACTIONAL LAPLACIAN
DIRECTION OF VORTICITY AND A REFINED BLOW-UP CRITERION FOR THE NAVIER-STOKES EQUATIONS WITH FRACTIONAL LAPLACIAN KENGO NAKAI Abstract. We give a refined blow-up criterion for solutions of the D Navier-
More informationLecture 4: Numerical solution of ordinary differential equations
Lecture 4: Numerical solution of ordinary differential equations Department of Mathematics, ETH Zürich General explicit one-step method: Consistency; Stability; Convergence. High-order methods: Taylor
More informationTHE L 2 -HODGE THEORY AND REPRESENTATION ON R n
THE L 2 -HODGE THEORY AND REPRESENTATION ON R n BAISHENG YAN Abstract. We present an elementary L 2 -Hodge theory on whole R n based on the minimization principle of the calculus of variations and some
More informationSPECTRAL PROPERTIES OF THE LAPLACIAN ON BOUNDED DOMAINS
SPECTRAL PROPERTIES OF THE LAPLACIAN ON BOUNDED DOMAINS TSOGTGEREL GANTUMUR Abstract. After establishing discrete spectra for a large class of elliptic operators, we present some fundamental spectral properties
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 informationDS-GA 1002 Lecture notes 0 Fall Linear Algebra. These notes provide a review of basic concepts in linear algebra.
DS-GA 1002 Lecture notes 0 Fall 2016 Linear Algebra These notes provide a review of basic concepts in linear algebra. 1 Vector spaces You are no doubt familiar with vectors in R 2 or R 3, i.e. [ ] 1.1
More information2 A Model, Harmonic Map, Problem
ELLIPTIC SYSTEMS JOHN E. HUTCHINSON Department of Mathematics School of Mathematical Sciences, A.N.U. 1 Introduction Elliptic equations model the behaviour of scalar quantities u, such as temperature or
More informationSYMMETRY IN REARRANGEMENT OPTIMIZATION PROBLEMS
Electronic Journal of Differential Equations, Vol. 2009(2009), No. 149, pp. 1 10. ISSN: 1072-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu ftp ejde.math.txstate.edu SYMMETRY IN REARRANGEMENT
More information1. Subspaces A subset M of Hilbert space H is a subspace of it is closed under the operation of forming linear combinations;i.e.,
Abstract Hilbert Space Results We have learned a little about the Hilbert spaces L U and and we have at least defined H 1 U and the scale of Hilbert spaces H p U. Now we are going to develop additional
More informationFinite-dimensional spaces. C n is the space of n-tuples x = (x 1,..., x n ) of complex numbers. It is a Hilbert space with the inner product
Chapter 4 Hilbert Spaces 4.1 Inner Product Spaces Inner Product Space. A complex vector space E is called an inner product space (or a pre-hilbert space, or a unitary space) if there is a mapping (, )
More informationGlobal Maxwellians over All Space and Their Relation to Conserved Quantites of Classical Kinetic Equations
Global Maxwellians over All Space and Their Relation to Conserved Quantites of Classical Kinetic Equations C. David Levermore Department of Mathematics and Institute for Physical Science and Technology
More information23 Elements of analytic ODE theory. Bessel s functions
23 Elements of analytic ODE theory. Bessel s functions Recall I am changing the variables) that we need to solve the so-called Bessel s equation 23. Elements of analytic ODE theory Let x 2 u + xu + x 2
More informationarxiv: v1 [math.ap] 5 Nov 2018
STRONG CONTINUITY FOR THE 2D EULER EQUATIONS GIANLUCA CRIPPA, ELIZAVETA SEMENOVA, AND STEFANO SPIRITO arxiv:1811.01553v1 [math.ap] 5 Nov 2018 Abstract. We prove two results of strong continuity with respect
More informationTakens embedding theorem for infinite-dimensional dynamical systems
Takens embedding theorem for infinite-dimensional dynamical systems James C. Robinson Mathematics Institute, University of Warwick, Coventry, CV4 7AL, U.K. E-mail: jcr@maths.warwick.ac.uk Abstract. Takens
More informationTHE HOT SPOTS CONJECTURE FOR NEARLY CIRCULAR PLANAR CONVEX DOMAINS
THE HOT SPOTS CONJECTURE FOR NEARLY CIRCULAR PLANAR CONVEX DOMAINS YASUHITO MIYAMOTO Abstract. We prove the hot spots conjecture of J. Rauch in the case that the domain Ω is a planar convex domain satisfying
More informationLinear Algebra Massoud Malek
CSUEB Linear Algebra Massoud Malek Inner Product and Normed Space In all that follows, the n n identity matrix is denoted by I n, the n n zero matrix by Z n, and the zero vector by θ n An inner product
More informationEigenvalues and Eigenfunctions of the Laplacian
The Waterloo Mathematics Review 23 Eigenvalues and Eigenfunctions of the Laplacian Mihai Nica University of Waterloo mcnica@uwaterloo.ca Abstract: The problem of determining the eigenvalues and eigenvectors
More informationTHE POINCARÉ RECURRENCE PROBLEM OF INVISCID INCOMPRESSIBLE FLUIDS
ASIAN J. MATH. c 2009 International Press Vol. 13, No. 1, pp. 007 014, March 2009 002 THE POINCARÉ RECURRENCE PROBLEM OF INVISCID INCOMPRESSIBLE FLUIDS Y. CHARLES LI Abstract. Nadirashvili presented a
More informationOn a Suitable Weak Solution of the Navier Stokes Equation with the Generalized Impermeability Boundary Conditions
Proceedings of the 3rd IASME/WSEAS Int. Conf. on FLUID DYNAMICS & AERODYNAMICS, Corfu, Greece, August -, 5 pp36-41 On a Suitable Weak Solution of the Navier Stokes Equation with the Generalized Impermeability
More informationRemarks on the inviscid limit for the compressible flows
Remarks on the inviscid limit for the compressible flows Claude Bardos Toan T. Nguyen Abstract We establish various criteria, which are known in the incompressible case, for the validity of the inviscid
More informationThe Divergence Theorem Stokes Theorem Applications of Vector Calculus. Calculus. Vector Calculus (III)
Calculus Vector Calculus (III) Outline 1 The Divergence Theorem 2 Stokes Theorem 3 Applications of Vector Calculus The Divergence Theorem (I) Recall that at the end of section 12.5, we had rewritten Green
More informationRecurrence Relations and Fast Algorithms
Recurrence Relations and Fast Algorithms Mark Tygert Research Report YALEU/DCS/RR-343 December 29, 2005 Abstract We construct fast algorithms for decomposing into and reconstructing from linear combinations
More informationhal , version 6-26 Dec 2012
ON THE UNIQUENESS IN THE 3D NAVIER-STOKES EQUATIONS ABDEHAFID YOUNSI Abstract. In this paper, we give a new regularity criterion on the uniqueness results of weak solutions for the 3D Navier-Stokes equations
More informationLinear Algebra. Session 12
Linear Algebra. Session 12 Dr. Marco A Roque Sol 08/01/2017 Example 12.1 Find the constant function that is the least squares fit to the following data x 0 1 2 3 f(x) 1 0 1 2 Solution c = 1 c = 0 f (x)
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 informationMORE NOTES FOR MATH 823, FALL 2007
MORE NOTES FOR MATH 83, FALL 007 Prop 1.1 Prop 1. Lemma 1.3 1. The Siegel upper half space 1.1. The Siegel upper half space and its Bergman kernel. The Siegel upper half space is the domain { U n+1 z C
More informationOn 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 informationPage 404. Lecture 22: Simple Harmonic Oscillator: Energy Basis Date Given: 2008/11/19 Date Revised: 2008/11/19
Page 404 Lecture : Simple Harmonic Oscillator: Energy Basis Date Given: 008/11/19 Date Revised: 008/11/19 Coordinate Basis Section 6. The One-Dimensional Simple Harmonic Oscillator: Coordinate Basis Page
More informationLECTURE 12 A MODEL OF LIQUIDITY PREFERENCE (CONCLUDED)
LECTURE 12 A MODEL OF LIQUIDITY PREFERENCE (CONCLUDED) JACOB T. SCHWARTZ EDITED BY KENNETH R. DRIESSEL Abstract. We prove Theorem 11.6 and Theorem 11.9 of the previous lecture. 1. Proof of Theorem 11.6
More informationLaplace 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 informationHomogenization and error estimates of free boundary velocities in periodic media
Homogenization and error estimates of free boundary velocities in periodic media Inwon C. Kim October 7, 2011 Abstract In this note I describe a recent result ([14]-[15]) on homogenization and error estimates
More informationEuler Equations: local existence
Euler Equations: local existence Mat 529, Lesson 2. 1 Active scalars formulation We start with a lemma. Lemma 1. Assume that w is a magnetization variable, i.e. t w + u w + ( u) w = 0. If u = Pw then u
More informationMath 302 Outcome Statements Winter 2013
Math 302 Outcome Statements Winter 2013 1 Rectangular Space Coordinates; Vectors in the Three-Dimensional Space (a) Cartesian coordinates of a point (b) sphere (c) symmetry about a point, a line, and a
More informationLocal null controllability of the N-dimensional Navier-Stokes system with N-1 scalar controls in an arbitrary control domain
Local null controllability of the N-dimensional Navier-Stokes system with N-1 scalar controls in an arbitrary control domain Nicolás Carreño Université Pierre et Marie Curie-Paris 6 UMR 7598 Laboratoire
More informationPattern generation, topology, and non-holonomic systems
Systems & Control Letters ( www.elsevier.com/locate/sysconle Pattern generation, topology, and non-holonomic systems Abdol-Reza Mansouri Division of Engineering and Applied Sciences, Harvard University,
More informationPartial Differential Equations
Part II Partial Differential Equations Year 2015 2014 2013 2012 2011 2010 2009 2008 2007 2006 2005 2015 Paper 4, Section II 29E Partial Differential Equations 72 (a) Show that the Cauchy problem for u(x,
More informationWEYL S LEMMA, ONE OF MANY. Daniel W. Stroock
WEYL S LEMMA, ONE OF MANY Daniel W Stroock Abstract This note is a brief, and somewhat biased, account of the evolution of what people working in PDE s call Weyl s Lemma about the regularity of solutions
More informationUniqueness of ground states for quasilinear elliptic equations in the exponential case
Uniqueness of ground states for quasilinear elliptic equations in the exponential case Patrizia Pucci & James Serrin We consider ground states of the quasilinear equation (.) div(a( Du )Du) + f(u) = 0
More informationIn this chapter we study elliptical PDEs. That is, PDEs of the form. 2 u = lots,
Chapter 8 Elliptic PDEs In this chapter we study elliptical PDEs. That is, PDEs of the form 2 u = lots, where lots means lower-order terms (u x, u y,..., u, f). Here are some ways to think about the physical
More informationRadial Symmetry of Minimizers for Some Translation and Rotation Invariant Functionals
journal of differential equations 124, 378388 (1996) article no. 0015 Radial Symmetry of Minimizers for Some Translation and Rotation Invariant Functionals Orlando Lopes IMECCUNICAMPCaixa Postal 1170 13081-970,
More informationLecture 6: Lies, Inner Product Spaces, and Symmetric Matrices
Math 108B Professor: Padraic Bartlett Lecture 6: Lies, Inner Product Spaces, and Symmetric Matrices Week 6 UCSB 2014 1 Lies Fun fact: I have deceived 1 you somewhat with these last few lectures! Let me
More informationBoundary layers for Navier-Stokes equations with slip boundary conditions
Boundary layers for Navier-Stokes equations with slip boundary conditions Matthew Paddick Laboratoire Jacques-Louis Lions Université Pierre et Marie Curie (Paris 6) Séminaire EDP, Université Paris-Est
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 Liapunov Method for Determining Stability (DRAFT)
44 The Liapunov Method for Determining Stability (DRAFT) 44.1 The Liapunov Method, Naively Developed In the last chapter, we discussed describing trajectories of a 2 2 autonomous system x = F(x) as level
More informationFourth Order Convergence of Compact Finite Difference Solver for 2D Incompressible Flow
Fourth Order Convergence of Compact Finite Difference Solver for D Incompressible Flow Cheng Wang 1 Institute for Scientific Computing and Applied Mathematics and Department of Mathematics, Indiana University,
More informationDegenerate Monge-Ampère equations and the smoothness of the eigenfunction
Degenerate Monge-Ampère equations and the smoothness of the eigenfunction Ovidiu Savin Columbia University November 2nd, 2015 Ovidiu Savin (Columbia University) Degenerate Monge-Ampère equations November
More informationEXISTENCE OF SOLUTIONS TO ASYMPTOTICALLY PERIODIC SCHRÖDINGER EQUATIONS
Electronic Journal of Differential Equations, Vol. 017 (017), No. 15, pp. 1 7. ISSN: 107-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu EXISTENCE OF SOLUTIONS TO ASYMPTOTICALLY PERIODIC
More informationNumerical Sequences and Series
Numerical Sequences and Series Written by Men-Gen Tsai email: b89902089@ntu.edu.tw. Prove that the convergence of {s n } implies convergence of { s n }. Is the converse true? Solution: Since {s n } is
More informationEXISTENCE AND REGULARITY RESULTS FOR SOME NONLINEAR PARABOLIC EQUATIONS
EXISTECE AD REGULARITY RESULTS FOR SOME OLIEAR PARABOLIC EUATIOS Lucio BOCCARDO 1 Andrea DALL AGLIO 2 Thierry GALLOUËT3 Luigi ORSIA 1 Abstract We prove summability results for the solutions of nonlinear
More informationVectors in Function Spaces
Jim Lambers MAT 66 Spring Semester 15-16 Lecture 18 Notes These notes correspond to Section 6.3 in the text. Vectors in Function Spaces We begin with some necessary terminology. A vector space V, also
More informationEXISTENCE AND REGULARITY OF SOLUTIONS FOR STOKES SYSTEMS WITH NON-SMOOTH BOUNDARY DATA IN A POLYHEDRON
Electronic Journal of Differential Equations, Vol. 2017 (2017), No. 147, pp. 1 10. ISSN: 1072-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu EXISTENCE AND REGULARITY OF SOLUTIONS FOR
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 informationI teach myself... Hilbert spaces
I teach myself... Hilbert spaces by F.J.Sayas, for MATH 806 November 4, 2015 This document will be growing with the semester. Every in red is for you to justify. Even if we start with the basic definition
More informationExistence of Weak Solutions to a Class of Non-Newtonian Flows
Existence of Weak Solutions to a Class of Non-Newtonian Flows 1. Introduction and statement of the result. Ladyzhenskaya [8]. Newtonian : Air, other gases, water, motor oil, alcohols, simple hydrocarbon
More informationLinear Ordinary Differential Equations
MTH.B402; Sect. 1 20180703) 2 Linear Ordinary Differential Equations Preliminaries: Matrix Norms. Denote by M n R) the set of n n matrix with real components, which can be identified the vector space R
More informationSpectral theory for magnetic Schrödinger operators and applicatio. (after Bauman-Calderer-Liu-Phillips, Pan, Helffer-Pan)
Spectral theory for magnetic Schrödinger operators and applications to liquid crystals (after Bauman-Calderer-Liu-Phillips, Pan, Helffer-Pan) Ryukoku (June 2008) In [P2], based on the de Gennes analogy
More informationA LOCALIZATION PROPERTY AT THE BOUNDARY FOR MONGE-AMPERE EQUATION
A LOCALIZATION PROPERTY AT THE BOUNDARY FOR MONGE-AMPERE EQUATION O. SAVIN. Introduction In this paper we study the geometry of the sections for solutions to the Monge- Ampere equation det D 2 u = f, u
More informationA Ginzburg-Landau Type Problem for Nematics with Highly Anisotropic Elastic Term
A Ginzburg-Landau Type Problem for Nematics with Highly Anisotropic Elastic Term Peter Sternberg In collaboration with Dmitry Golovaty (Akron) and Raghav Venkatraman (Indiana) Department of Mathematics
More informationOn Smoothness of Suitable Weak Solutions to the Navier-Stokes Equations
On Smoothness of Suitable Weak Solutions to the Navier-Stokes Equations G. Seregin, V. Šverák Dedicated to Vsevolod Alexeevich Solonnikov Abstract We prove two sufficient conditions for local regularity
More informationGAKUTO International Series
1 GAKUTO International Series Mathematical Sciences and Applications, Vol.**(****) xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx, pp. xxx-xxx NAVIER-STOKES SPACE TIME DECAY REVISITED In memory of Tetsuro Miyakawa,
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 informationand finally, any second order divergence form elliptic operator
Supporting Information: Mathematical proofs Preliminaries Let be an arbitrary bounded open set in R n and let L be any elliptic differential operator associated to a symmetric positive bilinear form B
More 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 informationA generalised Ladyzhenskaya inequality and a coupled parabolic-elliptic problem
A generalised Ladyzhenskaya inequality and a coupled parabolic-elliptic problem Dave McCormick joint work with James Robinson and José Rodrigo Mathematics and Statistics Centre for Doctoral Training University
More informationMULTIPLE SOLUTIONS FOR THE p-laplace EQUATION WITH NONLINEAR BOUNDARY CONDITIONS
Electronic Journal of Differential Equations, Vol. 2006(2006), No. 37, pp. 1 7. ISSN: 1072-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu ftp ejde.math.txstate.edu (login: ftp) MULTIPLE
More informationHARDY 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 informationTensors, and differential forms - Lecture 2
Tensors, and differential forms - Lecture 2 1 Introduction The concept of a tensor is derived from considering the properties of a function under a transformation of the coordinate system. A description
More information