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.

     

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.     

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.

     

Visual Investigation of Improvement in Extra-Heavy Oil Recovery by Borate-Assisted $$\hbox {CO}_{2}$$-Foam Injection

The significant reduction in heavy oil viscosity when mixed with \(\hbox {CO}_{2}\) is well documented. However, for \(\hbox {CO}_{2}\) injection to be an efficient method for improving heavy oil recovery, other mechanisms are required to improve the mobility ratio between the \(\hbox {CO}_{2}\) front and the resident heavy oil. In situ generation of \(\hbox {CO}_{2}\)-foam can improve \(\hbox {CO}_{2}\) injection performance by (a) increasing the effective viscosity of \(\hbox {CO}_{2}\) in the reservoir and (b) increasing the contact area between the heavy oil and injected \(\hbox {CO}_{2}\) and hence improving \(\hbox {CO}_{2}\) dissolution rate. However, in situ generation of stable \(\hbox {CO}_{2}\)-foam capable of travelling from the injection well to the production well is hard to achieve. We have previously published the results of a series of foam stability experiments using alkali and in the presence of heavy crude oil (Farzaneh and Sohrabi 2015). The results showed that stability of \(\hbox {CO}_{2}\)-foam decreased by addition of NaOH, while it increased by addition of \(\hbox {Na}_{2}\hbox {CO}_{3}\). However, the highest increase in \(\hbox {CO}_{2}\)-foam stability was achieved by adding borate to the surfactant solution. Borate is a mild alkaline with an excellent pH buffering ability. The previous study was performed in a foam column in the absence of a porous medium. In this paper, we present the results of a new series of experiments carried out in a high-pressure glass micromodel to visually investigate the performance of borate–surfactant \(\hbox {CO}_{2}\)-foam injection in an extra-heavy crude oil in a transparent porous medium. In the first part of the paper, the pore-scale interactions of \(\hbox {CO}_{2}\)-foam and extra-heavy oil and the mechanisms of oil displacement and hence oil recovery are presented through image analysis of micromodel images. The results show that very high oil recovery was achieved by co-injection of the borate–surfactant solution with \(\hbox {CO}_{2}\), due to in-situ formation of stable foam. Dissolution of \(\hbox {CO}_{2}\) in heavy oil resulted in significant reduction in its viscosity. \(\hbox {CO}_{2}\)-foam significantly increased the contact area between the oil and \(\hbox {CO}_{2}\) significantly and thus the efficiency of the process. The synergy effect between the borate and surfactant resulted in (1) alteration of the wettability of the porous medium towards water wet and (2) significant reduction of the oil–water IFT. As a result, a bank of oil-in-water (O/W) emulsion was formed in the porous medium and moved ahead of the \(\hbox {CO}_{2}\)-foam front. The in-situ generated O/W emulsion has a much lower viscosity than the original oil and plays a major role in the observed additional oil recovery in the range of performed experiments. Borate also made \(\hbox {CO}_{2}\)-foam more stable by changing the system to non-spreading oil and reducing coalescence of the foam bubbles. The results of these visual experiments suggest that borate can be a useful additive for improving heavy oil recovery in the range of the performed tests, by increasing \(\hbox {CO}_{2}\)-foam stability and producing O/W emulsions.     

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.     

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.     

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.

     

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

     

Comparative Study of Five Outcrop Chalks Flooded at Reservoir Conditions: Chemo-mechanical Behaviour and Profiles of Compositional Alteration

This study presents experimental results from a flooding test series performed at reservoir conditions for five high-porosity Cretaceous onshore chalks from Denmark, Belgium and the USA, analogous to North Sea reservoir chalk. The chalks are studied in regard to their chemo-mechanical behaviour when performing tri-axial compaction tests while injecting brines (0.219 mol/L \(\hbox {MgCl}_{2}\) or 0.657 mol/L NaCl) at reservoir conditions for 2–3 months (T = 130 \(^\circ \hbox {C}\); 1 PV/d). Each chalk type was examined in terms of its mineralogical and chemical composition before and after the mechanical flooding tests, using an extensive set of analysis methods, to evaluate the chalk- and brine-dependent chemical alterations. All \(\hbox {MgCl}_{2}\)-flooded cores showed precipitation of Mg-bearing minerals (mainly magnesite). The distribution of newly formed Mg-bearing minerals appears to be chalk-dependent with varying peaks of enrichment. The chalk samples from Aalborg originally contained abundant opal-CT, which was dissolved with both NaCl and \(\hbox {MgCl}_{2}\) and partly re-precipitated as Si-Mg-bearing minerals. The Aalborg core injected with \(\hbox {MgCl}_{2}\) indicated strongly increased specific surface area (from 4.9 \(\hbox {m}^{2}\hbox {/g}\) to within 7–9 \(\hbox {m}^{2}\hbox {/g}\)). Mineral precipitation effects were negligible in chalk samples flooded with NaCl compared to \(\hbox {MgCl}_{2}\). Silicates were the main mineralogical impurity in the studied chalk samples (0.3–6 wt%). The cores with higher \(\hbox {SiO}_{2}\) content showed less deformation when injecting NaCl brine, but more compaction when injecting \(\hbox {MgCl}_{2}\)-brine. The observations were successfully interpreted by mathematical geochemical modelling which suggests that the re-precipitation of Si-bearing minerals leads to enhanced calcite dissolution and mass loss (as seen experimentally) explaining the high compaction seen in \(\hbox {MgCl}_{2}\)-flooded Aalborg chalk. Our work demonstrates that the original mineralogy, together with the newly formed minerals, can control the chemo-mechanical interactions during flooding and should be taken into account when predicting reservoir behaviour from laboratory studies. This study improves the understanding of complex flow reaction mechanisms also relevant for field-scale dynamics seen during brine injection.     

A semi-locally scaled eddy viscosity formulation for LES wall models and flows at high speeds

We show that the mean wall-shear stresses in wall-modeled large-eddy simulations (WMLES) of high-speed flows can be off by up to \(\approx 100\%\) with respect to a DNS benchmark when using the van-Driest-based damping function, i.e., the conventional damping function. Errors in the WMLES-predicted wall-shear stresses are often attributed to the so-called log-layer mismatch, which, albeit also an error in wall-shear stresses \(\tau _\mathrm{w}\), is an error of about \(15\%\). The larger error identified here cannot be removed using the previously developed remedies for the log-layer mismatch. This error may be removed by using the semi-local scaling, i.e., \(l_\nu =\mu /\sqrt{\rho \tau _\mathrm{w}}\), in the damping function, where \(\mu \) and \(\rho \) are the local mean dynamic viscosity and density, respectively.     

A new type of functional chemical sensitizer $$\hbox {MgH}_{2}$$ for improving pressure desensitization resistance of emulsion explosives

In millisecond-delay blasting and deep water blasting projects, traditional emulsion explosives sensitized by the chemical sensitizer \(\hbox {NaNO}_{2}\) often encounter incomplete explosion or misfire problems because of the “pressure desensitization” phenomenon, which seriously affects blasting safety and construction progress. A \(\hbox {MgH}_{2}\)-sensitized emulsion explosive was invented to solve these problems. Experimental results show that \(\hbox {MgH}_{2}\) can effectively reduce the problem of pressure desensitization. In this paper, the factors which influence the pressure desensitization of two types of emulsion explosives are studied, and resistance to this phenomenon of \(\hbox {MgH}_{2}\)-sensitized emulsion explosives is discussed.     

A KAM Theorem for Higher Dimensional Forced Nonlinear Schrödinger Equations

In this paper, we prove an infinite dimensional KAM theorem. As an application, it is shown that there are many real-analytic small-amplitude linearly-stable quasi-periodic solutions for higher dimensional nonlinear Schrödinger equations with outer force

$$\begin{aligned} iu_t-\triangle u +M_\xi u+f(\bar{\omega }t)|u|^2u=0, \quad t\in {{\mathbb R}}, x\in {{\mathbb T}}^d \end{aligned}$$

where \(M_\xi \) is a real Fourier multiplier,\(f({\bar{\theta }})({\bar{\theta }}={\bar{\omega }} t)\) is real analytic and the forced frequencies \(\bar{\omega }\) are fixed Diophantine vectors.

     

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

Validation of numerical solvers for liquid metal flow in a complex geometry in the presence of a strong magnetic field

Following the magnetohydrodynamic (MHD) code validation and verification proposal by Smolentsev et al. (Fusion Eng Des 100:65–72, 2015), we perform code to code and code to experiment comparisons between two computational solvers, FLUIDYN and HIMAG, which are presently considered as two of the prospective CFD tools for fusion blanket applications. In such applications, an electrically conducting breeder/coolant circulates in the blanket ducts in the presence of a strong plasma-confining magnetic field at high Hartmann numbers, \(\textit{Ha}\) (\(\textit{Ha}^2\) is the ratio between electromagnetic and viscous forces) and high interaction parameters, \(\textit{N}\) (\(\textit{N}\) is the ratio of electromagnetic to inertial forces). The main objective of this paper is to provide the scientific and engineering community with common references to assist fusion researchers in the selection of adequate computational means to be used for blanket design and analysis. As an initial validation case, the two codes are applied to the classic problem of a laminar fully developed MHD flows in a rectangular duct. Both codes demonstrate a very good agreement with the analytical solution for \(\textit{Ha}\) up to 15, 000. To address the capabilities of the two codes to properly resolve complex geometry flows, we consider a case of three-dimensional developing MHD flow in a geometry comprising of a series of interconnected electrically conducting rectangular ducts. The computed electric potential distributions for two flows (Case A) \(\textit{Ha}=515\), \(\textit{N}=3.2\) and (Case B) \(\textit{Ha}=2059\), \(\textit{N}=63.8\) are in very good agreement with the experimental data, while the comparisons for the MHD pressure drop are still unsatisfactory. To better interpret the observed differences, the obtained numerical data are analyzed against earlier theoretical and experimental studies for flows that involve changes in the relative orientation between the flow and the magnetic field.     

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.

     

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.     

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.     

Steady States of Fokker–Planck Equations: III. Degenerate Diffusion

This is the third paper in a series concerning the study of steady states of a Fokker–Planck equation in a general domain in \({\mathbb R}^n\) with \(L^{p}_{loc}\) drift term and \(W^{1,p}_{loc}\) diffusion term for any \(p>n\). In this paper, we give some existence results of stationary measures of the Fokker–Planck equation under Lyapunov conditions which allow the degeneracy of 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\).