The Decay Rate and Higher Approximation of Mild Solutions to the Navier�CStokes Equations

In this paper, we consider v(t) = u(t) − e ^tΔ u ₀, where u(t) is the mild solution of the Navier–Stokes equations with the initial data u₀ ? L²(\mathbb Rⁿ)?Lⁿ(\mathbb Rⁿ){u_0\in L^2({\mathbb R}^n)\cap L^n({\mathbb R}^n)} . We shall show that the L ² norm of D ^β v(t) decays like t^{-\frac |b|-1 2-\frac n4}{t^{-\frac {|\beta|-1} {2}-\frac n4}} for |β| ≥ 0. Moreover, we will find the asymptotic profile u ₁(t) such that the L ² norm of D ^β (v(t) − u ₁(t)) decays faster for 3 ≤ n ≤ 5 and |β| ≥ 0. Besides, higher-order asymptotics of v(t) are deduced under some assumptions.     

On the Navier�CStokes equations with rotating effect and prescribed outflow velocity

We consider the equations of Navier–Stokes modeling viscous fluid flow past a moving or rotating obstacle in \mathbb R^d{\mathbb R^d} subject to a prescribed velocity condition at infinity. In contrast to previously known results, where the prescribed velocity vector is assumed to be parallel to the axis of rotation, in this paper we are interested in a general outflow velocity. In order to use L ^p -techniques we introduce a new coordinate system, in which we obtain a non-autonomous partial differential equation with an unbounded drift term. We prove that the linearized problem in \mathbb R^d{\mathbb R^d} is solved by an evolution system on L^p_s(\mathbb R^d){L^p_{\sigma}(\mathbb R^d)} for 1 < p < ∞. For this we use results about time-dependent Ornstein–Uhlenbeck operators. Finally, we prove, for p ≥ d and initial data u₀ ? L^p_s(\mathbb R^d){u_0\in L^p_{\sigma}(\mathbb R^d)}, the existence of a unique mild solution to the full Navier–Stokes system.     

Special Fast Diffusion with Slow Asymptotics: Entropy Method and Flow on a Riemannian Manifold

We consider the asymptotic behaviour of positive solutions u(t, x) of the fast diffusion equation ${u_t=\Delta (u^{m}/m)= {\rm div}\,(u^{m-1} \nabla u)}We consider the asymptotic behaviour of positive solutions u(t, x) of the fast diffusion equation u_t=D(u^m/m) = div (u^m-1 ?u){u_t=\Delta (u^{m}/m)= {\rm div}\,(u^{m-1} \nabla u)} posed for x ? \mathbb R^d{x\in\mathbb R^d}, t > 0, with a precise value for the exponent m = (d − 4)/(d − 2). The space dimension is d ≧ 3 so that m < 1, and even m = −1 for d = 3. This case had been left open in the general study (Blanchet et al. in Arch Rat Mech Anal 191:347–385, 2009) since it requires quite different functional analytic methods, due in particular to the absence of a spectral gap for the operator generating the linearized evolution. The linearization of this flow is interpreted here as the heat flow of the Laplace– Beltrami operator of a suitable Riemannian Manifold (\mathbb R^d,g){(\mathbb R^d,{\bf g})}, with a metric g which is conformal to the standard \mathbb R^d{\mathbb R^d} metric. Studying the pointwise heat kernel behaviour allows to prove suitable Gagliardo–Nirenberg inequalities associated with the generator. Such inequalities in turn allow one to study the nonlinear evolution as well, and to determine its asymptotics, which is identical to the one satisfied by the linearization. In terms of the rescaled representation, which is a nonlinear Fokker–Planck equation, the convergence rate turns out to be polynomial in time. This result is in contrast with the known exponential decay of such representation for all other values of m.     

Stability of Degenerate Stationary Waves for Viscous Gases

This paper is concerned with the asymptotic stability of degenerate stationary waves for viscous gases in the half space. We discuss the following two cases: (1) viscous conservation laws and (2) damped wave equations with nonlinear convection. In each case, we prove that the solution converges to the corresponding degenerate stationary wave at the rate t ^−α/4 as t → ∞, provided that the initial perturbation is in the weighted space L²_a=L²(\mathbb R₊; (1+x)^a dx){L^2_\alpha=L^2({\mathbb R}_+;\,(1+x)^\alpha dx)} . This convergence rate t ^−α/4 is weaker than the one for the non-degenerate case and requires the restriction α < α_*(q), where α_*(q) is the critical value depending only on the degeneracy exponent q. Such a restriction is reasonable because the corresponding linearized operator for viscous conservation laws cannot be dissipative in L²_a{L^2_\alpha} for α > α^*(q) with another critical value α^*(q). Our stability analysis is based on the space–time weighted energy method in which the spatial weight is chosen as a function of the degenerate stationary wave.     

Global Behavior of Compressible Navier-Stokes Equations with a Degenerate Viscosity Coefficient

In this paper, we study a free boundary problem for compressible Navier-Stokes equations with density-dependent viscosity. Precisely, the viscosity coefficient μ is proportional to ρ ^θ with , where ρ is the density, and γ > 1 is the physical constant of polytropic gas. Under certain assumptions imposed on the initial data, we obtain the global existence and uniqueness of the weak solution, give the uniform bounds (with respect to time) of the solution and show that it converges to a stationary one as time tends to infinity. Moreover, we estimate the stabilization rate in L ^∞ norm, (weighted) L ² norm and weighted H ¹ norm of the solution as time tends to infinity.     

Well-Posedness in Smooth Function Spaces for the Moving-Boundary Three-Dimensional Compressible Euler Equations in Physical Vacuum

We prove well-posedness for the three-dimensional compressible Euler equations with moving physical vacuum boundary, with an equation of state given by p(ρ) =  C _γ ρ ^γ for γ > 1. The physical vacuum singularity requires the sound speed c to go to zero as the square-root of the distance to the moving boundary, and thus creates a degenerate and characteristic hyperbolic free-boundary system wherein the density vanishes on the free-boundary, the uniform Kreiss–Lopatinskii condition is violated, and manifest derivative loss ensues. Nevertheless, we are able to establish the existence of unique solutions to this system on a short time-interval, which are smooth (in Sobolev spaces) all the way to the moving boundary, and our estimates have no derivative loss with respect to initial data. Our proof is founded on an approximation of the Euler equations by a degenerate parabolic regularization obtained from a specific choice of a degenerate artificial viscosity term, chosen to preserve as much of the geometric structure of the Euler equations as possible. We first construct solutions to this degenerate parabolic regularization using a higher-order version of Hardy’s inequality; we then establish estimates for solutions to this degenerate parabolic system which are independent of the artificial viscosity parameter. Solutions to the compressible Euler equations are found in the limit as the artificial viscosity tends to zero. Our regular solutions can be viewed as degenerate viscosity solutions. Our methodology can be applied to many other systems of degenerate and characteristic hyperbolic systems of conservation laws.     

Isentropic Approximation¶of the Compressible Euler System¶in One Space Dimension

Using the stability results of Bressan & Colombo [BC] for strictly hyperbolic 2 × 2 systems in one space dimension, we prove that the solutions of isentropic and non-isentropic Euler equations in one space dimension with the respective initial data (ρ₀, u ₀) and (ρ₀, u ₀, &\theta;₀=ρ₀ ^{γ− 1}) remain close as soon as the total variation of (ρ₀, u ₀) is sufficiently small. Accepted April 25, 2000?Published online November 24, 2000     

Billiards and Two-Dimensional Problems of Optimal Resistance

A body moves in a medium composed of noninteracting point particles; the interaction of the particles with the body is completely elastic. The problem is: find the body’s shape that minimizes or maximizes resistance of the medium to its motion. This is the general setting of the optimal resistance problem going back to Newton. Here, we restrict ourselves to the two-dimensional problems for rotating (generally non-convex) bodies. The main results of the paper are the following. First, to any compact connected set with piecewise smooth boundary B ì \mathbbR²{B \subset \mathbb{R}^2} we assign a measure ν _B on ∂(conv B)×[ − π/2, π/2] generated by the billiard in \mathbbR² \B{\mathbb{R}^2 \setminus B} and characterize the set of measures {ν _B }. Second, using this characterization, we solve various problems of minimal and maximal resistance of rotating bodies by reducing them to special Monge–Kantorovich problems.     

Bifurcations of Heteroclinic Orbits

The search for traveling wave solutions of a semilinear diffusion partial differential equation can be reduced to the search for heteroclinic solutions of the ordinary differential equation ü − cu̇ + f(u) = 0, where c is a positive constant and f is a nonlinear function. A heteroclinic orbit is a solution u(t) such that u(t) → γ ₁ as t → −∞ and u(t) → γ ₂ as t → ∞ where γ ₁, γ ₂ are zeros of f. We study the existence of heteroclinic orbits under various assumptions on the nonlinear function f and their bifurcations as c is varied. Our arguments are geometric in nature and so we make only minimal smoothness assumptions. We only assume that f is continuous and that the equation has a unique solution to the initial value problem. Under these weaker smoothness conditions we reprove the classical result that for large c there is a unique positive heteroclinic orbit from 0 to 1 when f(0) = f(1) = 0 and f(u) > 0 for 0 < u < 1. When there are more zeros of f, there is the possibility of bifurcations of the heteroclinic orbit as c varies. We give a detailed analysis of the bifurcation of the heteroclinic orbits when f is zero at the five points −1 < −θ < 0 < θ < 1 and f is odd. The heteroclinic orbit that tends to 1 as t → ∞ starts at one of the three zeros, −θ, 0, θ as t → −∞. It hops back and forth among these three zeros an infinite number of times in a predictable sequence as c is varied.     

Investigation of liquid–liquid drop coalescence using tomographic PIV

High-speed tomographic PIV was used to investigate the coalescence of drops placed on a liquid/liquid interface; the coalescence of a single drop and of a drop in the presence of an adjacent drop (side-by-side drops) was investigated. The viscosity ratio between the drop and surrounding fluids was 0.14, the Ohnesorge number (Oh = μ_d/(ρ_dσD)^1/2) was 0.011, and Bond numbers (Bo = (ρ _d  − ρ _s )gD ²/σ) were 3.1–7.5. Evolving volumetric velocity fields of the full coalescence process allowed for quantification of the velocity scales occurring over different time scales. For both single and side-by-side drops, the coalescence initiates with an off-axis film rupture and film retraction speeds an order of magnitude larger than the collapse speed of the drop fluid. This is followed by the formation and propagation of an outward surface wave along the coalescing interface with wavelength of approximately 2D. For side-by-side drops, the collapse of the first drop is asymmetric due to the presence of the second drop and associated interface deformation. Overall, tomographic PIV provides insight into the flow physics and inherent three-dimensionalities in the coalescence process that would not be achievable with flow visualization or planar PIV only.     

Radial Solutions and Phase Separation in a System of Two Coupled Schrödinger Equations

We consider the nonlinear elliptic system

Generalized Solution for 1-D Non-Newtonian Flow in a Porous Domain due to an Instantaneous Mass Injection

Non-Newtonian fluid flow through porous media is of considerable interest in several fields, ranging from environmental sciences to chemical and petroleum engineering. In this article, we consider an infinite porous domain of uniform permeability k and porosity f{\phi} , saturated by a weakly compressible non-Newtonian fluid, and analyze the dynamics of the pressure variation generated within the domain by an instantaneous mass injection in its origin. The pressure is taken initially to be constant in the porous domain. The fluid is described by a rheological power-law model of given consistency index H and flow behavior index n; n, < 1 describes shear-thinning behavior, n > 1 shear-thickening behavior; for n = 1, the Newtonian case is recovered. The law of motion for the fluid is a modified Darcy’s law based on the effective viscosity μ _ef , in turn a function of f, H, n{\phi, H, n} . Coupling the flow law with the mass balance equation yields the nonlinear partial differential equation governing the pressure field; an analytical solution is then derived as a function of a self-similar variable η = rt ^β (the exponent β being a suitable function of n), combining spatial coordinate r and time t. We revisit and expand the work in previous papers by providing a dimensionless general formulation and solution to the problem depending on a geometrical parameter d, valid for plane (d = 1), cylindrical (d = 2), and semi-spherical (d = 3) geometry. When a shear-thinning fluid is considered, the analytical solution exhibits traveling wave characteristics, in variance with Newtonian fluids; the front velocity is proportional to t ^(n-2)/2 in plane geometry, t ^{(2n-3)/(3−n)} in cylindrical geometry, and t ^{(3n-4)/[2(2−n)]} in semi-spherical geometry. To reflect the uncertainty inherent in the value of the problem parameters, we consider selected properties of fluid and matrix as independent random variables with an associated probability distribution. The influence of the uncertain parameters on the front position and the pressure field is investigated via a global sensitivity analysis evaluating the associated Sobol’ indices. The analysis reveals that compressibility coefficient and flow behavior index are the most influential variables affecting the front position; when the excess pressure is considered, compressibility and permeability coefficients contribute most to the total response variance. For both output variables the influence of the uncertainty in the porosity is decidedly lower.     

Positivity and Almost Positivity of Biharmonic Green’s Functions under Dirichlet Boundary Conditions

In general, for higher order elliptic equations and boundary value problems like the biharmonic equation and the linear clamped plate boundary value problem, neither a maximum principle nor a comparison principle or—equivalently—a positivity preserving property is available. The problem is rather involved since the clamped boundary conditions prevent the boundary value problem from being reasonably written as a system of second order boundary value problems. It is shown that, on the other hand, for bounded smooth domains W ì \mathbbRⁿ{\Omega \subset\mathbb{R}^n} , the negative part of the corresponding Green’s function is “small” when compared with its singular positive part, provided n\geqq 3{n\geqq 3} . Moreover, the biharmonic Green’s function in balls B ì \mathbbRⁿ{B\subset\mathbb{R}^n} under Dirichlet (that is, clamped) boundary conditions is known explicitly and is positive. It has been known for some time that positivity is preserved under small regular perturbations of the domain, if n = 2. In the present paper, such a stability result is proved for n\geqq 3{n\geqq 3} .     

The Dirichlet Problem for H-Systems with Small Boundary Data: BlowUp Phenomena and Nonexistence Results

Given H:ℝ³→ℝ of class C¹ and bounded, we consider a sequence (uⁿ) of solutions of the H-system in the unit open disc satisfying the boundary condition uⁿ=γⁿ on ∂. In the first part of this paper, assuming that (uⁿ) is bounded in H¹(,ℝ³) we study the behavior of (uⁿ) when the boundary data γⁿ shrink to zero. We show that either uⁿ→0 strongly in H¹(,ℝ³) or uⁿ blows up at least one H-bubble ω, namely a nonconstant, conformal solution of the H-system on ℝ². Under additional assumptions on H, we can obtain more precise information on the blow up. In the second part of this paper we investigate the multiplicity of solutions for the Dirichlet problem on the disc with small boundary datum. We detect a family of nonconstant functions H (even close to a nonzero constant in any reasonable topology) for which the Dirichlet problem cannot admit a ``large' solution at a mountain pass level when the boundary datum is small.     

Laminarisation and re-transition of a turbulent boundary layer subjected to favourable pressure gradient

 Experimental results are reported for the response of an initially turbulent boundary layer (Re _θ≈1700) to a favourable pressure gradient with a peak value of K≡(−υ/ρU ³ _E ) dp/dx equal to 4.4×10^-6. In the near-wall region of the boundary layer (y/δ<0.1) the turbulence intensity u′ scales roughly with the free-stream velocity U _E until close to the location where K is a maximum whereas in the outer region u′ remains essentially frozen. Once the pressure gradient is relaxed, the turbulence level increases throughout the boundary layer until K falls to zero when the near wall u′ levels show a significant decrease. The intermittency γ is the clearest indicator of a fundamental change in the turbulence structure: once K exceeds 3×10^-6, the value of γ in the immediate vicinity of the wall γ _s falls rapidly from unity, reaches zero at the location where K again falls below 3×10^-6 and then rises back to unity. Although γ is practically zero throughout the boundary layer in the vicinity of γ _s =0, the turbulence level remains high. The explanation for what appears to be a contradiction is that the turbulent frequencies are too low to induce turbulent mixing. The mean velocity profile changes shape abruptly where K exceeds 3×10^-6. Values for the skin friction coefficient, based upon hot-film measurements, peak at the same location as K and fall to a minimum close to the location where K drops back to zero. Received: 28 January 1998/Accepted: 8 April 1998     

Spectral Stability of Small-Amplitude Viscous Shock Waves in Several Space Dimensions

A planar viscous shock profile of a hyperbolic–parabolic system of conservation laws is a steady solution in a moving coordinate frame. The asymptotic stability of viscous profiles and the related vanishing-viscosity limit are delicate questions already in the well understood case of one space dimension and even more so in the case of several space dimensions. It is a natural idea to study the stability of viscous profiles by analyzing the spectrum of the linearization about the profile. The Evans function method provides a geometric dynamical-systems framework to study the eigenvalue problem. In this approach eigenvalues correspond to zeros of an essentially analytic function E(rl,rw){\mathcal{E}(\rho\lambda,\rho\omega)} which detects nontrivial intersections of the so-called stable and unstable spaces, that is, spaces of solutions that decay on one (“−∞”) or the other side (“ + ∞”) of the shock wave, respectively. In a series of pioneering papers, Kevin Zumbrun and collaborators have established in various contexts that spectral stability, that is, the non-vanishing of E(rl,rw){\mathcal{E}(\rho\lambda,\rho\omega)} and the non-vanishing of the Lopatinski–Kreiss–Majda function Δ(λ,ω), imply nonlinear stability of viscous shock profiles in several space dimensions. In this paper we show that these conditions hold true for small amplitude extreme shocks under natural assumptions. This is done by exploiting the slow-fast nature of the small-amplitude limit, which was used in a previous paper by the authors to prove spectral stability of small-amplitude shock waves in one space dimension. Geometric singular perturbation methods are applied to decompose the stable and unstable spaces into subbundles with good control over their limiting behavior. Three qualitatively different regimes are distinguished that relate the small strength e{\epsilon} of the shock wave to appropriate ranges of values of the spectral parameters (ρλ, ρ ω). Various rescalings are used to overcome apparent degeneracies in the problem caused by loss of hyperbolicity or lack of transversality.     

Influence of surfactant upon air entrainment hysteresis in curtain coating

The onset of air entrainment for curtain coating onto a pre-wetted substrate was studied experimentally in similar parameter regimes to commercial coating (Re = ρQ/μ = O(1), We = ρQ u _c /σ = O(10), Ca = μU/σ = O(1)). Impingement speed and viscosity were previously shown to be critical parameters in correlating air entrainment data with three qualitatively different regimes of hydrodynamic assist identified (Marston et al. in Exp Fluids 42(3):483–488, 2007a). The interaction of the impinging curtain with the pre-existing film also led to a significant hysteretic effect throughout the flow rate-substrate speed parameter space. For the first time, results considering the influence of surfactants are presented in attempt to elucidate the relative importance of surface tension in this inertia-dominated system. The results show quantitative and qualitative differences to previous results with much more complex hysteretic behaviour which has only been reported previously for rough surfaces.     

Well-posedness and Stability of a Multi-dimensional Tumor Growth Model

We study a moving boundary problem modeling the growth of in vitro tumors. This problem consists of two elliptic equations describing the distribution of the nutrient and the internal pressure, respectively, and a first-order partial differential equation describing the evolution of the moving boundary. An important feature is that the effect of surface tension on the moving boundary is taken into account. We show that this problem is locally well-posed for a large class of initial data by using analytic semi-group theory. We also prove that if the surface tension coefficient γ is larger than a threshold value γ _* then the unique flat equilibrium is asymptotically stable, whereas in the case γ  < γ _* this flat equilibrium is unstable.     

Prediction of CHF in concentric-tube open thermosiphon using artificial neural network and genetic algorithm

In this paper, an artificial neural network (ANN) for predicting critical heat flux (CHF) of concentric-tube open thermosiphon has been trained successfully based on the experimental data from the literature. The dimensionless input parameters of the ANN are density ratio, ρ _l/ρ _v; the ratio of the heated tube length to the inner diameter of the outer tube, L/D _i; the ratio of frictional area, d _i/(D _i + d _o); and the ratio of equivalent heated diameter to characteristic bubble size, D _he/[σ/g(ρ _l−ρ _v)]^0.5, the output is Kutateladze number, Ku. The predicted values of ANN are found to be in reasonable agreement with the actual values from the experiments with a mean relative error (MRE) of 8.46%. New correlations for predicting CHF were also proposed by using genetic algorithm (GA) and succeeded to correlate the existing CHF data with better accuracy than the existing empirical correlations.     

Optimal Decay Rate of the Compressible Navier–Stokes–Poisson System in {\mathbb {R}^3}

The compressible Navier–Stokes–Poisson (NSP) system is considered in ${\mathbb {R}^3}The compressible Navier–Stokes–Poisson (NSP) system is considered in \mathbb R³{\mathbb {R}^3} in the present paper, and the influences of the electric field of the internal electrostatic potential force governed by the self-consistent Poisson equation on the qualitative behaviors of solutions is analyzed. It is observed that the rotating effect of electric field affects the dispersion of fluids and reduces the time decay rate of solutions. Indeed, we show that the density of the NSP system converges to its equilibrium state at the same L ²-rate (1+t)^{-\frac 34}{(1+t)^{-\frac {3}{4}}} or L ^∞-rate (1 + t)^−3/2 respectively as the compressible Navier–Stokes system, but the momentum of the NSP system decays at the L ²-rate (1+t)^{-\frac 14}{(1+t)^{-\frac {1}{4}}} or L ^∞-rate (1 + t)⁻¹ respectively, which is slower than the L ²-rate (1+t)^{-\frac 34}{(1+t)^{-\frac {3}{4}}} or L ^∞-rate (1 + t)^−3/2 for compressible Navier–Stokes system [Duan et al., in Math Models Methods Appl Sci 17:737–758, 2007; Liu and Wang, in Comm Math Phys 196:145–173, 1998; Matsumura and Nishida, in J Math Kyoto Univ 20:67–104, 1980] and the L ^∞-rate (1 + t)^−p with p ? (1, 3/2){p \in (1, 3/2)} for irrotational Euler–Poisson system [Guo, in Comm Math Phys 195:249–265, 1998]. These convergence rates are shown to be optimal for the compressible NSP system.