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}$ .     

Hadamard Variational Formula for the Green’s Function of the Boundary Value Problem on the Stokes Equations

For every ${\varepsilon > 0}$ , we consider the Green’s matrix ${G_{\varepsilon}(x, y)}$ of the Stokes equations describing the motion of incompressible fluids in a bounded domain ${\Omega_{\varepsilon} \subset \mathbb{R}^d}$ , which is a family of perturbation of domains from ${\Omega\equiv \Omega_0}$ with the smooth boundary ${\partial\Omega}$ . Assuming the volume preserving property, that is, ${\mbox{vol.}\Omega_{\varepsilon} = \mbox{vol.}\Omega}$ for all ${\varepsilon > 0}$ , we give an explicit representation formula for ${\delta G(x, y) \equiv \lim_{\varepsilon\to +0}\varepsilon^{-1}(G_{\varepsilon}(x, y) - G_0(x, y))}$ in terms of the boundary integral on ${\partial \Omega}$ of ${G_0(x, y)}$ . Our result may be regarded as a classical Hadamard variational formula for the Green’s functions of the elliptic boundary value problems.     

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.     

Regularity and Asymptotic Behavior of Nonlinear Stefan Problems

We study the following nonlinear Stefan problem $$\left\{\begin{aligned}\!\!&u_t\,-\,d\Delta u = g(u) & &\quad{\rm for}\,x\,\in\,\Omega(t), t > 0, \\ & u = 0 \, {\rm and} u_t = \mu|\nabla_{x} u|^{2} &&\quad {\rm for}\,x\,\in\,\Gamma(t), t > 0, \\ &u(0, x) = u_{0}(x) &&\quad {\rm for}\,x\,\in\,\Omega_0,\end{aligned} \right.$$ where ${\Omega(t) \subset \mathbb{R}^{n}}$ ( ${n \geqq 2}$ ) is bounded by the free boundary ${\Gamma(t)}$ , with ${\Omega(0) = \Omega_0}$ , μ and d are given positive constants. The initial function u ₀ is positive in ${\Omega_0}$ and vanishes on ${\partial \Omega_0}$ . The class of nonlinear functions g(u) includes the standard monostable, bistable and combustion type nonlinearities. We show that the free boundary ${\Gamma(t)}$ is smooth outside the closed convex hull of ${\Omega_0}$ , and as ${t \to \infty}$ , either ${\Omega(t)}$ expands to the entire ${\mathbb{R}^n}$ , or it stays bounded. Moreover, in the former case, ${\Gamma(t)}$ converges to the unit sphere when normalized, and in the latter case, ${u \to 0}$ uniformly. When ${g(u) = au - bu^2}$ , we further prove that in the case ${\Omega(t)}$ expands to ${{\mathbb R}^n}$ , ${u \to a/b}$ as ${t \to \infty}$ , and the spreading speed of the free boundary converges to a positive constant; moreover, there exists ${\mu^* \geqq 0}$ such that ${\Omega(t)}$ expands to ${{\mathbb{R}}^n}$ exactly when ${\mu > \mu^*}$ .     

Spatiotemporal complexity of a three-species ratio-dependent food chain model

In this paper, we investigate the complex dynamics of a ratio-dependent spatially extended food chain model. Through a detailed analytical study of the reaction–diffusion model, we obtain some conditions for global stability. On the basis of bifurcation analysis, we present the evolutionary process of pattern formation near the coexistence equilibrium point $(N^*,P^*,Z^*)$ via numerical simulation. And the sequence cold spots $\rightarrow $ stripe–spots mixtures $\rightarrow $ stripes $\rightarrow $ hot stripe–spots mixtures $\rightarrow $ hot spots $\rightarrow $ chaotic wave patterns controlled by parameters $a_1$ or $c_1$ in the model are presented. These results indicate that the reaction–diffusion model is an appropriate tool for investigating fundamental mechanism of complex spatiotemporal dynamics.     

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.     

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.     

Incompressible 3D Navier–Stokes Equations as a Limit of a Nonlinear Parabolic System

In this paper, we consider the Cauchy problem for a nonlinear parabolic system ${u^\epsilon_t - \Delta u^\epsilon + u^\epsilon \cdot \nabla u^\epsilon + \frac{1}{2}u^\epsilon\, {\rm div}\, u^\epsilon - \frac{1}{\epsilon}\nabla\, {\rm div}\, u^\epsilon = 0}$ in ${\mathbb {R}^3 \times (0,\infty)}$ with initial data in Lebesgue spaces ${L^2(\mathbb {R}^3)}$ or ${L^3(\mathbb {R}^3)}$ . We analyze the convergence of its solutions to a solution of the incompressible Navier?CStokes system as ${\epsilon \to 0}$ .     

Residual Trapping Beneath Impermeable Barriers During Buoyant Migration of \mathrm{CO}_2

The presence of impermeable barriers in a reservoir can significantly impede the buoyant migration of $\mathrm{CO}_2$ injected deep into a heterogeneous geological formation. An important consequence of the presence of these impermeable barriers in terms of the long-term storage of $\mathrm{CO}_2$ is the residual trapping that takes place beneath the barriers, which acts to both increase the storage potential of the reservoir and improve the storage security of the $\mathrm{CO}_2$ . Analytical results for the total amount of $\mathrm{CO}_2$ trapped in a reservoir with an uncorrelated random distribution of impermeable barriers are obtained for both two and three-dimensional cases. In two dimensions, it is shown that the total amount of $\mathrm{CO}_2$ contained in this fashion scales as $n^{5/4}$ , where $n$ is the number of barriers in the vertical direction. In three dimensions, the trapped amount scales as $n^c$ , where $5/4 \le c \le 2$ depending on the aspect ratio of the barriers. The analytical two-dimensional results are compared with results of detailed numerical simulations, and good agreement is observed.     

Dynamical Parallelepipeds in Minimal Systems

For a topological dynamical system $(X,T)$ ( X , T ) and $d\in \mathbb N $ d ∈ N , the associated dynamical parallelepiped $\mathbf{Q}^{[d]}$ Q [ d ] was defined by Host–Kra–Maass. For a minimal distal system it was shown by them that the relation $\sim _{d-1}$ ～ d ? 1 defined on $\mathbf{Q}^{[d-1]}$ Q [ d ? 1 ] is an equivalence relation; the closing parallelepiped property holds, and for each $x\in X$ x ∈ X the collection of points in $\mathbf{Q}^{[d]}$ Q [ d ] with first coordinate $x$ x is a minimal subset under the face transformations. We give examples showing that the results do not extend to general minimal systems.     

Bounded and Unbounded Motions in Asymmetric Oscillators at Resonance

We consider the boundedness and unboundedness of solutions for the asymmetric oscillator $$\begin{aligned} x''+ax^+-bx^-+g(x)=p(t), \end{aligned}$$ where $x^+=\max \{x,0\},x^-=\max \{-x,0\}, a$ and $b$ are two positive constants, $ p(t)$ is a $2\pi $ -periodic smooth function and $g(x)$ satisfies $\lim _{|x|\rightarrow +\infty }x^{-1}g(x)=0$ . We establish some sharp sufficient conditions concerning the boundedness of all the solutions and the existence of unbounded solutions. It turns out that the boundedness of all the solutions and the existence of unbounded solutions have a close relation to the interaction of some well-defined functions $\Phi _p(\theta )$ and $\Lambda (h)$ . Some explicit conditions are given for the boundedness of all the solutions and the existence of unbounded solutions. Unlike many existing results in the literature where the function $g(x)$ is required to be a bounded function with asymptotic limits, here we allow $g(x)$ be unbounded or oscillatory without asymptotic limits.     

Γ-Convergence Analysis of Systems of Edge Dislocations: the Self Energy Regime

This paper deals with the elastic energy induced by systems of straight edge dislocations in the framework of linearized plane elasticity. The dislocations are introduced as point topological defects of the displacement-gradient fields. Following the core radius approach, we introduce a parameter ${\varepsilon > 0}$ representing the lattice spacing of the crystal, we remove a disc of radius ${\varepsilon}$ around each dislocation and compute the elastic energy stored outside the union of such discs, namely outside the core region. Then, we analyze the asymptotic behaviour of the elastic energy as ${\varepsilon \rightarrow 0}$ , in terms of Γ-convergence. We focus on the self energy regime of order ${\log\frac{1}{\varepsilon}}$ ; we show that configurations with logarithmic diverging energy converge, up to a subsequence, to a finite number of multiple dislocations and we compute the corresponding Γ-limit.     

A New Insight into the Convective Boundary Condition

The steady mixed convection boundary layer flows over a vertical surface adjacent to a Darcy porous medium and subject respectively to (i) a prescribed constant wall temperature, (ii) a prescribed variable heat flux, $q_\mathrm{w} =q_0 x^{-1/2}$ q w = q 0 x ? 1 / 2 , and (iii) a convective boundary condition are compared to each other in this article. It is shown that, in the characteristic plane spanned by the dimensionless flow velocity at the wall ${f}^{\prime }(0)\equiv \lambda $ f ′ ( 0 ) ≡ λ and the dimensionless wall shear stress $f^{\prime \prime }(0)\equiv S$ f ′ ′ ( 0 ) ≡ S , every solution $(\lambda , S)$ ( λ , S ) of one of these three flow problems at the same time is also a solution of the other two ones. There also turns out that with respect to the governing mixed convection and surface heat transfer parameters $\varepsilon $ ε and $\gamma $ γ , every solution $(\lambda , S)$ ( λ , S ) of the flow problem (iii) is infinitely degenerate. Specifically, to the very same flow solution $(\lambda , S)$ ( λ , S ) there corresponds a whole continuous set of values of $\varepsilon $ ε and $\gamma $ γ which satisfy the equation $S=-\gamma (1+\varepsilon -\lambda )$ S = ? γ ( 1 + ε ? λ ) . For the temperature solutions, however, the infinite degeneracy of the velocity solutions becomes lifted. These and further outstanding features of the convective problem (iii) are discussed in the article in some detail.     

Heat Convection of Compressible Viscous Fluids. II

We consider steady heat convections of compressible viscous fluids in the horizontal strip domain ${z_0 < z < z_0 + 1}$ under the gravity. Pattern formations are shown uniformly for ${z_0 \geq Z_0}$ . The limit of them as ${Z_0 \rightarrow + \infty}$ is that of Oberbeck-Boussinesq equations.     

Adjoint Methods for the Infinity Laplacian Partial Differential Equation

To study fine properties of certain smooth approximations ${u^\varepsilon}$ to a viscosity solution u of the infinity Laplacian partial differential equation (PDE), we introduce Green??s function ${\sigma^\varepsilon}$ for the linearization. We can then integrate by parts with respect to ${\sigma^\varepsilon}$ and derive various useful integral estimates. We are, in particular, able to use these estimates (i) to prove the everywhere differentiability of u and (ii) to rigorously justify interpreting the infinity Laplacian equation as a parabolic PDE.     

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.     

On Anisotropic Polynomial Relations for the Elasticity Tensor

In this paper, we explore new conditions for an elasticity tensor to belong to a given symmetry class. Our goal is to propose an alternative approach to the identification problem of the symmetry class, based on polynomial invariants and covariants of the elasticity tensor C, rather than on spectral properties of the Kelvin representation. We compute a set of algebraic relations which describe precisely the orthotropic ( $[\mathbb {D}_{2}]$ ), trigonal ( $[\mathbb {D}_{3}]$ ), tetragonal ( $[\mathbb {D}_{4}]$ ), transverse isotropic ([SO(2)]) and cubic ( $[\mathbb {O}]$ ) symmetry classes in $\mathbb {H}^{4}$ , the highest-order irreducible component in the decomposition of $\mathbb {E}\mathrm {la}$ . We provide a bifurcation diagram which describes how one “travels” in $\mathbb {H}^{4}$ from a given isotropy class to another. Finally, we study the link between these polynomial invariants and those obtained as the coefficients of the characteristic or the Betten polynomials. We show, in particular, that the Betten invariants do not separate the orbits of the elasticity tensors.     

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}$ .     

Mixed Convection Boundary-Layer Flow Over a Vertical Surface Embedded in a Porous Material Subject to a Convective Boundary Condition

The mixed convection boundary-layer flow on one face of a semi-infinite vertical surface embedded in a fluid-saturated porous medium is considered when the other face is taken to be in contact with a hot or cooled fluid maintaining that surface at a constant temperature $T_\mathrm{{f}}$ . The governing system of partial differential equations is transformed into a system of ordinary differential equations through an appropriate similarity transformation. These equations are solved numerically in terms of a dimensionless mixed convection parameter $\epsilon $ and a surface heat transfer parameter $\gamma $ . The results indicate that dual solutions exist for opposing flow, $\epsilon <0$ , with the dependence of the critical values $\epsilon _\mathrm{{c}}$ on $\gamma $ being determined, whereas for the assisting flow $\epsilon >0$ , the solution is unique. Limiting asymptotic forms for both $\gamma $ small and large and $\epsilon $ large are also discussed.