Further Dense Properties of the Space of Circle Diffeomorphisms with a Liouville Rotation Number

In continuation of Matsumoto’s paper (Nonlinearity 25:1495–1511, 2012) we show that various subspaces are \(C^{\infty }\)-dense in the space of orientation-preserving \(C^{\infty }\)-diffeomorphisms of the circle with rotation number \(\alpha \), where \(\alpha \in {\mathbb {S}}^1\) is any prescribed Liouville number. In particular, for every odometer \({\mathcal {O}}\) of product type we prove the denseness of the subspace of diffeomorphisms which are orbit-equivalent to \({\mathcal {O}}\).     

Resolvent Estimates for High-Contrast Elliptic Problems with Periodic Coefficients

We study the asymptotic behaviour of the resolvents \({(\mathcal{A}^\varepsilon+I)^{-1}}\) of elliptic second-order differential operators \({{\mathcal{A}}^\varepsilon}\) in \({\mathbb{R}^d}\) with periodic rapidly oscillating coefficients, as the period \({\varepsilon}\) goes to zero. The class of operators covered by our analysis includes both the “classical” case of uniformly elliptic families (where the ellipticity constant does not depend on \({\varepsilon}\)) and the “double-porosity” case of coefficients that take contrasting values of order one and of order \({\varepsilon^2}\) in different parts of the period cell. We provide a construction for the leading order term of the “operator asymptotics” of \({(\mathcal{A}^\varepsilon+I)^{-1}}\) in the sense of operator-norm convergence and prove order \({O(\varepsilon)}\) remainder estimates.     

Remarks on the Generalized Reflectionless Schrödinger Potentials

We discuss certain compact, translation-invariant subsets of the set \({\mathcal {R}}\) of the generalized reflectionless potentials for the one-dimensional Schrödinger operator. We determine a stationary ergodic subset of \({\mathcal {R}}\) whose Lyapunov exponent is discontinuous at a point. We also determine an almost automorphic, non-almost periodic minimal subset of \(\mathcal {R}\).     

The Unstable Set of a Periodic Orbit for Delayed Positive Feedback

In the paper [Large-amplitude periodic solutions for differential equations with delayed monotone positive feedback, JDDE 23 (2011), no. 4, 727–790], we have constructed large-amplitude periodic orbits for an equation with delayed monotone positive feedback. We have shown that the unstable sets of the large-amplitude periodic orbits constitute the global attractor besides spindle-like structures. In this paper we focus on a large-amplitude periodic orbit \({\mathcal {O}}_{p}\) with two Floquet multipliers outside the unit circle, and we intend to characterize the geometric structure of its unstable set \({\mathcal {W}}^{u}\left( {\mathcal {O}}_{p}\right) \). We prove that \({\mathcal {W}}^{u}\left( {\mathcal {O}}_{p}\right) \) is a three-dimensional \(C^{1}\)-submanifold of the phase space and admits a smooth global graph representation. Within \({\mathcal {W}}^{u}\left( {\mathcal {O}}_{p}\right) \), there exist heteroclinic connections from \({\mathcal {O}}_{p}\) to three different periodic orbits. These connecting sets are two-dimensional \(C^{1}\)-submanifolds of \({\mathcal {W}}^{u}\left( {\mathcal {O}}_{p}\right) \) and homeomorphic to the two-dimensional open annulus. They form \(C^{1}\)-smooth separatrices in the sense that they divide the points of \({\mathcal {W}}^{u}\left( {\mathcal {O}}_{p}\right) \) into three subsets according to their \(\omega \)-limit sets.     

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

Linearly Stable Quasi-Periodic Breathers in a Class of Random Hamiltonian Systems

In this paper, we construct linearly stable quasi-periodic breathers for the Hamiltonian systems in the form \({{\rm i} \dot{q}_n+v_n q_n+\delta|q_n|^2q_n+\varepsilon_n \left(q_{n+1}+q_{n-1} \right)=0,\quad n \in \mathbb{Z}}\) where \({\{v_n\}_{n \in \mathbb{Z}}}\) is a family of time independent identically distributed (i.i.d) random variables with common distribution \({g = dv_n, v_n \in [0,1]}\) and \({|\varepsilon_n| \leq \varepsilon e^{-\varrho |n|}}\) with \({\varepsilon,\varrho > 0}\) . We prove that for \({\varepsilon, \delta}\) sufficiently small, the equation admits a family of small-amplitude and linearly stable, time quasi-periodic solutions for most of the parameters \({\{v_n\}_{n \in \mathbb{Z}}}\) .     

Uniaxial experimental study of the acoustic emission and deformation behavior of composite rock based on 3D digital image correlation (DIC)

In this paper, uniaxial compression tests were carried out on a series of composite rock specimens with different dip angles, which were made from two types of rock-like material with different strength. The acoustic emission technique was used to monitor the acoustic signal characteristics of composite rock specimens during the entire loading process. At the same time, an optical non-contact 3 D digital image correlation technique was used to study the evolution of axial strain field and the maximal strain field before and after the peak strength at different stress levels during the loading process. The effect of bedding plane inclination on the deformation and strength during uniaxial loading was analyzed. The methods of solving the elastic constants of hard and weak rock were described. The damage evolution process, deformation and failure mechanism, and failure mode during uniaxial loading were fully determined. The experimental results show that the θ = 0?–45?specimens had obvious plastic deformation during loading, and the brittleness of the θ = 60?–90?specimens gradually increased during the loading process. When the anisotropic angle θincreased from 0?to 90?, the peak strength, peak strain,and apparent elastic modulus all decreased initially and then increased. The failure mode of the composite rock specimen during uniaxial loading can be divided into three categories:tensile fracture across the discontinuities(θ = 0?–30?), slid-ing failure along the discontinuities(θ = 45?–75?), and tensile-split along the discontinuities(θ = 90?). The axial strain of the weak and hard rock layers in the composite rock specimen during the loading process was significantly different from that of the θ = 0?–45?specimens and was almost the same as that of the θ = 60?–90?specimens. As for the strain localization highlighted in the maximum principal strain field, the θ = 0?–30?specimens appeared in the rock matrix approximately parallel to the loading direction,while in the θ = 45?–90?specimens it appeared at the hard and weak rock layer interface.     

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

     

A nilpotent algebra approach to Lagrangian mechanics and constrained motion

Lagrangian mechanics is extended to the so-called nilpotent Taylor algebra \({\mathbb {T}}\). It is shown that this extension yields a practical computational technique for the evaluation and analysis of the equations of motion of general constrained dynamical systems. The underlying \({\mathbb {T}}\)-algebra utilized herein permits the analysis of constrained dynamical systems without the need for analytical or symbolic differentiations. Instead, the algebra produces the necessary exact derivatives inherently through binary operations, thus permitting the numerical analysis of constrained dynamical systems using only the defining scalar functions (the Lagrangian \({\mathcal {L}}\) and the imposed constraints). The extension of the Lagrangian framework to the \({\mathbb {T}}\)-algebra is demonstrated analytically for a problem of constrained motion in a central field and numerically for the calculation of Lyapunov exponents of N-pendulum systems.     

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.

     

Asymmetric Domain Walls of Small Angle in Soft Ferromagnetic Films

We focus on a special type of domain wall appearing in the Landau–Lifshitz theory for soft ferromagnetic films. These domain walls are divergence-free \({\mathbb{S}^2}\)-valued transition layers that connect two directions \({m_\theta^\pm \in \mathbb{S}^2}\) (differing by an angle \({2\theta}\)) and minimize the Dirichlet energy. Our main result is the rigorous derivation of the asymptotic structure and energy of such “asymmetric” domain walls in the limit \({\theta \downarrow 0}\). As an application, we deduce that a supercritical bifurcation causes the transition from symmetric to asymmetric walls in the full micromagnetic model.     

On the Global Regularity of the Two-Dimensional Density Patch for Inhomogeneous Incompressible Viscous Flow

Regarding P.-L. Lions’ open question in Oxford Lecture Series in Mathematics and its Applications, Vol. 3 (1996) concerning the propagation of regularity for the density patch, we establish the global existence of solutions to the two-dimensional inhomogeneous incompressible Navier–Stokes system with initial density given by \({(1 - \eta){\bf 1}_{{\Omega}_{0}} + {\bf 1}_{{\Omega}_{0}^{c}}}\) for some small enough constant \({\eta}\) and some \({W^{k+2,p}}\) domain \({\Omega_{0}}\), with initial vorticity belonging to \({L^{1} \cap L^{p}}\) and with appropriate tangential regularities. Furthermore, we prove that the regularity of the domain \({\Omega_0}\) is preserved by time evolution.     

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.

     

Displacement and Dissolution Characteristics of CO$$_{}$$/Brine System in Unconsolidated Porous Media

This study investigated the dynamic displacement and dissolution of \(\hbox {CO}_{2}\) in porous media at 313 K and 6/8 MPa. Gaseous (\(\hbox {gCO}_{2}\)) at 6 MPa and supercritical \(\hbox {CO}_{2 }(\hbox {scCO}_{2}) \) at 8 MPa were injected downward into a glass bead pack at different flow rates, following upwards brine injection. The processes occurring during \(\hbox {CO}_{2}\) drainage and brine imbibition were visualized using magnetic resonance imaging. The drainage flow fronts were strongly influenced by the flow rates, resulting in different gas distributions. However, brine imbibition proceeded as a vertical compacted front due to the strong effect of gravity. Additionally, the effects of flow rate on distribution and saturation were analyzed. Then, the front movement of \(\hbox {CO}_{2}\) dissolution was visualized along different paths after imbibition. The determined \(\hbox {CO}_{2}\) concentrations implied that little \(\hbox {scCO}_{2}\) dissolved in brine after imbibition. The dissolution rate was from \(10^{-8}\) to \(10^{-9}\, \hbox {kg}\, \hbox {m}^{-3} \, \hbox {s}^{-1}\) and from \(10^{-6}\) to \(10^{-8}\, \hbox {kg}\, \hbox {m}^{-3} \, \hbox {s}^{-1}\) for \(\hbox {gCO}_{2}\) at 6 MPa and \(\hbox {scCO}_{2 }\) at 8 MPa, respectively. The total time for the \(\hbox {scCO}_{2}\) dissolution was short, indicating fast mass transfer between the \(\hbox {CO}_{2}\) and brine. Injection of \(\hbox {CO}_{2}\) under supercritical conditions resulted in a quick establishment of a steady state with high storage safety.     

A Spectral Criterion for Stability of a Steady Viscous Incompressible Flow Past an Obstacle

We show that the question of stability of a steady incompressible Navier-Stokes flow \({\mathrm{V}}\) in a 3D exterior domain \({\Omega}\) is essentially a finite-dimensional problem (Theorem 3.2). Although the associated linearized operator has an essential spectrum touching the imaginary axis, we show that certain assumptions on the eigenvalues of this operator guarantee the stability of flow \({\mathrm{V}}\) (Theorem 4.1). No assumption on the smallness of the steady flow \({\mathrm{V}}\) is required.     

Bifurcations of a creeping air–water flow in a conical container

This numerical study describes the eddy emergence and transformations in a slow steady axisymmetric air–water flow, driven by a rotating top disk in a vertical conical container. As water height \(H_{\mathrm{w}}\) and cone half-angle \(\beta \) vary, numerous flow metamorphoses occur. They are investigated for \(\beta =30^{\circ }, 45^{\circ }\), and \(60^{\circ }\). For small \(H_{\mathrm{w}}\), the air flow is multi-cellular with clockwise meridional circulation near the disk. The air flow becomes one cellular as \(H_{\mathrm{w}}\) exceeds a threshold depending on \(\beta \). For all \(\beta \), the water flow has an unbounded number of eddies whose size and strength diminish as the cone apex is approached. As the water level becomes close to the disk, the outmost water eddy with clockwise meridional circulation expands, reaches the interface, and induces a thin layer with anticlockwise circulation in the air. Then this layer expands and occupies the entire air domain. The physical reasons for the flow transformations are provided. The results are of fundamental interest and can be relevant for aerial bioreactors.     

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.