Delay-dependent stability analysis and \mathcal {H}_{\infty } control for LPV systems with parameter-varying state delays

This paper develops the stability analysis and delay-dependent \(\mathcal {H}_{\infty }\) control synthesis for linear parameter-varying (LPV) systems with time-varying state delays. On the basis of the Finsler’s lemma, sufficient conditions on \(\mathcal {H}_{\infty }\) performance analysis are formulated in terms of parameterized linear matrix inequalities. The interesting annihilator matrix is constituted by time-varying parameters of LPV systems to reduce the conservatism. A numerical example is presented to confirm the efficiency of the proposed method.     

Output feedback \varvec{\mathcal {H}}_{{\varvec{\infty }}} control of stochastic nonlinear time-delay systems with state and disturbance-dependent noises

This paper is concerned with the output feedback \(\mathcal {H}_\infty \) control problem for a class of stochastic nonlinear systems with time-varying state delays; the system dynamics is governed by the stochastic time-delay It \(\hat{o}\) -type differential equation with state and disturbance contaminated by white noises. The design of the output feedback \(\mathcal {H}_\infty \) control is based on the stochastic dissipative theory. By establishing the stochastic dissipation of the closed-loop system, the delay-dependent and delay-independent approaches are proposed for designing the output feedback \(\mathcal {H}_\infty \) controller. It is shown that the output feedback \(\mathcal {H}_\infty \) control problem for the stochastic nonlinear time-delay systems can be solved by two delay-involved Hamilton–Jacobi inequalities. A numerical example is provided to illustrate the effectiveness of the proposed methods.     

Multistage γ-level \mathcal{H}_{\infty} control for input-saturated systems with disturbances

For input-saturated systems with disturbances, states in the domain of attraction cannot converge to the origin, but only to neighborhood around it. In order to design the smallest possible target invariant set and the largest possible domain of attraction, in this paper, we introduce a multistage γ-level $\mathcal{H}_{\infty}$ control for achieving a smaller target invariant set within a given $\mathcal{H}_{\infty}$ performance level and a larger domain of attraction than results obtained in previous studies. In particular, for the case in which the disturbances satisfy a matched condition, this paper introduces an $\mathcal{H}_{\infty}$ control with an extra control part to perfectly reject these disturbances despite the uncertainties; the introduction of the $\mathcal{H}_{\infty}$ control with an extra control part causes the target invariant set to shrink to the origin and the $\mathcal{H}_{\infty}$ performance level to become zero.     

Miscible and Immiscible Foam Injection for Mobility Control and EOR in Fractured Oil-Wet Carbonate Rocks

Foam injection is a proven enhanced oil recovery (EOR) technique for heterogeneous reservoirs, but is less studied for EOR in fractured systems. We experimentally investigated tertiary \(\text {CO}_{2}\) injections, and \(\text {N}_{2}\) - and \(\text {CO}_{2}\) -foam injections for enhanced oil recovery in fractured, oil-wet limestone core plugs. Miscible \(\text {CO}_{2}\) and \(\text {CO}_{2}\) -foam was compared with immiscible \(\text {CO}_{2}\) - and \(\text {N}_{2}\) -foam as tertiary recovery techniques, subsequent to waterfloods, in fractured rocks with different wettability preferences. At water-wet conditions waterfloods produced approximately 40 % OOIP, by spontaneous imbibition. Waterflood oil recovery at oil-wet conditions was below 20 % OOIP, due to suppressed imbibition where water predominantly flowed through the fractures, unable to mobilize the oil trapped in the matrix. Tertiary, supercritical \(\text {CO}_{2}\) -mobilized oil trapped in the matrix, particularly at weakly oil-wet conditions, by diffusion. Recovery by diffusion was high due to small core samples, high initial oil saturation and a continuous oil phase at oil-wet conditions. Both immiscible \(\text {CO}_{2}\) - and \(\text {N}_{2}\) -foams and miscible, supercritical \(\text {CO}_{2}\) -foam demonstrated high ultimate oil recoveries, but immiscible foam was less efficient (30 pore volumes injected) compared to miscible foam (2 pore volumes injected) to reach ultimate recovery. This is explained by the capillary threshold pressure preventing the injected \(\text {N}_{2}\) gas from entering the matrix, verified by computed X-ray tomography, and the mobilized oil was displaced by the aqueous surfactant in the foam. At miscible conditions, there exists no capillary entry pressure between the oil-saturated matrix and the injected \(\text {CO}_{2}\) , allowing foam to invade the matrix for efficient oil recovery.     

Towards the Development of an Evolution Equation for Flame Turbulence Interaction in Premixed Turbulent Combustion

Flame turbulence interaction is one of the leading order terms in the scalar dissipation \(\left (\widetilde {\varepsilon }_{c}\right )\) transport equation [35] and is thus an important phenomenon in premixed turbulent combustion. Swaminathan and Grout [36] and Chakraborty and Swaminathan [15, 16] have shown that the effect of strain rate on the transport of \(\widetilde {\varepsilon }_{c}\) is dominated by the interaction between the fluctuating scalar gradients and the fluctuating strain rate, denoted here by \(\overline {\rho }\widetilde {\Delta }_{c}= \overline {\rho {\alpha }\nabla c^{\prime \prime }S_{ij}^{\prime \prime }\nabla c^{\prime \prime }}\) ; this represents the flame turbulence interaction. In order to obtain an accurate representation of this phenomenon, a new evolution equation for \(\widetilde {\Delta }_{c}\) has been proposed. This equation gives a detailed insight into flame turbulence interaction and provides an alternative approach to model the important physics represented by \(\widetilde {\Delta }_{c}\) . The \(\widetilde {\Delta }_{c}\) evolution equation is derived in detail and an order of magnitude analysis is carried out to determine the leading order terms in the \(\widetilde {\Delta }_{c}\) evolution equation. The leading order terms are then studied using a Direct Numerical Simulation (DNS) of premixed turbulent flames in the corrugated flamelet regime. It is found that the behaviour of \(\widetilde {\Delta }_{c}\) is determined by the competition between the source terms (pressure gradient and the reaction rate), diffusion/dissipation processes, turbulent strain rate and the dilatation rate. Closures for the leading order terms in \(\widetilde {\Delta }_{c}\) evolution equation have been proposed and compared with the DNS data.     

Flow Dynamics of \hbox {CO}_{2}/brine at the Interface Between the Storage Formation and Sealing Units in a Multi-layered Model

Pressure distribution and \(\hbox {CO}_{2}\) plume migration are two major interests in \(\hbox {CO}_{2}\) geologic storage as they determine the injectivity and storage capacity. In this study, we adopted a three-layer model comprising a storage formation and the over- and underlying seals and determined three distinct flow regions based on the vertical flux exchange of \(\hbox {CO}_{2}\) and native brine. Regions 1 and 2 showed \(\hbox {CO}_{2}\) flowing from the storage formation to adjacent seals with counter-flowing brine. The characteristics of these fluxes in Region 1 were governed by permeability change due to salt precipitation whereas buoyancy force controlled the flux pattern in Region 2. Region 3 showed brine flowing from storage formation toward the over- and underlying seals, which enabled the displaced brine to escape from the storage formation and make room for \(\hbox {CO}_{2}\) to store as well as reduce the pressure build-up. In the multi-layered model, the counter-flowing brine in flow Region 1 resulted in localized salt precipitation at the upper and lower boundary of storage formation. We assessed the bottom-hole pressure and \(\hbox {CO}_{2}\) mass in caprock with respect to reservoir size. While the formation thickness influenced the bottom-hole pressure in the early stage of injection, the horizontal extension of the reservoir was more influential to pressure build-up during the injection period, and to the stabilized pressure during the post-injection period. The \(\hbox {CO}_{2}\) mass in caprock gently increased during the injection period as well as during the post-injection period and reached about 4–5 % of injected \(\hbox {CO}_{2}\) . The percentage of escaped brine from the storage formation ranged from 80–100 % of the \(\hbox {CO}_{2}\) mass stored in the storage formation depending on the reservoir scale.     

Some Aspects of Stability for Semigroup Actions and Control Systems

The present paper introduces both the notions of Lagrange and Poisson stabilities for semigroup actions. Let \(S\) be a semigroup acting on a topological space \(X\) with mapping \(\sigma :S\times X\rightarrow X\) , and let \(\mathcal {F}\) be a family of subsets of \(S\) . For \(x\in X\) the motion \(\sigma _{x}:S\rightarrow X\) is said to be forward Lagrange stable if the orbit \(Sx\) has compact closure in \(X\) . The point \(x\) is forward \(\mathcal {F}\) -Poisson stable if and only if it belongs to the limit set \(\omega \left( x,\mathcal {F}\right) \) . The concept of prolongational limit set is also introduced and used to describe nonwandering points. It is shown that a point \(x\) is \( \mathcal {F}\) -nonwandering if and only if \(x\) lies in its forward \(\mathcal {F} \) -prolongational limit set \(J\left( x,\mathcal {F}\right) \) . The paper contains applications to control systems.     

The Effects of Brine Species on the Formation of Residual Water in a \hbox {CO}_{2}–Brine System

Geological storage of \(\hbox {CO}_{2}\) in deep saline aquifers is achieved by injecting \(\hbox {CO}_{2}\) into the aquifers and displacing the brine. Although most of the brine is displaced, some residual groundwater remains in the rock pores. We conducted experiments to investigate factors that influence how much of this residual water remains after \(\hbox {CO}_{2}\) is injected. A rock sample was saturated with brines of two different salts. Supercritical \(\hbox {CO}_{2}\) was injected into the samples at aquifer temperature and pressure, and the displaced water and water–gas mixtures were collected and measured. The results show that deionized water drains more completely than either of the two brines, and NaCl brine drains more completely than \(\hbox {CaCl}_{2}\) brine. The ranking of the irreducible water saturation at the end of the experiment is deionized \(\hbox {water}<\hbox {NaCl brine } <\hbox {CaCl}_{2}\) brine. The process of drainage can be divided into three stages according to the drainage flow rates; the Pushing Drainage, Portable Drainage, and Dissolved Drainage stages. This paper proposed a capillary model which is used to interpret the mechanisms that characterize these three stages.     

Stochastic Variational Inequalities and Applications to the Total Variation Flow Perturbed by Linear Multiplicative Noise

In this work, we introduce a new method to prove the existence and uniqueness of a variational solution to the stochastic nonlinear diffusion equation ${{\rm d}X(t) = {\rm div} \left[\frac{\nabla X(t)}{|\nabla X(t)|}\right]{\rm d}t + X(t){\rm d}W(t) {\rm in} (0, \infty) \times \mathcal{O},}$ where ${\mathcal{O}}$ is a bounded and open domain in ${\mathbb{R}^N, N \geqq 1}$ and W(t) is a Wiener process of the form ${W(t) = \sum^{\infty}_{k = 1}\mu_{k}e_{k}\beta_{k}(t), e_{k} \in C^{2}(\overline{\mathcal{O}}) \cap H^{1}_{0}(\mathcal{O}),}$ and ${\beta_{k}, k \in \mathbb{N}}$ are independent Brownian motions. This is a stochastic diffusion equation with a highly singular diffusivity term. One main result established here is that for all initial conditions in ${L^2(\mathcal{O})}$ , it is well posed in a class of continuous solutions to the corresponding stochastic variational inequality. Thus, one obtains a stochastic version of the (minimal) total variation flow. The new approach developed here also allows us to prove the finite time extinction of solutions in dimensions ${1\leqq N \leqq3}$ , which is another main result of this work.     

On the Characterization of the Navier–Stokes Flows with the Power-Like Energy Decay

We study the energy decay of the turbulent solutions to the Navier–Stokes equations in the whole three-dimensional space. We show as the main result that the solutions with the energy decreasing at the rate \({O(t^{-\alpha}), t \rightarrow \infty, \alpha \in [0, 5/2]}\) , are exactly characterized by their initial conditions belonging into the homogeneous Besov space \({\dot{B}^{-\alpha}_{2, \infty}}\) . Similarly, for a solution u and \({p \in [1, \infty]}\) the integral \({\int_{0}^{\infty} \|t^{\alpha/2} u(t)\|^p \frac{1}{t} dt}\) is finite if and only if the initial condition of u belongs to the homogeneous Besov space \({\dot{B}_{2, p}^{-\alpha}}\) . For the case \({\alpha \in (5/2, 9/2]}\) we present analogical results for some subclasses of turbulent solutions.     

Global Compensated Compactness Theorem for General Differential Operators of First Order

Let A ₁(x, D) and A ₂(x, D) be differential operators of the first order acting on l-vector functions ${u= (u_1, \ldots, u_l)}$ in a bounded domain ${\Omega \subset \mathbb{R}^{n}}$ with the smooth boundary ${\partial\Omega}$ . We assume that the H ¹-norm ${\|u\|_{H^{1}(\Omega)}}$ is equivalent to ${\sum_{i=1}^2\|A_iu\|_{L^2(\Omega)} + \|B_1u\|_{H^{\frac{1}{2}}(\partial\Omega)}}$ and ${\sum_{i=1}^2\|A_iu\|_{L^2(\Omega)} + \|B_2u\|_{H^{\frac{1}{2}}(\partial\Omega)}}$ , where B _i  = B _i (x, ν) is the trace operator onto ${\partial\Omega}$ associated with A _i (x, D) for i = 1, 2 which is determined by the Stokes integral formula (ν: unit outer normal to ${\partial\Omega}$ ). Furthermore, we impose on A ₁ and A ₂ a cancellation property such as ${A_1A_2^{\prime}=0}$ and ${A_2A_1^{\prime}=0}$ , where ${A^{\prime}_i}$ is the formal adjoint differential operator of A _i (i = 1, 2). Suppose that ${\{u_m\}_{m=1}^{\infty}}$ and ${\{v_m\}_{m=1}^{\infty}}$ converge to u and v weakly in ${L^2(\Omega)}$ , respectively. Assume also that ${\{A_{1}u_m\}_{m=1}^{\infty}}$ and ${\{A_{2}v_{m}\}_{m=1}^{\infty}}$ are bounded in ${L^{2}(\Omega)}$ . If either ${\{B_{1}u_m\}_{m=1}^{\infty}}$ or ${\{B_{2}v_m\}_{m=1}^{\infty}}$ is bounded in ${H^{\frac{1}{2}}(\partial\Omega)}$ , then it holds that ${\int_{\Omega}u_m\cdot v_m \,{\rm d}x \to \int_{\Omega}u\cdot v \,{\rm d}x}$ . We also discuss a corresponding result on compact Riemannian manifolds with boundary.     

Elucidating the Role of Interfacial Tension for Hydrological Properties of Two-Phase Flow in Natural Sandstone by an Improved Lattice Boltzmann Method

We investigated the interfacial tension (IFT) effect on fluid flow characteristics inside micro-scale, porous media by a highly efficient multi-phase lattice Boltzmann method using a graphics processing unit. IFT is one of the most important parameters for carbon capture and storage and enhanced oil recovery. Rock pores of Berea sandstone were reconstructed from micro-CT scanned images, and multi-phase flows were simulated for the digital rock model at extremely high resolution (3.2  \(\upmu \) m). Under different IFT conditions, numerical analyses were carried out first to investigate the variation in relative permeability, and then to clarify evolution of the saturation distribution of injected fluid. We confirmed that the relative permeability decreases with increasing IFT due to growing capillary trapping intensity. It was also observed that with certain pressure gradient \(\Delta P\) two crucial IFT values, \(\sigma _{1}\) and \(\sigma _{2}\) , exist, creating three zones in which the displacement process has totally different characteristics. When \(\sigma _{1}< \sigma < \sigma _{2}\) , the capillary fingering patterns are observed, while for \(\sigma < \sigma _{1}\) viscous fingering is dominant and most of the passable pore spaces were invaded. When \(\sigma > \sigma _{2}\) the invading fluid failed to break through. The pore-throat-size distribution estimated from these crucial IFT values ( \(\sigma _{1 }\) and \(\sigma _{2})\) agrees with that derived from mercury porosimetry measurements of Berea sandstone. This study demonstrates that the proposed numerical method is an efficient tool for investigating hydrological properties from pore structures.     

Semianalytical Solution for \text{ CO}_{2} Plume Shape and Pressure Evolution During \text{ CO}_{2} Injection in Deep Saline Formations

The injection of supercritical carbon dioxide ( $\text{ CO}_{2})$ in deep saline aquifers leads to the formation of a $\text{ CO}_{2}$ rich phase plume that tends to float over the resident brine. As pressure builds up, $\text{ CO}_{2}$ density will increase because of its high compressibility. Current analytical solutions do not account for $\text{ CO}_{2}$ compressibility and consider a volumetric injection rate that is uniformly distributed along the whole thickness of the aquifer, which is unrealistic. Furthermore, the slope of the $\text{ CO}_{2}$ pressure with respect to the logarithm of distance obtained from these solutions differs from that of numerical solutions. We develop a semianalytical solution for the $\text{ CO}_{2}$ plume geometry and fluid pressure evolution, accounting for $\text{ CO}_{2}$ compressibility and buoyancy effects in the injection well, so $\text{ CO}_{2}$ is not uniformly injected along the aquifer thickness. We formulate the problem in terms of a $\text{ CO}_{2}$ potential that facilitates solution in horizontal layers, with which we discretize the aquifer. Capillary pressure is considered at the interface between the $\text{ CO}_{2}$ rich phase and the aqueous phase. When a prescribed $\text{ CO}_{2}$ mass flow rate is injected, $\text{ CO}_{2}$ advances initially through the top portion of the aquifer. As $\text{ CO}_{2}$ is being injected, the $\text{ CO}_{2}$ plume advances not only laterally, but also vertically downwards. However, the $\text{ CO}_{2}$ plume does not necessarily occupy the whole thickness of the aquifer. We found that even in the cases in which the $\text{ CO}_{2}$ plume reaches the bottom of the aquifer, most of the injected $\text{ CO}_{2}$ enters the aquifer through the layers at the top. Both $\text{ CO}_{2}$ plume position and fluid pressure compare well with numerical simulations. This solution permits quick evaluations of the $\text{ CO}_{2}$ plume position and fluid pressure distribution when injecting supercritical $\text{ CO}_{2}$ in a deep saline aquifer.     

Influence of \text{ CF }_{3}\text{ H } and \text{ CCl }_{4} additives on acetylene detonation

The influence of $\text{ CF }_{3}\text{ H }$ and $\text{ CCl }_{4}$ admixtures (known as detonation suppressors for combustible mixtures) on the development of acetylene detonation was experimentally investigated in a shock tube. The time-resolved images of detonation wave development and propagation were registered using a high-speed streak camera. Shock wave velocity and pressure profiles were measured by five calibrated piezoelectric gauges and the formation of condensed particles was detected by laser light extinction. The induction time of detonation development was determined as the moment of a pressure rise at the end plate of the shock tube. It was shown that $\text{ CF }_{3}\text{ H }$ additive had no influence on the induction time. For $\text{ CCl }_{4}$ , a significant promoting effect was observed. A simplified kinetic model was suggested and characteristic rates of diacetylene $\text{ C }_{4}\text{ H }_{2}$ formation were estimated as the limiting stage of acetylene polymerisation. An analysis of the obtained data indicated that the promoting species is atomic chlorine formed by $\text{ CCl }_{4}$ pyrolysis, which interacts with acetylene and produces $\text{ C }_{2}\text{ H }$ radical, initiating a chain mechanism of acetylene decomposition. The results of kinetic modelling agree well with the experimental data.     

Elongational rheology of NIPAM-based hydrogels

Hydrogels of different composition based on the copolymerization of N-isopropyl acrylamide and surfmers of different chemical structure were tested in elongation using Hencky/real definitions for stress, strain, and strain rate, offering a more scientific insight into the effect of deformation on the properties. In a range between $\dot {\varepsilon }=10$ and 0.01 s $^{-1}$ , the material properties are independent of strain rate and show a very clear strain hardening with a “brittle” sudden fracture. The addition of surfmer increases the strain at break $\varepsilon _{\mathrm {H}}^{\max }$ and at the same time leads to a failure of hyperelastic models. The samples can be stretched up to Hencky strains $\varepsilon _{\mathrm {H}}^{\max }$ between 0.6 and 2.5, depending on the molecular structure, yielding linear Young’s moduli E $_{0}$ between 2,700 and 39,000 Pa. The strain-rate independence indicates an ideal rubberlike behavior and fracture in a brittle-like fashion. The resulting stress at break $\sigma _{\textrm max}$ can be correlated with $\varepsilon _{\mathrm {H}}^{\max } $ and $E_{0}$ as well as with the solid molar mass between the cross-linking points $M_{\mathrm {c}}^{\textrm {solids}} $ , derived from $E_{0}$ .     

On the Elastic Energy Density of Constrained Q-Tensor Models for Biaxial Nematics

Within the Landau–de Gennes theory, the order parameter describing a biaxial nematic liquid crystal assigns a symmetric traceless 3 × 3 matrix Q with three distinct eigenvalues to every point of the region Ω occupied by the system. In the constrained case of matrices Q with constant eigenvalues, the order parameter space is diffeomorphic to the eightfold quotient ${\mathbb{S}^3/\mathcal{H}}$ of the 3-sphere ${\mathbb{S}^3}$ , where ${\mathcal{H}}$ is the quaternion group, and a configuration of a biaxial nematic liquid crystal is described by a map from Ω to ${\mathbb{S}^3/\mathcal{H}}$ . We express the (simplest form of the) Landau–de Gennes elastic free-energy density as a density defined on maps ${q: \Omega \to \mathbb{S}^3}$ , whose functional dependence is restricted by the requirements that (1) it is well defined on the class of configuration maps from Ω to ${\mathbb{S}^3/\mathcal{H}}$ (residual symmetry) and (2) it is independent of arbitrary superposed rigid rotations (frame indifference). As an application of this representation, we then discuss some properties of the corresponding energy functional, including coercivity, lower semicontinuity and strong density of smooth maps. Other invariance properties are also considered. In the discussion, we take advantage of the identification of ${\mathbb{S}^3}$ with the Lie group of unit quaternions ${Sp(1) \cong SU(2)}$ and of the relations between quaternions and rotations in ${\mathbb{R}^3}$ and ${\mathbb{R}^4}$ .     

Generalized Resolvent Estimates of the Stokes Equations with First Order Boundary Condition in a General Domain

In this paper, we prove unique existence of solutions to the generalized resolvent problem of the Stokes operator with first order boundary condition in a general domain ${\Omega}$ of the N-dimensional Eulidean space ${\mathbb{R}^N, N \geq 2}$ . This type of problem arises in the mathematical study of the flow of a viscous incompressible one-phase fluid with free surface. Moreover, we prove uniform estimates of solutions with respect to resolvent parameter ${\lambda}$ varying in a sector ${\Sigma_{\sigma, \lambda_0} = \{\lambda \in \mathbb{C} \mid |\arg \lambda| < \pi-\sigma, \enskip |\lambda| \geq \lambda_0\}}$ , where ${0 < \sigma < \pi/2}$ and ${\lambda_0 \geq 1}$ . The essential assumption of this paper is the existence of a unique solution to a suitable weak Dirichlet problem, namely it is assumed the unique existence of solution ${p \in \hat{W}^1_{q, \Gamma}(\Omega)}$ to the variational problem: ${(\nabla p, \nabla \varphi) = (f, \nabla \varphi)}$ for any ${\varphi \in \hat W^1_{q', \Gamma}(\Omega)}$ . Here, ${1 < q < \infty, q' = q/(q-1), \hat W^1_{q, \Gamma}(\Omega)}$ is the closure of ${W^1_{q, \Gamma}(\Omega) = \{ p \in W^1_q(\Omega) \mid p|_\Gamma = 0\}}$ by the semi-norm ${\|\nabla \cdot \|_{L_q(\Omega)}}$ , and ${\Gamma}$ is the boundary of ${\Omega}$ . In fact, we show that the unique solvability of such a Dirichlet problem is necessary for the unique existence of a solution to the resolvent problem with uniform estimate with respect to resolvent parameter varying in ${(\lambda_0, \infty)}$ . Our assumption is satisfied for any ${q \in (1, \infty)}$ by the following domains: whole space, half space, layer, bounded domains, exterior domains, perturbed half space, perturbed layer, but for a general domain, we do not know any result about the unique existence of solutions to the weak Dirichlet problem except for q =  2.     

Maximal and Minimal Spreading Speeds for Reaction Diffusion Equations in Nonperiodic Slowly Varying Media

This paper investigates the asymptotic behavior of the solutions of the Fisher-KPP equation in a heterogeneous medium, $$\partial_t u = \partial_{xx} u + f(x,u),$$ associated with a compactly supported initial datum. A typical nonlinearity we consider is ${f(x,u) = \mu_0 (\phi (x)) u(1-u)}$ , where??? ₀ is a 1-periodic function and ${\phi}$ is a ${\mathcal{C}^1}$ increasing function that satisfies ${\lim_{x \to+\infty}\phi (x) = +\infty}$ and ${\lim_{x \to +\infty}\phi' (x) =0}$ . Although quite specific, the choice of such a reaction term is motivated by its highly heterogeneous nature. We exhibit two different behaviors for u for large times, depending on the speed of the convergence of ${\phi}$ at infinity. If ${\phi}$ grows sufficiently slowly, then we prove that the spreading speed of u oscillates between two distinct values. If ${\phi}$ grows rapidly, then we compute explicitly a unique and well determined speed of propagation w _??, arising from the limiting problem of an infinite period. We give a heuristic interpretation for these two behaviors.     

Liouville Type Property and Spreading Speeds of KPP Equations in Periodic Media with Localized Spatial Inhomogeneity

The current paper is devoted to the study of semilinear dispersal evolution equations of the form $$\begin{aligned} u_t(t,x)=(\mathcal {A}u)(t,x)+u(t,x)f(t,x,u(t,x)),\quad x\in \mathcal {H}, \end{aligned}$$ where $\mathcal {H}=\mathbb {R}^N$ or $\mathbb {Z}^N,\; \mathcal {A}$ is a random dispersal operator or nonlocal dispersal operator in the case $\mathcal {H}=\mathbb {R}^N$ and is a discrete dispersal operator in the case $\mathcal {H}=\mathbb {Z}^N$ , and $f$ is periodic in $t$ , asymptotically periodic in $x$ (i.e. $f(t,x,u)-f_0(t,x,u)$ converges to $0$ as $\Vert x\Vert \rightarrow \infty $ for some time and space periodic function $f_0(t,x,u)$ ), and is of KPP type in $u$ . It is proved that Liouville type property for such equations holds, that is, time periodic strictly positive solutions are unique. It is also proved that if $u\equiv 0$ is a linearly unstable solution to the time and space periodic limit equation of such an equation, then it has a unique stable time periodic strictly positive solution and has a spatial spreading speed in every direction.     

Stable disarrangement phases of elastic aggregates: a setting for the emergence of no-tension materials with non-linear response in compression

A recent field theory of elastic bodies undergoing non-smooth submacroscopic geometrical changes (disarrangements) provides a setting in which, for a given homogeneous macroscopic deformation \(F\) of the body, there are typically a number of different states \(G\) of smooth, submacroscopic deformation (disarrangement phases) available to the body. A tensorial consistency relation and the inequality \(\det G\le \det F\) that guarantees that \(F\) accommodates \(G\) determine the totality of disarrangement phases \(G\) corresponding to \(F\) , and it is natural to seek for a given \(F\) those disarrangement phases that minimize the Helmholtz free energy (stable disarrangement phases). We introduce these concepts in the particular context of continuous bodies comprised of many small elastic bodies (elastic aggregates) and in the context where disarrangements do not contribute to the Helmholtz free energy (purely dissipative disarrangements). In this setting, the Helmholtz free energy response \(G\longmapsto \varPsi (G)\) of the pieces of the aggregate determines the totality of disarrangement phases corresponding to \(F\) , which necessarily includes the phase \(G=F\) (compact phase) in which every piece of the aggregate undergoes the given macroscopic deformation \(F\) . When the response function \(\varPsi \) is isotropic and smooth, and when \(\varPsi \) possesses standard semiconvexity and growth properties, the body also admits phases of the form \(G=\zeta _{\min }R\) (loose phases) with \(R\) an arbitrary rotation, provided that \(\zeta _{\min }R \) satisfies the accommodation inequality \(\zeta _{\min }^{3}\le \det F\) . Loose phases, when available, achieve the global minimum \(\varPsi (\zeta _{\min }R)\) of the free energy and consequently are stable and stress-free. When \( \varPsi (G)\) has the specific form \(\varPsi _{\alpha \beta }(G)=(\alpha /2)(\det G)^{-2}+(\beta /2)tr(GG^{T})\) , with \(\alpha \) , \(\beta \) given elastic constants, we determine all of the disarrangement phases corresponding to \(F\) . These include not only the compact and loose phases, but also disarrangement phases \(G\) in which the stress \(D\varPsi (G)\) is uniaxial or planar. Our main result (“stability implies no-tension”) is the assertion that every stable disarrangement phase for \(\varPsi _{\alpha \beta }\) cannot support tensile tractions, and our treatment of elastic aggregates thus provides a natural setting for the emergence of no-tension materials whose response in compression is non-linear. Existing treatments of no-tension materials assume at the outset that the body cannot support tension and that the response in compression is linear.