Models of Elastic Shells in Contact with a Rigid Foundation: An Asymptotic Approach

We consider a family of linearly elastic shells with thickness \(2\varepsilon\) (where \(\varepsilon\) is a small parameter). The shells are clamped along a portion of their lateral face, all having the same middle surface \(S\), and may enter in contact with a rigid foundation along the bottom face.We are interested in studying the limit behavior of both the three-dimensional problems, given in curvilinear coordinates, and their solutions (displacements \(\boldsymbol{u}^{\varepsilon}\) of covariant components \(u_{i}^{\varepsilon}\)) when \(\varepsilon\) tends to zero. To do that, we use asymptotic analysis methods. On one hand, we find that if the applied body force density is \(O(1)\) with respect to \(\varepsilon\) and surface tractions density is \(O(\varepsilon)\), a suitable approximation of the variational formulation of the contact problem is a two-dimensional variational inequality which can be identified as the variational formulation of the obstacle problem for an elastic membrane. On the other hand, if the applied body force density is \(O(\varepsilon^{2})\) and surface tractions density is \(O(\varepsilon^{3})\), the corresponding approximation is a different two-dimensional inequality which can be identified as the variational formulation of the obstacle problem for an elastic flexural shell. We finally discuss the existence and uniqueness of solution for the limit two-dimensional variational problems found.     

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.     

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.     

Statistical Analysis of Scalar Dissipation Rate Transport in Turbulent Partially Premixed Flames: A Direct Numerical Simulation Study

Statistically planar turbulent premixed and partially premixed flames for different initial turbulence intensities are simulated for global equivalence ratios ??>?=?0.7 and ??>?=?1.0 using three-dimensional Direct Numerical Simulations (DNS) with simplified chemistry. For the simulations of partially premixed flames, a random distribution of equivalence ratio following a bimodal distribution of equivalence ratio is introduced in the unburned reactants ahead of the flame. The simulation parameters in all of the cases were chosen such that the combustion situation belongs to the thin reaction zones regime. The DNS data has been used to analyse the behaviour of the dissipation rate transports of both active and passive scalars (i.e. the fuel mass fraction Y _F and the mixture fraction ξ) in the context of Reynolds Averaged Navier–Stokes (RANS) simulations. The behaviours of the unclosed terms of the Favre averaged scalar dissipation rates of fuel mass fraction and mixture fraction (i.e. \(\widetilde {\varepsilon }_Y =\overline {\rho D\nabla Y_F^{\prime \prime } \cdot \nabla Y_F^{\prime \prime } } /\overline{\rho }\) and \(\widetilde {\varepsilon }_\xi =\overline {\rho D\nabla \xi ^{\prime \prime }\cdot \nabla \xi ^{\prime \prime }} /\overline {\rho })\) transport equations have been analysed in detail. In the case of the \(\widetilde {\varepsilon }_Y \) transport, it has been observed that the turbulent transport term of scalar dissipation rate remains small throughout the flame brush whereas the terms due to density variation, scalar–turbulence interaction, reaction rate and molecular dissipation remain the leading order contributors. The term arising due to density variation remains positive throughout the flame brush and the combined contribution of the reaction and molecular dissipation to the \(\widetilde {\varepsilon }_Y \) transport remains negative throughout the flame brush in all cases. However, the behaviour of scalar–turbulence interaction term of the \(\widetilde {\varepsilon }_Y \) transport equation is significantly affected by the relative strengths of turbulent straining and the straining due to chemical heat release. In the case of the \(\widetilde {\varepsilon }_\xi \) transport, the turbulent transport term remains small throughout the flame brush and the density variation term is found to be negligible in all cases, whilst the reaction rate term is exactly zero. The scalar–turbulence interaction term and molecular dissipation term remain the leading order contributors to the \(\widetilde {\varepsilon }_\xi \) transport throughout the flame brush in all cases that have been analysed in the present study. Performances of existing models for the unclosed terms of the transport equations of \(\widetilde {\varepsilon }_Y \) and \(\widetilde {\varepsilon }_\xi \) are assessed with respect to the corresponding quantities obtained from DNS data. Based on this exercise either suitable models have been identified or new models have been proposed for the accurate closure of the unclosed terms of both \(\widetilde {\varepsilon }_Y \) and \(\widetilde {\varepsilon }_\xi \) transport equations in the context of Reynolds Averaged Navier–Stokes (RANS) simulations.     

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.     

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.

     

The Mixed Convection Boundary Layer on a Vertical Melting Front in a Non-Darcian Porous Medium

The effect of surface melting on the dual solutions that can arise in the problem of the mixed convection boundary-layer flow past a vertical surface embedded in a non-Darcian porous medium is considered. The problem is described by M, melting parameter, \(\lambda \), mixed convection parameter, and \(\gamma \), the flow inertia coefficient, numerical results being obtained in terms of these three parameters. It is seen that the melting phenomenon reduces the heat transfer rate and enhances the boundary-layer separation at the solid–liquid interface. Asymptotic solutions for the forced convection, \(\lambda =0\), and free convection, large \(\lambda \), limits are derived.     

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

     

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.     

Spreading in a Shifting Environment Modeled by the Diffusive Logistic Equation with a Free Boundary

We investigate the influence of a shifting environment on the spreading of an invasive species through a model given by the diffusive logistic equation with a free boundary. When the environment is homogeneous and favourable, this model was first studied in Du and Lin (SIAM J Math Anal 42:377–405, 2010), where a spreading–vanishing dichotomy was established for the long-time dynamics of the species, and when spreading happens, it was shown that the species invades the new territory at some uniquely determined asymptotic speed \(c_0>0\). Here we consider the situation that part of such an environment becomes unfavourable, and the unfavourable range of the environment moves into the favourable part with speed \(c>0\). We prove that when \(c\ge c_0\), the species always dies out in the long-run, but when \(0<c<c_0\), the long-time behavior of the species is determined by a trichotomy described by (a) vanishing, (b) borderline spreading, or (c) spreading. If the initial population is written in the form \(u_0(x)=\sigma \phi (x)\) with \(\phi \) fixed and \(\sigma >0\) a parameter, then there exists \(\sigma _0>0\) such that vanishing happens when \(\sigma \in (0,\sigma _0)\), borderline spreading happens when \(\sigma =\sigma _0\), and spreading happens when \(\sigma >\sigma _0\).     

An investigation into the effect of strain rate on forming limit diagram using ductile fracture criteria

In this study, forming limit diagram (FLD) is experimentally acquired for aluminum alloy 3105 in usual velocities (Quasi-static condition). In addition, numerical simulation by commercially available finite element code ABAQUS/Explicit using ductile fracture criteria is performed. Simulation is done in quasi-static condition (\(\dot{\varepsilon} \le 0.01/s\)) and case of forming by low-impact (\(\dot{\varepsilon} \le 50/s\)).The results show that a substantial improvement in high-strain-rate formability of the aluminum sheet can be obtained.     

A practical synchronization approach for fractional-order chaotic systems

A practical synchronization approach is proposed for a class of fractional-order chaotic systems to realize perfect \(\delta \)-synchronization, and the nonlinear functions in the fractional-order chaotic systems are all polynomials. The \(\delta \)-synchronization scheme in this paper means that the origin in synchronization error system is stable. The reliability of \(\delta \)-synchronization has been confirmed on a class of fractional-order chaotic systems with detailed theoretical proof and discussion. Furthermore, the \(\delta \)-synchronization scheme for the fractional-order Lorenz chaotic system and the fractional-order Chua circuit is presented to demonstrate the effectiveness of the proposed method.     

Biorthogonal stretching and shearing of an impermeable surface in a uniformly rotating fluid system

The flow induced by an impermeable flat surface executing orthogonal stretching and orthogonal shearing in a rotating fluid system is investigated. Both the stretching and shearing are linear in the coordinates. An exact similarity reduction of the Navier–Stokes equations gives rise to a pair of nonlinearly-coupled ordinary differential equations governed by three parameters. In this study we set one parameter and analyze the problem which leads to flow for an impermeable surface with shearing and stretching due to velocity u along the x-axis of equal strength a while the shearing and stretching due to velocity v along the y-axis of equal strength b. These solutions depend on two parameters—a Coriolis (rotation) parameter \(\sigma = \Omega /a\) and a stretching/shearing ratio \(\lambda =b/a\). A symmetry in solutions is found for \(\lambda = 1\). The exact solution for \(\sigma = 0\) and the asymptotic behavior of solutions for \(|\sigma | \rightarrow \infty\) are determined and compared with numerical results. Oscillatory solutions are found whose strength increases with increasing values of \(|\sigma |\). It is shown that these solutions tend to the well-known Ekman solution as \(|\sigma | \rightarrow \infty\).     

Chirality Transitions in Frustrated Ferromagnetic Spin Chains: A Link with the Gradient Theory of Phase Transitions

We study chirality transitions in frustrated ferromagnetic spin chains, in view of a possible connection with the theory of Liquid Crystals. A variational approach to the study of these systems has been recently proposed by Cicalese and Solombrino, focusing close to the helimagnet/ferromagnet transition point corresponding to the critical value of the frustration parameter \(\alpha=4\). We reformulate this problem for any \(\alpha\geq0\) in the framework of surface energies in nonconvex discrete systems with nearest neighbours ferromagnetic and next-to-nearest neighbours antiferromagnetic interactions and we link it to the gradient theory of phase transitions, by showing a uniform equivalence by \(\varGamma \)-convergence on \([0,4]\) with Modica-Mortola type functionals.     

Exact solution and invariant for fractional Cattaneo anomalous diffusion of cells in two-dimensional comb framework

This paper presents an investigation on anomalous diffusion of cells in a two-dimensional comb framework with effects of fractional Cattaneo flux. Formulated governing equation is an evolution equation with the coexisting characteristics of parabolic (diffusion) and hyperbolic (wave) for \(\alpha \) in (0, 1). Exact solution is obtained by the special fractional integral transformations, and a novel invariant is established, i.e., \(\left\langle {x^{2}\left( t \right) } \right\rangle \cdot \left\langle P \right\rangle = 0.5\) (the mean square displacement multiplied by the total number of cells along the x-axis = 0.5). Moreover, the characteristics of cells distribution, the total number and the mean square displacement of cells along the x-axis with different involved parameters, especially with the fractional parameter evolution, are shown graphically and analyzed in detail. For the cells distribution versus x, it turns from parabolic and hyperbolic with the decrease in t or the increase in \(\alpha \) or \(\xi \). It is monotonically decreasing for the cells distribution versus \(\alpha \) with different x, t and \(\xi \). For the distribution versus t with different \(\alpha \) and \(\xi \) or versus \(\alpha \) with different t, it is monotonically decreasing for the distribution of total number while monotonically increasing for the distribution of mean square displacement. It is remarkable that the anomalous subdiffusion happens along the x-axis for arbitrary parameters which is different from the classical Cattaneo diffusion.     

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

Injectivity of the Cauchy-stress tensor along rank-one connected lines under strict rank-one convexity condition

In this note, we show that the Cauchy stress tensor \(\sigma\) in nonlinear elasticity is injective along rank-one connected lines provided that the constitutive law is strictly rank-one convex. This means that \(\sigma(F+\xi\otimes\eta)=\sigma(F)\) implies \(\xi \otimes\eta=0\) under strict rank-one convexity. As a consequence of this seemingly unnoticed observation, it follows that rank-one convexity and a homogeneous Cauchy stress imply that the left Cauchy-Green strain is homogeneous, as is shown in Mihai and Neff (Int. J. Non-Linear Mech., 2016, to appear).     

New Derived from Anosov Diffeomorphisms with Pathological Center Foliation

In this paper we focused our study on derived from Anosov diffeomorphisms (DA diffeomorphisms ) of the torus \(\mathbb {T}^3,\) it is, an absolute partially hyperbolic diffeomorphism on \(\mathbb {T}^3\) homotopic to a linear Anosov automorphism of the \(\mathbb {T}^3.\) We can prove that if \(f: \mathbb {T}^3 \rightarrow \mathbb {T}^3 \) is a volume preserving DA diffeomorphism homotopic to a linear Anosov A,  such that the center Lyapunov exponent satisfies \(\lambda ^c_f(x) > \lambda ^c_A > 0,\) with x belongs to a positive volume set, then the center foliation of f is non absolutely continuous. We construct a new open class U of non Anosov and volume preserving DA diffeomorphisms, satisfying the property \(\lambda ^c_f(x) > \lambda ^c_A > 0\) for \(m-\)almost everywhere \(x \in \mathbb {T}^3.\) Particularly for every \(f \in U,\) the center foliation of f is non absolutely continuous.     

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