In this paper, we investigate the boundary layer behavior of solutions to the one dimensional Broadwell model of the nonlinear Boltzmann equation for small mean free path. We consider the analogue of Maxwell's diffusive and the reflexive boundary conditions. It is found that even for such a simple model, there are boundary layers due to purely kinetic effects which cannot be detected by the corresponding Navier-Stokes system. It is also found numerically that a compressive boundary layer is not always stable in the sense that it may detach from the boundary and move into the interior of the gas as a shock layer.

In this paper, we study the fluid-dynamic limit for the one-dimensional Broadwell model of the nonlinear Boltzmann equation in the presence of boundaries. We consider an analogue of Maxwell's diffusive and reflective boundary conditions. The boundary layers can be classified as either compressive or expansive in terms of the associated characteristic fields. We show that both expansive and compressive boundary layers (before detachment) are nonlinearly stable and that the layer effects are localized so that the fluid dynamic approximation is valid away from the boundary. We also show that the same conclusion holds for short time without the structural conditions on the boundary layers. A rigorous estimate for the distance between the kinetic solution and the fluid-dynamic solution in terms of the mean-free path in theL∞-norm is obtained provided that the interior fluid flow is smooth. The rate of convergence is optimal.

We investigate, by numerical simulation, the shear layer instability associated with the outer layer of a spiral vortex formed behind an impulsively started thin ellipse. The unstable free shear layer undergoes a secondary instability. We connect this instability with the dynamics of corner vortices adjacent to the tip of the ellipse by observing that the typical turnover time of the corner vortex matches the period of the unstable mode in the shear layer. We suggest that the corner vortex acts as a signal generator, and produces periodic perturbation which triggers the instability.

We give an elementary proof of the convergence of the point vortex method (PVM) to a classical weak solution for the two-dimensional incompressible Euler equations with initial vorticity being a finite Radon measure of distinguished sign and the initial velocity of locally bounded energy. This includes the important example of vortex sheets, which exhibits the classical Kelvin-Helmholtz instability. A surprise fact is that although the velocity fields generated by the point vortex method do not have bounded local kinetic energy, the limiting velocity field is shown to have a bounded local kinetic energy.

Two nonlinear diffusion equations for thin film epitaxy, with or without slope selection, are studied in this work. The nonlinearity models the Ehrlich–Schwoebel effect – the kinetic asymmetry in attachment and detachment of adatoms to and from terrace boundaries. Both perturbation analysis and numerical simulation are presented to show that such an atomistic effect is the origin of a nonlinear morphological instability, in a rough-smooth-rough pattern, that has been experimentally observed as transient in an early stage of epitaxial growth on rough surfaces. Initial-boundary-value problems for both equations are proven to be well-posed, and the solution regularity is also obtained. Galerkin spectral approximations are studied to provide both a priori bounds for proving the well-posedness and numerical schemes for simulation. Numerical results are presented to confirm part of the analysis and to explore the difference between the two models on coarsening dynamics.

In this paper we study the effects of small viscosity term and the far-field boundary conditions for systems of convection-diffusion equations in the zero viscosity limit. The far-field boundary conditions are classified and the corresponding solution structures are analyzed. It is confirmed that the Neumann type of far-field boundary condition is preferred. On the other hand, we also identify a class of improperly coupled boundary conditions which lead to catastrophic reflection waves dominating the inlet in the zero viscosity limit. The analysis is performed on the linearized convection-diffusion model which well describes the behavior at the far field for many physical and engineering systems such as fluid dynamical equations and electro-magnetic equations. The results obtained here should provide some theoretical guidance for designing effective far field boundary conditions.

In this paper, we give a systematic study of the boundary layer behavior for linear convection-diffusion equation in the zero viscosity limit. We analyze the boundary layer structures in the viscous solution and derive the boundary condition satisfied by the viscosity limit as a solution of the inviscid equation. The results confirm that the Neumann type of far-field boundary condition is preferred in the outlet and characteristic boundary condition. Under some appropriate regularity and compatibility conditions on the initial and boundary data, we obtain optimal error estimates between the full viscous solution and the inviscid solution with suitable boundary layer corrections. These results hold in arbitrary space dimensions and similar statements also hold for the strip problem. This model well describes the behavior at the far-field for many physical and engineering systems such as fluid dynamical equation and electro-magnetic equation. The results obtained here should provide some theoretical guidance for designing effective far-field boundary conditions.

We study the detailed process of bifurcation in the flow’s topological structure for a two-dimensional (2-D) incompressible flow subject to no-slip boundary conditions and its connection with boundary-layer separation. The boundary-layer separation theory of M. Ghil, T. Ma and S. Wang, based on the structural-bifurcation concept, is translated into vorticity form. The vorticity formulation of the theory shows that structural bifurcation occurs whenever a degenerate singular point for the vorticity appears on the boundary; this singular point is characterized by nonzero tangential second-order derivative and nonzero time derivative of the vorticity. Furthermore, we prove the presence of an adverse pressure gradient at the critical point, due to reversal in the direction of the pressure force with respect to the basic shear flow at this point. A numerical example of 2-D driven-cavity flow, governed by the Navier Stokes equations, is presented; boundary-layer separation occurs, the bifurcation criterion is satisfied, and an adverse pressure gradient is shown to be present.

We study a continuum model for epitaxial growth of thin films in which the slope of mound structure of film surface increases. This model is a diffusion equation for the surface height profile h which is assumed to satisfy the periodic boundary condition. The equation happens to possess a Liapunov or �free-energy� functional. This functional consists of the term |? h|2, which represents the surface diffusion, and - log (1 + |? h|2), which describes the effect of kinetic asymmetry in the adatom attachment-detachment. We first prove for large time t that the interface width---the standard deviation of the height profile---is bounded above by O(t1/2), the averaged gradient is bounded above by O(t1/4), and the averaged energy is bounded below by O(- log t). We then consider a small coefficient e2 of |? h|2 with e = 1/L and L the linear size of the underlying system, and study the energy asymptotics in the large system limit e ? 0. We show that global minimizers of the free-energy functional exist for each e>0, the L2-norm of the gradient of any global minimizer scales as O(1/e), and the global minimum energy scales as O( log e). The existence of global energy minimizers and a scaling argument are used to construct a sequence of equilibrium solutions with different wavelengths. Finally, we apply our minimum energy estimates to derive bounds in terms of the linear system size L for the saturation interface width and the corresponding saturation time.

A reformulation of the planetary geostrophic equations (PGEs) with the inviscid balance equation is proposed and the existence of global weak solutions is established, provided that the mechanical force satisfies an integral constraint. There is only one prognostic equation for the temperature field, and the velocity field is statically determined by the planetary geostrophic balance combined with the incompressibility condition. Furthermore, the velocity profile can be accurately represented as a function of the temperature gradient. In particular, the vertical velocity depends only on the first-order derivative of the temperature. As a result, the bound for the L ∞ (0, t 1 ; L 2 ) ∩ L 2 (0, t 1 ; H 1 ) norm of the temperature field is sufficient to show the existence of the weak solution.

We show that in bounded domains with no-slip boundary conditions, the Navier-Stokes pressure can be determined in a such way that it is strictly dominated by viscosity. As a consequence, in a general domain we can treat the Navier-Stokes equations as a perturbed vector diffusion equation, instead of as a perturbed Stokes system. To illustrate the advantages of this view, we provide a simple proof of the unconditional stability of a difference scheme that is implicit only in viscosity and explicit in both pressure and convection terms, requiring no solution of stationary Stokes systems or inf-sup conditions.

This paper presents a general methodology to design macroscopic fluid models that take into account localized kinetic upscaling effects. The fluid models are solved in the whole domain together with a localized kinetic upscaling that corrects the fluid model wherever it is necessary. This upscaling is obtained by solving a kinetic equation on the nonequilibrium part of the distri-bution function. This equation is solved only locally and is related to the fluid equation through a downscaling effect. The method does not need to find an interface condition as do usual domain decomposition methods to match fluid and kinetic representations. We show our approach applies to problems that have a hydrodynamic time scale as well as to problems with diffusion time scale. Simple numerical schemes are proposed to discretize our models, and several numerical examples are used to validate the method.

We show that in a tubular domain with sufficiently small width, the normal and tangential gradients of a harmonic function have almost the same L 2 norm. This estimate yields a sharp estimate of the pressure in terms of the viscosity term in the Navier-Stokes equation with no-slip boundary condition. By consequence, one can analyze the Navier-Stokes equations simply as a perturbed vector diffusion equation instead of as a perturbed Stokes system. As an application, we describe a rather easy approach to establish a new isomorphism theorem for the non-homogeneous Stokes system.

The first part of this paper describes some important underlying themes in the mathematical theory of continuum mechanics that are distinct from formulating and analyzing governing equations. The main part of this paper is devoted to a survey of some concrete, conceptually simple, pretty problems that help illuminate the underlying themes. The paper concludes with a discussion of the crucial role of invariant constitutive equations in computation.

For strong solutions of the incompressible Navier-Stokes equations in bounded domains with velocity specified at the boundary, we establish the unconditional stability and convergence of discretization schemes that decouple the updates of pressure and velocity through explicit time stepping for pressure. These schemes require no solution of stationary Stokes systems, nor any compatibility between velocity and pressure spaces to ensure an inf-sup condition, and are representative of a class of highly efficient computational methods that have recently emerged. The proofs are simple, based upon a new, sharp estimate for the commutator of the Laplacian and Helmholtz projection operators. This allows us to treat an unconstrained formulation of the Navier-Stokes equations as a perturbed diffusion equation. © 2007 Wiley Periodicals, Inc.

This article studies the structure and its evolution of incompressible flows with the anti-symmetry using a combination of rigorous analysis and numerical simulations, with an application to an example of oceanic flow. In particular, necessary and sufficient conditions for 2D divergence-free vector fields with anti-symmetry are obtained, and a detailed numerical simulation for a simplified model of Marsigli oceanic flow is provided to explore and verify the structure and its transitions of the flow. It is expected that the study will lead to useful insights to the understanding of the flow dynamics from both the mathematical and physical points of view.

We consider the vorticity-stream formulation of axisymmetric incompressible ows and its equivalence with the primitive formulation. It is shown that, to characterize the regularity of a divergence free axisymmetric vector eld in terms of the swirling components, an extra set of pole condition is necessary to give a full description of the regularity. In addition, smooth solutions up to the axis of rotation gives rise to smooth solutions of primitive formulation in the case of Navier-Stokes equations, but not the Euler equations. We also establish proper weak formulations and show its equivalence to Leray's solutions.

We present a simple proof on the formation of flocking in the F. Cucker and S. Smale system [Jpn. J. Math. (3) 2, No. 1, 197–227 (2007; Zbl 1166.92323)] based on the explicit construction of a Lyapunov functional. Our results also provide a unified condition on the initial states in which the exponential convergence to the flocking state will occur. For large particle systems, we give a rigorous justification for the mean-field limit from the many particle Cucker-Smale system to the Vlasov equation with flocking dissipation as the number of particles goes to infinity.

We study a semi-implicit time-difference scheme for magnetohydrodynamics of a viscous and resistive incompressible fluid in a bounded smooth domain with a perfectly conducting boundary. In the scheme, the velocity and magnetic fields are updated by solving simple Helmholtz equations. Pressure is treated explicitly in time, by solving Poisson equations corresponding to a recently de- veloped formula for the Navier-Stokes pressure involving the commutator of Laplacian and Leray projection operators. We prove stability of the time-difference scheme, and deduce a local-time well- posedness theorem for MHD dynamics extended to ignore the divergence-free constraint on velocity and magnetic fields. These fields are divergence-free for all later time if they are initially so.

Nous considérons un modèle provenant de la biologie, composé d03équations de chimiotactisme couplées aux équations de fluide visqueux incompressible par le transport et le for04age externe. L03existence globale des solutions du problème de Cauchy est étudiée sous certaines conditions. Précisément, pour le système chimiotactisme–Navier–Stokes en deux dimensions d03espace, nous obtenons l03existence globale pour des données grandes. En trois dimensions d03espace, nous démontrons l03existence globale des solutions faibles pour le système chimiotactisme–Stokes avec une diffusion non-linéaire de la densité des cellules.

The purpose of this paper is to provide global existence and uniqueness results for a family of fluid transport equations by establishing convergence results for the particle method applied to these equations. The considered family of PDEs is a collection of strongly nonlinear equations which yield traveling wave solutions and can be used to model a variety of flows in fluid dynamics. We apply a particle method to the studied evolutionary equations and provide a new self-contained method for proving its convergence. The latter is accomplished by using the concept of space-time bounded variation and the associated compactness properties. From this result, we prove the existence of a unique global weak solution in some special cases and obtain stronger regularity properties of the solution than previously established.

The existence of solution for the 2D-Keller-Segel system in the subcritical case, i.e. when the initial mass is less than 8π, is reproved. Instead of using the entropy in the free energy and free energy dissipation, which was used in the proofs (Blanchet et al. in SIAM J. Numer. Anal. 46:691–721, 2008; Electron. J. Differ. Equ. Conf. 44:32, 2006 (electronic)), the potential energy term is fully utilized by adapting Delort’s theory on 2D incompressible Euler equation (Delort in J. Am. Math. Soc. 4:553–386, 1991).

Motivated by a phenomenon of phase transition in a model of alignment of self-propelled particles, we obtain a kinetic mean-field equation which is nothing else than the Doi equation (also called Smoluchowski equation) with dipolar potential. In a self-contained article, using only basic tools, we analyze the dynamics of this equation in any dimension. We first prove global well-posedness of this equation, starting with an initial condition in any Sobolev space. We then compute all possible steady-states. There is a threshold for the noise parameter: over this threshold, the only equilibrium is the uniform distribution, and under this threshold, there is also a family of non-isotropic equilibria. We give a rigorous prove of convergence of the solution to a steady-state as time goes to infinity. In particular we show that in the supercritical case, the only initial conditions leading to the uniform distribution in large time are those with vanishing momentum. For any positive value of the noise parameter, and any initial condition, we give rates of convergence towards equilibrium, exponentially for both supercritical and subcritical cases and algebraically for the critical case.

This paper deals with a degenerate diffusion Patlak–Keller–Segel system in $n\geq 3$ dimension. The main difference between the current work and many other recent studies on the same model is that we study the diffusion exponent $m=2n/(n+2)$, which is smaller than the usual exponent $m^{*}=2-2/n$ used in other studies. With the exponent $m=2n/(n+2)$, the associated free energy is conformal invariant, and there is a family of stationary solutions $U_{\lambda,x_0}(x)=C(\frac{\lambda} {\lambda^2+|x-x_0|^2})^{\frac{n+2}{2}}$ $\forall \lambda>0$, $x_0\in {\mathbb R}^n$. For radially symmetric solutions, we prove that if the initial data are strictly below $U_{\lambda,0}(x)$ for some $\lambda$, then the solution vanishes in $L^1_{loc}$ as $t\to\infty$; if the initial data are strictly above $U_{\lambda,0}(x)$ for some $\lambda$, then the solution either blows up at a finite time or has a mass concentration at $r=0$ as time goes to infinity. For general initial data, we prove that there is a global weak solution provided that the $L^m$ norm of initial density is less than a universal constant, and the weak solution vanishes as time goes to infinity. We also prove a finite time blow-up of the solution if the $L^m$ norm for initial data is larger than the $L^m$ norm of $U_{\lambda,x_0}(x)$, which is constant independent of $\lambda$ and $x_0$, and the free energy of initial data is smaller than that of $U_{\lambda,x_0}(x)$.

We establish two nonlinear compactness theorems in Lp(0,T;B) with hypothesis on time translations, which are nonlinear counterparts of two results by Simon (1987) [1]. The first theorem sharpens a result by Maitre (2003) [10] and is important in the study of doubly nonlinear elliptic–parabolic equations. Based on this theorem, we then obtain a time translation counterpart of a result by Dubinskii˘ (1965) [5], which is supposed to be useful in the study of some nonlinear kinetic equations (e.g. the FENE-type bead–spring chains model).