Flow field in the wake of a bluff body driven through a steady recirculating flow

The wake produced by a bluff body driven through a steady recirculating flow is studied experimentally in a water facility using particle image velocimetry. The bluff body has a rectangular cross section of height, \(H\), and width, \(D\), such that the aspect ratio, AR = H/D, is equal to 3. The motion of the bluff body is uniform and rectilinear, and corresponds to a Reynolds number based on width, Re _D  = 9,600. The recirculating flow is confined within a hemicylindrical enclosure and is generated by planar jets emanating from slots of width, \(h\), such that \(Re_h=500\). Under these conditions, experiments are performed in a closed-loop facility that enables complete optical access to the near-wake. Velocity fields are obtained up to a distance of \(13D\) downstream of the moving body. Data include a selection of phase-averaged velocity fields representative of the wake for a baseline case (no recirculation) and an interaction case (with recirculation). Results indicate that the transient downwash flow typically observed in wakes behind finite bodies of small aspect ratio is significantly perturbed by the recirculating flow. The wake is displaced from the ground plane and exhibits a shorter recirculation zone downstream of the body. In summary, it was found that the interaction between a bluff body wake and a recirculating flow pattern alters profoundly the dynamics of the wake, which has implications on scalar transport in the wake.     

Initial and Boundary Blow-U Problem for $$$$-Lalacian Parabolic Equation with General Absortion

In this article, we investigate the initial and boundary blow-up problem for the \(p\)-Laplacian parabolic equation \(u_t-\Delta _p u=-b(x,t)f(u)\) over a smooth bounded domain \(\Omega \) of \(\mathbb {R}^N\) with \(N\ge 2\), where \(\Delta _pu=\mathrm{div}(|\nabla u|^{p-2}\nabla u)\) with \(p>1\), and \(f(u)\) is a function of regular variation at infinity. We study the existence and uniqueness of positive solutions, and their asymptotic behaviors near the parabolic boundary.     

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

Generic Hamiltonian Dynamics

In this paper we contribute to the generic theory of Hamiltonians by proving that there is a \(C^2\)-residual \({\mathcal {R}}\) in the set of \(C^2\) Hamiltonians on a closed symplectic manifold \(M\), such that, for any \(H\in {\mathcal {R}}\), there is a full measure subset of energies \(e\) in \(H(M)\) such that the Hamiltonian level \((H,e)\) is topologically mixing; moreover these level sets are homoclinic classes.     

Rossby solitary waves excited by the unstable topography in weak shear flow

A new forced KdV equation including topography is derived and the numerical solutions are given. The topographic variable should be related with the temporal and spatial function, which is called unstable topography. The physical features of the solitary waves about the mass and energy are discussed by theoretical analysis. In further studies, the pseudo-spectral numerical methods are used to discuss the evolution of solitary wave generated by the topography when meridional wave number \(m=1\); in a similar way, we analyze the solitary wave when meridional wave number \(m=2\). At last, we make the comparison for the characteristics of waves between \(m=1\) and \(m=2\), the wave of meridional number \(m=1\) plays a leading role.     

Threshold Behavior and Non-quasiconvergent Solutions with Localized Initial Data for Bistable Reaction–Diffusion Equations

We consider bounded solutions of the semilinear heat equation \(u_t=u_{xx}+f(u)\) on \(R\), where \(f\) is of the unbalanced bistable type. We examine the \(\omega \)-limit sets of bounded solutions with respect to the locally uniform convergence. Our goal is to show that even for solutions whose initial data vanish at \(x=\pm \infty \), the \(\omega \)-limit sets may contain functions which are not steady states. Previously, such examples were known for balanced bistable nonlinearities. The novelty of the present result is that it applies to a robust class of nonlinearities. Our proof is based on an analysis of threshold solutions for ordered families of initial data whose limits at infinity are not necessarily zeros of \(f\).     

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.

     

Oriented Shadowing Property and $$\Omega $$-Stability for Vector Fields

We call that a vector field has the oriented shadowing property if for any \(\varepsilon >0\) there is \(d>0\) such that each \(d\)-pseudo orbit is \(\varepsilon \)-oriented shadowed by some real orbit. In this paper, we show that the \(C^1\)-interior of the set of vector fields with the oriented shadowing property is contained in the set of vector fields with the \(\Omega \)-stability.     

Insight into Heterogeneity Effects in Methane Hydrate Dissociation via Pore-Scale Modeling

A large number (1253) of high-quality streaming potential coefficient (\(C_\mathrm{sp})\) measurements have been carried out on Berea, Boise, Fontainebleau, and Lochaline sandstones (the latter two including both detrital and authigenic overgrowth forms), as a function of pore fluid salinity (\(C_\mathrm{f})\) and rock microstructure. All samples were saturated with fully equilibrated aqueous solutions of NaCl (10\(^{-5}\) and 4.5 mol/dm\(^{3})\) upon which accurate measurements of their electrical conductivity and pH were taken. These \(C_\mathrm{sp}\) measurements represent about a fivefold increase in streaming potential data available in the literature, are consistent with the pre-existing 266 measurements, and have lower experimental uncertainties. The \(C_\mathrm{sp}\) measurements follow a pH-sensitive power law behaviour with respect to \(C_\mathrm{f}\) at medium salinities (\(C_\mathrm{sp} =-\,1.44\times 10^{-9} C_\mathrm{f}^{-\,1.127} \), units: V/Pa and mol/dm\(^{3})\) and show the effect of rock microstructure on the low salinity \(C_\mathrm{sp}\) clearly, producing a smaller decrease in \(C_\mathrm{sp}\) per decade reduction in \(C_\mathrm{f}\) for samples with (i) lower porosity, (ii) larger cementation exponents, (iii) smaller grain sizes (and hence pore and pore throat sizes), and (iv) larger surface conduction. The \(C_\mathrm{sp}\) measurements include 313 made at \(C_\mathrm{f} > 1\) mol/dm\(^{3}\), which confirm the limiting high salinity \(C_\mathrm{sp}\) behaviour noted by Vinogradov et al., which has been ascribed to the attainment of maximum charge density in the electrical double layer occurring when the Debye length approximates to the size of the hydrated metal ion. The zeta potential (\(\zeta \)) was calculated from each \(C_\mathrm{sp}\) measurement. It was found that \(\zeta \) is highly sensitive to pH but not sensitive to rock microstructure. It exhibits a pH-dependent logarithmic behaviour with respect to \(C_\mathrm{f}\) at low to medium salinities (\(\zeta =0.01133 \log _{10} \left( {C_\mathrm{f} } \right) +0.003505\), units: V and mol/dm\(^{3})\) and a limiting zeta potential (zeta potential offset) at high salinities of \({\zeta }_\mathrm{o} = -\,17.36\pm 5.11\) mV in the pH range 6–8, which is also pH dependent. The sensitivity of both \(C_\mathrm{sp}\) and \(\zeta \) to pH and of \(C_\mathrm{sp}\) to rock microstructure indicates that \(C_\mathrm{sp}\) and \(\zeta \) measurements can only be interpreted together with accurate and equilibrated measurements of pore fluid conductivity and pH and supporting microstructural and surface conduction measurements for each sample.     

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.     

Blow-Up Rate Estimates and Liouville Type Theorems for a Semilinear Heat Equation with Weighted Source

We study the Liouville-type theorem for the semilinear parabolic equation \(u_t-\Delta u =|x|^a u^p\) with \(p>1\) and \(a\in {\mathbb R}\). Relying on the recent result of Quittner (Math Ann, doi: 10.1007/s00208-015-1219-7, 2015), we establish the optimal Liouville-type theorem in dimension \(N=2\), in the class of nonnegative bounded solutions. We also provide a partial result in dimension \(N\ge 3\). As applications of Liouville-type theorems, we derive the blow-up rate estimates for the corresponding Cauchy problem.     

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.     

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

     

Scalar and vector multipole and vortex solitons in the spatially modulated cubic–quintic nonlinear media

We derive scalar and vector multipole and vortex soliton solutions in the spatially modulated cubic–quintic nonlinear media, which is governed by a (3+1)-dimensional N-coupled cubic–quintic nonlinear Schrödinger equation with spatially modulated nonlinearity and transverse modulation. If the modulation depth \(q=1\), the vortex soliton is constructed, and if \(q=0\), the multipole soliton, including dipole, quadrupole, hexapole, octopole and dodecagon solitons, is constructed, respectively, when the topological charge \(k=1\)–5. If the topological charge \(k=0\), scalar solitons can be obtained. Moreover, the number of layers for the scalar and vector multipole and vortex solitons is decided by the value of the soliton order number n.     

One-Dimensional Super-Fast Diffusion: Persistence Versus Extinction Revisited—Extinction at Spatial Infinity

We consider positive classical solutions of

$$\begin{aligned} v_t=(v^{m-1}v_x)_x, \qquad x\in {\mathbb {R}}, \ t>0, \qquad (\star ) \end{aligned}$$

in the super-fast diffusion range \(m<-1\). Our main interest is in smooth positive initial data \(v_0=v(\cdot ,0)\) which decay as \(x\rightarrow +\infty \), but which are possibly unbounded as \(x\rightarrow -\infty \), having in mind monotonically decreasing data as prototypes. It is firstly proved that if \(v_0\) decays sufficiently fast only in one direction by satisfying

$$\begin{aligned} v_0(x) \le cx^{-\beta } \qquad \text{ for } \text{ all } ~x>0 \quad \hbox { with some }\quad \beta >\frac{2}{1-m} \end{aligned}$$

and some \(c>0\), then the so-called proper solution of (\(\star \)) vanishes identically in \({\mathbb {R}}\times (0,\infty )\), and accordingly no positive classical solution exists in any time interval in this case. Complemented by some sufficient criteria for solutions to remain positive either locally or globally in time, this condition for instantaneous extinction is shown to be optimal at least with respect to algebraic decay of the initial data. This partially extends some known nonexistence results for (\(\star \)) (Daskalopoulos and Del Pino in Arch Rat Mech Anal 137(4):363–380, 1997) in that it does not require any knowledge on the behavior of \(v_0(x)\) for \(x<0\). Next focusing on the phenomenon of extinction in finite time, we show that in this respect a mass influx from \(x=-\infty \) can interact with mass loss at \(x=+\infty \) in a nontrivial manner. Namely, we shall detect examples of monotone initial data, with critical decay as \(x\rightarrow +\infty \) and exponential growth as \(x\rightarrow -\infty \), that lead to solutions of (\(\star \)) which become extinct at a finite positive time, but which have empty extinction sets. This is in sharp contrast to known extinction mechanisms which are such that the corresponding extinction sets coincide with all of \({\mathbb {R}}\).

     

Numerical investigation of mixed convection heat transfer from two isothermal circular cylinders in tandem arrangement: buoyancy,spacing ratio,and confinement effects

This paper presents a two-dimensional numerical study for mixed convection in a laminar cross-flow with a pair of stationary equal-sized isothermal cylinders in tandem arrangement confined in a channel. The governing equations are solved using the control volume method on a nonuniform orthogonal Cartesian grid, and the immersed boundary method is employed to identify the cylinders placed in the flow field. The numerical scheme is first validated against standard cases of symmetrically confined isothermal circular cylinders in plane channels, and grid convergence tests were also examined. The objective of the present study was to investigate the influence of buoyancy and the blockage ratio constraint on the flow and heat transfer characteristics of the immersed cylinder array. Using a fixed Reynolds number based on cylinder diameter of \(Re_{D} = 200\), a fixed value of the Prandtl number of \(Pr = 7\), and a blockage ratio of \(D/H = 0.2\), all possible flow regimes are considered by setting the longitudinal spacing ratio (\(\sigma = L/D\)) between the cylinder axes to 2, 3, and 5 for values of the buoyancy parameter (Richardson number) in the range \(-1\le Ri\le 4\). The interference effects and complex flow features are presented in the form of mean and instantaneous velocity, vorticity, and temperature distributions. The results demonstrate how the buoyancy, spacing ratio, and wall confinement affect the wake structure and vortex dynamics. In addition, local and average heat transfer characteristics of both cylinders are comprehensively presented for a wide range in the parametric space.     

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

     

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.

     

On the Convexity of Nonlinear Elastic Energies in the Right Cauchy-Green Tensor

We present a sufficient condition under which a weak solution of the Euler-Lagrange equations in nonlinear elasticity is already a global minimizer of the corresponding elastic energy functional. This criterion is applicable to energies \(W(F)=\widehat{W}(F^{T}F)=\widehat{W}(C)\) which are convex with respect to the right Cauchy-Green tensor \(C=F^{T}F\), where \(F\) denotes the gradient of deformation. Examples of such energies exhibiting a blow up for \(\det F\to0\) are given.