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.     

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.     

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

     

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.     

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.

     

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.     

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.     

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 Note on the Relative Equilibria Bifurcations in the $$$$-Body Problem

Consider the planar Newtonian \((2N+1)\)-body problem, \(N\ge 1,\) with \(2N\) bodies of unit mass and one body of mass \(m\). Using the discrete symmetry due to the equal masses and reducing by the rotational symmetry, we show that solutions with the \(2N\) unit mass points at the vertices of two concentric regular \(N\)-gons and \(m\) at the centre at all times form invariant manifold. We study the regular \(2N\)-gon with central mass \(m\) relative equilibria within the dynamics on the invariant manifold described above. As \(m\) varies, we identify the bifurcations, relate our results to previous work and provide the spectral picture of the linearization at the relative equilibria.     

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

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.

     

Vortex solitons of the (3+1)-dimensional spatially modulated cubic–quintic nonlinear Schrödinger equation with the transverse modulation

Vortex solitons in the spatially modulated cubic–quintic nonlinear media are governed by a (3+1)-dimensional cubic–quintic nonlinear Schrödinger equation with spatially modulated nonlinearity and transverse modulation. Via the variable separation principle with the similarity transformation, we derive two families of vortex soliton solutions in the spatially modulated cubic–quintic nonlinear media. For the disappearing and parabolic transverse modulation, vortex solitons with different configurations are constructed. The similar configurations of vortex solitons exist for the same value of \(l-k\) with the topological charge k and degree number l. Moreover, the number of the inner layer structure of vortex solitons getting rid of the package covering layer is related to \((n-1)/2+1\) with the soliton order number n. For the disappearing transverse modulation, there exist phase azimuthal jumps around their cores of vortex solitons with \(2\pi \) phase change in every jump, and any two jumps one after another realize the change in \(\pi \). For the parabolic transverse modulation, all phases of vortex soliton exist k-jump, and every jump realizes the change in \(2\pi /k\); thus, k-jumps totally realize the azimuthal change in \(2\pi \) around their cores.     

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.     

The Effect of Grain Size Distribution on Nonlinear Flow Behavior in Sandy Porous Media

The current study provides new experimental data on nonlinear flow behavior in various uniformly graded granular materials (20 samples) ranging from medium sands (\(d_{50 }>0.39\) mm) to gravel (\(d_{50}=6.3\) mm). Generally, theoretical equations relate the Forchheimer parameters a and b to the porosity, as well as the characteristic pore length, which is assumed to be the median grain size \((d_{50})\) of the porous medium. However, numerical and experimental studies show that flow resistance in porous media is largely determined by the geometry of the pore structure. In this study, the effect of the grain size distribution was analyzed using subangular-subrounded sands and approximately equal compaction grades. We have used a reference dataset of 11 uniformly graded filter sands. Mixtures of filter sands were used to obtain a slightly more well-graded composite sand (increased \(C_{u}\) values by a factor of 1.19 up to 2.32) with respect to its associated reference sand at equal median grain size \((d_{50})\) and porosity. For all composite sands, the observed flow resistance was higher than in the corresponding reference sand at equal \(d_{50}\), resulting in increased a coefficients by factors up to 1.68, as well as increased b coefficients by factors up to 1.44. A modified Ergun relationship with Ergun constants of 139.1 for A and 2.2 for B, as well as the use of \(d_{m}-\sigma \) as characteristic pore length predicted the coefficients a and b accurately.     

Monte Carlo Assessment of the Impact of Oscillatory and Pulsating Boundary Conditions on the Flow Through Porous Media

The linear stability analysis of vertical throughflow of power law fluid for double-diffusive convection with Soret effect in a porous channel is investigated in this study. The upper and lower boundaries are assumed to be permeable, isothermal and isosolutal. The linear stability of vertical through flow is influenced by the interactions among the non-Newtonian Rayleigh number (Ra), Buoyancy ratio (N), Lewis number (Le), Péclet number (Pe), Soret parameter (Sr) and power law index (n). The results indicate that the Soret parameter has a significant influence on convective instability of power law fluid. It has also been noticed that buoyancy ratio has a dual effect on the instability of fluid flow. Further, it is noticed that the basic temperature and concentration profiles have singularities at \(Pe = 0\) and \(Le = 1\), the convective instability is looked into for the limiting case of \(Pe\rightarrow 0\) and \(Le \rightarrow 1\). For the case of pure thermal convection with no vertical throughflow, the present numerical results coincide with the solution of standard Horton–Rogers–Lapwood problem. The present results for critical Rayleigh number obtained using bvp4c and two-term Galerkin approximation are compared with those available in the literature and are tabulated.     

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.     

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.     

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.     

Existence of $$C^{k}$$-Invariant Foliations for Lorenz-Type Maps

Under conditions similar to those in Shashkov and Shil’nikov (Differ Uravn 30(4):586–595, 732, 1994) we show that a \(C^{k+1}\) Lorenz-type map T has a \(C^{k}\) codimension one foliation which is invariant under the action of T. This allows us to associate T to a \(C^{k}\) one-dimensional transformation.     

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.