Blow-up or No Blow-up? Fluid Dynamic Perspective of the Clay Millennium Problem

Size: px
Start display at page:

Download "Blow-up or No Blow-up? Fluid Dynamic Perspective of the Clay Millennium Problem"

Transcription

1 Blow-up or No Blow-up? Fluid Dynamic Perspective of the Clay Millennium Problem Thomas Y. Hou Applied and Comput. Mathematics, Caltech Research was supported by NSF Theodore Y. Wu Lecture in Aerospace, GALCIT, Caltech October 7, 2011

2 A tribute to Prof. Ted Wu I feel extremely honored to give this lecture in honor of Prof. Ted Wu. Ted is a great scientist, a powerful applied mathematician, an excellent teacher/mentor, and a real gentleman. His research covers a wide range of areas, and has made fundamental and highly original contributions in every area that he worked on. I very much admire his passion, his persistence and his drive to achieve scientific excellence. It is people like Prof. Ted Wu who have made Caltech such a special and unique university in the world. He is a role model for many of us. He has always been very generous in sharing his knowledge and wisdom with others, nurturing and supporting several generations of Ph.D. students, postdocs, and junior scientists.

3 Motivation The question of whether the 3D incompressible Navier-Stokes equations can develop a finite time singularity from smooth initial data is one of the seven Clay Millennium Problems. u t + (u )u = p + ν u, u = 0, (1) with initial condition u(x, 0) = u 0. Define vorticity ω = u, then ω is governed by ω t + (u )ω = u ω + ν ω. (2) Note that u is related to ω by a Riesz operator K of degree zero: u = K(ω), and we have ω L p u L p C ω L p for 1 < p <. Thus the vortex stretching term u ω is formally of the order ω 2. So far, most regularity analysis uses energy estimates and treats the nonlinear terms as a small perturbation to the diffusion term. The global regularity can be obtained only for small data.

4 A brief review Global existence for small data (Leray, Ladyzhenskaya, Kato, etc). If u 0 L p (p 3) or u 0 L 2 u 0 L 2 is small, then the 3D Navier-Stokes equations have a globally smooth solution. Non-blowup criteria due to G. Prodi 59, J. Serrin 63. A weak solution u of the 3D Navier-Stokes equations is smooth on [0, T ] R 3 provided that u L q t Lp x ([0,T ] R 3 ) < for some p, q satisfying 3 p + 2 q 1 with 3 < p and 2 q <. The critical case of p = 3, q = was proved by L. Escauriaza, G. Seregin, and V. Sverak in Partial regularity theory (Caffarelli-Kohn-Nirenberg 82, F. Lin 98) For any suitable weak solution of the 3D Navier-Stokes equations on an open set in space-time, the one-dimensional Hausdorff measure of the associated singular set is zero.

5 Review continud (Beale-Kato-Majda criterion, 1984) u ceases to be classical at T if and only if T 0 ω (t) dt =. Geometry regularity of direction field of ω: Constantin, Fefferman and Majda (1996). Let ω = ω ξ. The solution of 3D Euler is smooth solution up to T if (1) u L (Ω t) is bounded; and (2) t 0 ξ 2 L (Ω τ ) dτ is uniformly bounded for all t < T. Localized non-blow-up criteria, Deng-Hou-Yu, Denote by L(t) the arclength of a vortex line segment L t around the maximum vorticity. If (i) max Lt ( u ξ + u n ) C U (T t) A with A < 1, (ii) C L (T t) B L(t) C 0 / max Lt ( κ, ξ ) with B < 1 A, then the 3D Euler remains smooth up to T.

6 The role of convection and 3D Euler Equations Due to the incompressibility condition, the convection term does not contribute to the energy norm of velocity or the L p -norm of ω. ω t + (u )ω = u ω + ν ω. The convection term and the vortex stretching term can be reformulated as a commutator: ω t + (u )ω (ω )u = ν ω. When we consider the two terms together, we preserve the Lagrangian structure of the solution of the 3D Euler equations: ω(x (α, t), t) = X α (α, t)ω 0(α), det( X (α, t)) 1 α where X (α, t) is the flow map: X t = u(x, t), X (α, 0) = α. Convection tends to severely deform and flatten the support of maximum vorticity. Such deformation tends to weaken the nonlinearity of vortex stretching dynamicallty.

7 Numerical evidence of Euler singularity In 1993 (and 2005), R. Kerr [Phys. Fluids] presented numerical evidence of 3D Euler singularity for two anti-parallel vortex tubes: Pseudo-spectral in x and y, Chebyshev in z direction; Best resolution: ; Predicted singularity time T = 18.7, but his numerical solutions became under-resolved after T 1 = 17. ω L (T t) 1 ; u L (T t) 1/2 ; Anisotropic scaling: (T t) T t T t; Vortex lines: relatively straight, ξ (T t) 1/2 ; It is important to emphasize that T T 1 = 1.7 is not asymptotically small. Thus the extrapolated blow-up rate of ω L (T t) 1 is not valid asymptotically from T 1 to T. This makes the claim T T 1 ω dt = questionable. Kerr s computation falls into the critical case of Deng-Hou-Yu s non-blowup criterion with A = B = 1/2.

8 Computation of Hou and Li, J. Nonlinear Science, 2006 Figure : Two slightly perturbed antiparallel vortex tubes at t=0 and t=6

9 Figure : The local 3D vortex structures and vortex lines around the maximum vorticity at t = 17 with resolution

10 Resolution study for energy spectra: vs FS : the Fourier smoothing method 10 0 energy spectra comparison dashed:1024x768x2048, 2/3rd dealiasing dash dotted:1024x768x2048, FS solid:1536x1024x3072, FS

11 Maximum velocity in time 0.55 t [0,19], t [0,19], t [10,19], Figure : Maximum velocity u in time using different resolutions. With maximum velocity being bounded, the non-blowup criterion of Deng-Hou-Yu applies with A = 0 and B = 1/2, implying no blowup at least up to T = 19.

12 Dynamic depletion of vortex stretching 35 ξ u ω 30 c 1 ω log( ω ) c 2 ω Figure : Study of the vortex stretching term in time, resolution The fact ξ u ω c 1 ω log ω plus ω = ξ u ω implies ω bounded by doubly exponential. D Dt

13 Log log plot of maximum vorticity in time Figure : The plot of log log ω vs time, resolution

14 Log log plot of peak vorticity in time from Kerr s 93 paper 1 Double logarithm of peak vorticity in time from Kerr 93 paper time

15 The Effect of Convection Consider the 3D axi-symmetric incompressible Navier-Stokes equations ( ut θ + u r ur θ + u z uz θ = ν 2 1 ) r 2 u θ 1 r ur u θ, (3) ( ωt θ + u r ωr θ + u z ωz θ = ν 2 1 ) r 2 ω θ + 1 ( (u θ ) 2) r + 1 z r ur ω θ, (4) ( 2 1 r 2 ) ψ θ = ω θ, (5) where u θ, ω θ and ψ θ are the angular components of the velocity, vorticity and stream function respectively, and u r = (ψ θ ) z u z = 1 r (rψθ ) r. Note that equations (??)-(??) completely determine the evolution of the 3D axisymmetric Navier-Stokes equations.

16 Reformulation of axisymmetric Navier-Stokes equations In [CPAM 08], Hou and Li introduced the following new variables: u 1 = u θ /r, ω 1 = ω θ /r, ψ 1 = ψ θ /r, (6) and derived the following equivalent system that governs the dynamics of u 1, ω 1 and ψ 1 as follows: t u 1 + u r r u 1 + u z z u 1 = ν ( r r ) r + z 2 u1 + 2u 1 ψ 1z, ) z, t ω 1 + u r r ω 1 + u z z ω 1 = ν ( r r ) r + z 2 ω1 + ( u1 2 ( r r ) r + z 2 ψ1 = ω 1, where u r = rψ 1z, u z = 2ψ 1 + rψ 1r. Liu and Wang [SINUM07] showed that if u is a smooth velocity field, then u θ, ω θ and ψ θ must satisfy: u θ r=0 = ω θ r=0 = ψ θ r=0 = 0. Thus u 1, ψ 1 and ω 1 are well defiend. (7)

17 Study of dynamic stability through an exact 1D model In [Hou-Li, CPAM, 61 (2008), no. 5, ], we derived an exact 1D model along the z-axis for the Navier-Stokes equations: (u 1 ) t + 2ψ 1 (u 1 ) z = ν(u 1 ) zz + 2 (ψ 1 ) z u 1, (8) (ω 1 ) t + 2ψ 1 (ω 1 ) z = ν(ω 1 ) zz + ( u1 2 ) z, (9) (ψ 1 ) zz = ω 1. (10) Let ũ = u 1, ṽ = (ψ 1 ) z, and ψ = ψ 1. The above system becomes (ũ) t + 2 ψ(ũ) z = ν(ũ) zz 2ṽũ, (11) (ṽ) t + 2 ψ(ṽ) z = ν(ṽ) zz + (ũ) 2 (ṽ) 2 + c(t), (12) where ṽ = ( ψ) z, ṽ z = ω, and c(t) is an integration constant to enforce the mean of ṽ equal to zero.

18 The 1D model is exact! A surprising result is that the above 1D model is exact. Theorem 1. Let u 1, ψ 1 and ω 1 be the solution of the 1D model (??)-(??) and define u θ (r, z, t) = ru 1 (z, t), ω θ (r, z, t) = rω 1 (z, t), ψ θ (r, z, t) = rψ 1 (z, t). Then (u θ (r, z, t), ω θ (r, z, t), ψ θ (r, z, t)) is an exact solution of the 3D Navier-Stokes equations. Theorem 1 tells us that the 1D model (??)-(??) preserves some essential nonlinear structure of the 3D axisymmetric Navier-Stokes equations.

19 Global Well-Posedness of the full 1D Model Theorem 2. Assume that ũ(z, 0) and ṽ(z, 0) are in C m [0, 1] with m 1 and periodic with period 1. Then the solution (ũ, ṽ) of the 1D model will be in C m [0, 1] for all times and for ν 0. Proof. The key is to obtain a priori pointwise estimate for the Lyapunov function ũ 2 z + ṽ 2 z. Differentiating the ũ and ṽ-equations w.r.t z, we get (ũ z ) t + 2 ψ(ũ z ) z 2ṽũ z = 2ṽũ z 2ũṽ z + ν(ũ z ) zz, (ṽ z ) t + 2 ψ(ṽ z ) z 2ṽṽ z = 2ũũ z 2ṽṽ z + ν(ṽ z ) zz. Surprisingly, nonlinear terms cancel each other and we obtain (ũ2 z + ṽz 2 ) + 2 ψ ( ũ 2 t z + ṽz 2 ) = ν ( ũ 2 z z + ṽz 2 ) 2ν [ (ũ zz zz ) 2 + (ṽ zz ) 2]. Thus, (ũ 2 z + ṽ 2 z ) satisfies a maximum principle for all ν 0: ũ 2 z + ṽ 2 z L (ũ 0 ) 2 z + (ṽ 0 ) 2 z L.

20 Construction of a family of globally smooth solutions Theorem 3. Let φ(r) be a smooth cut-off function and u 1, ω 1 and ψ 1 be the solution of the 1D model. Define u θ (r, z, t) = ru 1 (z, t)φ(r) + ũ(r, z, t), ω θ (r, z, t) = rω 1 (z, t)φ(r) + ω(r, z, t), ψ θ (r, z, t) = rψ 1 (z, t)φ(r) + ψ(r, z, t). Then there exists a family of globally smooth functioons ũ, ω and ψ such that u θ, ω θ and ψ θ are globally smooth solutions of the 3D Navier-Stokes equations with finite energy.

21 A New 3D Model for NSE, [Hou-Lei, CPAM, 09] Recall the reformulated 3D Navier-Stokes equations: t u 1 + u r r u 1 + u z z u 1 = ν ( r r ) r + z 2 u1 + 2u 1 ψ 1z, t ω 1 + u r r ω 1 + u z z ω 1 = ν ( r r ) r + z 2 ω1 + ( ) u1 2 z, ( r r ) r + z 2 ψ1 = ω 1, (13) where u r = rψ 1z, u z = 2ψ 1 + rψ 1r. Our 3D model is derived by simply dropping the convective term from (??): t u 1 = ν ( r r ) r + z 2 u1 + 2u 1 ψ 1z, t ω 1 = ν ( r r ) r + z 2 ω1 + (u1 2) z, ( (14) r r ) r + z 2 ψ1 = ω 1. Note that (??) is already a closed system, and u 1 = u θ /r characterizes the axial vorticity near r = 0.

22 Properties of the 3D Model [Hou-Lei, CPAM, 09] This 3D model shares many important properties with the axisymmetric Navier-Stokes equations. First of all, one can define an incompressible velocity field in the model equations (??). u(t, x) = u r (t, r, z)e r + u θ (t, r, z)e θ + u z (t, r, z)e z, (15) u θ = ru 1, u r = rψ 1z, u z = 2ψ 1 + rψ 1r, (16) where x = (x 1, x 2, z), r = x1 2 + x 2 2. It is easy to check that u = r u r + z u z + ur r = 0, (17) which is the same as the Navier-Stokes equations.

23 Properties of the 3D Model continued Our model enjoys the following properties ([Hou-Lei, CPAM-09]): Theorem 4. Energy identity. The smooth solution of (??) satisfies 1 d ( u Dψ 1 2) ( Du1 r 3 drdz + ν D 2 ψ 1 2) r 3 drdz = 0, 2 dt which has been proved to be equivalent to that of the Navier-Stokes equations. Here D is the first order derivative operator defined in R 5. Theorem 5. A non-blowup criterion of Beale-Kato-Majda type. A smooth solution (u 1, ω 1, ψ 1 ) of the model (??) for 0 t < T blows up at time t = T if and only if T where u is defined in (??)-(??). 0 u BMO(R 3 )dt =,

24 Properties of the 3D Model continued Theorem 6. A non-blowup criterion of Prodi-Serrin type. A weak solution (u 1, ω 1, ψ 1 ) of the model (??) is smooth on [0, T ] R 3 provided that u θ L q t Lp x ([0,T ] R 3 ) < (18) for some p, q satisfying 3 p + 2 q 1 with 3 < p and 2 q <. Theorem 7. An analog of Caffarelli-Kohn-Nirenberg partial regularity result [Hou-Lei, CMP-09]. For any suitable weak solution of the 3D model equations (??) on an open set in space-time, the one-dimensional Hausdorff measure of the associated singular set is zero.

25 Numerical evidence for a potential finite time singularity Innitial condition we consider in our numerical computations is given by u 1 (z, r, 0) = (1 + sin(4πz))(r 2 1) 20 (r 2 1.2) 30, 0 r 1, ψ 1 (z, r, 0) = 0, ω 1 (z, r, 0) = 0. A second order finite difference discretization is used in space, and the classical fourth order Runge-Kutta method is used to discretize in time. We use the following coordinate transformation along the r-direction to achieve the adaptivity: r = f (α) α 0.9 sin(πα)/π. We use an effective resolution up to for the 3D problem.

26 A 3D view of u 1 at t = 0.02.

27 A 3D view of u 1 at t =

28 (T t) Asymptotic blowup fit: u 1 1 C, with limiting values T = and C = x 10 4 The inverse of u 1 in time, N=1024, 2048, 3072, 4096, ν=0.001 N=1024 N=2048 N=3072 N=4096 Asymptotic fit 1 u Time

29 Asymptotic blowup rate: u 1 T = and C = C (T t), with 10 x u (blue) vs C/(T t) (red) in time, T= , C= , N=4096, ν=0.001 u, N= Fitted solution u Time

30 Local alignment of u 1 and ψ 1z at t = Recall (u 1 ) t = 2u 1 ψ 1z + ν u 1, (ω 1 ) t = (u 2 1 ) z + ν ω u 1 vs (ψ 1 ) z at r=0, t=0.02, N z =4096, N r =400, t= , ν=0.001 u (ψ 1 ) z z

31 Local alignment of u 1 and ψ 1z at t = x u 1 vs (ψ 1 ) z at r=0, t=0.021, N z =4096, N r =400, t= , ν= u 1 8 (ψ 1 ) z z

32 The effect of convection. To study the effect of convection, we add the convection term back to the 3D model and solve the Navier-Stokes equations using the solution of the 3D model at t = 0.02 as the initial condition Maximum u 1 in time, N z =2048, N r =1024, t= , ν=0.001, full NSE u Time

33 Contours of initial data for u 1. Contour of u 1 at t=0.02, N z =4096, N r =400, t= , ν=0.001, 3D model r z

34 Contours of u 1 at t = 0.021, solution of full NSE. Contour of u 1 at t=0.021, N z =2048, N r =1024, t= , ν=0.001, full NSE r z

35 Contours of u 1 at t = , solution of full NSE. Contour of u 1 at t=0.0235, N z =2048, N r =1024, t= , ν=0.001, full NSE r z

36 Proof of finite time blow-up of the 3D model Theorem 8. [Hou-Shi-Wang, 09]. Consider the 3D inviscid model { u t = 2uψ z, ω t = (u 2 ) z, ( ) x 2 + y 2 + z 2 (19) ψ = ω, 0 x, y, z 1, with boundary condition ψ = 0 at x = 0, 1, y = 0, 1, z = 1, and (α ψ z + ψ) z=0 = 0 for some 0 < α < 1. If the initial conditions, u 0 and ψ 0, are smooth, satisfying u 0 = 0 at z = 0, 1, and log(u 0 )φ(x, y, z)dxdydz 0, (ψ 0 ) z φ(x, y, z)dxdydz > 0, [0,1] 3 [0,1] 3 where φ(x, y, z) = sin(x) sin(y) cosh(α(1 z)), then the 3D inviscid model must develop a finite time singularity. Moreover, if the ω-equation is viscous and ω satisfies the same boundary condition as ψ, then the 3D model with partial viscosity must develop a finite time singularity.

37 Recent progress on full Euler/NSE, joint with Dr. Guo Luo A recent study of the full 3D NSE and Euler equations reveals that for certain axi-symmetric initial data with swirl, the solution could develop a strong dynamic nonlinear alignment in the vortex stretching term. Our preliminary numerical study shows that the solution of the 3D Euler/NSE could experience a strong transitional growth in time with a quadratic nonlinearity, between u 1 and ψ 1z, u 1t + u r (u 1 ) r + u z (u 1 ) z = 2u 1 ψ 1z + ν u 1. An important question is whether the profile near the potentially singular region remains dynamically stable. The dynamic evolution of u 1 and w 1 on the (r, z)-plane can be visualized in the following contour plots.

38 Contour plots of u 1 at t = 0.02 Figure : Contour plot of u 1 at t=0.02.

39 Contour plots of u 1 at t = 0.03 Figure : Contour plot of u 1 at t=0.03.

40 Contour plots of u 1 at t = Figure : Contour plot of u 1 at t=0.036.

41 Contour plots of u 1 at t = Figure : Contour plot of u 1 at t=0.043.

42 Contour plots of u 1 at t = Figure : Contour plot of u 1 at t= Adaptive mesh is uded with smallest mesh proportional to 10 6.

43 Velocity field and streamlines at t = The solution can also be visualized in 3D. Figure : Velocity field (u r, u θ, u z) and streamlines at t =

44 Vorticity iso-surface at t = Figure : Iso-surface of vorticity ω (with ω 1 max ω ) at t = Note the nontrivial geometry of the iso-surface.

45 Vortex lines at t = Figure : Vortex lines at t = Note the high curvature of the vortex lines near the maximum vorticity region.

46 Log log plot of the maximum of u 1 and ω 1 The following figure shows the rapid, super-double-exponential growth of the maximum of u 1 and ω 1 in time. 1.9 log(log(max u1 )) as a function of t 2.6 log(log(max w1 )) as a function of t log(log(max u1 )) log(log(max w1 )) t t Figure : Log log plot of maximum u 1 and ω 1 as a function of time.

47 The quadratic nonlinearity, u 1t + u u 1 = 2u 1 ψ 1z We observe that the nonlinear u 1 -stretching term, 2ψ 1z, is comparable with u 1 at the maximum of u 1. This quadratic nonlinearity is the driving force behind the rapid growth of u 1 and ω Nonlinear stretching of u1 a function of t 2 phi1 z max u t Figure : Ratio of the u 1-stretching term, 2ψ 1, to u 1 at the maximum of u 1.

48 Approaching the axis r = 0 To see how fast the solution is approaching r = 0, we plot the location (r 1, z 1 ) of the maximum of u 1 as a function of time. 3.5 x 10 3 r1 2 as a function of t, where r1 is the r coordinate of maxloc u1 3.5 x 10 3 z1 2 as a function of t, where z1 is the z coordinate of maxloc u r1 2 z t t Figure : Plot of r1 2 and z1 2 as a function of time, where (r 1, z 1) is the point at which u 1 attains its maximum. This suggests that r 1(t) c 1(T t) 1/2 and z 1(t) c 2(T t) 1/2.

49 Front thickness as a function of time 0.07 Thickness of the front along r as a function of t lr t Figure : The front thickness as a function of time at the point at which u 1 attains its maximum. This seems to suggest that the thickness decreases like O(T t).

50 Asymptotic scaling for 1/ max u 1 Assuming that the solution blows up at a finite time T with a blowup rate O(T t) α, we can obtain a crude estimate of T and α using a line fitting procedure. 9 x /max u1 as a function of t /max u t Figure : Plot of 1/ max u 1 as a function of time. This suggests that u 1 L c/(t t).

51 Asymptotic scaling for 1/ max ω x /max w1 2/3 as a function of t 2 1/max w1 2/ t Figure : Plot of 1/ max ω 1 2/3 as a function of time. This seems to suggest that ω 1 L c/(t t) 3/2.

52 Asymptotic scaling for 1/ max ψ /max phi1 2 as a function of t /phi t Figure : Plot of 1/ max ψ 1 2 as a function of time. This seems to suggest that ψ 1 L c/(t t) 1/2. The estimated singularity time is T

53 Fluid dynamic instability seems to induce front instability Unfortunately (or fortunately), the fluid dynamics instability (Kelvin-Helmholtz) eventually kicks in when the thickness of the vortex sheet like structure is below some critical value. The vortex sheet develops roll-up and destroys the strong nonlinear alignment in the vortex stretching term. The flow subsequently develops many small scales dyanamically and transitions to turbulence.

54 Contour plots of u 1 at t = , Euler equation Figure : Contour plot of u 1 at t= , the effective adaptive mesh near the most singular region = 10 6.

55 Contour plots of u 1 at t = , Euler equation Figure : Contour plot of u 1 at t=

56 Contour plots of u 1 at t = , Euler equation Figure : Contour plot of u 1 at t=

57 Contour plots of u 1 at t = , NSE with ν = 10 5 Figure : Contour plot of u 1 at t= , the effective adaptive mesh near the most singular region =

58 Contour plots of u 1 at t = , NSE with ν = 10 5 Figure : Contour plot of u 1 at t=

59 Contour plots of u 1 at t = , NSE with ν = 10 5 Figure : Contour plot of u 1 at t=

60 Contour plots of u 1 at t = , NSE with ν = 10 5 Figure : Contour plot of u 1 at t=

61 Contour plots of u 1 at t = , ν = 10 5 Figure : Contour plot of u 1 at t=

62 Concluding Remarks Our study shows that convection could play a dual role in the dynamic stability or finite time singularity formation of 3D incompressible flow. One the one hand, we construct an exact 1D model, in which we show that convection plays an essential role to deplete the nonlinear vortex stretching, and the 1D model has global regularity. We also construct a new 3D model by neglecting the convection term. This 3D model shares almost all properties of the Navier-Stokes equations, but could develop finite time singularities. On the other hand, our recent study suggests that convection may not always stablize the nonlinear vortex stretching. For initial data with certain symmetry, the combined effect of convection and vortex stretching could lead to a strong nonlinear alignment. It remains to understand whether the nearly singular behavior we observed numerically is dynamically stable. If so, it may suggest a possible way to construct a singular solution analytically.

63 References J. Deng, T. Y. Hou, and X. Yu, Geometric Properties and the non-blow-up of the Three-Dimensional Euler Equation, Comm. PDEs, 30:1 (2005), T. Y. Hou and R. Li, Dynamic Depletion of Vortex Stretching and Non-Blowup of the 3-D Incompressible Euler Equations, J. Nonlinear Science, 16 (2006), T. Y. Hou and C. Li, Dynamic Stability of the 3D Axisymmetric Navier-Stokes Equations with Swirl, CPAM, 61 (2008), T. Y. Hou and Z. Lei, On the Stabilizing Effect of Convection in 3D Incompressible Flows, CPAM, 62, pp , T. Y. Hou, Z. Shi, and S. Wang, On Singularity Formation of a 3D Model for Incompressible Navier-Stokes Equations, 2009, arxiv: v1 [math.ap]. T. Y. Hou, Blow-up or no blowup? A unified computational and analytical approach to 3D incompressible Euler and Navier-Stokes equations, Acta Numerica, pp , 2009.

The Role of Convection and Nearly Singular Behavior of the 3D Navier-Stokes Equations

The Role of Convection and Nearly Singular Behavior of the 3D Navier-Stokes Equations The Role of Convection and Nearly Singular Behavior of the 3D Navier-Stokes Equations Thomas Y. Hou Applied and Comput. Mathematics, Caltech PDE Conference in honor of Blake Temple, University of Michigan

More information

Blow-up or No Blow-up? the Role of Convection in 3D Incompressible Navier-Stokes Equations

Blow-up or No Blow-up? the Role of Convection in 3D Incompressible Navier-Stokes Equations Blow-up or No Blow-up? the Role of Convection in 3D Incompressible Navier-Stokes Equations Thomas Y. Hou Applied and Comput. Mathematics, Caltech Joint work with Zhen Lei; Congming Li, Ruo Li, and Guo

More information

The interplay between theory and computation in the study of 3D Euler equations

The interplay between theory and computation in the study of 3D Euler equations The interplay between theory and computation in the study of 3D Euler equations Thomas Y. Hou Applied and Computational Mathematics California Institute of Technology Joint work with Guo Luo and Pengfei

More information

Potentially Singular Solutions of the 3D Axisymmetric Euler Equations

Potentially Singular Solutions of the 3D Axisymmetric Euler Equations Potentially Singular Solutions of the 3D Axisymmetric Euler Equations Thomas Y. Hou Applied and Computational Mathematics California Institute of Technology Joint work with Guo Luo IPAM Workshop on Mathematical

More information

Blow-up or no blow-up? A unified computational and analytic approach to 3D incompressible Euler and Navier Stokes equations

Blow-up or no blow-up? A unified computational and analytic approach to 3D incompressible Euler and Navier Stokes equations Acta Numerica (29), pp. 277 346 c Cambridge University Press, 29 doi: 1.117/S962492964218 Printed in the United Kingdom Blow-up or no blow-up? A unified computational and analytic approach to 3D incompressible

More information

Incompressible Navier-Stokes Equations in R 3

Incompressible Navier-Stokes Equations in R 3 Incompressible Navier-Stokes Equations in R 3 Zhen Lei ( ) School of Mathematical Sciences Fudan University Incompressible Navier-Stokes Equations in R 3 p. 1/5 Contents Fundamental Work by Jean Leray

More information

The Dual Role of Convection in 3D. Navier-Stokes equations.

The Dual Role of Convection in 3D. Navier-Stokes equations. 1 The Dual Role of Convection in 3D Navier-Stokes Equations Thomas Y. Hou a, Zuoqiang Shi b, Shu Wang c Abstract We investigate the dual role of convection on the large time behavior of the 3D incompressible

More information

The enigma of the equations of fluid motion: A survey of existence and regularity results

The enigma of the equations of fluid motion: A survey of existence and regularity results The enigma of the equations of fluid motion: A survey of existence and regularity results RTG summer school: Analysis, PDEs and Mathematical Physics The University of Texas at Austin Lecture 1 1 The review

More information

On Finite Time Singularity and Global Regularity of an Axisymmetric Model for the 3D Euler Equations

On Finite Time Singularity and Global Regularity of an Axisymmetric Model for the 3D Euler Equations On Finite Time Singularity and Global Regularity of an Axisymmetric Model for the 3D Euler Equations arxiv:23.298v [math.ap] 4 Mar 22 Thomas Y. Hou Zhen Lei Shu Wang Chen Zou March 5, 22 Abstract We investigate

More information

Topics in Fluid Dynamics: Classical physics and recent mathematics

Topics in Fluid Dynamics: Classical physics and recent mathematics Topics in Fluid Dynamics: Classical physics and recent mathematics Toan T. Nguyen 1,2 Penn State University Graduate Student Seminar @ PSU Jan 18th, 2018 1 Homepage: http://math.psu.edu/nguyen 2 Math blog:

More information

The 3D Euler and 2D surface quasi-geostrophic equations

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

On Global Well-Posedness of the Lagrangian Averaged Euler Equations

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

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

Vorticity and Dynamics

Vorticity and Dynamics Vorticity and Dynamics In Navier-Stokes equation Nonlinear term ω u the Lamb vector is related to the nonlinear term u 2 (u ) u = + ω u 2 Sort of Coriolis force in a rotation frame Viscous term ν u = ν

More information

of some mathematical models in physics

of some mathematical models in physics Regularity, singularity and well-posedness of some mathematical models in physics Saleh Tanveer (The Ohio State University) Background Mathematical models involves simplifying assumptions where "small"

More information

Global regularity of a modified Navier-Stokes equation

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

u t + u u = p (1.1) u = 0 (1.2)

u t + u u = p (1.1) u = 0 (1.2) METHODS AND APPLICATIONS OF ANALYSIS. c 2005 International Press Vol. 12, No. 4, pp. 427 440, December 2005 004 A LEVEL SET FORMULATION FOR THE 3D INCOMPRESSIBLE EULER EQUATIONS JIAN DENG, THOMAS Y. HOU,

More information

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

Level Set Dynamics and the Non-blowup of the 2D Quasi-geostrophic Equation

Level Set Dynamics and the Non-blowup of the 2D Quasi-geostrophic Equation Level Set Dynamics and the Non-blowup of the 2D Quasi-geostrophic Equation J. Deng, T. Y. Hou, R. Li, X. Yu May 31, 2006 Abstract In this article we apply the technique proposed in Deng-Hou-Yu [7] to study

More information

Multiscale Computation of Isotropic Homogeneous Turbulent Flow

Multiscale Computation of Isotropic Homogeneous Turbulent Flow Multiscale Computation of Isotropic Homogeneous Turbulent Flow Tom Hou, Danping Yang, and Hongyu Ran Abstract. In this article we perform a systematic multi-scale analysis and computation for incompressible

More information

arxiv: v2 [math.ap] 6 Sep 2007

arxiv: v2 [math.ap] 6 Sep 2007 ON THE REGULARITY OF WEAK SOLUTIONS OF THE 3D NAVIER-STOKES EQUATIONS IN B 1, arxiv:0708.3067v2 [math.ap] 6 Sep 2007 A. CHESKIDOV AND R. SHVYDKOY ABSTRACT. We show that if a Leray-Hopf solution u to the

More information

On the limiting behaviour of regularizations of the Euler Equations with vortex sheet initial data

On the limiting behaviour of regularizations of the Euler Equations with vortex sheet initial data On the limiting behaviour of regularizations of the Euler Equations with vortex sheet initial data Monika Nitsche Department of Mathematics and Statistics University of New Mexico Collaborators: Darryl

More information

arxiv: v2 [physics.flu-dyn] 8 Dec 2013

arxiv: v2 [physics.flu-dyn] 8 Dec 2013 POTENTIALLY SINGULAR SOLUTIONS OF THE 3D INCOMPRESSIBLE EULER EQUATIONS GUO LUO AND THOMAS Y. HOU arxiv:131.497v2 [physics.flu-dyn] 8 Dec 213 December 1, 213 Abstract Whether the 3D incompressible Euler

More information

Spatial Profiles in the Singular Solutions of the 3D Euler Equations and Simplified Models

Spatial Profiles in the Singular Solutions of the 3D Euler Equations and Simplified Models Spatial Profiles in the Singular Solutions of the 3D Euler Equations and Simplified Models Thesis by Pengfei Liu In Partial Fulfillment of the Requirements for the degree of Doctor of Philosophy CALIFORNIA

More information

On the well-posedness of the Prandtl boundary layer equation

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

Remarks on the blow-up criterion of the 3D Euler equations

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

arxiv: v1 [math.ap] 16 May 2007

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

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

Vortex stretching in incompressible and compressible fluids

Vortex stretching in incompressible and compressible fluids Vortex stretching in incompressible and compressible fluids Esteban G. Tabak, Fluid Dynamics II, Spring 00 1 Introduction The primitive form of the incompressible Euler equations is given by ( ) du P =

More information

FINITE TIME BLOW-UP FOR A DYADIC MODEL OF THE EULER EQUATIONS

FINITE TIME BLOW-UP FOR A DYADIC MODEL OF THE EULER EQUATIONS TRANSACTIONS OF THE AMERICAN MATHEMATICAL SOCIETY Volume 357, Number 2, Pages 695 708 S 0002-9947(04)03532-9 Article electronically published on March 12, 2004 FINITE TIME BLOW-UP FOR A DYADIC MODEL OF

More information

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

c 2014 Society for Industrial and Applied Mathematics

c 2014 Society for Industrial and Applied Mathematics MULTISCALE MODEL. SIMUL. Vol. 12, No. 4, pp. 1722 1776 c 214 Society for Industrial and Applied Mathematics TOWARD THE FINITE-TIME BLOWUP OF THE 3D AXISYMMETRIC EULER EQUATIONS: A NUMERICAL INVESTIGATION

More information

On Smoothness of Suitable Weak Solutions to the Navier-Stokes Equations

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

Nonuniqueness of weak solutions to the Navier-Stokes equation

Nonuniqueness of weak solutions to the Navier-Stokes equation Nonuniqueness of weak solutions to the Navier-Stokes equation Tristan Buckmaster (joint work with Vlad Vicol) Princeton University November 29, 2017 Tristan Buckmaster (Princeton University) Nonuniqueness

More information

Nonlinear Analysis. A regularity criterion for the 3D magneto-micropolar fluid equations in Triebel Lizorkin spaces

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

Lecture 3: The Navier-Stokes Equations: Topological aspects

Lecture 3: The Navier-Stokes Equations: Topological aspects Lecture 3: The Navier-Stokes Equations: Topological aspects September 9, 2015 1 Goal Topology is the branch of math wich studies shape-changing objects; objects which can transform one into another without

More information

Borel Summability in PDE initial value problems

Borel Summability in PDE initial value problems Borel Summability in PDE initial value problems Saleh Tanveer (Ohio State University) Collaborator Ovidiu Costin & Guo Luo Research supported in part by Institute for Math Sciences (IC), EPSRC & NSF. Main

More information

Journal of Differential Equations

Journal of Differential Equations J. Differential Equations 48 (1) 6 74 Contents lists available at ScienceDirect Journal of Differential Equations www.elsevier.com/locate/jde Two regularity criteria for the D MHD equations Chongsheng

More information

Applying Moser s Iteration to the 3D Axially Symmetric Navier Stokes Equations (ASNSE)

Applying Moser s Iteration to the 3D Axially Symmetric Navier Stokes Equations (ASNSE) Applying Moser s Iteration to the 3D Axially Symmetric Navier Stokes Equations (ASNSE) Advisor: Qi Zhang Department of Mathematics University of California, Riverside November 4, 2012 / Graduate Student

More information

α 2 )(k x ũ(t, η) + k z W (t, η)

α 2 )(k x ũ(t, η) + k z W (t, η) 1 3D Poiseuille Flow Over the next two lectures we will be going over the stabilization of the 3-D Poiseuille flow. For those of you who havent had fluids, that means we have a channel of fluid that is

More information

The Euler Equations and Non-Local Conservative Riccati Equations

The Euler Equations and Non-Local Conservative Riccati Equations The Euler Equations and Non-Local Conservative Riccati Equations Peter Constantin Department of Mathematics The University of Chicago November 8, 999 The purpose of this brief note is to present an infinite

More information

TOWARD THE FINITE-TIME BLOWUP OF THE 3D AXISYMMETRIC EULER EQUATIONS: A NUMERICAL INVESTIGATION

TOWARD THE FINITE-TIME BLOWUP OF THE 3D AXISYMMETRIC EULER EQUATIONS: A NUMERICAL INVESTIGATION TOWARD THE FINITE-TIME BLOWUP OF THE 3D AXISYMMETRIC EULER EQUATIONS: A NUMERICAL INVESTIGATION GUO LUO AND THOMAS Y. HOU Abstract. Whether the 3D incompressible Euler equations can develop a singularity

More information

A regularity criterion for the 3D NSE in a local version of the space of functions of bounded mean oscillations

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

Vortex stretching and anisotropic diffusion in the 3D NSE

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

3 Generation and diffusion of vorticity

3 Generation and diffusion of vorticity Version date: March 22, 21 1 3 Generation and diffusion of vorticity 3.1 The vorticity equation We start from Navier Stokes: u t + u u = 1 ρ p + ν 2 u 1) where we have not included a term describing a

More information

MULTISCALE ANALYSIS IN LAGRANGIAN FORMULATION FOR THE 2-D INCOMPRESSIBLE EULER EQUATION. Thomas Y. Hou. Danping Yang. Hongyu Ran

MULTISCALE ANALYSIS IN LAGRANGIAN FORMULATION FOR THE 2-D INCOMPRESSIBLE EULER EQUATION. Thomas Y. Hou. Danping Yang. Hongyu Ran DISCRETE AND CONTINUOUS Website: http://aimsciences.org DYNAMICAL SYSTEMS Volume 13, Number 5, December 2005 pp. 1153 1186 MULTISCALE ANALYSIS IN LAGRANGIAN FORMULATION FOR THE 2-D INCOMPRESSIBLE EULER

More information

Global well-posedness of the primitive equations of oceanic and atmospheric dynamics

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

ON THE REGULARITY OF WEAK SOLUTIONS OF THE 3D NAVIER-STOKES EQUATIONS IN B 1

ON THE REGULARITY OF WEAK SOLUTIONS OF THE 3D NAVIER-STOKES EQUATIONS IN B 1 ON THE REGULARITY OF WEAK SOLUTIONS OF THE 3D NAVIER-STOKES EQUATIONS IN B 1, A. CHESKIDOV AND R. SHVYDKOY ABSTRACT. We show that if a Leray-Hopf solution u to the 3D Navier- Stokes equation belongs to

More information

F11AE1 1. C = ρν r r. r u z r

F11AE1 1. C = ρν r r. r u z r F11AE1 1 Question 1 20 Marks) Consider an infinite horizontal pipe with circular cross-section of radius a, whose centre line is aligned along the z-axis; see Figure 1. Assume no-slip boundary conditions

More information

Fundamentals of Fluid Dynamics: Ideal Flow Theory & Basic Aerodynamics

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

A New Regularity Criterion for the 3D Navier-Stokes Equations via Two Entries of the Velocity Gradient

A New Regularity Criterion for the 3D Navier-Stokes Equations via Two Entries of the Velocity Gradient Acta Appl Math (014) 19:175 181 DOI 10.1007/s10440-013-9834-3 A New Regularity Criterion for the 3D Navier-Stokes Euations via Two Entries of the Velocity Gradient Tensor Zujin Zhang Dingxing Zhong Lin

More information

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

Multiscale Computation of Isotropic Homogeneous Turbulent Flow

Multiscale Computation of Isotropic Homogeneous Turbulent Flow Multiscale Computation of Isotropic Homogeneous Turbulent Flow Tom Hou, Danping Yang, and Hongyu Ran Abstract. In this article we perform a systematic multi-scale analysis and computation for incompressible

More information

Available online at J. Math. Comput. Sci. 4 (2014), No. 3, ISSN:

Available online at   J. Math. Comput. Sci. 4 (2014), No. 3, ISSN: Available online at http://scik.org J. Math. Comput. Sci. 4 (2014), No. 3, 587-593 ISSN: 1927-5307 A SMALLNESS REGULARITY CRITERION FOR THE 3D NAVIER-STOKES EQUATIONS IN THE LARGEST CLASS ZUJIN ZHANG School

More information

Frequency Localized Regularity Criteria for the 3D Navier Stokes Equations. Z. Bradshaw & Z. Grujić. Archive for Rational Mechanics and Analysis

Frequency Localized Regularity Criteria for the 3D Navier Stokes Equations. Z. Bradshaw & Z. Grujić. Archive for Rational Mechanics and Analysis Frequency Localized Regularity Criteria for the 3D Navier Stokes Equations Z. Bradshaw & Z. Gruić Archive for Rational Mechanics and Analysis ISSN 0003-9527 Arch Rational Mech Anal DOI 10.1007/s00205-016-1069-9

More information

OSGOOD TYPE REGULARITY CRITERION FOR THE 3D NEWTON-BOUSSINESQ EQUATION

OSGOOD TYPE REGULARITY CRITERION FOR THE 3D NEWTON-BOUSSINESQ EQUATION Electronic Journal of Differential Equations, Vol. 013 (013), No. 3, pp. 1 6. ISSN: 107-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu ftp ejde.math.txstate.edu OSGOOD TYPE REGULARITY

More information

arxiv: v1 [math.ap] 3 Nov 2017

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

Boundary conditions and estimates for the linearized Navier-Stokes equations on staggered grids

Boundary conditions and estimates for the linearized Navier-Stokes equations on staggered grids Boundary conditions and estimates for the linearized Navier-Stokes equations on staggered grids Wendy Kress and Jonas Nilsson Department of Scientific Computing Uppsala University Box S-75 4 Uppsala Sweden

More information

Citation 数理解析研究所講究録 (2016), 1995:

Citation 数理解析研究所講究録 (2016), 1995: Swirling flow of the axi-symmetric Titlenear a saddle point and no-slip bou New Developments for Elucidating In Author(s) 野津, 裕史 ; 米田, 剛 Citation 数理解析研究所講究録 (2016), 1995: 9-13 Issue Date 2016-04 URL http://hdl.handle.net/2433/224715

More information

Dissipative quasi-geostrophic equations with L p data

Dissipative quasi-geostrophic equations with L p data Electronic Journal of Differential Equations, Vol. (), No. 56, pp. 3. ISSN: 7-669. URL: http://ejde.math.swt.edu or http://ejde.math.unt.edu ftp ejde.math.swt.edu (login: ftp) Dissipative quasi-geostrophic

More information

Properties at potential blow-up times for the incompressible Navier-Stokes equations

Properties at potential blow-up times for the incompressible Navier-Stokes equations Properties at potential blow-up times for the incompressible Navier-Stokes equations Jens Lorenz Department of Mathematics and Statistics University of New Mexico Albuquerque, NM 87131, USA Paulo R. Zingano

More information

Euler Equations of Incompressible Ideal Fluids

Euler Equations of Incompressible Ideal Fluids Euler Equations of Incompressible Ideal Fluids Claude BARDOS, and Edriss S. TITI February 27, 27 Abstract This article is a survey concerning the state-of-the-art mathematical theory of the Euler equations

More information

In honor of Prof. Hokee Minn

In honor of Prof. Hokee Minn In honor of Prof. Hokee Minn It is my greatest honor to be here today to deliver the nd Minn Hokee memorial lecture. I would like to thank the Prof. Ha for the nomination and the selecting committee for

More information

Maximum Rate of Growth of Enstrophy in Solutions of the Fractional Burgers Equation

Maximum Rate of Growth of Enstrophy in Solutions of the Fractional Burgers Equation Maximum Rate of Growth of Enstrophy in Solutions of the Fractional Burgers Equation arxiv:161.9578v [physics.flu-dyn] 15 Jan 18 Dongfang Yun and Bartosz Protas Department of Mathematics and Statistics,

More information

0 = p. 2 x + 2 w. z +ν w

0 = p. 2 x + 2 w. z +ν w Solution (Elliptical pipe flow (a Using the Navier Stokes equations in three dimensional cartesian coordinates, given that u =, v = and w = w(x,y only, and assuming no body force, we are left with = p

More information

REGULARITY OF GENERALIZED NAVEIR-STOKES EQUATIONS IN TERMS OF DIRECTION OF THE VELOCITY

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

An Eulerian-Lagrangian Approach for Incompressible Fluids: Local Theory

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

UNIVERSITY OF EAST ANGLIA

UNIVERSITY OF EAST ANGLIA UNIVERSITY OF EAST ANGLIA School of Mathematics May/June UG Examination 2007 2008 FLUIDS DYNAMICS WITH ADVANCED TOPICS Time allowed: 3 hours Attempt question ONE and FOUR other questions. Candidates must

More information

Partial regularity for suitable weak solutions to Navier-Stokes equations

Partial regularity for suitable weak solutions to Navier-Stokes equations Partial regularity for suitable weak solutions to Navier-Stokes equations Yanqing Wang Capital Normal University Joint work with: Quansen Jiu, Gang Wu Contents 1 What is the partial regularity? 2 Review

More information

THE POINCARÉ RECURRENCE PROBLEM OF INVISCID INCOMPRESSIBLE FLUIDS

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

TRANSIENT FLOW AROUND A VORTEX RING BY A VORTEX METHOD

TRANSIENT FLOW AROUND A VORTEX RING BY A VORTEX METHOD Proceedings of The Second International Conference on Vortex Methods, September 26-28, 21, Istanbul, Turkey TRANSIENT FLOW AROUND A VORTEX RING BY A VORTEX METHOD Teruhiko Kida* Department of Energy Systems

More information

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

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

More information

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

PAPER 331 HYDRODYNAMIC STABILITY

PAPER 331 HYDRODYNAMIC STABILITY MATHEMATICAL TRIPOS Part III Thursday, 6 May, 016 1:30 pm to 4:30 pm PAPER 331 HYDRODYNAMIC STABILITY Attempt no more than THREE questions. There are FOUR questions in total. The questions carry equal

More information

On Some Variational Optimization Problems in Classical Fluids and Superfluids

On Some Variational Optimization Problems in Classical Fluids and Superfluids On Some Variational Optimization Problems in Classical Fluids and Superfluids Bartosz Protas Department of Mathematics & Statistics McMaster University, Hamilton, Ontario, Canada URL: http://www.math.mcmaster.ca/bprotas

More information

Multiscale method and pseudospectral simulations for linear viscoelastic incompressible flows

Multiscale method and pseudospectral simulations for linear viscoelastic incompressible flows Interaction and Multiscale Mechanics, Vol. 5, No. 1 (2012) 27-40 27 Multiscale method and pseudospectral simulations for linear viscoelastic incompressible flows Ling Zhang and Jie Ouyang* Department of

More information

Week 6 Notes, Math 865, Tanveer

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

More information

Lecture 1: Introduction to Linear and Non-Linear Waves

Lecture 1: Introduction to Linear and Non-Linear Waves Lecture 1: Introduction to Linear and Non-Linear Waves Lecturer: Harvey Segur. Write-up: Michael Bates June 15, 2009 1 Introduction to Water Waves 1.1 Motivation and Basic Properties There are many types

More information

Global well-posedness for the critical 2D dissipative quasi-geostrophic equation

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

This note presents an infinite-dimensional family of exact solutions of the incompressible three-dimensional Euler equations

This note presents an infinite-dimensional family of exact solutions of the incompressible three-dimensional Euler equations IMRN International Mathematics Research Notices 2000, No. 9 The Euler Equations and Nonlocal Conservative Riccati Equations Peter Constantin This note presents an infinite-dimensional family of exact solutions

More information

arxiv:math/ v1 [math.ap] 6 Mar 2007

arxiv:math/ v1 [math.ap] 6 Mar 2007 arxiv:math/73144v1 [math.ap] 6 Mar 27 ON THE GLOBAL EXISTENCE FOR THE AXISYMMETRIC EULER EQUATIONS HAMMADI ABIDI, TAOUFIK HMIDI, AND SAHBI KERAANI Abstract. This paper deals with the global well-posedness

More information

Blow Up of Incompressible Euler Solutions

Blow Up of Incompressible Euler Solutions BIT Numerical Mathematics (2006)46:000-000 c Springer 2006. DOI:10.1007/s10543-000-0000-x Blow Up of Incompressible Euler Solutions JOHAN HOFFMAN 1 and CLAES JOHNSON 1 April 1 2008 1 School of Computer

More information

A STOCHASTIC LAGRANGIAN REPRESENTATION OF THE 3-DIMENSIONAL INCOMPRESSIBLE NAVIER-STOKES EQUATIONS

A STOCHASTIC LAGRANGIAN REPRESENTATION OF THE 3-DIMENSIONAL INCOMPRESSIBLE NAVIER-STOKES EQUATIONS A STOCHASTIC LAGRANGIAN REPRESENTATION OF THE 3-DIMENSIONAL INCOMPRESSIBLE NAVIER-STOKES EQUATIONS PETER CONSTANTIN AND GAUTAM IYER Abstract. In this paper we derive a probabilistic representation of the

More information

REGULARITY CRITERIA FOR WEAK SOLUTIONS TO 3D INCOMPRESSIBLE MHD EQUATIONS WITH HALL TERM

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

TOPICS IN MATHEMATICAL AND COMPUTATIONAL FLUID DYNAMICS

TOPICS IN MATHEMATICAL AND COMPUTATIONAL FLUID DYNAMICS TOPICS IN MATHEMATICAL AND COMPUTATIONAL FLUID DYNAMICS Lecture Notes by Prof. Dr. Siddhartha Mishra and Dr. Franziska Weber October 4, 216 1 Contents 1 Fundamental Equations of Fluid Dynamics 4 1.1 Eulerian

More information

GFD 2012 Lecture 1: Dynamics of Coherent Structures and their Impact on Transport and Predictability

GFD 2012 Lecture 1: Dynamics of Coherent Structures and their Impact on Transport and Predictability GFD 2012 Lecture 1: Dynamics of Coherent Structures and their Impact on Transport and Predictability Jeffrey B. Weiss; notes by Duncan Hewitt and Pedram Hassanzadeh June 18, 2012 1 Introduction 1.1 What

More information

Development of high vorticity structures in incompressible 3D Euler equations.

Development of high vorticity structures in incompressible 3D Euler equations. Development of high vorticity structures in incompressible 3D Euler equations. D. Agafontsev (a), A. Mailybaev (b), (c) and E. Kuznetsov (d). (a) - P.P. Shirshov Institute of Oceanology of RAS, Moscow,

More information

Sublayer of Prandtl boundary layers

Sublayer of Prandtl boundary layers ublayer of Prandtl boundary layers Emmanuel Grenier Toan T. Nguyen May 12, 2017 Abstract The aim of this paper is to investigate the stability of Prandtl boundary layers in the vanishing viscosity limit:

More information

The Orchestra of Partial Differential Equations. Adam Larios

The Orchestra of Partial Differential Equations. Adam Larios The Orchestra of Partial Differential Equations Adam Larios 19 January 2017 Landscape Seminar Outline 1 Fourier Series 2 Some Easy Differential Equations 3 Some Not-So-Easy Differential Equations Outline

More information

Research Statement. 1 Overview. Zachary Bradshaw. October 20, 2016

Research Statement. 1 Overview. Zachary Bradshaw. October 20, 2016 Research Statement Zachary Bradshaw October 20, 2016 1 Overview My research is in the field of partial differential equations. I am primarily interested in the three dimensional non-stationary Navier-Stokes

More information

General introduction to Hydrodynamic Instabilities

General introduction to Hydrodynamic Instabilities KTH ROYAL INSTITUTE OF TECHNOLOGY General introduction to Hydrodynamic Instabilities L. Brandt & J.-Ch. Loiseau KTH Mechanics, November 2015 Luca Brandt Professor at KTH Mechanics Email: luca@mech.kth.se

More information

Good reasons to study Euler equation.

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

V (r,t) = i ˆ u( x, y,z,t) + ˆ j v( x, y,z,t) + k ˆ w( x, y, z,t)

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

Fourth Order Convergence of Compact Finite Difference Solver for 2D Incompressible Flow

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

Quasi-two Dimensional Hydrodynamics and Interaction of Vortex Tubes

Quasi-two Dimensional Hydrodynamics and Interaction of Vortex Tubes Quasi-two Dimensional Hydrodynamics and Interaction of Vortex Tubes Vladimir Zakharov 1 Introduction This paper is long overdue. The most of the results presented here were obtained in 1986-87. Just a

More information

3.5 Vorticity Equation

3.5 Vorticity Equation .0 - Marine Hydrodynamics, Spring 005 Lecture 9.0 - Marine Hydrodynamics Lecture 9 Lecture 9 is structured as follows: In paragraph 3.5 we return to the full Navier-Stokes equations (unsteady, viscous

More information

Evolution of sharp fronts for the surface quasi-geostrophic equation

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

Shear Turbulence. Fabian Waleffe. Depts. of Mathematics and Engineering Physics. Wisconsin

Shear Turbulence. Fabian Waleffe. Depts. of Mathematics and Engineering Physics. Wisconsin Shear Turbulence Fabian Waleffe Depts. of Mathematics and Engineering Physics, Madison University of Wisconsin Mini-Symposium on Subcritical Flow Instability for training of early stage researchers Mini-Symposium

More information

Extreme Vortex States and the Growth of Enstrophy in 3D Incompressible Flows

Extreme Vortex States and the Growth of Enstrophy in 3D Incompressible Flows arxiv:1605.05742v2 [physics.flu-dyn] 28 Feb 2017 Extreme Vortex States and the Growth of Enstrophy in 3D Incompressible Flows Diego Ayala 1,2 and Bartosz Protas 2, 1 Department of Mathematics, University

More information

Anisotropic partial regularity criteria for the Navier-Stokes equations

Anisotropic partial regularity criteria for the Navier-Stokes equations Anisotropic partial regularity criteria for the Navier-Stokes equations Walter Rusin Department of Mathematics Mathflows 205 Porquerolles September 7, 205 The question of regularity of the weak solutions

More information