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.     

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.     

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.     

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.     

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

     

FEM solution to natural convection flow of a micropolar nanofluid in the presence of a magnetic field

The two-dimensional, laminar, unsteady natural convection flow in a square enclosure filled with aluminum oxide (\(\hbox {Al}_{2} \hbox {O}_{3}\))–water nanofluid under the influence of a magnetic field, is considered numerically. The nanofluid is considered as Newtonian and incompressible, the nanoparticles and water are assumed to be in thermal equilibrium. The mathematical modelling results in a coupled nonlinear system of partial differential equations. The equations are solved using finite element method (FEM) in space, whereas, the implicit backward difference scheme is used in time direction. The results are obtained for Rayleigh (Ra), Hartmann (Ha) numbers, and nanoparticles volume fractions (\(\phi\)), in the ranges of \(10^3 \le Ra \le 10^7\), \(0\le Ha \le 500\) and \(0 \le \phi \le 0.2\), respectively. The streamlines and microrotation contours are observed to show similar behaviors with altering magnitudes. For low Ra values, when \(Ha=0\), symmetric vortices near the walls and a central vortex in opposite direction are observed in vorticity. As Ra increases, the central vortex splits into two due to the circulation in the effect of the buoyant flow. Boundary layer formation is observed when Ha increases for almost all Rayleigh numbers in both streamlines and vorticity. The isotherms have horizontal profiles for high Ra values owing to convective dominance over conduction. As Ha is increased, the convection effect is reduced, and isotherms tend to have vertical profiles. This study presents the first FEM application for solving highly nonlinear PDEs defining micropolar nanofluid flow especially for large values of Rayleigh and Hartmann numbers.     

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.     

On the formation of hidden chaotic attractors and nested invariant tori in the Sprott A system

We consider the well-known Sprott A system, which depends on a single real parameter a and, for \(a=1\), was shown to present a hidden chaotic attractor. We study the formation of hidden chaotic attractors as well as the formation of nested invariant tori in this system, performing a bifurcation analysis by varying the parameter a. We prove that, for \(a=0\), the Sprott A system has a line of equilibria in the z-axis, the phase space is foliated by concentric invariant spheres with two equilibrium points located at the south and north poles, and each one of these spheres is filled by heteroclinic orbits of south pole–north pole type. For \(a\ne 0\), the spheres are no longer invariant algebraic surfaces and the heteroclinic orbits are destroyed. We do a detailed numerical study for \(a>0\) small, showing that small nested invariant tori and a limit set, which encompasses these tori and is the \(\alpha \)- and \(\omega \)-limit set of almost all orbits in the phase space, are formed in a neighborhood of the origin. As the parameter a increases, this limit set evolves into a hidden chaotic attractor, which coexists with the nested invariant tori. In particular, we find hidden chaotic attractors for \(a<1\). Furthermore, we make a global analysis of Sprott A system, including the dynamics at infinity via the Poincaré compactification, showing that for \(a>0\), the only orbit which escapes to infinity is the one contained in the z-axis and all other orbits are either homoclinic to a limit set (or to a hidden chaotic attractor, depending on the value of a), or contained on an invariant torus, depending on the initial condition considered.     

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.

     

Topological Pressure of Generic Points for Amenable Group Actions

Let (X, G) be a G-action topological dynamical system (t.d.s. for short), where G is a countably infinite discrete amenable group. In this paper, we study the topological pressure of the sets of generic points. We show that when the system satisfies the almost specification property, for any G-invariant measure \(\mu \) and any continuous map \(\varphi \),

$$\begin{aligned} P\left( X_{\mu },\varphi ,\{F_n\}\right) = h_{\mu }(X)+\int \varphi d\mu , \end{aligned}$$

where \(\{F_n\}\) is a Følner sequence, \(X_{\mu }\) is the set of generic points of \(\mu \) with respect to (w.r.t. for short) \(\{F_n\}\), \(P(X_{\mu },\varphi ,\{F_n\})\) is the topological pressure of \(X_{\mu }\) for \(\varphi \) w.r.t. \(\{F_n\}\) and \(h_{\mu }(X)\) is the measure-theoretic entropy.

     

Linear Stability of the Double-Diffusive Convection in a Horizontal Porous Layer with Open Top: Soret and Viscous Dissipation Effects

The linear stability of the double-diffusive convection in a horizontal porous layer is studied considering the upper boundary to be open. A horizontal temperature gradient is applied along the upper boundary. It is assumed that the viscous dissipation and Soret effect are significant in the medium. The governing parameters are horizontal Rayleigh number (\(Ra_\mathrm{H}\)), solutal Rayleigh number (\(Ra_\mathrm{S}\)), Lewis number (Le), Gebhart number (Ge) and Soret parameter (Sr). The Rayleigh number (Ra) corresponding to the applied heat flux at the bottom boundary is considered as the eigenvalue. The influence of the solutal gradient caused due to the thermal diffusion on the double-diffusive instability is investigated by varying the Soret parameter. A horizontal basic flow is induced by the applied horizontal temperature gradient. The stability of this basic flow is analyzed by calculating the critical Rayleigh number (\(Ra_\mathrm{cr}\)) using the Runge–Kutta scheme accompanied by the Shooting method. The longitudinal rolls are more unstable except for some special cases. The Soret parameter has a significant effect on the stability of the flow when the upper boundary is at constant pressure. The critical Rayleigh number is decreasing in the presence of viscous dissipation except for some positive values of the Soret parameter. How a change in Soret parameter is attributing to the convective rolls is presented.     

Global Attractors of Sixth Order PDEs Describing the Faceting of Growing Surfaces

A spatially two-dimensional sixth order PDE describing the evolution of a growing crystalline surface h(x, y, t) that undergoes faceting is considered with periodic boundary conditions, as well as its reduced one-dimensional version. These equations are expressed in terms of the slopes \(u_1=h_{x}\) and \(u_2=h_y\) to establish the existence of global, connected attractors for both equations. Since unique solutions are guaranteed for initial conditions in \(\dot{H}^2_{per}\), we consider the solution operator \(S(t): \dot{H}^2_{per} \rightarrow \dot{H}^2_{per}\), to gain our results. We prove the necessary continuity, dissipation and compactness properties.     

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.

     

Delta Shock Waves in the Shallow Water System

We consider a Riemann problem for the shallow water system \(u_{t} +\big (v+\textstyle \frac{1}{2}u^{2}\big )_{x}=0\), \(v_{t}+\big (u+uv\big )_{x}=0\) and evaluate all singular solutions of the form \(u(x,t)=l(t)+b(t)H\big (x-\gamma (t)\big )+a(t)\delta \big (x-\gamma (t)\big )\), \(v(x,t)=k(t)+c(t)H\big (x-\gamma (t)\big )\), where \(l,b,a,k,c,\gamma :\mathbb {R}\rightarrow \mathbb {R}\) are \(C^{1}\)-functions of time t, H is the Heaviside function, and \(\delta \) stands for the Dirac measure with support at the origin. A product of distributions, not constructed by approximation processes, is used to define a solution concept, that is a consistent extension of the classical solution concept. Results showing the advantage of this framework are briefly presented in the introduction.     

Biorthogonal stretching and shearing of an impermeable surface in a uniformly rotating fluid system

The flow induced by an impermeable flat surface executing orthogonal stretching and orthogonal shearing in a rotating fluid system is investigated. Both the stretching and shearing are linear in the coordinates. An exact similarity reduction of the Navier–Stokes equations gives rise to a pair of nonlinearly-coupled ordinary differential equations governed by three parameters. In this study we set one parameter and analyze the problem which leads to flow for an impermeable surface with shearing and stretching due to velocity u along the x-axis of equal strength a while the shearing and stretching due to velocity v along the y-axis of equal strength b. These solutions depend on two parameters—a Coriolis (rotation) parameter \(\sigma = \Omega /a\) and a stretching/shearing ratio \(\lambda =b/a\). A symmetry in solutions is found for \(\lambda = 1\). The exact solution for \(\sigma = 0\) and the asymptotic behavior of solutions for \(|\sigma | \rightarrow \infty\) are determined and compared with numerical results. Oscillatory solutions are found whose strength increases with increasing values of \(|\sigma |\). It is shown that these solutions tend to the well-known Ekman solution as \(|\sigma | \rightarrow \infty\).     

On Entropy of Graph Maps That Give Hereditarily Indecomposable Inverse Limits

We prove that if \(f:G\rightarrow G\) is a map on a topological graph G such that the inverse limit \(\varprojlim (G,f)\) is hereditarily indecomposable, and entropy of f is positive, then there exists an entropy set with infinite topological entropy. When G is the circle and the degree of f is positive then the entropy is always infinite and the rotation set of f is nondegenerate. This shows that the Anosov-Katok type constructions of the pseudo-circle as a minimal set in volume-preserving smooth dynamical systems, or in complex dynamics, obtained previously by Handel, Herman and Chéritat cannot be modeled on inverse limits. This also extends a previous result of Mouron who proved that if \(G=[0,1]\), then \(h(f)\in \{0,\infty \}\), and combined with a result of Ito shows that certain dynamical systems on compact finite-dimensional Riemannian manifolds must either have zero entropy on their invariant sets or be non-differentiable.     

A Reaction–Diffusion–Advection Equation with Mixed and Free Boundary Conditions

We investigate a reaction–diffusion–advection equation of the form \(u_t-u_{xx}+\beta u_x=f(u)\) \((t>0,\,0<x<h(t))\) with mixed boundary condition at \(x=0\) and Stefan free boundary condition at \(x=h(t)\). Such a model may be applied to describe the dynamical process of a new or invasive species adopting a combination of random movement and advection upward or downward along the resource gradient, with the free boundary representing the expanding front. The goal of this paper is to understand the effect of advection environment and no flux across the left boundary on the dynamics of this species. For the case \(|\beta |<c_0\), we first derive the spreading–vanishing dichotomy and sharp threshold for spreading and vanishing, and then provide a much sharper estimate for the spreading speed of h(t) and the uniform convergence of u(t, x) when spreading happens. For the case \(|\beta |\ge c_0\), some results concerning virtual spreading, vanishing and virtual vanishing are obtained. Here \(c_0\) is the minimal speed of traveling waves of the differential equation.     

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.     

A Shilnikov Phenomenon Due to State-Dependent Delay,by Means of the Fixed Point Index

The first part of this paper is a general approach towards chaotic dynamics for a continuous map \(f:X\supset M\rightarrow X\) which employs the fixed point index and continuation. The second part deals with the differential equation

$$\begin{aligned} x'(t)=-\alpha \,x(t-d_{{\varDelta }}(x_t)). \end{aligned}$$

with state-dependent delay. For a suitable parameter \(\alpha \) close to \(5\pi /2\) we construct a delay functional \(d_{{\varDelta }}\), constant near the origin, so that the previous equation has a homoclinic solution, \(h(t)\rightarrow 0\) as \(t\rightarrow \pm \infty \), with certain regularity properties of the linearization of the semiflow along the flowline \(t\mapsto h_t\). The third part applies the method from the beginning to a return map which describes solution behaviour close to the homoclinic loop, and yields the existence of chaotic motion.