Regularity and Asymptotic Behavior of Nonlinear Stefan Problems

We study the following nonlinear Stefan problem $$\left\{\begin{aligned}\!\!&u_t\,-\,d\Delta u = g(u) & &\quad{\rm for}\,x\,\in\,\Omega(t), t > 0, \\ & u = 0 \, {\rm and} u_t = \mu|\nabla_{x} u|^{2} &&\quad {\rm for}\,x\,\in\,\Gamma(t), t > 0, \\ &u(0, x) = u_{0}(x) &&\quad {\rm for}\,x\,\in\,\Omega_0,\end{aligned} \right.$$ where ${\Omega(t) \subset \mathbb{R}^{n}}$ ( ${n \geqq 2}$ ) is bounded by the free boundary ${\Gamma(t)}$ , with ${\Omega(0) = \Omega_0}$ , μ and d are given positive constants. The initial function u ₀ is positive in ${\Omega_0}$ and vanishes on ${\partial \Omega_0}$ . The class of nonlinear functions g(u) includes the standard monostable, bistable and combustion type nonlinearities. We show that the free boundary ${\Gamma(t)}$ is smooth outside the closed convex hull of ${\Omega_0}$ , and as ${t \to \infty}$ , either ${\Omega(t)}$ expands to the entire ${\mathbb{R}^n}$ , or it stays bounded. Moreover, in the former case, ${\Gamma(t)}$ converges to the unit sphere when normalized, and in the latter case, ${u \to 0}$ uniformly. When ${g(u) = au - bu^2}$ , we further prove that in the case ${\Omega(t)}$ expands to ${{\mathbb R}^n}$ , ${u \to a/b}$ as ${t \to \infty}$ , and the spreading speed of the free boundary converges to a positive constant; moreover, there exists ${\mu^* \geqq 0}$ such that ${\Omega(t)}$ expands to ${{\mathbb{R}}^n}$ exactly when ${\mu > \mu^*}$ .     

Hadamard Variational Formula for the Green’s Function of the Boundary Value Problem on the Stokes Equations

For every ${\varepsilon > 0}$ , we consider the Green’s matrix ${G_{\varepsilon}(x, y)}$ of the Stokes equations describing the motion of incompressible fluids in a bounded domain ${\Omega_{\varepsilon} \subset \mathbb{R}^d}$ , which is a family of perturbation of domains from ${\Omega\equiv \Omega_0}$ with the smooth boundary ${\partial\Omega}$ . Assuming the volume preserving property, that is, ${\mbox{vol.}\Omega_{\varepsilon} = \mbox{vol.}\Omega}$ for all ${\varepsilon > 0}$ , we give an explicit representation formula for ${\delta G(x, y) \equiv \lim_{\varepsilon\to +0}\varepsilon^{-1}(G_{\varepsilon}(x, y) - G_0(x, y))}$ in terms of the boundary integral on ${\partial \Omega}$ of ${G_0(x, y)}$ . Our result may be regarded as a classical Hadamard variational formula for the Green’s functions of the elliptic boundary value problems.     

Global Compensated Compactness Theorem for General Differential Operators of First Order

Let A ₁(x, D) and A ₂(x, D) be differential operators of the first order acting on l-vector functions ${u= (u_1, \ldots, u_l)}$ in a bounded domain ${\Omega \subset \mathbb{R}^{n}}$ with the smooth boundary ${\partial\Omega}$ . We assume that the H ¹-norm ${\|u\|_{H^{1}(\Omega)}}$ is equivalent to ${\sum_{i=1}^2\|A_iu\|_{L^2(\Omega)} + \|B_1u\|_{H^{\frac{1}{2}}(\partial\Omega)}}$ and ${\sum_{i=1}^2\|A_iu\|_{L^2(\Omega)} + \|B_2u\|_{H^{\frac{1}{2}}(\partial\Omega)}}$ , where B _i  = B _i (x, ν) is the trace operator onto ${\partial\Omega}$ associated with A _i (x, D) for i = 1, 2 which is determined by the Stokes integral formula (ν: unit outer normal to ${\partial\Omega}$ ). Furthermore, we impose on A ₁ and A ₂ a cancellation property such as ${A_1A_2^{\prime}=0}$ and ${A_2A_1^{\prime}=0}$ , where ${A^{\prime}_i}$ is the formal adjoint differential operator of A _i (i = 1, 2). Suppose that ${\{u_m\}_{m=1}^{\infty}}$ and ${\{v_m\}_{m=1}^{\infty}}$ converge to u and v weakly in ${L^2(\Omega)}$ , respectively. Assume also that ${\{A_{1}u_m\}_{m=1}^{\infty}}$ and ${\{A_{2}v_{m}\}_{m=1}^{\infty}}$ are bounded in ${L^{2}(\Omega)}$ . If either ${\{B_{1}u_m\}_{m=1}^{\infty}}$ or ${\{B_{2}v_m\}_{m=1}^{\infty}}$ is bounded in ${H^{\frac{1}{2}}(\partial\Omega)}$ , then it holds that ${\int_{\Omega}u_m\cdot v_m \,{\rm d}x \to \int_{\Omega}u\cdot v \,{\rm d}x}$ . We also discuss a corresponding result on compact Riemannian manifolds with boundary.     

Maximal and Minimal Spreading Speeds for Reaction Diffusion Equations in Nonperiodic Slowly Varying Media

This paper investigates the asymptotic behavior of the solutions of the Fisher-KPP equation in a heterogeneous medium, $$\partial_t u = \partial_{xx} u + f(x,u),$$ associated with a compactly supported initial datum. A typical nonlinearity we consider is ${f(x,u) = \mu_0 (\phi (x)) u(1-u)}$ , where??? ₀ is a 1-periodic function and ${\phi}$ is a ${\mathcal{C}^1}$ increasing function that satisfies ${\lim_{x \to+\infty}\phi (x) = +\infty}$ and ${\lim_{x \to +\infty}\phi' (x) =0}$ . Although quite specific, the choice of such a reaction term is motivated by its highly heterogeneous nature. We exhibit two different behaviors for u for large times, depending on the speed of the convergence of ${\phi}$ at infinity. If ${\phi}$ grows sufficiently slowly, then we prove that the spreading speed of u oscillates between two distinct values. If ${\phi}$ grows rapidly, then we compute explicitly a unique and well determined speed of propagation w _??, arising from the limiting problem of an infinite period. We give a heuristic interpretation for these two behaviors.     

On Formation of a Locally Self-Similar Collapse in the Incompressible Euler Equations

The paper addresses the question of the existence of a locally self-similar blow-up for the incompressible Euler equations. Several exclusion results are proved based on the L ^p -condition for velocity or vorticity and for a range of scaling exponents. In particular, in N dimensions if in self-similar variables ${u \in L^p}$ and ${u \sim \frac{1}{t^{\alpha/(1+\alpha)}}}$ , then the blow-up does not occur, provided ${\alpha > N/2}$ or ${-1 < \alpha \leq N\,/p}$ . This includes the L ³ case natural for the Navier–Stokes equations. For ${\alpha = N\,/2}$ we exclude profiles with asymptotic power bounds of the form ${ |y|^{-N-1+\delta} \lesssim |u(y)| \lesssim |y|^{1-\delta}}$ . Solutions homogeneous near infinity are eliminated, as well, except when homogeneity is scaling invariant.     

Regularization of Point Vortices Pairs for the Euler Equation in Dimension Two

In this paper, we construct stationary classical solutions of the incompressible Euler equation approximating singular stationary solutions of this equation. This procedure is carried out by constructing solutions to the following elliptic problem $$\left\{\begin{array}{l@{\quad}l} -\varepsilon^2 \Delta u = \sum\limits_{i=1}^m \chi_{\Omega_i^{+}} \left(u - q - \frac{\kappa_i^{+}}{2\pi} {\rm ln} \frac{1}{\varepsilon}\right)_+^p\\ \quad - \sum_{j=1}^n \chi_{\Omega_j^{-}} \left(q - \frac{\kappa_j^{-}}{2\pi} {\rm \ln} \frac{1}{\varepsilon} - u\right)_+^p , \quad \quad x \in \Omega,\\ u = 0, \quad \quad \quad \quad \quad \quad \quad \quad \quad \quad \quad \quad \quad \quad x \in \partial \Omega,\end{array}\right.$$ where p > 1, ${\Omega \subset \mathbb{R}^2}$ is a bounded domain, ${\Omega_i^{+}}$ and ${\Omega_j^{-}}$ are mutually disjoint subdomains of Ω and ${\chi_{\Omega_i^{+}} ({\rm resp}.\; \chi_{\Omega_j^{-}})}$ are characteristic functions of ${\Omega_i^{+}({\rm resp}. \;\Omega_j^{-}})$ , q is a harmonic function. We show that if Ω is a simply-connected smooth domain, then for any given C ¹-stable critical point of Kirchhoff–Routh function ${\mathcal{W}\;(x_1^{+},\ldots, x_m^{+}, x_1^{-}, \ldots, x_n^{-})}$ with ${\kappa^{+}_i > 0\,(i = 1,\ldots, m)}$ and ${\kappa^{-}_j > 0\,(j = 1,\ldots,n)}$ , there is a stationary classical solution approximating stationary m + n points vortex solution of incompressible Euler equations with total vorticity ${\sum_{i=1}^m \kappa^{+}_i -\sum_{j=1}^n \kappa_j^{-}}$ . The case that n = 0 can be dealt with in the same way as well by taking each ${\Omega_j^{-}}$ as an empty set and set ${\chi_{\Omega_j^{-}} \equiv 0,\,\kappa^{-}_j=0}$ .     

Resistance Coefficients of Polymer Membrane with Concentration Polarization

The Kedem-Katchalsky equations, modified by means of symmetric transformations of Peusner thermodynamic networks, were applied to interpret the membrane transport in concentration polarization conditions. The results from the study demonstrate that the resistance coefficients counted for membrane transport of aqueous solutions of glucose through Nephrophan membrane in horizontal plane are nonlinearly dependent on mean concentration of glucose in the membrane ${(\bar{C})}$ . It was also shown that the threshold value of concentration ${(\bar{C}_{cr})}$ existed, and for ${\bar{C} > \bar{C}_{cr}}$ , the resistance coefficients depend, while for ${\bar{C} < \bar{C}_{cr}}$ , they do not depend on the membrane system configuration. Increase of mean glucose concentration in the membrane (in the range ${\bar{C} > \bar{C}_{cr})}$ causes decrease of difference between resistance coefficients of the membrane system in homogeneous conditions (solutions mechanically stirred) and in conditions with hydrodynamic instabilities (configuration B). Besides increase of mean glucose concentration in the membrane (in the range ${\bar{C} > \bar{C}_{cr})}$ causes increase of the difference between resistance coefficients for membrane system with concentration polarization without hydrodynamic instabilities (configuration A) and membrane system in homogeneous conditions.     

Elongational rheology of NIPAM-based hydrogels

Hydrogels of different composition based on the copolymerization of N-isopropyl acrylamide and surfmers of different chemical structure were tested in elongation using Hencky/real definitions for stress, strain, and strain rate, offering a more scientific insight into the effect of deformation on the properties. In a range between $\dot {\varepsilon }=10$ and 0.01 s $^{-1}$ , the material properties are independent of strain rate and show a very clear strain hardening with a “brittle” sudden fracture. The addition of surfmer increases the strain at break $\varepsilon _{\mathrm {H}}^{\max }$ and at the same time leads to a failure of hyperelastic models. The samples can be stretched up to Hencky strains $\varepsilon _{\mathrm {H}}^{\max }$ between 0.6 and 2.5, depending on the molecular structure, yielding linear Young’s moduli E $_{0}$ between 2,700 and 39,000 Pa. The strain-rate independence indicates an ideal rubberlike behavior and fracture in a brittle-like fashion. The resulting stress at break $\sigma _{\textrm max}$ can be correlated with $\varepsilon _{\mathrm {H}}^{\max } $ and $E_{0}$ as well as with the solid molar mass between the cross-linking points $M_{\mathrm {c}}^{\textrm {solids}} $ , derived from $E_{0}$ .     

Stochastic Variational Inequalities and Applications to the Total Variation Flow Perturbed by Linear Multiplicative Noise

In this work, we introduce a new method to prove the existence and uniqueness of a variational solution to the stochastic nonlinear diffusion equation ${{\rm d}X(t) = {\rm div} \left[\frac{\nabla X(t)}{|\nabla X(t)|}\right]{\rm d}t + X(t){\rm d}W(t) {\rm in} (0, \infty) \times \mathcal{O},}$ where ${\mathcal{O}}$ is a bounded and open domain in ${\mathbb{R}^N, N \geqq 1}$ and W(t) is a Wiener process of the form ${W(t) = \sum^{\infty}_{k = 1}\mu_{k}e_{k}\beta_{k}(t), e_{k} \in C^{2}(\overline{\mathcal{O}}) \cap H^{1}_{0}(\mathcal{O}),}$ and ${\beta_{k}, k \in \mathbb{N}}$ are independent Brownian motions. This is a stochastic diffusion equation with a highly singular diffusivity term. One main result established here is that for all initial conditions in ${L^2(\mathcal{O})}$ , it is well posed in a class of continuous solutions to the corresponding stochastic variational inequality. Thus, one obtains a stochastic version of the (minimal) total variation flow. The new approach developed here also allows us to prove the finite time extinction of solutions in dimensions ${1\leqq N \leqq3}$ , which is another main result of this work.     

On the effect of non-smooth Coulomb damping on flutter-type self-excitation in a non-gyroscopic circulatory 2-DoF-system

This article deals with self-excited vibrations, attractivity of stationary solutions, and the corresponding bifurcation behavior of two-dimensional differential inclusions of the type $\mathbf{M}\mathbf{q}'' + \mathbf{D}\mathbf{q}' + (\mathbf{K} + \bar{\mu}\mathbf{N})\mathbf{q} \in-\mathbf{R}\operatorname{Sign}(\mathbf{q}')$ . For the smooth case R=0, the equilibrium may become unstable due to non-conservative positional forces stemming from the circulatory matrix N. This type of instability is usually referred to as flutter instability and the loss of stability is related to a Hopf bifurcation of the steady state, which occurs for a critical parameter $\bar{\mu}= \bar{\mu}_{\mathrm{crit}}$ . For R≠0, the steady state is a set of equilibria, which turns out to be attractive for all values of the bifurcation parameter $\bar{\mu}$ . Depending on $\bar{\mu}$ , the basin of attraction of the equilibrium set can be infinite or finite. The transition from an infinite to a finite basin of attraction occurs at the stability threshold $\bar{\mu}_{\mathrm{crit}}$ of the underlying smooth problem. For the finite basin of attraction, its size is proportional to the Coulomb friction and inverse-proportional to $(\bar{\mu}- \bar{\mu}_{\mathrm{crit}})$ . By adding Coulomb damping the notion of steady state stability for the smooth problem is replaced by the question whether the basin of attraction of the steady state is infinite or finite. Simultaneously, the local Hopf-bifurcation is replaced by a global bifurcation. This implies that in the presence of Coulomb damping the occurrence of self-excited vibrations can only be investigated with regard to the perturbation level.     

Generalized Resolvent Estimates of the Stokes Equations with First Order Boundary Condition in a General Domain

In this paper, we prove unique existence of solutions to the generalized resolvent problem of the Stokes operator with first order boundary condition in a general domain ${\Omega}$ of the N-dimensional Eulidean space ${\mathbb{R}^N, N \geq 2}$ . This type of problem arises in the mathematical study of the flow of a viscous incompressible one-phase fluid with free surface. Moreover, we prove uniform estimates of solutions with respect to resolvent parameter ${\lambda}$ varying in a sector ${\Sigma_{\sigma, \lambda_0} = \{\lambda \in \mathbb{C} \mid |\arg \lambda| < \pi-\sigma, \enskip |\lambda| \geq \lambda_0\}}$ , where ${0 < \sigma < \pi/2}$ and ${\lambda_0 \geq 1}$ . The essential assumption of this paper is the existence of a unique solution to a suitable weak Dirichlet problem, namely it is assumed the unique existence of solution ${p \in \hat{W}^1_{q, \Gamma}(\Omega)}$ to the variational problem: ${(\nabla p, \nabla \varphi) = (f, \nabla \varphi)}$ for any ${\varphi \in \hat W^1_{q', \Gamma}(\Omega)}$ . Here, ${1 < q < \infty, q' = q/(q-1), \hat W^1_{q, \Gamma}(\Omega)}$ is the closure of ${W^1_{q, \Gamma}(\Omega) = \{ p \in W^1_q(\Omega) \mid p|_\Gamma = 0\}}$ by the semi-norm ${\|\nabla \cdot \|_{L_q(\Omega)}}$ , and ${\Gamma}$ is the boundary of ${\Omega}$ . In fact, we show that the unique solvability of such a Dirichlet problem is necessary for the unique existence of a solution to the resolvent problem with uniform estimate with respect to resolvent parameter varying in ${(\lambda_0, \infty)}$ . Our assumption is satisfied for any ${q \in (1, \infty)}$ by the following domains: whole space, half space, layer, bounded domains, exterior domains, perturbed half space, perturbed layer, but for a general domain, we do not know any result about the unique existence of solutions to the weak Dirichlet problem except for q =  2.     

Attractors and Chain Recurrence for Semigroup Actions

Let ${\mathcal{S}}$ be a semigroup acting on a topological space M. We define attractors for the action of ${\mathcal{S}}$ on M. This concept depends on a family ${\mathcal{F}}$ of subsets of ${\mathcal{S}}$ . For certain semigroups and families it recovers the concept of attractors for flows or semiflows. We define and study the complementary repeller of an attractor. We also characterize the set of chain recurrent points in terms of attractors.     

Random Data Cauchy Theory for the Generalized Incompressible Navier–Stokes Equations

In this paper, we consider the generalized Navier?CStokes equations where the space domain is ${\mathbb{T}^N}$ or ${\mathbb{R}^N, N\geq3}$ . The generalized Navier?CStokes equations here refer to the equations obtained by replacing the Laplacian in the classical Navier?CStokes equations by the more general operator (???) ^?? with ${\alpha\in (\frac{1}{2},\frac{N+2}{4})}$ . After a suitable randomization, we obtain the existence and uniqueness of the local mild solution for a large set of the initial data in ${H^s, s\in[-\alpha,0]}$ , if ${1 < \alpha < \frac{N+2}{4}, s\in(1-2\alpha,0]}$ , if ${\frac{1}{2} < \alpha\leq 1}$ . Furthermore, we obtain the probability for the global existence and uniqueness of the solution. Specially, our result shows that, in some sense, the Cauchy problem of the classical Navier?CStokes equation is local well-posed for a large set of the initial data in H ^?1+, exhibiting a gain of ${\frac{N}{2}-}$ derivatives with respect to the critical Hilbert space ${H^{\frac{N}{2}-1}}$ .     

On the Stokes and Oseen Problems with Singular Data

We consider the steady Stokes and Oseen problems in bounded and exterior domains of ${\mathbb{R}^n}$ of class C ^k-1,1 (n = 2, 3; k ≥ 2). We prove existence and uniqueness of a very weak solution for boundary data a in ${W^{2-k-1/q,q} (\partial\Omega)}$ . If ${\Omega}$ is of class ${C^\infty}$ , we can assume a to be a distribution on ${\partial\Omega}$ .     

Incompressible 3D Navier–Stokes Equations as a Limit of a Nonlinear Parabolic System

In this paper, we consider the Cauchy problem for a nonlinear parabolic system ${u^\epsilon_t - \Delta u^\epsilon + u^\epsilon \cdot \nabla u^\epsilon + \frac{1}{2}u^\epsilon\, {\rm div}\, u^\epsilon - \frac{1}{\epsilon}\nabla\, {\rm div}\, u^\epsilon = 0}$ in ${\mathbb {R}^3 \times (0,\infty)}$ with initial data in Lebesgue spaces ${L^2(\mathbb {R}^3)}$ or ${L^3(\mathbb {R}^3)}$ . We analyze the convergence of its solutions to a solution of the incompressible Navier?CStokes system as ${\epsilon \to 0}$ .     

Extensions of Serrin’s Uniqueness and Regularity Conditions for the Navier–Stokes Equations

Consider a smooth bounded domain ${\Omega \subseteq {\mathbb{R}}^3}$ , a time interval [0, T), 0?<?T?≤?∞, and a weak solution u of the Navier–Stokes system. Our aim is to develop several new sufficient conditions on u yielding uniqueness and/or regularity. Based on semigroup properties of the Stokes operator we obtain that the local left-hand Serrin condition for each ${t\in (0,T)}$ is sufficient for the regularity of u. Somehow optimal conditions are obtained in terms of Besov spaces. In particular we obtain such properties under the limiting Serrin condition ${u \in L_{\rm loc}^\infty([0,T);L^3(\Omega))}$ . The complete regularity under this condition has been shown recently for bounded domains using some additional assumptions in particular on the pressure. Our result avoids such assumptions but yields global uniqueness and the right-hand regularity at each time when ${u \in L_{\rm loc}^\infty([0,T);L^3(\Omega))}$ or when ${u(t)\in L^3(\Omega)}$ pointwise and u satisfies the energy equality. In the last section we obtain uniqueness and right-hand regularity for completely general domains.     

Drainage of Capillary-Trapped Oil by an Immiscible Gas: Impact of Transient and Steady-State Water Displacement on Three-Phase Oil Permeability

In a previous paper (Dehghanpour et al., Phys Rev E 83:065302, 2011a), we showed that relative permeability of mobilized oil, $k_\mathrm{ro}$ , measured during tertiary gravity drainage, is significantly higher than that of the same oil saturation in other tests where oil is initially a continuous phase. We also showed that tertiary $k_\mathrm{ro}$ strongly correlates to both water saturation, $S_\mathrm{w}$ , water flux (water relative permeability), $k_\mathrm{rw}$ , and the change in water saturation with time, $\mathrm{d}S_\mathrm{w}/\mathrm{d}t$ . To develop a model and understanding of the enhanced oil transport, identifying which of these parameters ( $S_\mathrm{w},\,k_{\mathrm{rw}}$ , or $\mathrm{d}S_\mathrm{w}/\mathrm{d}t$ ) plays the controlling role is necessary, but in the previous experiments these could not be deconvolved. To answer the remaining question, we conduct specific three-phase displacement experiments in which $k_{\mathrm{rw}}$ is controlled by applying a fixed water influx, and $S_\mathrm{w}$ develops naturally. We obtain $k_{\mathrm{ro}}$ by using the saturation data measured in time and space. The results suggest that steady-state water influx, in contrast to transient water displacement, does not enhance $k_{\mathrm{ro}}$ . Instead, reducing water influx rate results in excess oil flow. Furthermore, according to our pore scale hydraulic conductivity calculations, viscous coupling and fluid positioning do not sufficiently explain the observed correlation between $k_{\mathrm{ro}}$ and $S_{\mathrm{w}}$ . We conclude that tertiary $k_{\mathrm{ro}}$ is controlled by the oil mobilization rate, which in turn is linked to the rate of water saturation decrease with time, $\mathrm{d}S_\mathrm{w}/\mathrm{d}t$ . Finally, we develop a simple model which relates tertiary $k_{\mathrm{ro}}$ to transient two-phase gas/water relative permeability.     

On the Elastic Energy Density of Constrained Q-Tensor Models for Biaxial Nematics

Within the Landau–de Gennes theory, the order parameter describing a biaxial nematic liquid crystal assigns a symmetric traceless 3 × 3 matrix Q with three distinct eigenvalues to every point of the region Ω occupied by the system. In the constrained case of matrices Q with constant eigenvalues, the order parameter space is diffeomorphic to the eightfold quotient ${\mathbb{S}^3/\mathcal{H}}$ of the 3-sphere ${\mathbb{S}^3}$ , where ${\mathcal{H}}$ is the quaternion group, and a configuration of a biaxial nematic liquid crystal is described by a map from Ω to ${\mathbb{S}^3/\mathcal{H}}$ . We express the (simplest form of the) Landau–de Gennes elastic free-energy density as a density defined on maps ${q: \Omega \to \mathbb{S}^3}$ , whose functional dependence is restricted by the requirements that (1) it is well defined on the class of configuration maps from Ω to ${\mathbb{S}^3/\mathcal{H}}$ (residual symmetry) and (2) it is independent of arbitrary superposed rigid rotations (frame indifference). As an application of this representation, we then discuss some properties of the corresponding energy functional, including coercivity, lower semicontinuity and strong density of smooth maps. Other invariance properties are also considered. In the discussion, we take advantage of the identification of ${\mathbb{S}^3}$ with the Lie group of unit quaternions ${Sp(1) \cong SU(2)}$ and of the relations between quaternions and rotations in ${\mathbb{R}^3}$ and ${\mathbb{R}^4}$ .     

Global Regularity for the Two-Dimensional Anisotropic Boussinesq Equations with Vertical Dissipation

This paper establishes the global in time existence of classical solutions to the two-dimensional anisotropic Boussinesq equations with vertical dissipation. When only vertical dissipation is present, there is no direct control on the horizontal derivatives and the global regularity problem is very challenging. To solve this problem, we bound the derivatives in terms of the ${L^\infty}$ -norm of the vertical velocity v and prove that ${\|v\|_{L^{r}}}$ with ${2\leqq r < \infty}$ does not grow faster than ${\sqrt{r \log r}}$ at any time as r increases. A delicate interpolation inequality connecting ${\|v\|_{L^\infty}}$ and ${\|v\|_{L^r}}$ then yields the desired global regularity.