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.     

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.     

Besov Space Regularity Conditions for Weak Solutions of the Navier–Stokes Equations

Consider a bounded domain ${{\Omega \subseteq \mathbb{R}^3}}$ with smooth boundary, some initial value ${{u_0 \in L^2_{\sigma}(\Omega )}}$ , and a weak solution u of the Navier–Stokes system in ${{[0,T) \times\Omega,\,0 < T \le \infty}}$ . Our aim is to develop regularity and uniqueness conditions for u which are based on the Besov space $$B^{q,s}(\Omega ):=\left\{v\in L^2_{\sigma}(\Omega ); \|v\|_{B^{q,s}(\Omega )} := \left(\int\limits^{\infty}_0 \left\|e^{-\tau A}v\right\|^s_q {\rm d} \tau\right)^{1/s}<\infty \right\}$$ with ${{2 < s < \infty,\,3 < q <\infty,\,\frac2{s}+\frac{3}{q} = 1}}$ ; here A denotes the Stokes operator. This space, introduced by Farwig et al. (Ann. Univ. Ferrara 55:89–110, 2009 and J. Math. Fluid Mech. 14: 529–540, 2012), is a subspace of the well known Besov space ${{{\mathbb{B}}^{-2/s}_{q,s}(\Omega )}}$ , see Amann (Nonhomogeneous Navier–Stokes Equations with Integrable Low-Regularity Data. Int. Math. Ser. pp. 1–28. Kluwer/Plenum, New York, 2002). Our main results on the regularity of u exploits a variant of the space ${{B^{q,s}(\Omega )}}$ in which the integral in time has to be considered only on finite intervals (0, δ ) with ${{\delta \to 0}}$ . Further we discuss several criteria for uniqueness and local right-hand regularity, in particular, if u satisfies Serrin’s limit condition ${{u\in L^{\infty}_{\text{loc}}([0,T);L^3_{\sigma}(\Omega ))}}$ . Finally, we obtain a large class of regular weak solutions u defined by a smallness condition ${{\|u_0\|_{B^{q,s}(\Omega )} \le K}}$ with some constant ${{K=K(\Omega, q)>0}}$ .     

Multistage γ-level \mathcal{H}_{\infty} control for input-saturated systems with disturbances

For input-saturated systems with disturbances, states in the domain of attraction cannot converge to the origin, but only to neighborhood around it. In order to design the smallest possible target invariant set and the largest possible domain of attraction, in this paper, we introduce a multistage γ-level $\mathcal{H}_{\infty}$ control for achieving a smaller target invariant set within a given $\mathcal{H}_{\infty}$ performance level and a larger domain of attraction than results obtained in previous studies. In particular, for the case in which the disturbances satisfy a matched condition, this paper introduces an $\mathcal{H}_{\infty}$ control with an extra control part to perfectly reject these disturbances despite the uncertainties; the introduction of the $\mathcal{H}_{\infty}$ control with an extra control part causes the target invariant set to shrink to the origin and the $\mathcal{H}_{\infty}$ performance level to become zero.     

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.     

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^*}$ .     

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.     

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.     

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 the Stability in Weak Topology of the Set of Global Solutions to the Navier–Stokes Equations

Let X be a suitable function space and let ${\mathcal{G} \subset X}$ be the set of divergence free vector fields generating a global, smooth solution to the incompressible, homogeneous three-dimensional Navier–Stokes equations. We prove that a sequence of divergence free vector fields converging in the sense of distributions to an element of ${\mathcal{G}}$ belongs to ${\mathcal{G}}$ if n is large enough, provided the convergence holds “anisotropically” in frequency space. Typically, this excludes self-similar type convergence. Anisotropy appears as an important qualitative feature in the analysis of the Navier–Stokes equations; it is also shown that initial data which do not belong to ${\mathcal{G}}$ (hence which produce a solution blowing up in finite time) cannot have a strong anisotropy in their frequency support.     

Well-Posedness of Nematic Liquid Crystal Flow in {L^{3}_{\rm uloc}(\mathbb{R}^3)}

In this paper, we establish the local well-posedness for the Cauchy problem of a simplified version of hydrodynamic flow of nematic liquid crystals in ${\mathbb{R}^3}$ for any initial data (u ₀, d ₀) having small ${L^{3}_{\rm uloc}}$ -norm of ${(u_{0}, \nabla d_{0})}$ . Here ${L^{3}_{\rm uloc}(\mathbb{R}^3)}$ is the space of uniformly locally L ³-integrable functions. For any initial data (u ₀, d ₀) with small ${\|(u_0, \nabla d_0)\|_{L^{3}(\mathbb{R}^3)}}$ , we show that there exists a unique, global solution to the problem under consideration which is smooth for t > 0 and has monotone deceasing L ³-energy for ${t \geqq 0}$ .     

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

Experimental Characterization of Porosity Structure and Transport Property Changes in Limestone Undergoing Different Dissolution Regimes

Limestone dissolution by $\hbox {CO}_2$ -rich brine induces critical changes of the pore network geometrical parameters such as the pore size distribution, the connectivity, and the tortuosity which govern the macroscopic transport properties (permeability and dispersivity) that are required to parameterize the models, simulating the injection and the fate of $\hbox {CO}_2$ . A set of four reactive core-flood experiments reproducing underground conditions ( $T = 100\,^{\circ }\hbox {C}$ and $P = 12$ MPa) has been conducted for different $\hbox {CO}_2$ partial pressures $(0.034 < P_{\mathrm{CO}_2}< 3.4\; \hbox {MPa})$ in order to study the different dissolution regimes. X-ray microtomographic images have been used to characterize the changes in the structural properties from pore scale to Darcy scale, while time-resolved pressure loss and chemical fluxes enabled the determination of the sample-scale change in porosity and permeability. The results show the growth of localized dissolution features associated with high permeability increase for the highest $P_{\mathrm{CO}_2}$ , whereas dissolution tends to be more homogeneously distributed for lower values of $P_{\mathrm{CO}_2}$ . For the latter, the higher the $P_{\mathrm{CO}_2}$ , the more the dissolution patterns display ramified structures and permeability increase. For the lowest value of $P_{\mathrm{CO}_2}$ , the preferential dissolution of the calcite cement associated with the low dissolution kinetics triggers the transport that may locally accumulate and form a microporous material that alters permeability and produces an anti-correlated porosity–permeability relationship. The combined analysis of the pore network geometry and the macroscopic measurements shows that $P_{\mathrm{CO}_2}$ regulates the tortuosity change during dissolution. Conversely, the increase of the exponent value of the observed power law permeability–porosity trend while $P_{\mathrm{CO}_2}$ increases, which appears to be strongly linked to the increase of the effective hydraulic diameter, depends on the initial rock structure.     

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}$ .     

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}$ .