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.

     

Real Analytic Quasi-Periodic Solutions with More Diophantine Frequencies for Perturbed KdV Equations

In this paper, we consider the perturbed KdV equation with Fourier multiplier

$$\begin{aligned} u_{t} =- u_{xxx} + \big (M_{\xi }u+u^3 \big )_{x},\quad u(t,x+2\pi )=u(t,x),\quad \int _0^{2\pi }u(t,x)dx=0, \end{aligned}$$

with analytic data of size \(\varepsilon \). We prove that the equation admits a Whitney smooth family of small amplitude, real analytic quasi-periodic solutions with \(\tilde{J}\) Diophantine frequencies, where the order of \(\tilde{J}\) is \(O(\frac{1}{\varepsilon })\). The proof is based on a conserved quantity \(\int _0^{2\pi } u^2 dx\), Töplitz–Lipschitz property and an abstract infinite dimensional KAM theorem. By taking advantage of the conserved quantity \(\int _0^{2\pi } u^2 dx\) and Töplitz–Lipschitz property, our normal form part is independent of angle variables in spite of the unbounded perturbation.

     

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.     

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.     

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.     

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.     

SIF-based fracture criterion for interface cracks

The complex stress intensity factor K governing the stress field of an interface crack tip may be split into two parts, i.e.,■ and s~(-iε), so that K = ■ s~(-iε), s is a characteristic length and ε is the oscillatory index. ■ has the same dimension as the classical stress intensity factor and characterizes the interface crack tip field. That means a criterion for interface cracks may be formulated directly with■, as Irwin(ASME J. Appl. Mech. 24:361–364, 1957) did in 1957 for the classical fracture mechanics. Then, for an interface crack,it is demonstrated that the quasi Mode I and Mode II tip fields can be defined and distinguished from the coupled mode tip fields. Built upon SIF-based fracture criteria for quasi Mode I and Mode II, the stress intensity factor(SIF)-based fracture criterion for mixed mode interface cracks is proposed and validated against existing experimental results.     

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

     

Information entropies of multi-qubit Rabi model beyond the rotating wave approximation

In this paper, we present an analytical solution for the multi-qubit Rabi model for the general case without the rotating wave approximation. We discuss three different cases, one-qubit, two-qubit and three-qubit Rabi model, as applications to the presented solution. In addition, we study the effect of the qubit-photon coupling constants and the detuning parameters on the information entropies \(H(\sigma _{X})\), \(H(\sigma _{Y})\) and \(H(\sigma _{Z}),\) by which we note that, in the case of two-qubit Rabi model, there are collapse and revivals in the curves of the information entropies \(H(\sigma _{Z})\) and \(H(\sigma _{X})\) that vary in their numbers and amplitudes depending on the detuning parameter. The effects of the coupling constants and the detuning parameters on other cases are discussed. Besides, we explain that the dynamic nonlinear properties of the system have strong proportional relationship with the number of atoms.     

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.     

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

     

Infinitely Many Stability Switches in a Problem with Sublinear Oscillatory Boundary Conditions

We consider the elliptic equation \(-\Delta u +u =0\) with nonlinear boundary condition \(\frac{\partial u}{\partial n}= \lambda u + g(\lambda ,x,u), \) where \(\frac{g(\lambda ,x,s)}{s} \rightarrow 0, \hbox { as }|s|\rightarrow \infty \) and g is oscillatory. We provide sufficient conditions on g for the existence of unbounded sequences of stable solutions, unstable solutions, and turning points, even in the absence of resonant solutions.     

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

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

     

L.A.S. and Negative Schwarzian Derivative Do Not Imply G.A.S. in Clark’s Equation

It has been conjectured that a locally attracting equilibrium of Clark’s equation \(x_{n+1}=\alpha x_n + (1-\alpha )h(x_{n-k})\) is also globally attracting whenever h is a unimodal or decreasing map with negative Schwarzian derivative. The main aim of this paper is to show that the conjecture is false when \(k\ge 3\). This is done by studying the (Neimark–Sacker) bifurcation at the parameter \(\alpha \) where the locally attracting equilibrium of the equation becomes unstable. Our results, on the other hand, reinforce the validity of the conjecture in the cases \(k=1,2\).     

Chaotic characteristics analysis of the sintering process system with unknown dynamic functions based on phase space reconstruction and chaotic invariables

We develop a local discontinuous Galerkin finite element method for the distributed-order time and Riesz space-fractional convection–diffusion and Schrödinger-type equations. The stability of the presented schemes is proved and optimal order of convergence \(\mathcal {O}(h^{N+1}+(\Delta t)^{1+\frac{\theta }{2}}+\theta ^{2})\) for the Riesz space-fractional diffusion and Schrödinger-type equations with distributed order in time, an order of convergence of \(\mathcal {O}(h^{N+\frac{1}{2}}+(\Delta t)^{1+\frac{\theta }{2}}\) \(+\theta ^{2})\) is provided for the Riesz space-fractional convection–diffusion equations with distributed order in time where h, \(\theta \) and \(\Delta t\) are space step size, the distributed-order variables and the step sizes in time, respectively. Finally, the performed numerical examples confirm the optimal convergence order and illustrate the effectiveness of the method.     

Whirl Mappings on Generalised Annuli and the Incompressible Symmetric Equilibria of the Dirichlet Energy

In this paper we show a striking contrast in the symmetries of equilibria and extremisers of the total elastic energy of a hyperelastic incompressible annulus subject to pure displacement boundary conditions. Indeed upon considering the equilibrium equations, here, the nonlinear second order elliptic system formulated for the deformation \(u=(u_{1}, \ldots, u_{N})\):

$$ {\mathbb{E}} {\mathbb{L}}[u, {\mathbf {X}}] = \left \{ \textstyle\begin{array}{l@{\quad}l} \Delta u = \operatorname{div}(\mathscr{P} (x) \operatorname{cof} \nabla u) & \textrm{in }{\mathbf {X}},\\ \det\nabla u = 1 & \textrm{in }{\mathbf {X}},\\ u \equiv\varphi& \textrm{on }\partial{\mathbf {X}}, \end{array}\displaystyle \right . $$

where \({\mathbf {X}}\) is a finite, open, symmetric \(N\)-annulus (with \(N \ge2\)), \(\mathscr{P}=\mathscr{P}(x)\) is an unknown hydrostatic pressure field and \(\varphi\) is the identity mapping, we prove that, despite the inherent rotational symmetry in the system, when \(N=3\), the problem possesses no non-trivial symmetric equilibria whereas in sharp contrast, when \(N=2\), the problem possesses an infinite family of symmetric and topologically distinct equilibria. We extend and prove the counterparts of these results in higher dimensions by way of showing that a similar dichotomy persists between all odd vs. even dimensions \(N \ge4\) and discuss a number of closely related issues.

     

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.

     

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.     

On the Damped Harmonic Oscillator with Time Dependent Damping Coefficient

Conditions guaranteeing asymptotic stability for the differential equation

$$\begin{aligned} x''+h(t)x'+\omega ^2x=0 \qquad (x\in \mathbb {R}) \end{aligned}$$

are studied, where the damping coefficient \(h:[0,\infty )\rightarrow [0,\infty )\) is a locally integrable function, and the frequency \(\omega >0\) is constant. Our conditions need neither the requirement \(h(t)\le \overline{h}<\infty \) (\(t\in [0,\infty )\); \(\overline{h}\) is constant) (“small damping”), nor \(0< \underline{h}\le h(t)\) (\(t\in [0,\infty )\); \(\underline{h}\) is constant) (“large damping”); in other words, they can be applied to the general case \(0\le h(t)<\infty \) (\(t\in [0,\infty \))). We establish a condition which combines weak integral positivity with Smith’s growth condition

$$\begin{aligned} \int ^\infty _0 \exp [-H(t)]\int _0^t \exp [H(s)]\,\mathrm{{d}}s\,\mathrm{{d}}t=\infty \qquad \left( H(t):=\int _0^t h(\tau )\,\mathrm{{d}}\tau \right) , \end{aligned}$$

so it is able to control both the small and the large values of the damping coefficient simultaneously.