Contour dynamics of incompressible 3-D fluids in a porous medium with different densities
|
|
- Jewel Martin
- 5 years ago
- Views:
Transcription
1 Contour dynamics of incompressible 3-D fluids in a porous medium with different densities Diego Córdoba Francisco Gancedo Abstract We consider the problem of the evolution of the interface given by two incompressible fluids through a porous medium, which is known as the Muskat problem in two dimensions it is mathematically analogous to the two-phase Hele-Shaw cell. We focus on a fluid interface given by a jump of densities, being the equation of the evolution obtained using Darcy s law. We prove local well-posedness when the smaller density is above (stable case) in the unstable case we show ill-posedness. Introduction The evolution of a fluid in a porous medium is an important interesting topic of fluid mechanics (see [3]). This phenomena is based on an experimental physical principle given by H. Darcy in 856. Darcy s law for a 3-D fluid is given by the momentum equation µ v = p (,, g ρ), κ where v is the incompressible velocity, p is the pressure, µ is the dynamic viscosity, κ is the permeability of the medium, ρ is the liquid density g is the acceleration due to gravity. A different problem is the motion of a -D fluid in a Hele Shaw cell (see []). In this case the fluid is set between two fixed parallel plates. These plates are close enough in such a way that the mean velocity is described by µ v = p (, g ρ), b where b denotes the distance between the plates. Considering that the fluid in the porous medium only moves in two directions suppressing one of the variables in the horizontal plane, these two different physical phenomena of fluid dynamics become nevertheless mathematically analogous if we identify the permeability of the medium κ the constant b /. The Muskat problem (see [4]) the two-phase Hele-Shaw flow (see [7]) model the evolution of an interface between two fluids (in a porous medium in a Hele Shaw cell respectively) with different viscosities densities. A lot of information can be found in the literature about both problems (see references in [9] [3]). These free boundary problems are considered with surface tension using the Laplace Young condition also without surface tension in which case the pressures are equal on the interface. With surface tension, in the two dimensional case has been proven that the problems have classical solutions (see []). Without surface tension, Siegel, Caflisch Howison [8] proved ill-posedness in an unstable -D case, namely when the higher-viscosity fluid contracts, they show global-in-time existence of small initial data in the stable case when the higher-viscosity fluid exps. The results rely on the assumption that the Atwood number A µ = µ µ µ + µ
2 Ω (t) X 3 Ω (t) X X is nonzero where µ µ are the viscosities of the fluids. In the same year, Ambrose [] treated the -D problem with an initial data fulfilling the following condition (ρ ρ )g cos(θ(α, )) + A µ U(α, ) >, (x(α, ) x(α, )) + (y(α, ) y(α, )) (α α ) >, () where the interface is the curve (x(α, t), y(α, t)), ρ ρ are the densities of the fluids, θ is the angle that the tangent to the curve forms with the horizontal U is the normal velocity (given by the Birkhoff-Rott integral). We are interested in the case A µ = that presents the evolution of the interface for different densities. This case, for example, models moist dry regions in a porous medium. Meanwhile the work of Ambrose is based on the arclength the tangent angle formulation used by Hou, Lowengrub Shelley [3], due to the particular form of the vorticity in the case A µ =, we get to parameterize the curve in the two dimensional problem getting the condition () for any time (see equation (6)). The free boundary problems given by fluids with different densities are been intensely studied. Notice the classical paper of Taylor [] the works of Wu [3] [4] where the full water wave problem is solved considering the water with positive density the air with zero density. A study of the two-dimensional case can be found in [] due to Ambrose Masmoudi. In order to simplify the notation, we consider µ/k = µ/b = g =. Thus, the 3-D system is written as v(x, x, x 3, t) = p(x, x, x 3, t) (,, ρ(x, x, x 3, t)), () where (x, x, x 3 ) R 3 are the spatial variables t denotes the time. Here ρ is defined by { ρ in Ω ρ(x, x, x 3, t) = (t) ρ in Ω (t), with ρ, ρ constants ρ ρ. We show in section that in this case it is not necessary to assume any condition on the pressure along the interface to obtain the contour equation. Furthermore, we illustrate below that the solutions to this model are weak solutions to the following conservation of mass equation Dρ Dt = ρ t + v ρ =, (3)
3 where div v =. We notice the similarity with the -D vortex patch problem given by the two-dimensional Euler equation where the vorticity is conserved along trajectories in a weak sense. The vorticity is considered to be a characteristic function of a domain. Chemin [8] proved global-in-time existence using paradifferential calculus. A simpler proof can be found in [4] due to Bertozzi Constantin. In section we show that due to () the velocity can be determined from the density by singular integral operators (see []). It makes the equation more singular than the -D vortex patch problem where the velocity is given by the Biot-Savart law. A singular problem, more analogous to (3), is the evolution of the -D quasi-geostrophic (QG) equation for sharp fronts. The QG equation models the dynamics of cold hot air the formation of fronts. Here the temperature θ is conserved along particle trajectories the velocity is given by singular integral operators in the following form v = ( R θ, R θ), where R R are the Riesz transforms (see [] for more details of the QG equation). Rodrigo [9] proposed the contour equation of the sharp fronts where the temperature is concentrated in a domain proved local existence uniqueness. The paper is organized as follows. In section we derive the contour equation. We show that this equation fulfills the conservation of mass equation in section 3. In section 4 we prove local existence uniqueness of the stable case. In section 5 we get a family of global solutions of the -D stable case with small initial data. Finally, as a consequence of the previous section, in 6 we prove ill-posedness for the 3-D unstable case. The Contour Equation We consider the equation with (x, x, x 3 ) R 3 where the fluid has different densities, that is ρ is represented by { ρ, {x ρ(x, x, x 3, t) = 3 > f(x, x, t)} (4) ρ, {x 3 < f(x, x, t)}, being f the interface. Using Darcy s Law we get curl curl v = ( x x3 ρ, x x3 ρ, ( x + x )ρ). Since div v = we have curl curl v = v, therefore it follows v = ( x x3 ρ, x x3 ρ, ( x + x ) ρ). (5) The integral operators x x are given by the kernels K (x, x, x 3 ) = 4π respectively, thus the velocity can be expressed by Since ρ satisfies (4) we have x (x + x + x 3 )3/, K (x, x, x 3 ) = 4π x (x + x + x 3 )3/, v = (K x3 ρ, K x3 ρ, K x ρ K x ρ). (6) ρ = (ρ ρ )( x f(x, x, t), x f(x, x, t), )δ(x 3 f(x, x, t)), (7) 3
4 where δ is the Dirac distribution. Using (6) we obtain v(x, x, x 3, t) = ρ ρ (y, y, f(x y, t) y) P V dy, (8) 4π R [ y + (x 3 f(x y, t)) ] 3/ where we note x = (x, x ), y = (y, y ) f(x y, t) y = x f(x y, t)y + x f(x y, t)y. In (8) x 3 f(x, t) the principal value is taken at infinity (see []). When x 3 approaches f(x, t) in the normal direction, we get a discontinuity on the velocity due to the fact that the vorticity is concentrated on the interface. Thus, for ε > we define v (x, f(x, t), t) = lim ε v(x ε x f(x, t), x ε x f(x, t), f(x, t) + ε, t), v (x, f(x, t), t) = lim ε v(x + ε x f(x, t), x + ε x f(x, t), f(x, t) ε, t). It follows v (x, f(x, t), t) = ρ ρ (y, y, f(x y, t) y) P V 4π R [ y + (f(x, t) f(x y, t)) ] + ρ ρ + ρ ρ x f(x, t)(,, x f(x, t)) + ( x f(x, t)) + ( x f(x, t)) x f(x, t)(,, x f(x, t)) + ( x f(x, t)) + ( x f(x, t)), 3/ dy v (x, f(x, t), t) = ρ ρ (y, y, f(x y, t) y) P V dy 4π R [ y + (f(x, t) f(x y, t)) ] 3/ ρ ρ ρ ρ x f(x, t)(,, x f(x, t)) + ( x f(x, t)) + ( x f(x, t)) x f(x, t)(,, x f(x, t)) + ( x f(x, t)) + ( x f(x, t)). (9) () The velocity in the tangential directions only moves the particles on the surface f(x, t); i.e., if we rewrite the velocity in the tangential directions, we only make a change on the parametrization do not alter the shape of the interface. Thus, it follows that due to the fact that the terms v(x, f(x, t), t) = ρ ρ P V 4π ± ρ ρ ± ρ ρ (y, y, f(x y, t) y) dy, () R [ y + (f(x, t) f(x y, t)) ] 3/ x f(x, t)(,, x f(x, t)) + ( x f(x, t)) + ( x f(x, t)), x f(x, t)(,, x f(x, t)) + ( x f(x, t)) + ( x f(x, t)), are in the tangential directions. Moreover, if we add the following tangential terms to () ρ ρ y P V 4π R [ y + (f(x, t) f(x y, t)) ] dy(,, 3/ x f(x, t)), 4
5 we obtain ρ ρ y P V 4π R [ y + (f(x, t) f(x y, t)) ] dy(,, 3/ x f(x, t)), v(x, f(x, t), t) = ρ ρ ( f(x, t) f(x y, t)) y (,, P V dy). () 4π R [ y + (f(x, t) f(x y, t)) ] 3/ Finally we have the contour equation given by df dt (x, t) = ρ ρ ( f(x, t) f(x y, t)) y P V dy, 4π [ y + (f(x, t) f(x y, t)) ] 3/ f(x, ) = f (x). R In the periodic case, we can obtain an equivalent equation to (3) due to the integral operators x x can be presented by the kernels K p (x, x, x 3 ) = 4π ( x (x + x + x 3 )3/ L(x, x, x 3 ) + M(x, x, x 3 )), K p (x, x, x 3 ) = 4π ( x (x + x + x 3 )3/ L(x, x, x 3 ) + M(x, x, x 3 )), respectively for (x, x, x 3 ) T R with T = [ π, π] the functions L, M C (T R) (see [] for the kernel of the Riesz potentials on the n-torus). Adding an appropiate function to the singular part of K p KP, we can choose L C c (T R), L, supp L {x + x + x 3 4}, L = in {x + x + x 3 } L( x, x, x 3 ) = L(x, x, x 3 ). (3) (4) The function M belongs to C b (T R) M(,, ) =. The velocity can be expressed by v = (K p x 3 ρ, K p x 3 ρ, K p x ρ K p x ρ), due to (7) it follows (suppressing the dependence on t) v(x, x, x 3 ) = ρ ρ (y, y, f(x y) y) P V 4π T [ y + (x 3 f(x y)) ] L(y, x 3/ 3 f(x y))dy ρ ρ (,, f(x y) (, ))M(y, x 3 f(x y))dy, 4π T if x 3 f(x). Adding a term in the tangential direction we obtain v(x, f(x)) = ρ ρ ( f(x) f(x y)) y (,, L(y, f(x) f(x y))dy 4π T [ y + (f(x) f(x y)) ] 3/ + ( f(x) f(x y)) (, ))M(y, f(x) f(x y))dy). T Finally we have the contour equation in the periodic case given by df dt (x, t) =ρ ρ ( f(x, t) f(x y, t)) y L(y, f(x, t) f(x y, t))dy 4π T [ y + (f(x, t) f(x y, t)) ] 3/ + ρ ρ ( f(x, t) f(x y, t)) (, ))M(y, f(x, t) f(x y, t))dy, 4π T f(x, ) =f (x). (5) 5
6 We use both formulations throughout the paper. Suppose that the function f(x) only depends on x in equation (3). Then the contour equation in the -D case (with a -D interface) follows df dt (x, t) = ρ ρ P V π f(x, ) = f (x); x R. R ( x f(x, t) x f(x α, t))α α + (f(x, t) f(x α, t)) dα, This equation can be obtained in a similar way that (3) using the stream function. Performing a two-dimensional analysis using the stream function, we obtain an equivalent equation to (6) in the two dimensional periodic case as follows df dt (x, t) =ρ ρ π T + ρ ρ π f(x, ) =f (x), ( x f(x, t) x f(x α, t))α α P (α, f(x, t) f(x α, t))dα + (f(x, t) f(x α, t)) ( x f(x, t) x f(x α, t))q(α, f(x, t) f(x α, t))dα, T (6) (7) with P (x, x ) C c (T R), P, supp P {x + x 4}, P = in {x + x } P ( x, x ) = P (x, x ). The function Q(x, x ) belongs to Cb (T R) Q(, ) =. If we consider the linearized equation of the motion, we obtain a dissipative equation when ρ < ρ (the greater density is below) an unstable equation when ρ > ρ. The unstable linearized equation presents an instability similar to the Kelvin-Helmholtz s (see [6]). As usual, we note the Riesz transforms in R (see []) R f(x) = π P.V. y f(x y)dy, R y 3 R f(x) = π P.V. y f(x y)dy, R y 3 the operator Λ s f defined by the fourier transform Λ s f(ξ) = ξ s f(ξ). Suppose that f(x) is uniformly small we can neglect the terms of order grater than one in (3), then it reduces to the following linear equation f t = ρ ρ (R x f + R x f) = ρ ρ Λf, f(x, ) = f (x). (8) Applying the fourier transform we get ˆf(ξ) = ˆf (ξ)e ρ ρ ξ t, therefore (8) is a dissipative equation when ρ < ρ an ill posed problem in the case ρ > ρ with a general initial data in the Schwartz class. We need an analytic initial data in order to get a well posed problem for ρ > ρ. 6
7 3 The conservation of mass equation We show that if ρ is defined by (4) f(x, t) is convected by the velocity () then ρ is a weak solution of the conservation of mass equation (3) conversely. From now on, Ω is equal to R or T x = (x, x, x 3 ). Definition 3. The density ρ is a weak solution of the conservation of mass equation if for any ϕ C (Ω R (, T )), ϕ with compact support in the real case periodic in (x, x ) otherwise, we have T (ρ( x, t) t ϕ( x, t) + v( x, t)ρ( x, t) ϕ( x, t))d xdt =, (9) Ω R where the incompressible velocity v is given by Darcy s law. Then Proposition 3. If f(x, t) satisfies (3) ρ( x, t) is defined by (4), then ρ is a weak solution of the conservation of mass equation. Furthermore, if ρ is a weak solution of the conservation of mass equation given by (4), then f(x, t) satisfies (3). Proof: Let ρ be a weak solution of (3) defined by (4). Integrating by parts we have T T T I = ρ t ϕd xdt = ρ t ϕd xdt + ρ t ϕd xdt Ω = (ρ ρ ) R T Ω {x 3 >f} ϕ(x, f(x, t), t) t f(x, t)dxdt. On the other h, due to (9) () we obtain T T T J = ρv ϕ d xdt = ρ v ϕd xdt + ρ = T Ω Ω = (ρ ρ ) R {x 3 >f} {x 3 <f} {x 3 <f} v ϕd xdt ϕ(x, f(x, t), t)(ρ v (x, f(x, t), t) ρ v (x, f(x, t), t)) ( x f(x, t), x f(x, t), )dxdt T Ω ϕ(x, f(x, t), t)v(x, f(x, t), t) ( x f(x, t), x f(x, t), )dxdt, where v(x, f(x, t), t) is given by (). We get J = (ρ ρ ) T ϕ(x, f(x, t), t)p V 4π Ω R ( f(x, t) f(x y, t)) y dt. [ y + (f(x, t) f(x y, t)) ] 3/ Then I + J = due to (9). Thus, if we choose ϕ( x, t) = ϕ(x, t) for x 3 [ f L, f L ] follows that f(x, t) fulfils (3). Following the same arguments it is easy to check that if f(x, t) satisfies (3), then ρ is a weak solution given by (4). Remark 3.3 Note that due to (5), the velocity satisfies v = (R (R 3 ρ), R (R 3 ρ), (R + R )(ρ)), where the operators R, R R 3 are the Riesz transforms in three dimensions (see []). Since ρ L (Ω R) then v belongs to BMO (bounded mean oscillation) therefore v is in L (Ω R) locally (see [] for the definitions properties of the BMO space). 7
8 4 Local well-posedness for the stable case In this section we prove local existence uniqueness for the stable case using energy estimates. First we study the case Ω = R at the end of the section we give the main differences with the periodic domain. Denote the Sobolev spaces by H k, the Holder spaces by C k,δ with δ < the Holder continuity the hessian matrix of a function f(x) by f(x). The norms of H k C k,δ are defined as follows f H k = f L + Λk f L, 4. Case Ω = R f C k,δ The main theorem in this section is the following = f C k + max max x i x j f(x) x i x j f(y) i+j=k x y x y δ. Theorem 4. Let f (x) H k (R ) for k 4 ρ > ρ. Then there exists a time T > so that there is a unique solution to (3) in C ([, T ]; H k (R )) with f(x, ) = f (x). We choose ρ ρ = 4π without loss of generality, then df ( f(x, t) f(x y, t)) y (x, t) = P V dy, dt [ y + (f(x, t) f(x y, t)) ] 3/ f(x, ) = f (x). R () We show the proof with k = 4 being analogous for k > 4. We apply energy methods (see [5] for more details). Then d dt f L (t) = ( f(x) f(x y)) y f(x)p V R R [ y + (f(x) f(x y)) ] 3/ ( f(x) f(x y)) y = f(x) R y < [ y + (f(x) f(x y)) ] 3/ f(x) y + f(x)p V R y > [ y + (f(x) f(x y)) ] 3/ f(x y) y f(x)p V R y > [ y + (f(x) f(x y)) ] 3/ = I + I + I 3. The identity yields I C C xi f(x) xi f(x y) = ds ds y < y < xi f(x + (s )y) y ds, y f(x) f(x + (s )y) dxdy R [ + ((f(x) f(x y)) y ] 3/ y dy f L x i x j f L C f H. i+j= 8
9 Integrating by parts, the term I is written I = 3 (f(x) f(x y))( f(x) f(x y)) y f(x) dxdy y > R [ y + ((f(x) f(x y)) ] 5/ C y 3 f(x) f(x) f(x y) y f(x) f(x y) dxdy y > R [ + ((f(x) f(x y)) y ] 5/ C f L f H. Integrating by parts in I 3, it follows R I 3 = f(x)f(x y) y (f(x) f(x y)) dxdy y > [ y + (f(x) f(x y)) ] 5/ (f(x) f(x y)) f(x y) y + 3 f(x)f(x y) y > R [ y + (f(x) f(x y)) ] dxdy 5/ y f(x)f(x y) dxdσ(y) y = R [ y + (f(x) f(x y)) ] 3/ C( f L + ) f H. Using Sobolev inequalities, we get finally We consider the quantity where d dt f L (t) C( f 3 H(t) + ). () d dt 4 x f L (t) = I 4 + I 5 + I 6 + I 7 + I 8, R I 4 = 4 ( x 4 x f(x)p V f(x) x 4 f(x y)) y R [ y + (f(x) f(x y)) ], 3/ I 5 = 4 x 4 f(x) ( x 3 f(x) x 3 f(x y)) y x A(x, y), R R I 6 = 6 x 4 f(x) ( x f(x) x f(x y)) y x A(x, y), R R I 7 = 4 x 4 f(x) ( x f(x) x f(x y)) y x 3 A(x, y), R R I 8 = x 4 f(x) ( f(x) f(x y)) y x 4 A(x, y), R R A(x, y) = [ y + (f(x) f(x y)) ] 3/. The most singular term is I 4. In order to estimate it R I 4 = 4 x 4 x f(x)p V f(x, t) y R [ y + (f(x, t) f(x y, t)) ] 3/ R 4 x 4 x f(x)p V f(y, t) (x y) R [ x y + (f(x, t) f(y, t)) ] 3/ = J + J. 9
10 Integrating by parts J = 3 x 4 f(x) P V R = 3 x 4 f(x) ( R R y > (f(x) f(x y))( f(x) f(x y)) y [ y + ((f(x) f(x y)) ] 5/ dy + P V dy)dx y < () 3 f C 4 x f L + 3 M(f) 4 x f L, where M(f) = max x P V y < (f(x) f(x y))( f(x) f(x y)) y dy. [ y + ((f(x) f(x y)) ] 5/ We estimate this maximum in the following form (f(x) f(x y) f(x) y)( f(x) f(x y)) y M(f) max dy x y < [ y + ((f(x) f(x y)) ] 5/ ( f(x) y)(( f(x) f(x y)) y y f(x) y) + max dy x y < [ y + ((f(x) f(x y)) ] 5/ + max ( f(x) y)(y f(x) y)(b(x, y) C(x, y))dy x y < ( f(x) y)(y f(x) y) + max P V dy, x [ y + ( f(x) y) ] 5/ y < (3) where B(x, y) = [ y + ((f(x) f(x y)) ] 5/, C(x, y) = [ y + ( f(x) y) ] 5/. Making the change of variables y = z, we obtain that the last integral in (3) is null, then we can estimate M(f) by y < + f C f C,δ M(f) f C max x with < δ <, having finally y [ + ((f(x) f(x y)) y ) ] 5/ dy y < y +δ dy + f C f C C( f C + f C f C,δ + f C f C ), In order to estimate J, we integrate by parts getting y < y dy J C( f 4 C +,δ ) 4 x f L. (4) R J = 4 y ( x 4 x f(x)p V f(x) x 4 f(y)) (x y) R [ x y + (f(x) f(y)) ] 3/ = K + K, with R K = 4 x 4 x f(x)p V f(x) x 4 f(y), R [ x y + (f(x) f(y)) ] 3/
11 K = x 4 f(x) ( x 4 f(x) 4 3(f(x) f(y))(f(x) f(y) f(y) (x y)) x f(y)). R R [ x y + (f(x) f(y)) ] 5/ Making a change of variables we can obtain K = P V 4 x 4 x f(x) f(x) x 4 f(y) R R [ x y + (f(x) f(y)) ] 3/ = P V 4 x 4 x f(y) f(x) x 4 f(y) R R [ x y + (f(x) f(y)) ] 3/ = R ( x 4 f(x) x 4 f(y)) [ x y + (f(x) f(y)) ] 3/. R Here we observe the main difference with the unstable case in which we obtain the opposite sign. Now we consider K = L + L + L 3, being L = 3 x 4 f(x) P V R L = 3 4 P V with R R (f(x) f(y))(f(x) f(y) f(y) (x y)), [ x y + (f(x) f(y)) ] 5/ x 4 f(x) x 4 f(y) (f(x) f(y))(x y) ( f(x) + f(y)) (x y), R [ x y + (f(x) f(y)) ] 5/ L 3 = 3 x 4 f(x) x 4 f(y)(f(x) f(y))d(x, y), R R D(x, y) = (f(x) f(y) f(y) (x y) 4 (x y) ( f(x) + f(y)) (x y)) [ x y + (f(x) f(y)) ] 5/. The term L can be estimated like J in () one finds that L C( + f 4 C,δ) 4 x f L. Exchanging x for y we obtain that L =. For the last term L 3, it follows L 3 C x 4 f(x) f(x) f(y) D(x, y) R R + C x 4 f(y) f(x) f(y) D(x, y) dxdy R R C x 4 f(x) f(x) f(x y) D(x, x y) R R + C x 4 f(y) f(x + y) f(y) D(x + y, y) dxdy R R C x 4 f(x) dx( dy + dy) + x 4 f(y) dy( R y < y > R C f C f C,δ 4 x f L. x < dx + dx) x >
12 Finally J = K + K K = L + L + L 3 = L + L 3 C( f 4 C,δ + ) 4 x f L, due to (4) we obtain I 4 C( f 4 C +,δ ) 4 x f L. (5) Now we estimate the integral I 5. We have I 5 = J 3 + J 4 where J 3 = 4 x 4 f(x) x 3 f(x) P V y (f(x) f(x y))( x f(x) x f(x y)), R R [ y + (f(x) f(x y)) ] 5/ J 4 = 4P V = 4P V R R x 4 f(x) x 3 f(y) (x y) (f(x) f(y))( x f(x) x f(y)) R [ x y + (f(x) f(y)) ] 5/ x 4 f(x) x x 3 f(y)(x y ) (f(x) f(y))( x f(x) x f(y)). R [ x y + (f(x) f(y)) ] 5/ The way to estimate J 3 is similar to the term J in (), we find that J 3 C( + f 4 C,δ) f H 4. We decompose the term J 4 = K 3 + K 4 + K 5 + K 6 as follows K 3 = 4P V R x 4 f(x) x x 3 f(x y)y E(x, y), R K 4 = 4P V R x 4 f(x) x x 3 f(x y)y F (x, y), R K 5 = 4P V x 4 f(x) x x 3 f(x y)y ( f(x) y)( x f(x) y)(b(x, y) C(x, y)) where R y < K 6 = 4P V R y < x 4 f(x) x x 3 ( f(x) y)( x f(x) y) f(x y)y [ y + ( f(x) y) ], 5/ E(x, y) = (f(x) f(x y) f(x) y)( x f(x) x f(x y)) [ y + (f(x) f(x y)) ] 5/, F (x, y) = ( f(x) y)( x f(x) x f(x y) x f(x) yx { y <} ) [ y + (f(x) f(x y)) ] 5/. The terms K 3, K 4, K 5 are estimated as J. Then we find that K 3 C( + f C ) f H 4, We rewrite K 6 we get K 4 C( + f C f C,δ) f H 4 K 5 C( + f C f C ) f H 4. K 6 = 4P V x 4 f(x)s( x x 3 f)(x)dx, R
13 with the operator S defined by S(g)(x) = P V The function Σ(x, y) satisfies that y < Σ(x, y) = y ( f(x) y Σ(x, y) y g(x y)dy, y y )( x f(x) [ + ( f(x) (i) Σ(x, λy) = Σ(x, y), λ >, (ii) (iii) Σ(x, y) = Σ(x, y), sup Σ(x, y) x f L, x y y ) y. (6) y ) ] 5/ therefore S is a bounded linear map on L p (R ) for < p < S p C x f L (see [] references therein for more details). Then K 6 C f C 4 x f L x 3 x f L. We obtain finally I 5 C( + f 4 C,δ) f H 4. In order to estimate the term I 6 we take I 6 = 6 ds dy x 4 f(x) y ( x f(x + (s )y)) y x A(x, y)dx R R ds( dy + dy) x 4 f(x) x f(x + (s )y) x A(x, y) y dx y < y > R C( y +δ dy + y 3 dy)( + f 4 C,δ) f H 4. y < y > The most singular term of I 7 is K 7 K 7 = x 4 f(x)( x f(x) x f(x y)) y G(x, y), R R where G(x, y) = (f(x) f(x y))( 3 x f(x) 3 x f(x y)) [ y + (f(x) f(x y)) ] 5/. Due to x f(x) x f(x y) f C,δ y δ writing 3 x f(x) 3 x f(x y) = 3 x f(x + (s )y) y ds, we obtain K 7 C( f C + ) f,δ H I 4 7 C( + f 4 C ) f,δ H. 4 The most singular term of I 8 is K 8 K 8 = x 4 f(x)( x 4 f(x) x 4 f(x y)) H(x, y), R R where H(x, y) = (f(x) f(x y))( f(x) f(x y)) y [ y + (f(x) f(x y)) ] 5/. 3
14 Then K 8 = x 4 f(x) P V R R (f(x) f(x y))( f(x) f(x y)) y, [ y + (f(x) f(x y)) ] 5/ is controlled as before. We obtain K 8 C( + f 4 C,δ ) 4 x f L I 8 C( + f 4 C,δ ) f H 4. Finally, we have d dt 4 x f L (t) C( + f 4 C,δ(t)) f H 4(t), using Sobolev inequalities we get In a similar way we obtain d dt 4 x f L (t) C( f 6 H4(t) + ). (7) d dt 4 x f L (t) C( f 6 H4(t) + ), (8) since we can define f H = f 4 L + 4 x f L + 4 x f L, due to (), (7) (8) it follows d dt f H 4(t) C( f 5 H4(t) + ). Using Gronwall s inequality we get that the quantity f H 4 is bounded up to a time T = T ( f H 4). Then, applying energy methods the local existence result follows. Let the functions f (x, t), f (x, t) be two solutions of equation (3) with f (x, ) = f (x, ) = f (x), f = f f. Then d dt f L (t) = I 9 + I + I, with I 9 = f(x) f(x) P V R y[ y + (f (x) f (x y)) ] 3/, R I = f(x)p V R f(y, t) (x y)[ x y + (f (x) f (y)) ] 3/ R I = f(x)p V R ( f (x) f (x y)) yn(x, y), R with N(x, y) = [ y + (f (x) f (x y)) ] 3/ [ y + (f (x) f (x y)) ] 3/. Integrating by part in I 9 we have I 9 C( f H 4) f L, I C( f H 4, f H 4) f L. The term I = f(x)p V y (f(y) f(x)) (x y)[ x y + (f (x) f (y)) ] 3/ R R = P V f(x)(f(x) f(y))[ x y + (f (x) f (y)) ] 3/ R R + P V f(x)(f(x) f(y)) 3(f (x) f (y))(f (x) f (y) f (x)(x y)). R R [ x y + (f (x) f (y)) ] 5/ Then we have that I J 5 + J 6 where J 5 = f(x) 3(f (x) f (y))(f (x) f (y) f (x)(x y)) P V R R [ x y + (f (x) f (y)) ] 5/ 4
15 J 6 = f(x)p V f(x y) 3(f (x) f (x y))(f (x) f (x y) f (x) y). R R [ y + (f (x) f (x y)) ] 5/ The term J 5 is estimated as J obtaining J 5 C( f H 4) f L. The term J 6 can be expressed as J 6 = K 9 + K with (f (x) f (x y))g(x, y) K 9 = 3 f(x) f(x y), R R [ y + (f (x) f (x y)) ] 5/ where the function G(x, y) is given by G(x, y) = f (x) f (x y) f (x) y 4 y ( f (x) + f (x y)) y. One finds that the principal value K is null. Therefore, we obtain finally J 6 C( f H 4) f L. Applying Gronwall s inequality we get uniqueness. 4. Case Ω = T In the periodic case we give the theorem of local well-posedness the differences with Ω = R. Theorem 4. Let f (x) H k (T ) for k 4 ρ > ρ. Then there exists a time T > so that there is a unique solution to (5) in C ([, T ]; H k (T )) with f(x, ) = f (x). The proof is similar to the theorem 4. but we must use the properties of the function L in (4). We consider without loss of generality ρ ρ = 4π. In order to control the evolution of the quantity 4 x f L, the most singular term is T I = 4 ( x 4 x f(x) f(x) x 4 f(x y)) y L(y, f(x) f(x y)) T [ y + (f(x) f(x y)) ] 3/ = x 4 f(x) 4 L(y, f(x) f(x y)) x f(x) P V y T T [ y + (f(x) f(x y)) ] 3/ T 4 x 4 x f(x)p V f(y) (x y) L(x y, f(x) f(y)) T [ x y + (f(x) f(y)) ] 3/ = J + J. Integrating by parts J = 3 x 4 f(x) P V A(x, y)l(y, f(x) f(x y)) T T x 4 f(x) L x3 (y, f(x) f(x y))( f(x) f(x y)) y P V T T [ y + (f(x) f(x y)) ] 3/ = K + K, (9) where A(x, y) = (f(x) f(x y))( f(x) f(x y)) y [ y + ((f(x) f(x y)) ] 5/ L x3 (x, x, x 3 ) = x3 L(x, x, x 3 ). 5
16 Due to L(x, x, x 3 ) C (x, x, x 3 ) we have K C( + f 4 C,δ) f H 4 + C 4 x f L max x C( + f 4 C,δ) f H 4. P V ( f(x) y)(y f(x) y) T [ y + ( f(x) y)) ] 5/ Using that f(x) f(x y) f C y L x3 = in {x + x + x 3 4}, we have that K = x 4 f(x) L x3 (y, f(x) f(x y))( f(x) f(x y)) y T [ y + (f(x) f(x y)) ] 3/ C( f 4 C,δ + ) f H 4. y > + f C In order to estimates J, we integrate by parts getting T J = 4 y ( x 4 x f(x)p V f(y) x 4 f(x)) (x y) L(x y, f(x) f(y)) T [ x y + (f(x) f(y)) ] 3/ = K 3 + K 4 + K 5 + K 6 with K 3 = x 4 f(x)p V T x 4 f(x) x 4 f(y) L(x y, f(x) f(y)), [ x y + (f(x) f(y)) ] 3/ K 4 = x 4 f(x) T ( x 4 f(x) x 4 f(y))b(x, y)l(x y, f(x) f(y)), T K 5 = x 4 f(x)( x 4 f(y) x 4 f(x))c(x, y)l x3 (x y, f(x) f(y)), T T K 6 = T x 4 f(x)( x 4 f(y) x 4 f(x))d(x, y), T B(x, y) = 3(f(x) f(y))(f(x) f(y) f(y) (x y)) [ x y + (f(x) f(y)) ] 5/, f(y) (x y) C(x, y) = [ x y + (f(x) f(y)) ], 3/ D(x, y) = L x (x y, f(x) f(y))(x y ) + L x (x y, f(x) f(y))(x y ) [ x y + (f(x) f(y)) ] 3/, L x (x, x, x 3 ) = x L(x, x, x 3 ), L x (x, x, x 3 ) = x L(x, x, x 3 ). Exchanging the variables x y we can obtain K 3. The terms K 4, K 5 K 6 can be estimated in a similar way as K. Therefore, we obtain d dt 4 x f L (t) C( f 5 H4(t) + ), analogously d dt 4 x f L (t) C( f 5 H4(t) + ). This proves local existence. The proof of the uniqueness is similar to the case Ω = R. dy 6
17 4.3 -D case Using the equation (6) in R the equation (7) in the periodic case we obtain the following theorem Theorem 4.3 Let f (x) H k for k 3 ρ > ρ. Then there exists a time T > so that there is a unique solution to (6) in C ([, T ]; H k ) with f(x, ) = f (x). The proof is similar to the theorem s Global solution for the -D stable case In this section we obtain a family of global solutions for the -D stable case with a small initial data with respect to a fixed norm. Indeed, we can get the result with an initial data with the property f H s = for s > 3/. In this section we consider x R f a = ˆf(k) e a k. For a >, if f a <, then the function f can be extended analytically on the strip Iz < a. Furthermore x f a C f b b a, (3) for b > a. The main result of this section is Theorem 5. Let f (x) be a function such that T f (x) dx =, x f ε for ε small enough xf b(t) εe b(t) ( + b(t) γ ), (3) with < γ <, b(t) = a (ρ ρ )t/, ρ > ρ a (ρ ρ )t/. Then, there exists a unique solution of (6) with f(x, ) = f (x) ρ > ρ satisfying x f a (t) C(ε) exp((σa (ρ ρ )t)/4), (3) xf a (t) C(ε)( + σa ρ ρ t γ ) exp((σa (ρ ρ )t)/4), (33) for a ρ ρ σ t, σ = + δ < δ <. The condition (3) can be satisfied for example if Λ +γ f < ε ˆf () = ˆf () = ˆf ( ) = since xf b(t) e b(t) Λ +γ f max k k γ e b(t)( k ). In order to prove the theorem, we use the Cauchy-Kowalewski method (see [5] [6]) in a similar way as Caflisch Orellana [7] Siegel, Caflisch Howison [8]. We show the proof with ρ ρ = without loss of generality. Let g(x, t) h(x, t) be functions satisfying g t = Λg, g(x, ) = f (x), h t = Λh + T (g + h), h(x, ) =, (34) 7
18 with T (f) = π R x f(x) x f(x α) α Then the function f(x, t) = g(x, t) + h(x, t) is a solution of (6). First, we show some properties of the nonlinear operator T. Lemma 5. If x f a, x g a < for a then ( f(x) f(x α) ) α + ( f(x) f(x α) α ) dα. (35) T (f)() =, (36) x T (f) a C xf a x f a, (37) x T (f) x T (g) a C ( xf a + xg a ) x f x g a + C ( x f a + x g a ) xf xg a, (38) with C = 4( x f a) C = 4( x f a) + ( x g a). Proof of the Lemma: Due to the inequality x f(x) x f a < by (35) we obtain T (f) = π n ( ) n R x f(x) x f(x α) α ( ) n f(x) f(x α) dα, (39) α T (f) = π x n ( ) n n + R ( ) n+ f(x) f(x α) dα. α Thus T (f)() =. Using (39) T (f)(k) = π n = ( ) n n n δ( R k,...,k n j= n δ( k,...,k n j= ( ) n k j, k)ik n j= ˆf(k j ) e iαkj dα α n k j, k) M n (k,..., k n ) ik ˆf(k j ), j= where M n (k,..., k n ) = π n R j= e iαkj dα. (4) α We get M n (k,..., k n ) = ( ) n m n (k,..., k n ) n j= k j, 8
19 with It follows m n (k,..., k n ) = π ds... with m n (k,..., k n ). We have therefore = π ds... = i π ds... A = ds... T (f)(k) = n n j= e a k k T (f)(k) k k ds n R ds n P V e iαk R ds n P V α ( n exp iα (s j )k j )dα j= ( n ) dα exp iα (s j )k j α R ds n (sing A sing B), (s j )k j, B = k + n δ( k,...,k n j= n n (n + ) j= ( exp iαk + iα n j= k j, k) m n (k,..., k n ) k,...,k n e a k k δ( x T (f) a xf a (n + ) x f n n n j= k j, k) s j k j. n j= n k,...,k n e a k k ˆf(k ) We get (37) for x f a <. In a similar way we obtain (38). From (34) g can be expressed as follows by the hypothesis of the initial data we have ĝ(k, t) = e k t ˆf (k), a j= n n j= k j ˆf(kj ), k j ˆf(k j ) j= (s j )k j ) dα α e a k j k j ˆf(k j ), = xf a 3 x f 3 a x f 4 a ( x f a). x g a (t) εe a t, (4) xg a (t) εe a t ( + (t a) γ ), (4) for t a. We will prove the existence of h by an induction argument on the iterative equation: t h n+ = Λh n+ + T (g + h n ), h n+ (x, ) =, h =, 9
20 or ĥ n+ (k, t) = h =. t For h we obtain the following estimates e k (t s) (T (g + h n )) b (k, s)ds, x h a (t) Using (37), (4) (4) we get By (4) (4) we have t I C t a t T (g) a+s t (s)ds = t a t + = I + I. t a t a t a I e a t e s x T (g) (s)ds Ce a t e s xg (s) x g (s)ds t a Cε e a t e s ( + s γ )ds Cε ( + γ) e a t. γ xg a+s t (s) x g a+s t (s)ds Cε e (a t) a( + (t a) γ ) Cε e a t, δ due to the inequalities (aδ) γ > (t a) γ ae a t δ for σa < t. Then Choosing b = a + s t + t a we have x h a (t) 5Cε δγ ea t. where xh a (t) t xt (g) a+s t (s)ds ( t a t t ) + = I 3 + I 4, t a x T (g) b (s) b (a + s t) ds t x T (g) b t a I 3 a t Ce t a t a Cε γ e a t ( + (t a) γ ), e s xg (s) x g (s)ds Cε e a t t a t a e s ( + s γ )ds Therefore I 4 C t t a xg b (s) x g b (s)ds Cε t a t a ea t ( + ( t a )γ )( t + a ) 3Cε e a t ( + (t a) γ ). δ
21 Define r n+ = h n+ h n, R n = M n = x h a (t) 5Cε δγ ea t, (43) xh a (t) 5Cε δγ e a t ( + (t a) γ ). (44) sup a < σa < t sup a < σa < t ( x r n a + xr n ) a + (t σa) γ e t σa, ( x h n a + xh n ) a + (t σa) γ e t σa. Take M = R 5Cε δγ ε suppose that M j, R j ε for any j =,..., n, then x r n+ a Using (38) we have t x T (g + h n ) x T (g + h n ) a+s t (s)ds = t a t + = I 7 + I 8. t a t a I 7 Ce a t e s ( x r n (s)( xg + xh n (s) + xg + xh n (s)) t a + Ce a t e s ( xr n (s)( x g + x h n (s) + x g + x h n (s))ds t a Cε R n e a t ( + s γ )ds Cε γ R ne σa t, t I 8 C + C t a t ( x r n a+s t (s)( xg + xh n a+s t (s) + xg + xh n a+s t (s))ds t a ( xr n a+s t (s)( x g + x h n a+s t (s) + x g + x h n a+s t (s))ds t Cε R n e δs σ(t a) ( + (σ(t a) δs) γ )ds t a Cε t a R n e x ( + x γ )dx 6Cε δ t σa γδ R ne σa t. We obtain for b = a + s t + σ(t a) δs σ xr n+ a (t) t t x(t (g + h n ) T (g + h n )) a+s t (s)ds x (T (g + h n ) T (g + h n ) b (s) ds b (a + s t) t x (T (g + h n ) T (g + h n ) ( σ b σ σ(t a) δs σ+ (t a) t + σ σ+ (t a) ) = I 9 + I.
22 We have σ(t a) δs > σ σ σ+ (t a) for s σ+ (t a) therefore we obtain I 9 (σ + )C t a + (σ + )C t a (σ + )Cε t a σ σ+ (t a) σ σ+ (t a) R n e b ( x r n (s)( xg + xh n (s) + xg + xh n (s)) σ σ+ (t a) Using (38) the induction hypothesis we get e b ( xr n (s)( x g + x h n (s) + x g + x h n (s))ds t I 4σCε R n e σb s ( + (s σb)γ ) ds σ σ+ (t a) σ(t a) δs t 4σCε R n e δs σ(t a) + ( σ(t a) δs ) γ σ σ+ (t a) σ(t a) δs σ σ+ (t a) R n 4σCε δ t σa Due to the estimates for I 7, I 8, I 9 I we obtain e b s ( + s γ )ds 4σCε R n e a t ( + (t a) γ ). γ ds e x (x + x γ )dx 8σCε δ( γ) R ne σa t ( + (σa t) γ ). Choosing ε small enough we get R n+ Cσε δγ(γ ) R n. (45) R n+ R n... n R ε n+, n+ M n+ R n+ ε. j= Therefore, we obtain the function h = lim n hn satisfying x h a (t) n R n e σa t ε e σa t. Taking f(x, t) = g(x, t) + h(x, t), we get (3) for ρ ρ =. In order to show the uniqueness, we write the equation (3) for ρ ρ = in the following form f t = Λf + T (f), f(x, ) = f (x). Suppose that there exist two solutions f f with f (x, ) = f (x, ). Define R by It follows that R R = sup a < σa < t ( x f x f a + xf xf ) a + (t σa) γ e t σa. C(ε)σ C(ε)σ δγ(γ ) R for ε small enough it yields δγ(γ ) < therefore f = f.
23 6 Ill-posedness for the unstable case Here we show ill-posedness for the unstable case ρ > ρ. We use the global solution for the -D stable case f(x, t) satisfying (3) with Λ +γ f < C Λ +γ+ζ f = for γ, ζ >. Making a change of variables, we define f λ (x, t) = λ f(λx, λt + λ / ) obtaining {f λ } λ> a family of solutions to the unstable case. Using (3) follows f λ H s() = λ s 3 f H s(λ / ) C λ s 3 f (λ / ) C λ s 3 e ρ ρ 4 λ /, f λ H s(λ / ) = λ s 3 f H s() λ s 3 C k +γ+ζ ˆf (k) =, for s > 3/ γ, ζ small enough. We obtain an ill posed problem for s > 3/. Theorem 6. Let s > 3/, then for any ε > there exists a solution f of (6) with ρ > ρ < δ < ε such that f H s() ε f H s(δ) =. Remark 6. If one considers a solution of the 3-D problem satisfying f(x, x, t) = f(x, t), from the equation (3) one obtains a solution of (6). This shows that solutions of the -D case are solutions of the 3-D problem therefore, using the above theorem, one obtains ill-posedness for the 3-D case with ρ > ρ. k References [] D. Ambrose. Well-posedness of Two-phase Hele-Shaw Flow without Surface Tension. Euro. Jnl. of Applied Mathematics , 4. [] D. Ambrose N. Masmoudi. The zero surface tension limit of two-dimensional water waves. Comm. Pure Appl. Math , 5. [3] J. Bear, Dynamics of Fluids in Porous Media, American Elsevier, New York, 97. [4] A. L. Bertozzi P. Constantin. Global regularity for vortex patches. Comm. Math. Phys. 5 (): 9 8, 993. [5] A. L. Bertozzi A. J. Majda. Vorticity the Mathematical Theory of Incompresible Fluid Flow. Cambridge Press,. [6] G. Birkhoff. Helmholtz Taylor instability. Hydrodynamics Instability, Proc. Symp. Appl. Math. XII A.M.S., 55 76, 96. [7] R. Caflisch O. Orellana. Singular solutions ill-posedness for the evolution of vortex sheets. SIAM J. Math. Anal. (): 93 37, 989. [8] J.Y. Chemin. Persistence of geometric structures in two-dimensional incompressible fluids. Ann. Sci. Ecole. Norm. Sup. 6 (4): 57 54, 993. [9] P. Constantin, T.F. Dupont, R.E. Goldstein, L.P. Kadanoff, M.J. Shelley S.M. Zhou. Droplet breakup in a model of the Hele-Shaw cell. Physical Review E, 47, , 993. [] P. Constantin, A. J. Majda, E. Tabak. Formation of strong fronts in the -D quasigeostrophic thermal active scalar. Nonlinearity, 7: ,
24 [] J. Escher G. Simonett. Classical solutions for Hele-Shaw models with surface tension. Adv. Differential Equations, :69-64, 997. [] Hele-Shaw. Nature 58, 34, 898. [3] T.Y. Hou, J.S. Lowengrub M.J. Shelley. Removing the Stiffness from Interfacial Flows with Surface Tension. J. Comput. Phys., 4: 3-338, 994. [4] M. Muskat. The flow of homogeneous fluids through porous media. New York, 937. [5] L. Nirenberg. An abstract form of the nonlinear Cauchy-Kowalewski theorem. J. Differential Geometry, , 97. [6] T. Nishida. A note on a theorem of Nirenberg. J. Differential Geometry, , 977. [7] P.G. Saffman Taylor. The penetration of a fluid into a porous medium or Hele-Shaw cell containing a more viscous liquid. Proc. R. Soc. London, Ser. A 45, 3-39, 958. [8] M. Siegel, R. Caflisch S. Howison. Global Existence, Singular Solutions, Ill-Posedness for the Muskat Problem. Comm. Pure Appl. Math., 57: 374-4, 4. [9] J.L. Rodrigo. On the Evolution of Sharp Fronts for the Quasi-Geostrophic Equation. Comm. Pure Appl. Math., 58: 8-866, 5. [] E. Stein G. Weiss. Introduction to Fourier Analysis on Euclidean spaces. Princeton University Press. Princeton, NJ, 97. [] E. Stein. Harmonic Analysis. Princeton University Press. Princeton, NJ, 993. [] G. Taylor. The instability of liquid surfaces when accelerated in a direction perpendicular to their planes. I. Proc. Roy. Soc. London. Ser. A., [3] S. Wu. Well-posedness in Sobolev spaces of the full water wave problem in -D. Invent. math. 3, [4] S. Wu. Well-posedness in Sobolev spaces of the full water wave problem in 3-D. J. Amer. Math. Soc.,
Contour dynamics of incompressible 3-D fluids in a porous medium with different densities
Contour dynamics of incompressible 3-D fluids in a porous medium with different densities Diego Córdoba Francisco Gancedo September 4, 6 Abstract We consider the problem of the evolution of the interface
More informationThe Rayleigh-Taylor condition for the evolution of irrotational fluid interfaces.
The Rayleigh-Taylor condition for the evolution of irrotational fluid interfaces. Antonio Cordoba Universidad Autonoma de Madrid ctra. de Colmenar Viejo, Km. 15 28049 Madrid e-mail: antonio.cordoba@uam.es
More informationDiego Córdoba. Interface dynamics for incompressible flows in 2D
Interface dynamics for incompressible flows in 2D Diego Córdoba Equation ρ t + u ρ = 0, u = 0, { ρ ρ(x 1, x 2, t) = 1, x Ω 1 (t) ρ 2, x Ω 2 (t), with ρ(x, t) an active scalar, (x, t) R 2 R +, ρ 1 ρ 2 are
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 informationBlow up of solutions for a 1D transport equation with nonlocal velocity and supercritical dissipation
Blow up of solutions for a 1D transport equation with nonlocal velocity and supercritical dissipation Dong Li a,1 a School of Mathematics, Institute for Advanced Study, Einstein Drive, Princeton, NJ 854,
More informationA non-stiff boundary integral method for 3D porous media flow with surface tension
A non-stiff boundary integral method for 3D porous media flow with surface tension D. M. Ambrose 1 and M. Siegel 1 Department of Mathematics, Drexel University, Philadelphia, PA 191 Department of Mathematical
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 informationJournées ÉQUATIONS AUX DÉRIVÉES PARTIELLES. Port d Albret, 7 juin 11 juin Angel Castro, Diego Córdoba and Francisco Gancedo
Journées ÉQUATIONS AUX DÉIVÉES PATIELLES Port d Albret, 7 juin 11 juin 010 Angel Castro, Diego Córdoba and Francisco Gancedo Some recent results on the Muskat problem J. É. D. P. (010), Eposé n o V, 14
More informationANALYTICAL BEHAVIOR OF 2-D INCOMPRESSIBLE FLOW IN POROUS MEDIA
ANALYTICAL BEHAVIOR OF -D INCOMPRESSIBLE FLOW IN POROUS MEDIA DIEGO CÓRDOBA, FRANCISCO GANCEDO AND RAFAEL ORIVE Abstract In this paper we study the analytic structure of a two-dimensional mass balance
More informationfor all subintervals I J. If the same is true for the dyadic subintervals I D J only, we will write ϕ BMO d (J). In fact, the following is true
3 ohn Nirenberg inequality, Part I A function ϕ L () belongs to the space BMO() if sup ϕ(s) ϕ I I I < for all subintervals I If the same is true for the dyadic subintervals I D only, we will write ϕ BMO
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 informationComputing Time-Periodic Solutions of Nonlinear Systems of Partial Differential Equations
1 Computing ime-periodic Solutions of Nonlinear Systems of Partial Differential Equations David M. Ambrose Department of Mathematics, Drexel University, Philadelphia, PA 191 USA E-mail: ambrose@math.drexel.edu
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 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 Global Well-Posedness of the Lagrangian Averaged Euler Equations
On Global Well-Posedness of the Lagrangian Averaged Euler Equations Thomas Y. Hou Congming Li Abstract We study the global well-posedness of the Lagrangian averaged Euler equations in three dimensions.
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 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 informationMy doctoral research concentrated on elliptic quasilinear partial differential equations of the form
1 Introduction I am interested in applied analysis and computation of non-linear partial differential equations, primarily in systems that are motivated by the physics of fluid dynamics. My Ph.D. thesis
More informationContents (I): The model: The physical scenario and the derivation of the model: Muskat and the confined Hele-Shaw cell problems.
Contents (I): The model: The physical scenario and the derivation of the model: Muskat and the confined Hele-Shaw cell problems. First ideas and results: Existence and uniqueness But, are the boundaries
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 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 informationREGULARITY OF GENERALIZED NAVEIR-STOKES EQUATIONS IN TERMS OF DIRECTION OF THE VELOCITY
Electronic Journal of Differential Equations, Vol. 00(00), No. 05, pp. 5. ISSN: 07-669. UR: http://ejde.math.txstate.edu or http://ejde.math.unt.edu ftp ejde.math.txstate.edu REGUARITY OF GENERAIZED NAVEIR-STOKES
More informationThe 3D Euler and 2D surface quasi-geostrophic equations
The 3D Euler and 2D surface quasi-geostrophic equations organized by Peter Constantin, Diego Cordoba, and Jiahong Wu Workshop Summary This workshop focused on the fundamental problem of whether classical
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 informationEvolution of sharp fronts for the surface quasi-geostrophic equation
Evolution of sharp fronts for the surface quasi-geostrophic equation José Luis Rodrigo University of Warwick supported by ERC Consolidator 616797 UIMP - Frontiers of Mathematics IV July 23rd 2015 José
More informationarxiv: v1 [math.ap] 18 May 2017
Littlewood-Paley-Stein functions for Schrödinger operators arxiv:175.6794v1 [math.ap] 18 May 217 El Maati Ouhabaz Dedicated to the memory of Abdelghani Bellouquid (2/2/1966 8/31/215) Abstract We study
More informationThe 2D Magnetohydrodynamic Equations with Partial Dissipation. Oklahoma State University
The 2D Magnetohydrodynamic Equations with Partial Dissipation Jiahong Wu Oklahoma State University IPAM Workshop Mathematical Analysis of Turbulence IPAM, UCLA, September 29-October 3, 2014 1 / 112 Outline
More 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 information4 Riesz Kernels. Since the functions k i (ξ) = ξ i. are bounded functions it is clear that R
4 Riesz Kernels. A natural generalization of the Hilbert transform to higher dimension is mutiplication of the Fourier Transform by homogeneous functions of degree 0, the simplest ones being R i f(ξ) =
More informationarxiv:math/ v1 [math.ap] 28 Oct 2005
arxiv:math/050643v [math.ap] 28 Oct 2005 A remark on asymptotic completeness for the critical nonlinear Klein-Gordon equation Hans Lindblad and Avy Soffer University of California at San Diego and Rutgers
More informationV (r,t) = i ˆ u( x, y,z,t) + ˆ j v( x, y,z,t) + k ˆ w( x, y, z,t)
IV. DIFFERENTIAL RELATIONS FOR A FLUID PARTICLE This chapter presents the development and application of the basic differential equations of fluid motion. Simplifications in the general equations and common
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 informationg t + Λ α g + θg x + g 2 = 0.
NONLINEAR MAXIMUM PRINCIPLES FOR DISSIPATIVE LINEAR NONLOCAL OPERATORS AND APPLICATIONS PETER CONSTANTIN AND VLAD VICOL ABSTRACT. We obtain a family of nonlinear maximum principles for linear dissipative
More informationOn the displacement of two immiscible Stokes fluids in a 3D Hele-Shaw cell
On the displacement of two immiscible Stokes fluids in a 3D Hele-Shaw cell Gelu Paşa Abstract. In this paper we study the linear stability of the displacement of two incompressible Stokes fluids in a 3D
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 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 informationFree boundaries in fractional filtration equations
Free boundaries in fractional filtration equations Fernando Quirós Universidad Autónoma de Madrid Joint work with Arturo de Pablo, Ana Rodríguez and Juan Luis Vázquez International Conference on Free Boundary
More informationExponential Energy Decay for the Kadomtsev-Petviashvili (KP-II) equation
São Paulo Journal of Mathematical Sciences 5, (11), 135 148 Exponential Energy Decay for the Kadomtsev-Petviashvili (KP-II) equation Diogo A. Gomes Department of Mathematics, CAMGSD, IST 149 1 Av. Rovisco
More informationFundamentals of Fluid Dynamics: Ideal Flow Theory & Basic Aerodynamics
Fundamentals of Fluid Dynamics: Ideal Flow Theory & Basic Aerodynamics Introductory Course on Multiphysics Modelling TOMASZ G. ZIELIŃSKI (after: D.J. ACHESON s Elementary Fluid Dynamics ) bluebox.ippt.pan.pl/
More informationINSTITUTE OF MATHEMATICS THE CZECH ACADEMY OF SCIENCES. Note on the fast decay property of steady Navier-Stokes flows in the whole space
INSTITUTE OF MATHEMATICS THE CZECH ACADEMY OF SCIENCES Note on the fast decay property of stea Navier-Stokes flows in the whole space Tomoyuki Nakatsuka Preprint No. 15-017 PRAHA 017 Note on the fast
More informationarxiv: v2 [math.ap] 3 Oct 2011
Splash singularity for water waves arxiv:1106.2120v2 [math.ap] 3 Oct 2011 Angel Castro, Diego Córdoba, Charles Fefferman, Francisco Gancedo and Javier Gómez-Serrano Abstract We exhibit smooth initial data
More information1. Nonlinear Equations. This lecture note excerpted parts from Michael Heath and Max Gunzburger. f(x) = 0
Numerical Analysis 1 1. Nonlinear Equations This lecture note excerpted parts from Michael Heath and Max Gunzburger. Given function f, we seek value x for which where f : D R n R n is nonlinear. f(x) =
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 informationFinite-time singularity formation for Euler vortex sheet
Finite-time singularity formation for Euler vortex sheet Daniel Coutand Maxwell Institute Heriot-Watt University Oxbridge PDE conference, 20-21 March 2017 Free surface Euler equations Picture n N x Ω Γ=
More informationREGULARITY CRITERIA FOR WEAK SOLUTIONS TO 3D INCOMPRESSIBLE MHD EQUATIONS WITH HALL TERM
Electronic Journal of Differential Equations, Vol. 2018 (2018), No. 10, pp. 1 12. ISSN: 1072-6691. UL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu EGULAITY CITEIA FO WEAK SOLUTIONS TO D INCOMPESSIBLE
More informationGlobal well-posedness of the primitive equations of oceanic and atmospheric dynamics
Global well-posedness of the primitive equations of oceanic and atmospheric dynamics Jinkai Li Department of Mathematics The Chinese University of Hong Kong Dynamics of Small Scales in Fluids ICERM, Feb
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 informationGlobal well-posedness and decay for the viscous surface wave problem without surface tension
Global well-posedness and decay for the viscous surface wave problem without surface tension Ian Tice (joint work with Yan Guo) Université Paris-Est Créteil Laboratoire d Analyse et de Mathématiques Appliquées
More informationRandom Variables. Random variables. A numerically valued map X of an outcome ω from a sample space Ω to the real line R
In probabilistic models, a random variable is a variable whose possible values are numerical outcomes of a random phenomenon. As a function or a map, it maps from an element (or an outcome) of a sample
More informationFinite time singularity for the modified SQG patch equation
Finite time singularity for the modified SQG patch equation Alexander Kiselev Lenya Ryzhik Yao Yao Andrej Zlatoš September 4, 2015 Abstract It is well known that the incompressible Euler equations in two
More informationA regularity criterion for the 3D NSE in a local version of the space of functions of bounded mean oscillations
Ann. I. H. Poincaré AN 27 (2010) 773 778 www.elsevier.com/locate/anihpc A regularity criterion for the 3D NSE in a local version of the space of functions of bounded mean oscillations Zoran Grujić a,,
More informationin Bounded Domains Ariane Trescases CMLA, ENS Cachan
CMLA, ENS Cachan Joint work with Yan GUO, Chanwoo KIM and Daniela TONON International Conference on Nonlinear Analysis: Boundary Phenomena for Evolutionnary PDE Academia Sinica December 21, 214 Outline
More informationVortex stretching and anisotropic diffusion in the 3D NSE
University of Virginia SIAM APDE 2013 ; Lake Buena Vista, December 7 3D Navier-Stokes equations (NSE) describing a flow of 3D incompressible viscous fluid read u t + (u )u = p + ν u, supplemented with
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 informationThe Navier Stokes Equations for Incompressible Flows: Solution Properties at Potential Blow Up Times
The Navier Stokes Equations for Incompressible Flows: Solution Properties at Potential Blow Up Times Jens Lorenz Department of Mathematics and Statistics, UNM, Albuquerque, NM 873 Paulo Zingano Dept. De
More informationWELL-POSEDNESS OF THE TWO-DIMENSIONAL GENERALIZED BENJAMIN-BONA-MAHONY EQUATION ON THE UPPER HALF PLANE
WELL-POSEDNESS OF THE TWO-DIMENSIONAL GENERALIZED BENJAMIN-BONA-MAHONY EQUATION ON THE UPPER HALF PLANE YING-CHIEH LIN, C. H. ARTHUR CHENG, JOHN M. HONG, JIAHONG WU, AND JUAN-MING YUAN Abstract. This paper
More informationMAE210C: Fluid Mechanics III Spring Quarter sgls/mae210c 2013/ Solution II
MAE210C: Fluid Mechanics III Spring Quarter 2013 http://web.eng.ucsd.edu/ sgls/mae210c 2013/ Solution II D 4.1 The equations are exactly the same as before, with the difference that the pressure in the
More informationGlobal well-posedness for the critical 2D dissipative quasi-geostrophic equation
Invent. math. 167, 445 453 (7) DOI: 1.17/s-6--3 Global well-posedness for the critical D dissipative quasi-geostrophic equation A. Kiselev 1, F. Nazarov, A. Volberg 1 Department of Mathematics, University
More informationII. FOURIER TRANSFORM ON L 1 (R)
II. FOURIER TRANSFORM ON L 1 (R) In this chapter we will discuss the Fourier transform of Lebesgue integrable functions defined on R. To fix the notation, we denote L 1 (R) = {f : R C f(t) dt < }. The
More informationProceedings of the 5th International Conference on Inverse Problems in Engineering: Theory and Practice, Cambridge, UK, 11-15th July 2005
Proceedings of the 5th International Conference on Inverse Problems in Engineering: Theory and Practice, Cambridge, UK, 11-15th July 2005 SOME INVERSE SCATTERING PROBLEMS FOR TWO-DIMENSIONAL SCHRÖDINGER
More informationNonlinear Analysis. A regularity criterion for the 3D magneto-micropolar fluid equations in Triebel Lizorkin spaces
Nonlinear Analysis 74 (11) 5 Contents lists available at ScienceDirect Nonlinear Analysis journal homepage: www.elsevier.com/locate/na A regularity criterion for the 3D magneto-micropolar fluid equations
More informationGLOBAL REGULARITY OF LOGARITHMICALLY SUPERCRITICAL 3-D LAMHD-ALPHA SYSTEM WITH ZERO DIFFUSION
GLOBAL REGULARITY OF LOGARITHMICALLY SUPERCRITICAL 3-D LAMHD-ALPHA SYSTEM WITH ZERO DIFFUSION KAZUO YAMAZAKI Abstract. We study the three-dimensional Lagrangian-averaged magnetohydrodynamicsalpha system
More informationExistence and Continuation for Euler Equations
Existence and Continuation for Euler Equations David Driver Zhuo Min Harold Lim 22 March, 214 Contents 1 Introduction 1 1.1 Physical Background................................... 1 1.2 Notation and Conventions................................
More informationHOMEOMORPHISMS OF BOUNDED VARIATION
HOMEOMORPHISMS OF BOUNDED VARIATION STANISLAV HENCL, PEKKA KOSKELA AND JANI ONNINEN Abstract. We show that the inverse of a planar homeomorphism of bounded variation is also of bounded variation. In higher
More informationMultivariable Calculus
2 Multivariable Calculus 2.1 Limits and Continuity Problem 2.1.1 (Fa94) Let the function f : R n R n satisfy the following two conditions: (i) f (K ) is compact whenever K is a compact subset of R n. (ii)
More informationRegularity and compactness for the DiPerna Lions flow
Regularity and compactness for the DiPerna Lions flow Gianluca Crippa 1 and Camillo De Lellis 2 1 Scuola Normale Superiore, Piazza dei Cavalieri 7, 56126 Pisa, Italy. g.crippa@sns.it 2 Institut für Mathematik,
More informationA small-scale decomposition for 3D boundary integral computations with surface tension
A small-scale decomposition for 3D boundary integral computations with surface tension D. M. Ambrose b, M. Siegel a,, S. Tlupova c a Department of Mathematical Sciences and Center for Applied Mathematics
More informationOPTIMAL CONVERGENCE RATES FOR THE COMPRESSIBLE NAVIER-STOKES EQUATIONS WITH POTENTIAL FORCES
OPTIMAL CONVERGENCE RATES FOR THE COMPRESSIBLE NAVIER-STOKES EQUATIONS WITH POTENTIAL FORCES RENJUN DUAN Department of Mathematics, City University of Hong Kong 83 Tat Chee Avenue, Kowloon, Hong Kong,
More informationKelvin Helmholtz Instability
Kelvin Helmholtz Instability A. Salih Department of Aerospace Engineering Indian Institute of Space Science and Technology, Thiruvananthapuram November 00 One of the most well known instabilities in fluid
More informationarxiv: v1 [math.ap] 3 Nov 2017
BLOW UP OF A HYPERBOLIC SQG MODEL HANG YANG arxiv:7.255v [math.ap] 3 Nov 27 Abstract. This paper studies of a variation of the hyperbolic blow up scenario suggested by Hou and Luo s recent numerical simulation
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 informationFormulation of the problem
TOPICAL PROBLEMS OF FLUID MECHANICS DOI: https://doi.org/.43/tpfm.27. NOTE ON THE PROBLEM OF DISSIPATIVE MEASURE-VALUED SOLUTIONS TO THE COMPRESSIBLE NON-NEWTONIAN SYSTEM H. Al Baba, 2, M. Caggio, B. Ducomet
More informationImplicit Functions, Curves and Surfaces
Chapter 11 Implicit Functions, Curves and Surfaces 11.1 Implicit Function Theorem Motivation. In many problems, objects or quantities of interest can only be described indirectly or implicitly. It is then
More informationRegularity and Decay Estimates of the Navier-Stokes Equations
Regularity and Decay Estimates of the Navier-Stokes Equations Hantaek Bae Ulsan National Institute of Science and Technology (UNIST), Korea Recent Advances in Hydrodynamics, 216.6.9 Joint work with Eitan
More informationStarting from Heat Equation
Department of Applied Mathematics National Chiao Tung University Hsin-Chu 30010, TAIWAN 20th August 2009 Analytical Theory of Heat The differential equations of the propagation of heat express the most
More informationSobolev Spaces. Chapter 10
Chapter 1 Sobolev Spaces We now define spaces H 1,p (R n ), known as Sobolev spaces. For u to belong to H 1,p (R n ), we require that u L p (R n ) and that u have weak derivatives of first order in L p
More informationTEST CODE: MMA (Objective type) 2015 SYLLABUS
TEST CODE: MMA (Objective type) 2015 SYLLABUS Analytical Reasoning Algebra Arithmetic, geometric and harmonic progression. Continued fractions. Elementary combinatorics: Permutations and combinations,
More information1 Assignment 1: Nonlinear dynamics (due September
Assignment : Nonlinear dynamics (due September 4, 28). Consider the ordinary differential equation du/dt = cos(u). Sketch the equilibria and indicate by arrows the increase or decrease of the solutions.
More informationAn Eulerian-Lagrangian Approach for Incompressible Fluids: Local Theory
An Eulerian-Lagrangian Approach for Incompressible Fluids: Local Theory Peter Constantin Department of Mathematics The University of Chicago December 22, 2000 Abstract We study a formulation of the incompressible
More informationScientific Computing
Lecture on Scientific Computing Dr. Kersten Schmidt Lecture 4 Technische Universität Berlin Institut für Mathematik Wintersemester 2014/2015 Syllabus Linear Regression Fast Fourier transform Modelling
More informationMATH 6337 Second Midterm April 1, 2014
You can use your book and notes. No laptop or wireless devices allowed. Write clearly and try to make your arguments as linear and simple as possible. The complete solution of one exercise will be considered
More informationChapter 1 Mathematical Foundations
Computational Electromagnetics; Chapter 1 1 Chapter 1 Mathematical Foundations 1.1 Maxwell s Equations Electromagnetic phenomena can be described by the electric field E, the electric induction D, the
More informationThe oblique derivative problem for general elliptic systems in Lipschitz domains
M. MITREA The oblique derivative problem for general elliptic systems in Lipschitz domains Let M be a smooth, oriented, connected, compact, boundaryless manifold of real dimension m, and let T M and T
More informationThe incompressible Navier-Stokes equations in vacuum
The incompressible, Université Paris-Est Créteil Piotr Bogus law Mucha, Warsaw University Journées Jeunes EDPistes 218, Institut Elie Cartan, Université de Lorraine March 23th, 218 Incompressible Navier-Stokes
More informationThe h-principle for the Euler equations
The h-principle for the Euler equations László Székelyhidi University of Leipzig IPAM Workshop on Turbulent Dissipation, Mixing and Predictability 09 January 2017 Incompressible Euler equations @ t v +
More informationBLOWUP THEORY FOR THE CRITICAL NONLINEAR SCHRÖDINGER EQUATIONS REVISITED
BLOWUP THEORY FOR THE CRITICAL NONLINEAR SCHRÖDINGER EQUATIONS REVISITED TAOUFIK HMIDI AND SAHBI KERAANI Abstract. In this note we prove a refined version of compactness lemma adapted to the blowup analysis
More informationarxiv: v1 [math.ap] 16 May 2007
arxiv:0705.446v1 [math.ap] 16 May 007 Regularity criterion for 3D Navier-Stokes equations in terms of the direction of the velocity Alexis Vasseur October 3, 018 Abstract In this short note, we give a
More informationExplosive Solution of the Nonlinear Equation of a Parabolic Type
Int. J. Contemp. Math. Sciences, Vol. 4, 2009, no. 5, 233-239 Explosive Solution of the Nonlinear Equation of a Parabolic Type T. S. Hajiev Institute of Mathematics and Mechanics, Acad. of Sciences Baku,
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 informationDecay in Time of Incompressible Flows
J. math. fluid mech. 5 (23) 231 244 1422-6928/3/3231-14 c 23 Birkhäuser Verlag, Basel DOI 1.17/s21-3-79-1 Journal of Mathematical Fluid Mechanics Decay in Time of Incompressible Flows Heinz-Otto Kreiss,
More informationOn Some Estimates of the Remainder in Taylor s Formula
Journal of Mathematical Analysis and Applications 263, 246 263 (2) doi:.6/jmaa.2.7622, available online at http://www.idealibrary.com on On Some Estimates of the Remainder in Taylor s Formula G. A. Anastassiou
More informationMATH 220 solution to homework 4
MATH 22 solution to homework 4 Problem. Define v(t, x) = u(t, x + bt), then v t (t, x) a(x + u bt) 2 (t, x) =, t >, x R, x2 v(, x) = f(x). It suffices to show that v(t, x) F = max y R f(y). We consider
More informationGood reasons to study Euler equation.
Good reasons to study Euler equation. Applications often correspond to very large Reynolds number =ratio between the strenght of the non linear effects and the strenght of the linear viscous effects. R
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 informationESTIMATES FOR ELLIPTIC HOMOGENIZATION PROBLEMS IN NONSMOOTH DOMAINS. Zhongwei Shen
W,p ESTIMATES FOR ELLIPTIC HOMOGENIZATION PROBLEMS IN NONSMOOTH DOMAINS Zhongwei Shen Abstract. Let L = div`a` x, > be a family of second order elliptic operators with real, symmetric coefficients on a
More informationLagrangian analyticity for vortex patches
Lagrangian analyticity for vortex patches Matthew Hernandez It is a known result that the incompressible Euler equations have the property that by imposing only a small amount of regularity in the initial
More informationA Product Property of Sobolev Spaces with Application to Elliptic Estimates
Rend. Sem. Mat. Univ. Padova Manoscritto in corso di stampa pervenuto il 23 luglio 2012 accettato l 1 ottobre 2012 A Product Property of Sobolev Spaces with Application to Elliptic Estimates by Henry C.
More informationMaxima and Minima. (a, b) of R if
Maxima and Minima Definition Let R be any region on the xy-plane, a function f (x, y) attains its absolute or global, maximum value M on R at the point (a, b) of R if (i) f (x, y) M for all points (x,
More informationNonlinear stabilization via a linear observability
via a linear observability Kaïs Ammari Department of Mathematics University of Monastir Joint work with Fathia Alabau-Boussouira Collocated feedback stabilization Outline 1 Introduction and main result
More informationOne-sided operators in grand variable exponent Lebesgue spaces
One-sided operators in grand variable exponent Lebesgue spaces ALEXANDER MESKHI A. Razmadze Mathematical Institute of I. Javakhishvili Tbilisi State University, Georgia Porto, June 10, 2015 One-sided operators
More information