Unsteady MHD double-diffusive convection boundary-layer flow past a radiate hot vertical surface in porous media in the presence of chemical reaction and heat sink

This paper presents an analytical study of the unsteady MHD free convective heat and mass transfer flow of a viscous, incompressible, gray, absorbing-emitting but non-scattering, optically-thick and electrically conducting fluid occupying a semi-infinite porous regime adjacent to an infinite moving hot vertical plate with constant velocity. We employ a Darcian viscous flow model for the porous medium the Rosseland diffusion approximation is used to describe the radiative heat flux in the energy equation. The homogeneous chemical reaction of first order is accounted in mass diffusion equation. The governing equations are solved in closed form by Laplace-transform technique. A parametric study of all involved parameters is conducted and representative set of numerical results for the velocity, temperature, concentration, shear stress function $\frac{\partial u}{\partial y} \vert_{y=0}$ , temperature gradient $\frac{\partial \theta }{ \partial y}\vert_{y=0}$ , and concentration gradient $\frac{ \partial \phi }{\partial y}\vert_{y=0}$ is illustrated graphically and physical aspects of the problem are discussed.     

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

An Updated Portrait of Transition to Turbulence in Laminar Pipe Flows with Periodic Time Dependence (A Correlation Study)

Transition to turbulence in axially symmetrical laminar pipe flows with periodic time dependence classified as pure oscillating and pulsatile (pulsating) ones is the concern of the paper. The current state of art on the transitional characteristics of pulsatile and oscillating pipe flows is introduced with a particular attention to the utilized terminology and methodology. Transition from laminar to turbulent regime is usually described by the presence of the disturbed flow with small amplitude perturbations followed by the growth of turbulent bursts. The visual treatment of velocity waveforms is therefore a preferred inspection method. The observation of turbulent bursts first in the decelerating phase and covering the whole cycle of oscillation are used to define the critical states of the start and end of transition, respectively. A correlation study referring to the available experimental data of the literature particularly at the start of transition are presented in terms of the governing periodic flow parameters. In this respect critical oscillating and time averaged Reynolds numbers at the start of transition; Re _os,crit and Re _ta,crit are expressed as a major function of Womersley number, $\sqrt {\omega ^\prime } $ defined as dimensionless frequency of oscillation, f. The correlation study indicates that in oscillating flows, an increase in Re _os,crit with increasing magnitudes of $\sqrt {\omega ^\prime } $ is observed in the covered range of $1<\sqrt {\omega ^\prime } <72$ . The proposed equation (Eq. 7), ${\rm{Re}}_{os,crit} ={\rm{Re}}_{os,crit} \left( {\sqrt {\omega ^\prime } } \right)$ , can be utilized to estimate the critical magnitude of $\sqrt {\omega ^\prime }$ at the start of transition with an accuracy of ±12?% in the range of $\sqrt {\omega ^\prime } <41$ . However in pulsatile flows, the influence of $\sqrt {\omega ^\prime }$ on Re _ta,crit seems to be different in the ranges of $\sqrt {\omega ^\prime } <8$ and $\sqrt {\omega ^\prime } >8$ . Furthermore there is rather insufficient experimental data in pulsatile flows considering interactive influences of $\sqrt {\omega ^\prime } $ and velocity amplitude ratio, A ₁. For the purpose, the measurements conducted at the start of transition of a laminar sinusoidal pulsatile pipe flow test case covering the range of 0.21<?A ₁?<0.95 with $\sqrt {\omega ^\prime } <8$ are evaluated. In conformity with the literature, the start of transition corresponds to the observation of first turbulent bursts in the decelerating phase of oscillation. The measured data indicate that increase in $\sqrt {\omega ^\prime } $ is associated with an increase in Re _ta,crit up to $\sqrt {\omega ^\prime } =3.85$ while a decrease in Re _ta,crit is observed with an increase in $\sqrt {\omega ^\prime } $ for $\sqrt {{\omega }'} >3.85$ . Eventually updated portrait is pointing out the need for further measurements on i) the end of transition both in oscillating and pulsatile flows with the ranges of $\sqrt {\omega ^\prime } <8$ and $\sqrt {\omega ^\prime } >8$ , and ii) the interactive influences of $\sqrt {\omega ^\prime } $ and A ₁ on Re _ta,crit in pulsatile flows with the range of $\sqrt {\omega ^\prime } >8$ .     

Generalized Flows Satisfying Spatial Boundary Conditions

In a region D in ${\mathbb{R}^2}$ or ${\mathbb{R}^3}$ , the classical Euler equation for the regular motion of an inviscid and incompressible fluid of constant density is given by $$\partial_t v+(v\cdot \nabla_x)v=-\nabla_x p, {\rm div}_x v=0,$$ where v(t, x) is the velocity of the particle located at ${x\in D}$ at time t and ${p(t,x)\in\mathbb{R}}$ is the pressure. Solutions v and p to the Euler equation can be obtained by solving $$\left\{\begin{array}{l} \nabla_x\left\{\partial_t\phi(t,x,a) + p(t,x)+(1/2)|\nabla_x\phi(t,x,a)|^2 \right\}=0\,{\rm at}\,a=\kappa(t,x),\\ v(t,x)=\nabla_x \phi(t,x,a)\,{\rm at}\,a=\kappa(t,x), \\ \partial_t\kappa(t,x)+(v\cdot\nabla_x)\kappa(t,x)=0, \\ {\rm div}_x v(t,x)=0, \end{array}\right. \quad\quad\quad\quad\quad(0.1)$$ where $$\phi:\mathbb{R}\times D\times \mathbb{R}^l\rightarrow\mathbb{R}\,{\rm and}\, \kappa:\mathbb{R}\times D \rightarrow \mathbb{R}^l$$ are additional unknown mappings (l?≥ 1 is prescribed). The third equation in the system says that ${\kappa\in\mathbb{R}^l}$ is convected by the flow and the second one that ${\phi}$ can be interpreted as some kind of velocity potential. However vorticity is not precluded thanks to the dependence on a. With the additional condition κ(0, x)?=?x on D (and thus l?=?2 or 3), this formulation was developed by Brenier (Commun Pure Appl Math 52:411–452, 1999) in his Eulerian–Lagrangian variational approach to the Euler equation. He considered generalized flows that do not cross ${\partial D}$ and that carry each “particle” at time t?=?0 at a prescribed location at time t?=?T?>?0, that is, κ(T, x) is prescribed in D for all ${x\in D}$ . We are concerned with flows that are periodic in time and with prescribed flux through each point of the boundary ${\partial D}$ of the bounded region D (a two- or three-dimensional straight pipe). More precisely, the boundary condition is on the flux through ${\partial D}$ of particles labelled by each value of κ at each point of ${\partial D}$ . One of the main novelties is the introduction of a prescribed “generalized” Bernoulli’s function ${H:\mathbb{R}^l\rightarrow \mathbb{R}}$ , namely, we add to (0.1) the requirement that $$\partial_t\phi(t,x,a) +p(t,x)+(1/2)|\nabla_x\phi(t,x,a)|^2=H(a)\,{\rm at}\,a=\kappa(t,x)\quad\quad\quad\quad\quad(0.2)$$ with ${\phi,p,\kappa}$ periodic in time of prescribed period T?>?0. Equations (0.1) and (0.2) have a geometrical interpretation that is related to the notions of “Lamb’s surfaces” and “isotropic manifolds” in symplectic geometry. They may lead to flows with vorticity. An important advantage of Brenier’s formulation and its present adaptation consists in the fact that, under natural hypotheses, a solution in some weak sense always exists (if the boundary conditions are not contradictory). It is found by considering the functional $$(\kappa,v)\rightarrow \int\limits_{0}^T \int\limits_D\left\{\frac 1 2 |v(t,x)|^2+H(\kappa(t,x))\right\}dt\, dx$$ defined for κ and v that are T-periodic in t, such that $$\partial_t\kappa(t,x)+(v\cdot\nabla_x)\kappa(t,x)=0, {\rm div}_x v(t,x)=0,$$ and such that they satisfy the boundary conditions. The domain of this functional is enlarged to some set of vector measures and then a minimizer can be obtained. For stationary planar flows, the approach is compared with the following standard minimization method: to minimize $$\int\limits_{]0,L[\times]0,1[} \{(1/2)|\nabla \psi|^2+H(\psi)\}dx\,{\rm for}\,\psi\in W^{1,2}(]0,L[\times]0,1[)$$ under appropriate boundary conditions, where ψ is the stream function. For a minimizer, corresponding functions ${\phi}$ and κ are given in terms of the stream function ψ.     

Mixed Convection Boundary-Layer Flow Along a Vertical Cylinder Embedded in a Porous Medium Filled by a Nanofluid

The steady mixed convection boundary-layer flow on a vertical circular cylinder embedded in a porous medium filled by a nanofluid is studied for both cases of a heated and a cooled cylinder. The governing system of partial differential equations is reduced to ordinary differential equations by assuming that the surface temperature of the cylinder and the velocity of the external (inviscid) flow vary linearly with the axial distance x measured from the leading edge. Solutions of the resulting ordinary differential equations for the flow and heat transfer characteristics are evaluated numerically for various values of the governing parameters, namely the nanoparticle volume fraction ${\phi}$ , the mixed convection or buoyancy parameter ?? and the curvature parameter ??. Results are presented for the specific case of copper nanoparticles. A critical value ?? _c of ?? with ?? _c <?0 is found, with the values of | ?? _c| increasing as the curvature parameter ?? or nanoparticle volume fraction ${\phi}$ is increased. Dual solutions are seen for all values of ?? >??? _c for both aiding, ?? >?0 and opposing, ?? <?0, flows. Asymptotic solutions are also determined for both the free convection limit ${(\lambda \gg 1)}$ and for large curvature parameter ${(\gamma \gg 1)}$ .     

A well posed problem in singular Fickian diffusion

We prove that the problem of solving $$u_t = (u^{m - 1} u_x )_x {\text{ for }} - 1< m \leqq 0$$ with initial conditionu(x, 0)=φ(x) and flux conditions at infinity \(\mathop {\lim }\limits_{x \to \infty } u^{m - 1} u_x = - f(t),\mathop {\lim }\limits_{x \to - \infty } u^{m - 1} u_x = g(t)\) , admits a unique solution \(u \in C^\infty \{ - \infty< x< \infty ,0< t< T\} \) for every φεL¹(R), φ≧0, φ≡0 and every pair of nonnegative flux functionsf, g ε L _loc ^∞ [0, ∞) The maximal existence time is given by $$T = \sup \left\{ {t:\smallint \phi (x)dx > \int\limits_0^t {[f} (s) + g(s)]ds} \right\}$$ This mixed problem is ill posed for anym outside the above specified range.     

Advective Transport in Heterogeneous Formations: The Impact of Spatial Anisotropy on the Breakthrough Curve

Water flow and solute transport take place in formations of spatially variable conductivity K. The logconductivity Y?= ln K is modeled as a random stationary space function, of normal univariate pdf (of mean In K _G and variance ${\sigma_{Y}^{2}}$ ) and of axisymmetric autocorrelation of integral scales I _h,I _v (anisotropy ratio f?=?I _v/I _h?<?1). The head gradient and the velocity are uniform in the mean, parallel to bedding, and of constant and given as J and U, respectively. Transport is ruled by advection, which typically overwhelms pore scale dispersion in the breakthrough curve (BTC) determination. In the present study we analyze the impact of anisotropy f on the BTC of a passive solute, which is related to the mass flux??? (t, x) at a control plane at x. While a considerable body of literature dealt with weakly heterogeneous formations ( ${\sigma _{Y}^{2} <1 }$ ), the present study addresses the case of ${\sigma _{Y}^{2} >1 }$ , which is of interest for many aquifers and is more difficult to solve either numerically or by approximations. We approach the three dimensional problem by modeling the structure as an ensemble of densely packed oblate spheroids of semi-major and semi-minor axis R and f R, respectively, and independent lognormal K, submerged in a matrix of uniform conductivity K _ef, the effective conductivity of the ensemble. The detailed numerical simulations of transport show that the BTC is insensitive to the value of the anisotropy ratio f, i.e.,??? (t, x) I _h/U depends only on ${\sigma _{Y}^{2}}$ (except for small differences in the tail). This important result implies that transport, as quantified by BTCs or spatial longitudinal mass distributions, can be modeled accurately by the much simpler solutions developed in the past for isotropic media, like e.g., the semi-analytical self-consistent approximation.     

The approach of solutions of nonlinear diffusion equations to travelling front solutions

The paper is concerned with the asymptotic behavior as t → ∞ of solutions u(x, t) of the equation u_t—u_xx—∞;(u)=O, x∈(—∞, ∞) , in the case ∞(0)=∞(1)=0, ∞′(0)<0, ∞′(1)<0. Commonly, a travelling front solution u=U(x-ct), U(-∞)=0, U(∞)=1, exists. The following types of global stability results for fronts and various combinations of them will be given.

1. Let u(x, 0)=u ₀(x) satisfy 0≦u ₀≦1. Let \(a\_ = \mathop {\lim \sup u0}\limits_{x \to - \infty } {\text{(}}x{\text{), }}\mathop {\lim \inf u0}\limits_{x \to \infty } {\text{(}}x{\text{)}}\) . Then u approaches a translate of U uniformly in x and exponentially in time, if a_? is not too far from 0, and a₊ not too far from 1.
2. Suppose \(\int\limits_{\text{0}}^{\text{1}} {f{\text{(}}u{\text{)}}du} > {\text{0}}\) . If a _? and a ₊ are not too far from 0, but u₀ exceeds a certain threshold level for a sufficiently large x-interval, then u approaches a pair of diverging travelling fronts.
3. Under certain circumstances, u approaches a “stacked” combination of wave fronts, with differing ranges.

     

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.     

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.     

Eulerian–Lagrangian bridge for the energy and dissipation spectra in isotropic turbulence

We study, numerically and analytically, the relationship between the Eulerian spectrum of kinetic energy, E _E(k, t), in isotropic turbulence and the corresponding Lagrangian frequency energy spectrum, E _L(ω, t), for which we derive an evolution equation. Our DNS results show that not only E _L(ω, t) but also the Lagrangian frequency spectrum of the dissipation rate ${\varepsilon_{\rm L} (\omega, t)}$ has its maximum at low frequencies (about the turnover frequency of energy-containing eddies) and decays exponentially at large frequencies ω (about a half of the Kolmogorov microscale frequency) for both stationary and decaying isotropic turbulence. Our main analytical result is the derivation of equations that bridge the Eulerian and Lagrangian spectra and allow the determination of the Lagrangian spectrum, E _L (ω) for a given Eulerian spectrum, E _E (k), as well as the Lagrangian dissipation, ${\varepsilon_{\rm L}(\omega)}$ , for a given Eulerian counterpart, ${\varepsilon_{\rm E} (k)=2\nu k^2 E_{\rm E}(k)}$ . These equations were derived from the Navier–Stokes equations in the sweeping-free coordinate system (intermediate between the Eulerian and Lagrangian frameworks) which eliminates the effect of the kinematic sweeping of the small eddies by the larger eddies. We show that both analytical relationships between E _L (ω) and E _E (k) and between ${\varepsilon_{\rm L} (\omega)}$ and ${\varepsilon_{\rm E} (k)}$ are in very good quantitative agreement with our DNS results and explain how ${\varepsilon_{\rm L} (\omega, t)}$ has its maximum at low frequencies and decays exponentially at large frequencies.     

A General Class of Free Boundary Problems for Fully Nonlinear Elliptic Equations

In this paper we study the fully nonlinear free boundary problem $$\left\{\begin{array}{ll}F(D^{2}u) = 1 & {\rm almost \, everywhere \, in}\, B_{1} \cap \Omega\\ |D^{2} u| \leqq K & {\rm almost \, everywhere \, in} \, B_{1} \setminus \Omega,\end{array}\right.$$ where K > 0, and Ω is an unknown open set. Our main result is the optimal regularity for solutions to this problem: namely, we prove that W ^2,n solutions are locally C ^1,1 inside B ₁. Under the extra condition that ${\Omega \supset \{D{u} \neq 0 \}}$ and a uniform thickness assumption on the coincidence set {D u = 0}, we also show local regularity for the free boundary ${\partial \Omega \cap B_1}$ .     

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

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.     

Visualization and acoustic sounding of the fine structure of a stratified flow behind a vertical plate

High-resolution shadow visualization and high-frequency sonar detection are applied to separate out the density wake and a fine streaky structure in the vicinity of a vertical plate in motion in salt-stratified water. The length of the sounding acoustic wave is taken to be approximately equal to the universal microscale $\delta _N^v = \sqrt {{v \mathord{\left/ {\vphantom {v N}} \right. \kern-0em} N}}$ , where ν and N are the kinematic viscosity and the buoyance frequency. In the spectra of the vertical oscillations of the acoustic contrast some characteristic frequencies ω are separated out and used to calculate the local Stokes microscales $\delta _\omega ^v = \sqrt {{v \mathord{\left/ {\vphantom {v \omega }} \right. \kern-0em} \omega }}$ in the density wake region. The scales determined from the data of independent optical and acoustic measurements are in agreement with each other.     

Integrability of the constrained rigid body

The integrability theory for the differential equations, which describe the motion of an unconstrained rigid body around a fixed point is well known. When there are constraints the theory of integrability is incomplete. The main objective of this paper is to analyze the integrability of the equations of motion of a constrained rigid body around a fixed point in a force field with potential U(γ)=U(γ ₁,γ ₂,γ ₃). This motion subject to the constraint 〈ν,ω〉=0 with ν is a constant vector is known as the Suslov problem, and when ν=γ is the known Veselova problem, here ω=(ω ₁,ω ₂,ω ₃) is the angular velocity and 〈?,?〉 is the inner product of $\mathbb{R}^{3}$ . We provide the following new integrable cases. (i) The Suslov’s problem is integrable under the assumption that ν is an eigenvector of the inertial tensor I and the potential is such that $$U=-\frac{1}{2I_1I_2}\bigl(I_1\mu^2_1+I_2 \mu^2_2\bigr), $$ where I ₁,I ₂, and I ₃ are the principal moments of inertia of the body, μ ₁ and μ ₂ are solutions of the first-order partial differential equation $$\gamma_3 \biggl(\frac{\partial\mu_1}{\partial\gamma_2}- \frac{\partial\mu_2}{\partial \gamma_1} \biggr)- \gamma_2\frac{\partial \mu_1}{\partial\gamma_3}+\gamma_1\frac{\partial\mu_2}{\partial \gamma_3}=0. $$ (ii) The Veselova problem is integrable for the potential $$U=-\frac{\varPsi^2_1+\varPsi^2_2}{2(I_1\gamma^2_2+I_2\gamma^2_1)}, $$ where Ψ ₁ and Ψ ₂ are the solutions of the first-order partial differential equation where $p=\sqrt{I_{1}I_{2}I_{3} (\frac{\gamma^{2}_{1}}{I_{1}}+\frac{\gamma^{2}_{2}}{I_{2}}+ \frac{\gamma^{2}_{3}}{I_{3}} )}$ . Also it is integrable when the potential U is a solution of the second-order partial differential equation where $\tau_{2}=I_{1}\gamma^{2}_{1}+I_{2}\gamma^{2}_{2}+I_{3}\gamma^{2}_{3}$ and $\tau_{3}=\frac{\gamma^{2}_{1}}{I_{1}}+\frac{\gamma^{2}_{2}}{I_{2}}+ \frac{\gamma^{2}_{3}}{I_{3}}$ . Moreover, we show that these integrable cases contain as a particular case the previous known results.     

A Widder’s Type Theorem for the Heat Equation with Nonlocal Diffusion

The main goal of this work is to prove that every non-negative strong solution u(x, t) to the problem $$u_t + (-\Delta)^{\alpha/2}{u} = 0 \,\, {\rm for} (x, t) \in {\mathbb{R}^n} \times (0, T ), \, 0 < \alpha < 2,$$ can be written as $$u(x, t) = \int_{\mathbb{R}^n} P_t (x - y)u(y, 0) dy,$$ where $$P_t (x) = \frac{1}{t^{n/ \alpha}}P \left(\frac{x}{t^{1/ \alpha}}\right),$$ and $$P(x) := \int_{\mathbb{R}^n} e^{i x\cdot\xi-|\xi |^\alpha} d\xi.$$ This result shows uniqueness in the setting of non-negative solutions and extends some classical results for the heat equation by Widder in [15] to the nonlocal diffusion framework.     

An experimental study on turbulence modification in the near-wall boundary layer of a dilute gas-particle channel flow

Turbulence modifications of a dilute gas-particle flow are experimentally investigated in the lower boundary layer of a horizontal channel by means of a simultaneous two-phase PIV measurement technique. The measurements are conducted in the near-wall region with y ⁺?<?250 at Re _τ (based on the wall friction velocity u _τ and half channel height h)?=?430. High spatial resolution and small interrogation window are used to minimize the PIV measurement uncertainty due to the velocity gradient near the wall. Polythene beads with the diameter of 60?μm (d _p ⁺ ?=?1.71, normalized by the fluid kinematic viscosity ν and u _τ) are used as dispersed phase, and three low mass loading ratios (Φ _m ) ranging from 10^?4 to 10^?3 are tested. It is found that the addition of the particles noticeably modifies the mean velocity and turbulent intensities of the gas-phase, as well as the turbulence coherent structures, even at Φ _m ?=?0.025?%. Particle inertia changes the viscous sublayer of the gas turbulence with a smaller thickness and a larger streamwise velocity gradient, which increases the peak value of the streamwise fluctuation velocity ( $ u_{\text{rms}}^{ + } $ ) of the gas-phase with its location shifting to the wall. Particle sedimentation increases the roughness of the bottom wall, which significantly increases the wall-normal fluctuation velocity ( $ v_{\text{rms}}^{ + } $ ) and Reynolds shear stress ( $ - \langle u^{ \prime } v^{\prime } \rangle^{ + } $ ) of the gas-phase in the inner region of the boundary layer (y ⁺?<?10). Under effect of particle–wall collision, the Q2 events (ejections) of the gas-phase are slightly increased by particles, while the Q4 events (sweeps) are obviously decreased. The spatial scale of the coherent structures near the wall shrinks remarkably with the presence of the particles, which may be attributed to the intensified crossing-trajectory effects due to particle saltation near the bottom wall. Meanwhile, the $ v_{\text{rms}}^{ + } $ and $ - \langle u^{ \prime } v^{\prime } \rangle^{ + } $ of the gas-phase are significantly reduced in the outer region of the boundary layer (y ⁺?>?20).