Long-time Dynamics of Resonant Weakly Nonlinear CGL Equations

Consider a weakly nonlinear CGL equation on the torus \(\mathbb {T}^d\):

$$\begin{aligned} u_t+i\Delta u=\epsilon [\mu (-1)^{m-1}\Delta ^{m} u+b|u|^{2p}u+ ic|u|^{2q}u]. \end{aligned}$$

(*)

Here \(u=u(t,x)\), \(x\in \mathbb {T}^d\), \(0<\epsilon <<1\), \(\mu \geqslant 0\), \(b,c\in \mathbb {R}\) and \(m,p,q\in \mathbb {N}\). Define \(I(u)=(I_{\mathbf {k}},\mathbf {k}\in \mathbb {Z}^d)\), where \(I_{\mathbf {k}}=v_{\mathbf {k}}\bar{v}_{\mathbf {k}}/2\) and \(v_{\mathbf {k}}\), \(\mathbf {k}\in \mathbb {Z}^d\), are the Fourier coefficients of the function \(u\) we give. Assume that the equation \((*)\) is well posed on time intervals of order \(\epsilon ^{-1}\) and its solutions have there a-priori bounds, independent of the small parameter. Let \(u(t,x)\) solve the equation \((*)\). If \(\epsilon \) is small enough, then for \(t\lesssim {\epsilon ^{-1}}\), the quantity \(I(u(t,x))\) can be well described by solutions of an effective equation:

$$\begin{aligned} u_t=\epsilon [\mu (-1)^{m-1}\Delta ^m u+ F(u)], \end{aligned}$$

where the term \(F(u)\) can be constructed through a kind of resonant averaging of the nonlinearity \(b|u|^{2p}+ ic|u|^{2q}u\).

     

Convergence of Radially Symmetric Solutions for (<Emphasis Type="Italic">p</Emphasis>, <Emphasis Type="Italic">q</Emphasis>)-Laplacian Elliptic Equations with a Damping Term

This study considers the quasilinear elliptic equation with a damping term,

$$\begin{aligned} \text {div}(D(u)\nabla u) + \frac{k(|{\mathbf {x}}|)}{|{\mathbf {x}}|}\,{\mathbf {x}}\cdot (D(u)\nabla u) + \omega ^2\big (|u|^{p-2}u + |u|^{q-2}u\big ) = 0, \end{aligned}$$

where \({\mathbf {x}}\) is an N-dimensional vector in \(\big \{{\mathbf {x}} \in \mathbb {R}^N: |{\mathbf {x}}| \ge \alpha \big \}\) for some \(\alpha > 0\) and \(N \in {\mathbb {N}}\setminus \{1\}\); \(D(u) = |\nabla u|^{p-2} + |\nabla u|^{q-2}\) with \(1 < q \le p\); k is a nonnegative and locally integrable function on \([\alpha ,\infty )\); and \(\omega \) is a positive constant. A necessary and sufficient condition is given for all radially symmetric solutions to converge to zero as \(|{\mathbf {x}}|\rightarrow \infty \). Our necessary and sufficient condition is expressed by an improper integral related to the damping coefficient k. The case that k is a power function is explained in detail.

     

Boundedness for Some Doubly Nonlinear Parabolic Equations in Measure Spaces

In the context of measure spaces equipped with a doubling non-trivial Borel measure supporting a Poincaré inequality, we derive local and global sup bounds of the nonnegative weak subsolutions of

$$\begin{aligned} (u^{q})_t-\nabla \cdot {(|\nabla u|^{p-2}\nabla u)}=0, \quad \mathrm {in} \ U_\tau = U \times (\tau _1, \tau _2] , \quad p>1,\quad q>1 \end{aligned}$$

and of its associated Dirichlet problem, respectively. For particular ranges of the exponents p and q, we show that any locally nonnegative weak subsolution, taken in \(Q (\subset \bar{Q}\subset U_\tau )\), is controlled from above by the \(L^\alpha (\bar{Q}) \)-norm, for \(\alpha = \max \{p, q+1\}\). As for the global setting, under suitable assumptions on the boundary datum g and on the initial datum, we obtain sup bounds for u, in \(U \times \{ t\}\), which depend on the \(\sup g\) and on the \(L^{q+1}(U \times (\tau _1, \tau _1+t])\)-norm of \((u-\sup g)_+\), for all \(t \in (0, \tau _2-\tau _1]\). On the critical ranges of p and q, a priori local and global \(L^\infty \) estimates require extra qualitative information on u.

     

A Note on the Convergence of Singularly Perturbed Second Order Potential-Type Equations

In this paper we study the limit as \(\varepsilon \rightarrow 0\) of the singularly perturbed second order equation \(\varepsilon ^2 \ddot{u}_\varepsilon + \nabla _{\!x} V(t,u_\varepsilon (t))=0\), where V(t, x) is a potential. We assume that \(u_0(t)\) is one of its equilibrium points such that \(\nabla _{\!x}V(t,u_0(t))=0\) and \(\nabla _{\!x}^2V(t,u_0(t))>0\). We find that, under suitable initial data, the solutions \(u_\varepsilon \) converge uniformly to \(u_0\), by imposing mild hypotheses on V. A counterexample shows that they cannot be weakened.     

Synchronization transitions induced by partial time delay in a excitatory–inhibitory coupled neuronal network

Information transmission delays are an inherent factor of neuronal systems as a consequence of the finite propagation speeds and time lapses occurring by both dendritic and synaptic processes. In real neuronal systems, some delay between two neurons is too small and can be ignored, which results in partial time delay. In this paper, we focus on investigating influences of partial time delay on synchronization transitions in a excitatory–inhibitory (E–I) coupled neuronal networks. Here, we suppose time delay between two neurons equals to \(\tau \) with probability \(p_{\mathrm{delay}}\) and investigate effect of partial time delay on synchronization transitions of the neuronal networks by controlling \(\tau \) and \(p_{\mathrm{delay}}\) under three cases. In these three cases, excitatory synapses are always considered to delayed with probability \(p_{\mathrm{delay}}\), while inhibitory synapses are considered to be without delays (case I), delayed with probability \(p_{\mathrm{delay}}\) (case II), and always delayed (case III), respectively. It is revealed that, in the first two cases, partial time delay has little influences on synchronization of the neuronal network for small \(p_{\mathrm{delay}}\), while it could induce synchronization transitions at \(\tau \) around integer multiples of the period of individual neuron T when \(p_{\mathrm{delay}}\) is large enough, while in the case III, partial time delay could induce synchronization transitions at \(\tau \) being around odd integer multiples of T / 2 for small \(p_{\mathrm{delay}}\) and at \(\tau \) being around integer multiples of T for large \(p_{\mathrm{delay}}\). Most interesting observation is that partial time delay could induce frequent synchronization transitions at \(\tau \) being around integer multiples of T / 2 for intermediate \(p_{\mathrm{delay}}\). Moreover, effect of rewiring probability on synchronization transitions induced by partial time delay has been discussed. It is found that synchronization transitions induced by partial time delay are robust to rewiring probability for large \(p_{\mathrm{delay}}\) under the three cases.     

Whirl Mappings on Generalised Annuli and the Incompressible Symmetric Equilibria of the Dirichlet Energy

In this paper we show a striking contrast in the symmetries of equilibria and extremisers of the total elastic energy of a hyperelastic incompressible annulus subject to pure displacement boundary conditions. Indeed upon considering the equilibrium equations, here, the nonlinear second order elliptic system formulated for the deformation \(u=(u_{1}, \ldots, u_{N})\):

$$ {\mathbb{E}} {\mathbb{L}}[u, {\mathbf {X}}] = \left \{ \textstyle\begin{array}{l@{\quad}l} \Delta u = \operatorname{div}(\mathscr{P} (x) \operatorname{cof} \nabla u) & \textrm{in }{\mathbf {X}},\\ \det\nabla u = 1 & \textrm{in }{\mathbf {X}},\\ u \equiv\varphi& \textrm{on }\partial{\mathbf {X}}, \end{array}\displaystyle \right . $$

where \({\mathbf {X}}\) is a finite, open, symmetric \(N\)-annulus (with \(N \ge2\)), \(\mathscr{P}=\mathscr{P}(x)\) is an unknown hydrostatic pressure field and \(\varphi\) is the identity mapping, we prove that, despite the inherent rotational symmetry in the system, when \(N=3\), the problem possesses no non-trivial symmetric equilibria whereas in sharp contrast, when \(N=2\), the problem possesses an infinite family of symmetric and topologically distinct equilibria. We extend and prove the counterparts of these results in higher dimensions by way of showing that a similar dichotomy persists between all odd vs. even dimensions \(N \ge4\) and discuss a number of closely related issues.

     

Infinitely Many Stability Switches in a Problem with Sublinear Oscillatory Boundary Conditions

We consider the elliptic equation \(-\Delta u +u =0\) with nonlinear boundary condition \(\frac{\partial u}{\partial n}= \lambda u + g(\lambda ,x,u), \) where \(\frac{g(\lambda ,x,s)}{s} \rightarrow 0, \hbox { as }|s|\rightarrow \infty \) and g is oscillatory. We provide sufficient conditions on g for the existence of unbounded sequences of stable solutions, unstable solutions, and turning points, even in the absence of resonant solutions.     

Multiple Solutions to Neumann Problems with Indefinite Weight and Bounded Nonlinearities

We study the Neumann boundary value problem for the second order ODE

$$\begin{aligned} u^{\prime \prime } + (a^+(t)-\mu a^-(t))g(u) = 0, \qquad t \in [0,T], \end{aligned}$$

(1)

where \(g \in {\mathcal {C}}^1({\mathbb {R}})\) is a bounded function of constant sign, \(a^+,a^-: [0,T] \rightarrow {\mathbb {R}}^+\) are the positive/negative part of a sign-changing weight \(a(t)\) and \(\mu > 0\) is a real parameter. Depending on the sign of \(g^{\prime }(u)\) at infinity, we find existence/multiplicity of solutions for \(\mu \) in a “small” interval near the value

$$\begin{aligned} \mu _c = \frac{\int _0^T a^+(t) \, dt}{\int _0^T a^-(t) \, dt}\,. \end{aligned}$$

The proof exploits a change of variables, transforming the sign-indefinite Eq. (1) into a forced perturbation of an autonomous planar system, and a shooting argument. Nonexistence results for \(\mu \rightarrow 0^+\) and \(\mu \rightarrow +\infty \) are given, as well.

     

New Derived from Anosov Diffeomorphisms with Pathological Center Foliation

In this paper we focused our study on derived from Anosov diffeomorphisms (DA diffeomorphisms ) of the torus \(\mathbb {T}^3,\) it is, an absolute partially hyperbolic diffeomorphism on \(\mathbb {T}^3\) homotopic to a linear Anosov automorphism of the \(\mathbb {T}^3.\) We can prove that if \(f: \mathbb {T}^3 \rightarrow \mathbb {T}^3 \) is a volume preserving DA diffeomorphism homotopic to a linear Anosov A,  such that the center Lyapunov exponent satisfies \(\lambda ^c_f(x) > \lambda ^c_A > 0,\) with x belongs to a positive volume set, then the center foliation of f is non absolutely continuous. We construct a new open class U of non Anosov and volume preserving DA diffeomorphisms, satisfying the property \(\lambda ^c_f(x) > \lambda ^c_A > 0\) for \(m-\)almost everywhere \(x \in \mathbb {T}^3.\) Particularly for every \(f \in U,\) the center foliation of f is non absolutely continuous.     

On the Justification of Viscoelastic Elliptic Membrane Shell Equations

We consider a family of linearly viscoelastic shells with thickness \(2\varepsilon\), clamped along their entire lateral face, all having the same middle surface \(S=\boldsymbol{\theta}(\bar{\omega})\subset \mathbb{R}^{3}\), where \(\omega\subset\mathbb{R}^{2}\) is a bounded and connected open set with a Lipschitz-continuous boundary \(\gamma\). We make an essential geometrical assumption on the middle surface \(S\), which is satisfied if \(\gamma\) and \(\boldsymbol{\theta}\) are smooth enough and \(S\) is uniformly elliptic. We show that, if the applied body force density is \(O(1)\) with respect to \(\varepsilon\) and surface tractions density is \(O(\varepsilon)\), the solution of the scaled variational problem in curvilinear coordinates, \(\boldsymbol{u}( \varepsilon)\), defined over the fixed domain \(\varOmega=\omega\times (-1,1)\) for each \(t\in[0,T]\), converges to a limit \(\boldsymbol{u}\) with \(u_{\alpha}(\varepsilon)\rightarrow u_{\alpha}\) in \(W^{1,2}(0,T,H ^{1}(\varOmega))\) and \(u_{3}(\varepsilon)\rightarrow u_{3}\) in \(W^{1,2}(0,T,L^{2}(\varOmega))\) as \(\varepsilon\to0\). Moreover, we prove that this limit is independent of the transverse variable. Furthermore, the average \(\bar{\boldsymbol{u}}= \frac{1}{2}\int_{-1}^{1} \boldsymbol{u}dx_{3}\), which belongs to the space \(W^{1,2}(0,T, V_{M}( \omega))\), where

$$V_{M}(\omega)=H^{1}_{0}(\omega)\times H^{1}_{0}(\omega)\times L ^{2}(\omega), $$

satisfies what we have identified as (scaled) two-dimensional equations of a viscoelastic membrane elliptic shell, which includes a long-term memory that takes into account previous deformations. We finally provide convergence results which justify those equations.

     

A semi-locally scaled eddy viscosity formulation for LES wall models and flows at high speeds

We show that the mean wall-shear stresses in wall-modeled large-eddy simulations (WMLES) of high-speed flows can be off by up to \(\approx 100\%\) with respect to a DNS benchmark when using the van-Driest-based damping function, i.e., the conventional damping function. Errors in the WMLES-predicted wall-shear stresses are often attributed to the so-called log-layer mismatch, which, albeit also an error in wall-shear stresses \(\tau _\mathrm{w}\), is an error of about \(15\%\). The larger error identified here cannot be removed using the previously developed remedies for the log-layer mismatch. This error may be removed by using the semi-local scaling, i.e., \(l_\nu =\mu /\sqrt{\rho \tau _\mathrm{w}}\), in the damping function, where \(\mu \) and \(\rho \) are the local mean dynamic viscosity and density, respectively.     

Spreading in a Shifting Environment Modeled by the Diffusive Logistic Equation with a Free Boundary

We investigate the influence of a shifting environment on the spreading of an invasive species through a model given by the diffusive logistic equation with a free boundary. When the environment is homogeneous and favourable, this model was first studied in Du and Lin (SIAM J Math Anal 42:377–405, 2010), where a spreading–vanishing dichotomy was established for the long-time dynamics of the species, and when spreading happens, it was shown that the species invades the new territory at some uniquely determined asymptotic speed \(c_0>0\). Here we consider the situation that part of such an environment becomes unfavourable, and the unfavourable range of the environment moves into the favourable part with speed \(c>0\). We prove that when \(c\ge c_0\), the species always dies out in the long-run, but when \(0<c<c_0\), the long-time behavior of the species is determined by a trichotomy described by (a) vanishing, (b) borderline spreading, or (c) spreading. If the initial population is written in the form \(u_0(x)=\sigma \phi (x)\) with \(\phi \) fixed and \(\sigma >0\) a parameter, then there exists \(\sigma _0>0\) such that vanishing happens when \(\sigma \in (0,\sigma _0)\), borderline spreading happens when \(\sigma =\sigma _0\), and spreading happens when \(\sigma >\sigma _0\).     

Influence of hydrodynamic slip on convective transport in flow past a circular cylinder

The presence of a finite tangential velocity on a hydrodynamically slipping surface is known to reduce vorticity production in bluff body flows substantially while at the same time enhancing its convection downstream and into the wake. Here, we investigate the effect of hydrodynamic slippage on the convective heat transfer (scalar transport) from a heated isothermal circular cylinder placed in a uniform cross-flow of an incompressible fluid through analytical and simulation techniques. At low Reynolds (\({\textit{Re}}\ll 1\)) and high Péclet (\({\textit{Pe}}\gg 1\)) numbers, our theoretical analysis based on Oseen and thermal boundary layer equations allows for an explicit determination of the dependence of the thermal transport on the non-dimensional slip length \(l_s\). In this case, the surface-averaged Nusselt number, Nu transitions gradually between the asymptotic limits of \(Nu \sim {\textit{Pe}}^{1/3}\) and \(Nu \sim {\textit{Pe}}^{1/2}\) for no-slip (\(l_s \rightarrow 0\)) and shear-free (\(l_s \rightarrow \infty \)) boundaries, respectively. Boundary layer analysis also shows that the scaling \(Nu \sim {\textit{Pe}}^{1/2}\) holds for a shear-free cylinder surface in the asymptotic limit of \({\textit{Re}}\gg 1\) so that the corresponding heat transfer rate becomes independent of the fluid viscosity. At finite \({\textit{Re}}\), results from our two-dimensional simulations confirm the scaling \(Nu \sim {\textit{Pe}}^{1/2}\) for a shear-free boundary over the range \(0.1 \le {\textit{Re}}\le 10^3\) and \(0.1\le {\textit{Pr}}\le 10\). A gradual transition from the lower asymptotic limit corresponding to a no-slip surface, to the upper limit for a shear-free boundary, with \(l_s\), is observed in both the maximum slip velocity and the Nu. The local time-averaged Nusselt number \(Nu_{\theta }\) for a shear-free surface exceeds the one for a no-slip surface all along the cylinder boundary except over the downstream portion where unsteady separation and flow reversal lead to an appreciable rise in the local heat transfer rates, especially at high \({\textit{Re}}\) and Pr. At a Reynolds number of \(10^3\), the formation of secondary recirculating eddy pairs results in appearance of additional local maxima in \(Nu_{\theta }\) at locations that are in close proximity to the mean secondary stagnation points. As a consequence, Nu exhibits a non-monotonic variation with \(l_s\) increasing initially from its lowermost value for a no-slip surface and then decreasing before rising gradually toward the upper asymptotic limit for a shear-free cylinder. A non-monotonic dependence of the spanwise-averaged Nu on \(l_s\) is observed in three dimensions as well with the three-dimensional wake instabilities that appear at sufficiently low \(l_s\), strongly influencing the convective thermal transport from the cylinder. The analogy between heat transfer and single-component mass transfer implies that our results can directly be applied to determine the dependency of convective mass transfer of a single solute on hydrodynamic slip length in similar configurations through straightforward replacement of Nu and \({\textit{Pr}}\) with Sherwood and Schmidt numbers, respectively.     

Glimpses of Flow Development and Degradation During Type B Drag Reduction by Aqueous Solutions of Polyacrylamide B1120

Flow development and degradation during Type B turbulent drag reduction by 0.10 to 10 wppm solutions of a partially-hydrolysed polyacrylamide B1120 of MW \(=\) 18x10⁶ was studied in a smooth pipe of ID \(=\) 4.60 mm and L/D \(=\) 210 at Reynolds numbers from 10000 to 80000 and wall shear stresses Tw from 8 to 600 Pa. B1120 solutions exhibited facets of a Type B ladder, including segments roughly parallel to, but displaced upward from, the P-K line; those that attained asymptotic maximum drag reduction at low Re√ f but departed downwards into the polymeric regime at a higher retro-onset Re√ f; and segments at MDR for all Re√ f. Axial flow enhancement profiles of S\(^{\prime }\) vs L/D reflected a superposition of flow development and polymer degradation effects, the former increasing and the latter diminishing S\(^{\prime }\) with increasing distance downstream. Solutions that induced normalized flow enhancements S\(^{\prime }\)/S\(^{\prime }_{\mathrm {m}} <\) 0.4 developed akin to solvent, with L_e,p/D \(=\) L_e,n/D \(<\) 42.3, while those at maximum drag reduction showed entrance lengths L_e,m/D \(\sim \) 117, roughly 3 times the solvent L_e,n/D. Degradation kinetics were inferred by first detecting a falloff point (Re√f^∧, S\(^{{\prime }\wedge }\)), of maximum observed flow enhancement, for each polymer solution. A plot of S\(^{{\prime }\wedge }\)vs C revealed S\(^{{\prime }\wedge }\)linear in C at low C, with lower bound [S\(^{\prime }\)] \(=\) 5.0 wppm^??1, and S\(^{{\prime }\wedge }\) independent of C at high C, with upper bound S\(^{\prime }_{\mathrm {m}} =\) 15.9. The ratio S\(^{\prime }\)/S\(^{{\prime }\wedge }\) in any pipe section was interpreted to be the undegraded fraction of original polymer therein. Semi-log plots of (S\(^{\prime }\)/S\(^{{\prime }\wedge }\)) at a section vs transit time from pipe entrance thereto revealed first order kinetics, from which apparent degradation rate constants k_deg s^??1 and entrance severities ?ln(S\(^{\prime }\)/S\(^{{\prime }\wedge }\))₀ were extracted. At constant C, k_deg increased linearly with increasing wall shear stress Tw, and at constant Tw, k_deg was independent of C, providing a B1120 degradation modulus (k_deg/Tw) \(=\) (0.012 \(\pm \) 0.001) (Pa s)^??1 for 8 \(<\) Tw Pa \(<\) 600, 0.30 \(<\) C wppm \(<\) 10. Entrance severities were negligible below a threshold Twe \(\sim \) 30 Pa and increased linearly with increasing Tw for Tw \(>\) Twe. The foregoing methods were applied to Type A drag reduction by 0.10 to 10 wppm solutions of a polyethyleneoxide PEO P309, MW \(=\) 11x10⁶, in a smooth pipe of ID \(=\) 7.77 mm and L/D \(=\) 220 at Re from 4000 to 115000. P309 solutions that induced S\(^{\prime }\)/S\(^{\prime }_{\mathrm {m}} <\) 0.4 developed akin to solvent, with L_e,p/D \(=\) L_e,n/D \(<\) 23, while those at MDR had entrance lengths L_e,m/D \(\sim \) 93, roughly 4 times the solvent L_e,n/D. P309 solutions described a Type A fan distorted by polymer degradation. A typical trajectory departed the P-K line at an onset point Re√ f* followed by ascending and descending polymeric regime segments separated by a falloff point Re√f^∧, of maximum flow enhancement; for all P309 solutions, onset Re√ f* = 550 \(\pm \) 100 and falloff Re√f^∧ = 2550 \(\pm \) 250, the interval between them delineating Type A drag reduction unaffected by degradation. A plot of falloff S\(^{{\prime }\wedge }\) vs C for PEO P309 solutions bore a striking resemblance to the analogous S\(^{{\prime }\wedge }\) vs C plot for solutions of PAMH B1120, indicating that the initial Type A drag reduction by P309 after onset at Re√ f* had evolved to Type B drag reduction by falloff at Re√f^∧. Presuming that Type B behaviour persisted past falloff permitted inference of P309 degradation kinetics; k_deg was found to increase linearly with increasing Tw at constant C and was independent of C at constant Tw, providing a P309 degradation modulus (k_deg/Tw) \(=\) (0.011 \(\pm \) 0.002) (Pa s)^??1 for 4 \(<\) Tw Pa \(<\) 400, 0.10 \(<\) C wppm < 5.0. Comparisons between the present degradation kinetics and previous literature showed (k_deg/Tw) data from laboratory pipes of D \(\sim \) 0.01 m to lie on a simple extension of (k_deg/Tw) data from pipelines of D \(\sim \) 0.1 m and 1.0 m, along a power-law relation (k_deg/Tw) \(=\) 10^??5.4.D^??1.6. Intrinsic slips derived from PAMH B1120 and PEO P309-at-falloff experiments were compared with previous examples from Type B drag reduction by polymers with vinylic and glycosidic backbones, showing: (i) For a given polymer, [S\(^{\prime }\)] was independent of Re√ f and pipe ID, implying insensitivity to both micro- and macro-scales of turbulence; and (ii) [S\(^{\prime }\)] increased linearly with increasing polymer chain contour length Lc, the proportionality constant \(\beta =\) 0.053 \(\pm \) 0.036 enabling estimation of flow enhancement S\(^{\prime } =\) C.Lc.β for all Type B drag reduction by polymers.     

Complex Dynamics in a ODE Model Related to Phase Transition

Motivated by some recent studies on the Allen–Cahn phase transition model with a periodic nonautonomous term, we prove the existence of complex dynamics for the second order equation

$$\begin{aligned} -\ddot{x} + \left( 1 + \varepsilon ^{-1} A(t)\right) G'(x) = 0, \end{aligned}$$

where A(t) is a nonnegative T-periodic function and \(\varepsilon > 0\) is sufficiently small. More precisely, we find a full symbolic dynamics made by solutions which oscillate between any two different strict local minima \(x_0\) and \(x_1\) of G(x). Such solutions stay close to \(x_0\) or \(x_1\) in some fixed intervals, according to any prescribed coin tossing sequence. For convenience in the exposition we consider (without loss of generality) the case \(x_0 =0\) and \(x_1 = 1\).

     

Stationary Solutions and Connecting Orbits for <Emphasis Type="Italic">p</Emphasis>-Laplace Equation

We deal with one dimensional p-Laplace equation of the form

$$\begin{aligned} u_t = (|u_x|^{p-2} u_x )_x + f(x,u), \ x\in (0,l), \ t>0, \end{aligned}$$

under Dirichlet boundary condition, where \(p>2\) and \(f:[0,l]\times {\mathbb {R}}\rightarrow {\mathbb {R}}\) is a continuous function with \(f(x,0)=0\). We will prove that if there is at least one eigenvalue of the p-Laplace operator between \(\lim _{u\rightarrow 0} f(x,u)/|u|^{p-2}u\) and \(\lim _{|u|\rightarrow +\infty } f(x,u)/|u|^{p-2}u\), then there exists a nontrivial stationary solution. Moreover we show the existence of a connecting orbit between stationary solutions. The results are based on Conley index and detect stationary states even when those based on fixed point theory do not apply. In order to compute the Conley index for nonlinear semiflows deformation along p is used.

     

Existence of Traveling Waves of General Gray-Scott Models

This work gives a rigorous proof of the existence of propagating traveling waves of a nonlinear reaction–diffusion system which is a general Gray-Scott model of the pre-mixed isothermal autocatalytic chemical reaction of order m (\(m > 1\)) between two chemical species, a reactant A and an auto-catalyst B, \( A + m B \rightarrow (m+1) B\), and a super-linear decay of order \( n > 1\), \( B \rightarrow C\), where \( 1< n < m\). Here C is an inert product. Moreover, we establish that the speed set for existence must lie in a bounded interval for a given initial value \(u_0\) at \( - \infty \). The explicit bound is also derived in terms of \(u_0\) and other parameters. The same system also appears in a mathematical model of SIR type in infectious diseases.     

Statistical Analysis of Scalar Dissipation Rate Transport in Turbulent Partially Premixed Flames: A Direct Numerical Simulation Study

Statistically planar turbulent premixed and partially premixed flames for different initial turbulence intensities are simulated for global equivalence ratios ??>?=?0.7 and ??>?=?1.0 using three-dimensional Direct Numerical Simulations (DNS) with simplified chemistry. For the simulations of partially premixed flames, a random distribution of equivalence ratio following a bimodal distribution of equivalence ratio is introduced in the unburned reactants ahead of the flame. The simulation parameters in all of the cases were chosen such that the combustion situation belongs to the thin reaction zones regime. The DNS data has been used to analyse the behaviour of the dissipation rate transports of both active and passive scalars (i.e. the fuel mass fraction Y _F and the mixture fraction ξ) in the context of Reynolds Averaged Navier–Stokes (RANS) simulations. The behaviours of the unclosed terms of the Favre averaged scalar dissipation rates of fuel mass fraction and mixture fraction (i.e. \(\widetilde {\varepsilon }_Y =\overline {\rho D\nabla Y_F^{\prime \prime } \cdot \nabla Y_F^{\prime \prime } } /\overline{\rho }\) and \(\widetilde {\varepsilon }_\xi =\overline {\rho D\nabla \xi ^{\prime \prime }\cdot \nabla \xi ^{\prime \prime }} /\overline {\rho })\) transport equations have been analysed in detail. In the case of the \(\widetilde {\varepsilon }_Y \) transport, it has been observed that the turbulent transport term of scalar dissipation rate remains small throughout the flame brush whereas the terms due to density variation, scalar–turbulence interaction, reaction rate and molecular dissipation remain the leading order contributors. The term arising due to density variation remains positive throughout the flame brush and the combined contribution of the reaction and molecular dissipation to the \(\widetilde {\varepsilon }_Y \) transport remains negative throughout the flame brush in all cases. However, the behaviour of scalar–turbulence interaction term of the \(\widetilde {\varepsilon }_Y \) transport equation is significantly affected by the relative strengths of turbulent straining and the straining due to chemical heat release. In the case of the \(\widetilde {\varepsilon }_\xi \) transport, the turbulent transport term remains small throughout the flame brush and the density variation term is found to be negligible in all cases, whilst the reaction rate term is exactly zero. The scalar–turbulence interaction term and molecular dissipation term remain the leading order contributors to the \(\widetilde {\varepsilon }_\xi \) transport throughout the flame brush in all cases that have been analysed in the present study. Performances of existing models for the unclosed terms of the transport equations of \(\widetilde {\varepsilon }_Y \) and \(\widetilde {\varepsilon }_\xi \) are assessed with respect to the corresponding quantities obtained from DNS data. Based on this exercise either suitable models have been identified or new models have been proposed for the accurate closure of the unclosed terms of both \(\widetilde {\varepsilon }_Y \) and \(\widetilde {\varepsilon }_\xi \) transport equations in the context of Reynolds Averaged Navier–Stokes (RANS) simulations.     

Exact solution and invariant for fractional Cattaneo anomalous diffusion of cells in two-dimensional comb framework

This paper presents an investigation on anomalous diffusion of cells in a two-dimensional comb framework with effects of fractional Cattaneo flux. Formulated governing equation is an evolution equation with the coexisting characteristics of parabolic (diffusion) and hyperbolic (wave) for \(\alpha \) in (0, 1). Exact solution is obtained by the special fractional integral transformations, and a novel invariant is established, i.e., \(\left\langle {x^{2}\left( t \right) } \right\rangle \cdot \left\langle P \right\rangle = 0.5\) (the mean square displacement multiplied by the total number of cells along the x-axis = 0.5). Moreover, the characteristics of cells distribution, the total number and the mean square displacement of cells along the x-axis with different involved parameters, especially with the fractional parameter evolution, are shown graphically and analyzed in detail. For the cells distribution versus x, it turns from parabolic and hyperbolic with the decrease in t or the increase in \(\alpha \) or \(\xi \). It is monotonically decreasing for the cells distribution versus \(\alpha \) with different x, t and \(\xi \). For the distribution versus t with different \(\alpha \) and \(\xi \) or versus \(\alpha \) with different t, it is monotonically decreasing for the distribution of total number while monotonically increasing for the distribution of mean square displacement. It is remarkable that the anomalous subdiffusion happens along the x-axis for arbitrary parameters which is different from the classical Cattaneo diffusion.     

Analysis on nonlinear turning motion of multi-spherical soft robots

A modular multi-spherical soft robot, which consists of five deformable spherical cells, two friction feet, the electromagnetic valves and the control systems, is constructed. According to the deflating action and the inflating action of the spherical cells, the size and the shape of each spherical cell can be changed. With two friction feet sticking with the ground in turn, the soft robot can move forwards, make a turning motion and avoid the obstacle. This paper creates a nonlinear relation between the pressure P and the inflation radius \(\left( r \right) \) at different original radii \(\left( {r_0 } \right) \) and obtains the inflation or deflation velocity \(v_r \). Six inflating and deflating steps to finish the turning motion are presented. Based on the geometric relationship between the inflation radius (r) and the original radius \((r_0 )\) of each cell, the nonlinear turning process is described to control the center positions (x, y, z) of the spherical cell. Last, a simulation and an experiment of five spherical cells are shown to emulate the turning process. Experiment results show that the robot has a maximum turning capability of \(20{^{\circ }}\) in one period.