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.     

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.     

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.     

Physical Clogging of Uniformly Graded Porous Media Under Constant Flow rates

This study investigated the physical clogging of uniformly graded porous media under constant flow rates using natural porous media and suspensions. The porous media selected for this experimental study was a fine-to-medium sandy soil fractioned into thirteen uniformly graded beds: seven unisize beds and six uniform beds. The physical clogging of the beds was studied using two types of silt suspensions as along with two suspension concentrations and three water discharges. It was found that the permeability reduction due to physical clogging \([(K_\mathrm{i} - K_\mathrm{t})/K_\mathrm{i}]\) increased with decreasing \({D}_{15}/{d}_{85}\) ratios until a critical value of \({D}_{15}/{d}_{85}\), after which a surface mat of suspension was formed on the porous media. It was also found that the value of reduced permeability at any time (at any number of pore volumes of injected suspension-laden water), \(K_\mathrm{t}\), is directly proportional to square of \({D}_{15}\) and inversely proportional to \({C}_{\mathrm{u}}\) of the porous media and \({d}_{85}\) of suspensions. The effects of suspension type and flow rates on physical clogging seemed to depend on the size of the pores in the porous media.     

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.     

Approximation of Flat <Emphasis Type="Italic">W</Emphasis><Superscript>2,2</Superscript> Isometric Immersions by Smooth Ones

Let \({S\subset\mathbb{R}^2}\) be a bounded Lipschitz domain and denote by \({W^{2,2}_{\text{iso}}(S; \mathbb{R}^3)}\) the set of mappings \({u\in W^{2,2}(S;\mathbb{R}^3)}\) which satisfy \({(\nabla u)^T(\nabla u) = Id}\) almost everywhere. Under an additional regularity condition on the boundary \({\partial S}\) (which is satisfied if \({\partial S}\) is piecewise continuously differentiable), we prove that the strong W ^2,2 closure of \({W^{2,2}_{\text{iso}}(S; \mathbb{R}^3)\cap C^{\infty}(\overline{S};\mathbb{R}^3)}\) agrees with \({W^{2,2}_{\text{iso}}(S; \mathbb{R}^3)}\).     

A Simple Method for the Determination of Adsorption Kinetic Parameters Using Recycle Fixed-Bed Adsorber

This study focuses on the development of a novel analysis technique for determining the intraparticle diffusivity \((D_{\mathrm{S}})\) and fluid film mass transfer coefficient \((k_\mathrm{F})\) from a concentration history curve of a recycle fixed-bed reactor without using the linear driving force approximation and empirical equations for the estimation of \(k_\mathrm{F}\). The recycle fixed-bed method requires lesser amounts of working fluid for experiment purposes, which is an advantage over the usual fixed-bed method. Based on the characterization results of the concentration history curves, simulated by rigorous numerical calculations, a novel analysis technique was established. The \(D_{\mathrm{S}}\) value can be determined from the experimentally obtained time at which the concentration of the curve is minimal \((t_{C\mathrm{min}})\). The \(k_\mathrm{F}\) value can be determined from the \(D_{\mathrm{S}}\) value and Biot number (Bi), which can be estimated from the experimental ratio of the maximum concentration to the minimum concentration \((c_{\max }/c_{\min })\). The \(D_{\mathrm{S}}\) and \(k_\mathrm{F}\) values of phenol adsorption on an activated carbon material were determined experimentally using the proposed analysis method. This method enables the determination of reliable adsorption kinetic parameters through a simple and economical experiment. However, appropriate experimental data must be acquired under regulated experimental conditions, especially in the case of fluctuation of the concentration history curve.     

Exponential Decay of the Vorticity in the Steady-State Flow of a Viscous Liquid Past a Rotating Body

We construct a Sobolev homeomorphism in dimension \({n \geqq 4,\,f \in W^{1,1}((0, 1)^n,\mathbb{R}^n)}\) such that \({J_f = {\rm det} Df > 0}\) on a set of positive measure and J _f  < 0 on a set of positive measure. It follows that there are no diffeomorphisms (or piecewise affine homeomorphisms) f _k such that \({f_k\to f}\) in \({W^{1,1}_{\rm loc}}\).     

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.     

Exact solution of 3D Timoshenko beam problem using linked interpolation of arbitrary order

For arbitrary polynomial loading and a sufficient finite number of nodal points N, the solution for the 3D Timoshenko beam differential equations is polynomial and given as \({{\varvec \theta} = \sum_{i=1}^N I_i {\varvec \theta}_i}\) for the rotation field and \({{\bf u} = \sum_{i=1}^{N+1} J_i {\bf u}_i}\) for the displacement field, where I _i and J _i are the Lagrangian polynomials of order N?1 and N, respectively. It has been demonstrated in this work that the exact solution for the displacement field may be also written in a number of alternative ways involving contributions of the nodal rotations including \({{\bf u} = \sum_{i=1}^N I_i \left[ {\bf u}_i + \frac 1 N ( {\varvec \theta} - {\varvec \theta}_i ) \times {\bf R}_i \right]}\), where R _i are the beam nodal positions.     

Fine Level Set Structure of Flat Isometric Immersions

A result by Pogorelov asserts that C ¹ isometric immersions u of a bounded domain \({S \subset \mathbb R^2}\) into \({\mathbb {R}^3}\) whose normal takes values in a set of zero area enjoy the following regularity property: the gradient \({f := \nabla u}\) is ‘developable’ in the sense that the nondegenerate level sets of f consist of straight line segments intersecting the boundary of S at both endpoints. Motivated by applications in nonlinear elasticity, we study the level set structure of such f when S is an arbitrary bounded Lipschitz domain. We show that f can be approximated by uniformly bounded maps with a simplified level set structure. We also show that the domain S can be decomposed (up to a controlled remainder) into finitely many subdomains, each of which admits a global line of curvature parametrization.     

Uniqueness of Positive Ground State Solutions of the Logarithmic Schrödinger Equation

We prove the uniqueness of positive ground state solutions of the problem \({ {\frac {d^{2}u}{dr^{2}}} + {\frac {n-1}{r}}{\frac {du}{dr}} + u \ln(|u|) = 0}\), \({u(r) > 0~\forall r \ge 0}\), and \({(u(r),u'(r)) \to (0, 0)}\) as \({r \to \infty}\). This equation is derived from the logarithmic Schrödinger equation \({{\rm i}\psi_{t} = {\Delta} \psi + u \ln \left(|u|^{2}\right)}\), and also from the classical equation \({{\frac {\partial u}{\partial t}} = {\Delta} u +u \left(|u|^{p-1}\right) -u}\). For each \({n \ge 1}\), a positive ground state solution is \({ u_{0}(r) = \exp \left(-{\frac{r^2}{4}} + {\frac{n}{2}}\right),~0 \le r < \infty}\). We combine \({u_{0}(r)}\) with energy estimates and associated Ricatti equation estimates to prove that, for each \({n \in \left[1, 9 \right]}\), \({u_{0}(r)}\) is the only positive ground state. We also investigate the stability of \({u_{0}(r)}\). Several open problems are stated.     

Entire Solutions to Equivariant Elliptic Systems with Variational Structure

We consider the system Δu ? W _u (u) = 0, where \({u : \mathbb{R}^n \to \mathbb{R}^n}\) , for a class of potentials \({W : \mathbb{R}^n \to \mathbb{R}}\) that possess several global minima and are invariant under a general finite reflection group G. We establish existence of nontrivial G-equivariant entire solutions connecting the global minima of W along certain directions at infinity.     

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.

     

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.

     

Temporal Homogenization of Linear ODEs,with Applications to Parametric Super-Resonance and Energy Harvest

We consider the temporal homogenization of linear ODEs of the form \({\dot{x}=Ax+\epsilon P(t)x+f(t)}\), where P(t) is periodic and \({\epsilon}\) is small. Using a 2-scale expansion approach, we obtain the long-time approximation \({x(t)\approx {\rm exp}(At) \left( \Omega(t)+\int_0^t {\rm exp}(-A \tau) f(\tau) {\rm d}\tau \right)}\), where \({\Omega}\) solves the cell problem \({\dot{\Omega}=\epsilon B \Omega + \epsilon F(t)}\) with an effective matrix B and an explicitly-known F(t). We provide necessary and sufficient conditions for the accuracy of the approximation (over a \({{\mathcal{O}}(\epsilon^{-1})}\) time-scale), and show how B can be computed (at a cost independent of \({\epsilon}\)). As a direct application, we investigate the possibility of using RLC circuits to harvest the energy contained in small scale oscillations of ambient electromagnetic fields (such as Schumann resonances). Although a RLC circuit parametrically coupled to the field may achieve such energy extraction via parametric resonance, its resistance R needs to be smaller than a threshold \({\kappa}\) proportional to the fluctuations of the field, thereby limiting practical applications. We show that if n RLC circuits are appropriately coupled via mutual capacitances or inductances, then energy extraction can be achieved when the resistance of each circuit is smaller than \({n\kappa}\). Hence, if the resistance of each circuit has a non-zero fixed value, energy extraction can be made possible through the coupling of a sufficiently large number n of circuits (\({n\approx 1000}\) for the first mode of Schumann resonances and contemporary values of capacitances, inductances and resistances). The theory is also applied to the control of the oscillation amplitude of a (damped) oscillator.     

Stability of Traveling Waves of Nonlinear Schrödinger Equation with Nonzero Condition at Infinity

We study effective elastic behavior of the incompatibly prestrained thin plates, where the prestrain is independent of thickness and uniform through the plate’s thickness h. We model such plates as three-dimensional elastic bodies with a prescribed pointwise stress-free state characterized by a Riemannian metric G, and seek the limiting behavior as \({h \to 0}\). We first establish that when the energy per volume scales as the second power of h, the resulting \({\Gamma}\) -limit is a Kirchhoff-type bending theory. We then show the somewhat surprising result that there exist non-immersible metrics G for whom the infimum energy (per volume) scales smaller than h². This implies that the minimizing sequence of deformations carries nontrivial residual three-dimensional energy but it has zero bending energy as seen from the limit Kirchhoff theory perspective. Another implication is that other asymptotic scenarios are valid in appropriate smaller scaling regimes of energy. We characterize the metrics G with the above property, showing that the zero bending energy in the Kirchhoff limit occurs if and only if the Riemann curvatures R₁₂₁₃, R₁₂₂₃ and R₁₂₁₂ of G vanish identically. We illustrate our findings with examples; of particular interest is an example where \({G_{2 \times 2}}\), the two-dimensional restriction of G, is flat but the plate still exhibits the energy scaling of the Föppl–von Kármán type. Finally, we apply these results to a model of nematic glass, including a characterization of the condition when the metric is immersible, for \({G = Id_{3} + \gamma n \otimes n}\) given in terms of the inhomogeneous unit director field distribution \({ n \in \mathbb{R}^3}\).     

Local Fixed Point Indices of Iterations of Planar Maps

Let \({f: U\rightarrow {\mathbb R}^2}\) be a continuous map, where U is an open subset of \({{\mathbb R}^2}\). We consider a fixed point p of f which is neither a sink nor a source and such that {p} is an isolated invariant set. Under these assumption we prove, using Conley index methods and Nielsen theory, that the sequence of fixed point indices of iterations \({\{{\rm ind}(f^n,p)\}_{n=1}^\infty}\) is periodic, bounded from above by 1, and has infinitely many non-positive terms, which is a generalization of Le Calvez and Yoccoz theorem (Annals of Math., 146, 241–293 (1997)) onto the class of non-injective maps. We apply our result to study the dynamics of continuous maps on 2-dimensional sphere.