Drainage of Capillary-Trapped Oil by an Immiscible Gas: Impact of Transient and Steady-State Water Displacement on Three-Phase Oil Permeability

In a previous paper (Dehghanpour et al., Phys Rev E 83:065302, 2011a), we showed that relative permeability of mobilized oil, $k_\mathrm{ro}$ , measured during tertiary gravity drainage, is significantly higher than that of the same oil saturation in other tests where oil is initially a continuous phase. We also showed that tertiary $k_\mathrm{ro}$ strongly correlates to both water saturation, $S_\mathrm{w}$ , water flux (water relative permeability), $k_\mathrm{rw}$ , and the change in water saturation with time, $\mathrm{d}S_\mathrm{w}/\mathrm{d}t$ . To develop a model and understanding of the enhanced oil transport, identifying which of these parameters ( $S_\mathrm{w},\,k_{\mathrm{rw}}$ , or $\mathrm{d}S_\mathrm{w}/\mathrm{d}t$ ) plays the controlling role is necessary, but in the previous experiments these could not be deconvolved. To answer the remaining question, we conduct specific three-phase displacement experiments in which $k_{\mathrm{rw}}$ is controlled by applying a fixed water influx, and $S_\mathrm{w}$ develops naturally. We obtain $k_{\mathrm{ro}}$ by using the saturation data measured in time and space. The results suggest that steady-state water influx, in contrast to transient water displacement, does not enhance $k_{\mathrm{ro}}$ . Instead, reducing water influx rate results in excess oil flow. Furthermore, according to our pore scale hydraulic conductivity calculations, viscous coupling and fluid positioning do not sufficiently explain the observed correlation between $k_{\mathrm{ro}}$ and $S_{\mathrm{w}}$ . We conclude that tertiary $k_{\mathrm{ro}}$ is controlled by the oil mobilization rate, which in turn is linked to the rate of water saturation decrease with time, $\mathrm{d}S_\mathrm{w}/\mathrm{d}t$ . Finally, we develop a simple model which relates tertiary $k_{\mathrm{ro}}$ to transient two-phase gas/water relative permeability.     

Predominant Gas Transport in Microporous Hydrotalcite–Silica Membrane

The prepared microporous hydrotalcite (HT)–silica membrane was found to exhibit the molecular sieving characteristic of pristine silica material and high $\mathrm{CO}_{2}$ adsorption capacity of HT. The combined properties made enhanced $\mathrm{CO}_{2}$ permeability and separability from $\mathrm{CH}_{4}$ possible. The gas transport in the membrane was predominantly surface adsorption. The porous membrane overcame the Knudsen limitation and yielded the highest separation selectivity of 120 at 40 % $\mathrm{CO}_{2}$ feed concentration, $30\,^{\circ }\mathrm{C}$ operating temperature, and 100 kPa pressure difference.     

Experimental Characterization of Porosity Structure and Transport Property Changes in Limestone Undergoing Different Dissolution Regimes

Limestone dissolution by $\hbox {CO}_2$ -rich brine induces critical changes of the pore network geometrical parameters such as the pore size distribution, the connectivity, and the tortuosity which govern the macroscopic transport properties (permeability and dispersivity) that are required to parameterize the models, simulating the injection and the fate of $\hbox {CO}_2$ . A set of four reactive core-flood experiments reproducing underground conditions ( $T = 100\,^{\circ }\hbox {C}$ and $P = 12$ MPa) has been conducted for different $\hbox {CO}_2$ partial pressures $(0.034 < P_{\mathrm{CO}_2}< 3.4\; \hbox {MPa})$ in order to study the different dissolution regimes. X-ray microtomographic images have been used to characterize the changes in the structural properties from pore scale to Darcy scale, while time-resolved pressure loss and chemical fluxes enabled the determination of the sample-scale change in porosity and permeability. The results show the growth of localized dissolution features associated with high permeability increase for the highest $P_{\mathrm{CO}_2}$ , whereas dissolution tends to be more homogeneously distributed for lower values of $P_{\mathrm{CO}_2}$ . For the latter, the higher the $P_{\mathrm{CO}_2}$ , the more the dissolution patterns display ramified structures and permeability increase. For the lowest value of $P_{\mathrm{CO}_2}$ , the preferential dissolution of the calcite cement associated with the low dissolution kinetics triggers the transport that may locally accumulate and form a microporous material that alters permeability and produces an anti-correlated porosity–permeability relationship. The combined analysis of the pore network geometry and the macroscopic measurements shows that $P_{\mathrm{CO}_2}$ regulates the tortuosity change during dissolution. Conversely, the increase of the exponent value of the observed power law permeability–porosity trend while $P_{\mathrm{CO}_2}$ increases, which appears to be strongly linked to the increase of the effective hydraulic diameter, depends on the initial rock structure.     

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

Double-Diffusive Natural Convection in Anisotropic Porous Medium Bounded by Finite Thickness Walls: Validity of Local Thermal Equilibrium Assumption

Double-diffusive natural convection in fluid-saturated porous medium inside a vertical enclosure bounded by finite thickness walls with opposing temperature, concentration gradients on vertical walls as well as adiabatic and impermeable horizontal ones has been performed numerically. The Darcy model was used to predict fluid flow inside the porous material, while thermal fields are simulated based on two-energy equations for fluid and solid phases on the basis of a local thermal non-equilibrium model. Computations have been performed for different controlling parameters such as the buoyancy ratio $N$ , the Lewis number Le, the anisotropic permeability ratio $R_\mathrm{p}$ , the fluid-to-solid thermal conductivity ratio $R_\mathrm{c}$ , the interphase heat transfer coefficient $\mathcal{H}$ , the ratio of the wall thickness to its height $D$ , the wall-to-porous medium thermal diffusivity ratio $R_\mathrm{w}$ , and the solid-to-fluid heat capacity ratio $\gamma $ . Thus, the effects of the controlling parameters on heat and mass transfer characteristics are discussed in detail. Moreover, the validity domain of the local thermal equilibrium (LTE) assumption has been delimited for different set of the governing parameters. It has been shown that Le has a noticeable significant effect on fluid temperature profiles and that higher $N$ values lead to a significant enhancement in heat and mass transfer rates. Moreover, for higher $\mathcal{H}, R_\mathrm{c}$ , $R_\mathrm{p}, R_\mathrm{w}$ , or $D$ values and/or lower $\gamma $ values, the solid and fluid phases tend toward LTE.     

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.     

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.     

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.     

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.     

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

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.     

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

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.     

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.     

Estimability Analysis and Optimisation of Soil Hydraulic Parameters from Field Lysimeter Data

Modelling the water-flow in the vadose zone requires accurate hydraulic parameters to be obtained at the relevant scale. Weighable lysimeters enable us to monitor hydraulic data at an intermediate scale between lab and field scales and they can be used to optimise these parameters. Parameter optimisation using inverse methods may be limited by the non-uniqueness of the solution. In this contribution, an estimability method has been used to assess the estimability of the van Genuchten–Mualem parameters, to evaluate the information content of the data collected from a bare field lysimeter and to optimise the estimable model parameters. Daily data were monitored from a 2  $\text{ m }^{3}$ lysimeter, filled with the soil of a former coking plant: pressure heads and water contents were measured at three depths (50, 100, 150 cm), cumulative boundary water fluxes. Water-flow was represented using the one-dimensional single-porosity model implemented in HYDRUS-1D code. The estimability of the 5 van Genuchten–Mualem hydraulic parameters and the information content of different data were evaluated by sequentially calculating a sensitivity coefficient matrix. Optimisation was achieved by the Levenberg-Marquardt algorithm. The estimability analysis revealed that estimability of the soil hydraulic parameters, based on the combination of daily pressure heads and water contents, was higher than those based on these data separately. In case of 2.4 being considered as a cut-off criterion for this study, all the parameters were considered estimable from daily data in the decreasing order: $\theta _\mathrm{s},\, \textit{n}, K_\mathrm{s},\, \alpha ,\, \theta _\mathrm{r}$ . Hydraulic parameters were optimised in four scenarios: $\theta _\mathrm{s}$ and n were estimated with reliability while $\alpha ,\, K_\mathrm{s}$ and $\theta _\mathrm{r}$ were uncertain. However, the narrow variations in measured data restricted parameter optimisation.     

On the effect of non-smooth Coulomb damping on flutter-type self-excitation in a non-gyroscopic circulatory 2-DoF-system

This article deals with self-excited vibrations, attractivity of stationary solutions, and the corresponding bifurcation behavior of two-dimensional differential inclusions of the type $\mathbf{M}\mathbf{q}'' + \mathbf{D}\mathbf{q}' + (\mathbf{K} + \bar{\mu}\mathbf{N})\mathbf{q} \in-\mathbf{R}\operatorname{Sign}(\mathbf{q}')$ . For the smooth case R=0, the equilibrium may become unstable due to non-conservative positional forces stemming from the circulatory matrix N. This type of instability is usually referred to as flutter instability and the loss of stability is related to a Hopf bifurcation of the steady state, which occurs for a critical parameter $\bar{\mu}= \bar{\mu}_{\mathrm{crit}}$ . For R≠0, the steady state is a set of equilibria, which turns out to be attractive for all values of the bifurcation parameter $\bar{\mu}$ . Depending on $\bar{\mu}$ , the basin of attraction of the equilibrium set can be infinite or finite. The transition from an infinite to a finite basin of attraction occurs at the stability threshold $\bar{\mu}_{\mathrm{crit}}$ of the underlying smooth problem. For the finite basin of attraction, its size is proportional to the Coulomb friction and inverse-proportional to $(\bar{\mu}- \bar{\mu}_{\mathrm{crit}})$ . By adding Coulomb damping the notion of steady state stability for the smooth problem is replaced by the question whether the basin of attraction of the steady state is infinite or finite. Simultaneously, the local Hopf-bifurcation is replaced by a global bifurcation. This implies that in the presence of Coulomb damping the occurrence of self-excited vibrations can only be investigated with regard to the perturbation level.     

Random Data Cauchy Theory for the Generalized Incompressible Navier–Stokes Equations

In this paper, we consider the generalized Navier?CStokes equations where the space domain is ${\mathbb{T}^N}$ or ${\mathbb{R}^N, N\geq3}$ . The generalized Navier?CStokes equations here refer to the equations obtained by replacing the Laplacian in the classical Navier?CStokes equations by the more general operator (???) ^?? with ${\alpha\in (\frac{1}{2},\frac{N+2}{4})}$ . After a suitable randomization, we obtain the existence and uniqueness of the local mild solution for a large set of the initial data in ${H^s, s\in[-\alpha,0]}$ , if ${1 < \alpha < \frac{N+2}{4}, s\in(1-2\alpha,0]}$ , if ${\frac{1}{2} < \alpha\leq 1}$ . Furthermore, we obtain the probability for the global existence and uniqueness of the solution. Specially, our result shows that, in some sense, the Cauchy problem of the classical Navier?CStokes equation is local well-posed for a large set of the initial data in H ^?1+, exhibiting a gain of ${\frac{N}{2}-}$ derivatives with respect to the critical Hilbert space ${H^{\frac{N}{2}-1}}$ .     

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.