An experimental—Numerical evaluation of the<Emphasis Type="Italic">T</Emphasis><Stack><Subscript>ε</Subscript><Superscript>*</Superscript></Stack> integral for a three-dimensional crack front

TheT _ε ^* integral was calculated on the surface of single edge notched, three-point bend (SE(B)) specimens using experimentally obtained displacements. Comparison was made withT _ε ^* calculated with the measured surface displacements andT _ε ^* calculated at several points through the thickness of a finite element (FE) model of the SE(B) specimen. Good comparison was found between the surfaceT _ε ^* calculated from displacements extracted from the FE model and the surfaceT _ε ^* calculated from experimentally obtained displacements. The computedT _ε ^* integral was also observed to decrease as the crack front was traversed from the surface to the mid-plane of the specimen. Mid-planeT _ε ^* values tend to be approximately 10% of the surface values.     

<Emphasis Type="Italic">L</Emphasis><Superscript>1</Superscript> Stability of the Boltzmann Equation for the Hard-Sphere Model

We study the L¹ stability of classical solutions to the Boltzmann equation for a hard-sphere model, when initial datum is a small perturbation of a vacuum, and tends to zero exponentially fast at infinity in the phase space. For this, we introduce nonlinear functionals measuring potential interactions between particles with different velocities and L¹ distance between classical solutions. We use pointwise estimates for a solution and the gain term of a collision operator to control the time-evolution of nonlinear functionals.Dedicated to Marshall Slemrod on the occasion of his 60th birthday     

Convergence Rates in <Emphasis Type="Italic">L</Emphasis><Superscript>2</Superscript> for Elliptic Homogenization Problems

We study rates of convergence of solutions in L ² and H ^1/2 for a family of elliptic systems {L_e}{\{\mathcal{L}_\varepsilon\}} with rapidly oscillating coefficients in Lipschitz domains with Dirichlet or Neumann boundary conditions. As a consequence, we obtain convergence rates for Dirichlet, Neumann, and Steklov eigenvalues of {L_e}{\{\mathcal{L}_\varepsilon\}} . Most of our results, which rely on the recently established uniform estimates for the L ² Dirichlet and Neumann problems in Kenig and Shen (Math Ann 350:867–917, 2011; Commun Pure Appl Math 64:1–44, 2011) are new even for smooth domains.     

A Theory of <Emphasis Type="Italic">L</Emphasis><Superscript>1</Superscript>-Dissipative Solvers for Scalar Conservation Laws with Discontinuous Flux

We propose a general framework for the study of L ¹ contractive semigroups of solutions to conservation laws with discontinuous flux:

A new zig-zag theory and C<Superscript>0</Superscript> plate bending element for composite and sandwich plates

A higher-order zig-zag theory for laminated composite and sandwich structures is proposed. The proposed theory satisfies the interlaminar continuity conditions and free surface conditions of transverse shear stresses. Moreover, the number of unknown variables involved in present model is independent of the number of layers. Compared to the zig-zag theory available in literature, the merit of present theory is that the first derivatives of transverse displacement have been taken out from the in-plane displacement fields, so that the C⁰ interpolation functions is only required during its finite element implementation. To obtain accurately transverse shear stresses by integrating three-dimensional equilibrium equations within one element, a six-node triangular element is employed to model the present zig-zag theory. Numerical results show that the present zig-zag theory can predict more accurate in-plane displacements and stresses in comparison with other zig-zag theories. Moreover, it is convenient to obtain transverse shear stresses by integration of equilibrium equations, as the C⁰ shape functions is only used when implemented in a finite element.     

<Emphasis Type="Italic">L</Emphasis><Superscript>2</Superscript> Formulation of Multidimensional Scalar Conservation Laws

We show that Kruzhkov’s theory of entropy solutions to multidimensional scalar conservation laws (Kruzhkov in Mat Sb (N.S.), 81(123), 228–255, 1970) can be entirely recast in L ² and fits into the general theory of maximal monotone operators in Hilbert spaces. Our approach is based on a combination of level-set, kinetic and transport-collapse approximations, in the spirit of previous works by Brenier (in C R Acad Sci Paris Ser I Math, 292, 563–566, 1981; in J Diff Equ, 50, 375–390, 1983; in SIAM J Numer Anal, 21, 1013–1037; in Methods Appl Anal, 11, 515–532, 2004), Giga and Miyakawa (in Duke Math J, 50, 505–515, 1983), and Tsai et al. (in Math Comp, 72, 159–181, 2003).     

A Remark on <Emphasis Type="Italic">L</Emphasis><Superscript>∞</Superscript> Solutions to the 2-D Navier–Stokes Equations

A single-exponential growth estimate of the solutions to the 2-dimensional Navier–Stokes equations in the whole space for nondecaying initial velocity is established. The crucial idea is to decompose velocity into high and low frequency parts. Moreover, if the linear term and the vorticity decay exponentially, the velocity is bounded uniformly in time.      

Solutions for Quasilinear Nonsmooth Evolution Systems in <Emphasis Type="Italic">L</Emphasis><Superscript><Emphasis Type="Italic">p</Emphasis></Superscript>

We prove that nonsmooth quasilinear parabolic systems admit a local solution in L ^p strongly differentiable with respect to time over a bounded three-dimensional polyhedral space domain. The proof rests essentially on new elliptic regularity results for polyhedral Laplace interface problems for anisotropic materials. These results are based on sharp pointwise estimates for Greens function, which are also of independent interest. To treat the nonlinear problem, we then apply a classical theorem of Sobolevskii for abstract parabolic equations and recently obtained resolvent estimates for elliptic operators and interpolation results. As applications we have in mind primarily reaction-diffusion systems. The treatment of such equations in an L ^p context seems to be new and allows (by Gauss theorem) the proper definition of the normal component of currents across the boundary.     

On the <Emphasis Type="Italic">H</Emphasis><Superscript><Emphasis Type="Italic">s</Emphasis></Superscript> Theory of Hydrostatic Euler Equations

In this paper we study the two-dimensional hydrostatic Euler equations in a periodic channel. We prove the local existence and uniqueness of H ^s solutions under the local Rayleigh condition. This extends Brenier’s (Nonlinearity 12(3):495–512, 1999) existence result by removing an artificial condition and proving uniqueness. In addition, we prove weak–strong uniqueness, mathematical justification of the formal derivation and stability of the hydrostatic Euler equations. These results are based on weighted H ^s a priori estimates, which come from a new type of nonlinear cancellation between velocity and vorticity.     

Regularity and finite dimensionality of attractor for plate equation on R<Superscript><Emphasis Type="Italic">n</Emphasis></Superscript>

This paper addresses the regularity and finite dimensionality of the global attractor for the plate equation on the unbounded domain. The existence of the attractor in the phase space has been established in an earlier work of the author. It is shown that the attractor is actually a bounded set of the phase space and has finite fractal dimensionality.     

Approximation of Flat <Emphasis Type="Italic">W</Emphasis><Superscript>2,2</Superscript> Isometric Immersions by Smooth Ones

Let \({S\subset\mathbb{R}^2}\) be a bounded Lipschitz domain and denote by \({W^{2,2}_{\text{iso}}(S; \mathbb{R}^3)}\) the set of mappings \({u\in W^{2,2}(S;\mathbb{R}^3)}\) which satisfy \({(\nabla u)^T(\nabla u) = Id}\) almost everywhere. Under an additional regularity condition on the boundary \({\partial S}\) (which is satisfied if \({\partial S}\) is piecewise continuously differentiable), we prove that the strong W ^2,2 closure of \({W^{2,2}_{\text{iso}}(S; \mathbb{R}^3)\cap C^{\infty}(\overline{S};\mathbb{R}^3)}\) agrees with \({W^{2,2}_{\text{iso}}(S; \mathbb{R}^3)}\).     

<Emphasis Type="Italic">L</Emphasis><Superscript>2</Superscript> Stability Estimates for Shock Solutions of Scalar Conservation Laws Using the Relative Entropy Method

We consider scalar nonviscous conservation laws with strictly convex flux in one spatial dimension, and we investigate the behavior of bounded L ² perturbations of shock wave solutions to the Riemann problem using the relative entropy method. We show that up to a time-dependent translation of the shock, the L ² norm of a perturbed solution relative to the shock wave is bounded above by the L ² norm of the initial perturbation.     

HAMILTON OPERATORS AND HOMOTHETIC MOTIONS IN R^3

Quaternion is a division ring. It is shown that planes passing through the origin can be made a field with the quaternion product in R3. The Hamiltonian operators help us define the homothetic motions on these planes. New characterizations for these motions are investigated.     

On solutions continuously differentiable on ℝ<Superscript>+</Superscript> for systems of linear functional differential equations with linearly transformed argument

We obtain conditions for the existence of solutions continuous on ℝ⁺ for a system of linear functional differential equations with linearly transformed argument. __________ Translated from Neliniini Kolyvannya, Vol. 10, No. 3, pp. 322–327, July–September, 2007.     

Three-Dimensional Simulation of Turbulent Hot-Jet Ignition for Air-CH<Subscript>4</Subscript>-H<Subscript>2</Subscript> Deflagration in a Confined Volume

This work describes essential aspects of the ignition and deflagration process initiated by the injection of a hot transient gas jet into a narrowly confined volume containing air-CH₄-H₂ mixture. Driven by the pressure difference between a prechamber and a long narrow constant-volume-combustion (CVC) chamber, the developing jet or puff involves complex processes of turbulent jet penetration and evolution of multi-scale vortices in the shear layer, jet tip, and adjacent confined spaces. The CVC chamber contains stoichiometric mixtures of air with gaseous fuel initially at atmospheric conditions. Fuel reactivity is varied using two different CH₄/H₂ blends. Jet momentum is varied using different pre-chamber pressures at jet initiation. The jet initiation and the subsequent ignition events generate pressure waves that interact with the mixing region and the propagating flame, depositing baroclinic vorticity. Transient three-dimensional flow simulations with detailed chemical kinetics are used to model CVC mixture ignition. Pre-ignition gas properties are then examined to develop and verify criteria to predict ignition delay time using lower-cost non-reacting flow simulations for this particular case of study.     

An Experimental Study of CO<Subscript>2</Subscript> Exsolution and Relative Permeability Measurements During CO<Subscript>2</Subscript> Saturated Water Depressurization

Dissolution of CO₂ into brine is an important and favorable trapping mechanism for geologic storage of CO₂. There are scenarios, however, where dissolved CO₂ may migrate out of the storage reservoir. Under these conditions, CO₂ will exsolve from solution during depressurization of the brine, leading to the formation of separate phase CO₂. For example, a CO₂ sequestration system with a brine-permeable caprock may be favored to allow for pressure relief in the sequestration reservoir. In this case, CO₂-rich brine may be transported upwards along a pressure gradient caused by CO₂ injection. Here we conduct an experimental study of CO₂ exsolution to observe the behavior of exsolved gas under a wide range of depressurization. Exsolution experiments in highly permeable Berea sandstones and low permeability Mount Simon sandstones are presented. Using X-ray CT scanning, the evolution of gas phase CO₂ and its spatial distribution is observed. In addition, we measure relative permeability for exsolved CO₂ and water in sandstone rocks based on mass balances and continuous observation of the pressure drop across the core from 12.41 to 2.76 MPa. The results show that the minimum CO₂ saturation at which the exsolved CO₂ phase mobilization occurs is from 11.7 to 15.5%. Exsolved CO₂ is distributed uniformly in homogeneous rock samples with no statistical correlation between porosity and CO₂ saturation observed. No gravitational redistribution of exsolved CO₂ was observed after depressurization, even in the high permeability core. Significant differences exist between the exsolved CO₂ and water relative permeabilities, compared to relative permeabilities derived from steady-state drainage relative permeability measurements in the same cores. Specifically, very low CO₂ and water relative permeabilities are measured in the exsolution experiments, even when the CO₂ saturation is as high as 40%. The large relative permeability reduction in both the water and CO₂ phases is hypothesized to result from the presence of disconnected gas bubbles in this two-phase flow system. This feature is also thought to be favorable for storage security after CO₂ injection.     

Transport of Lead (Pb<Superscript>2+</Superscript>) Ions Through Silty-Clayey Soils Under Acidic Conditions

This study aimed to identify effects of pH on the transport of Pb²⁺ ions through a saturated silty-clayey soil layer by using advection–dispersion equation (ADE). The predictive accuracy of the solution of ADE depends on the proper determination of the retardation by adsorption and, therefore, the adsorption mechanism of lead onto silty-clayey soil was investigated first by performing batch equilibrium experiments. These results showed that the sorption mechanism of lead onto silty-clayey soil depended on pH and could be best described by the Langmuir isotherm. Based on the results of the sequential experiments, it was also concluded that the pH dependent charges in silty-clayey soil were mainly associated with the surfaces of carbonates and the specific adsorption of lead ions. The numerical solutions of the combined form of ADE with the Langmuir isotherm indicated that the migration profiles of lead in silty-clayey soil were a strong function of the parameters of the Langmuir isotherm rather than the infiltration velocity.     

Injection and Storage of CO<Subscript>2</Subscript> in Deep Saline Aquifers: Analytical Solution for CO<Subscript>2</Subscript> Plume Evolution During Injection

Injection of fluids into deep saline aquifers is practiced in several industrial activities, and is being considered as part of a possible mitigation strategy to reduce anthropogenic emissions of carbon dioxide into the atmosphere. Injection of CO₂ into deep saline aquifers involves CO₂ as a supercritical fluid that is less dense and less viscous than the resident formation water. These fluid properties lead to gravity override and possible viscous fingering. With relatively mild assumptions regarding fluid properties and displacement patterns, an analytical solution may be derived to describe the space–time evolution of the CO₂ plume. The solution uses arguments of energy minimization, and reduces to a simple radial form of the Buckley–Leverett solution for conditions of viscous domination. In order to test the applicability of the analytical solution to the CO₂ injection problem, we consider a wide range of subsurface conditions, characteristic of sedimentary basins around the world, that are expected to apply to possible CO₂ injection scenarios. For comparison, we run numerical simulations with an industry standard simulator, and show that the new analytical solution matches a full numerical solution for the entire range of CO₂ injection scenarios considered. The analytical solution provides a tool to estimate practical quantities associated with CO₂ injection, including maximum spatial extent of a plume and the shape of the overriding less-dense CO₂ front.