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.     

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.     

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.     

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.     

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

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.     

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.     

Impact of Effective Stress and Matrix Deformation on the Coal Fracture Permeability

The permeability of coal is an important parameter in mine methane control and coal bed methane exploitation because it determines the practicability of methane extraction. We developed a new coal permeability model under tri-axial stress conditions. In our model, the coal matrix is compressible and Biot’s coefficient, which is considered to be 1 in existing models, varies between 0 and 1. Only a portion of the matrix deformation, which is represented by the effective coal matrix deformation factor $f_\mathrm{m}$ , contributes to fracture deformation. The factor $f_\mathrm{m}$ is a parameter of the coal structure and is a constant between 0 and 1 for a specific coal. Laboratory tests indicate that the Sulcis coal sample has an $f_\mathrm{m}$ value of 0.1794 for $\hbox {N}_{2}$ and $\hbox {CO}_{2}$ . The proposed permeability model was evaluated using published data for the Sulcis coal sample and is compared to three popular permeability models. The proposed model agrees well with the observed permeability changes and can predict the permeability of coal better than the other models. The sensitivity of the new model to changes in the physical, mechanical and adsorption deformation parameters of the coal was investigated. Biot’s coefficient and the bulk modulus mainly affect the effective stress term in the proposed model. The sorption deformation parameters and the factor $f_\mathrm{m}$ affect the coal matrix deformation term.     

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.     

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.     

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.     

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.     

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.     

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.     

Conjugate Natural Convection in a Porous Enclosure Sandwiched by Finite Walls Under the Influence of Non-uniform Heat Generation and Radiation

Conjugate natural convection in a square porous enclosure sandwiched by finite walls under the influence of non-uniform heat generation and radiation is studied numerically in the present article. The horizontal heating is considered, where the vertical walls heated isothermally at different temperatures, while the horizontal walls are kept adiabatic. The Darcy model is used in the mathematical formulation for the porous layer and finite difference method is applied to solve the dimensionless governing equations. The governing parameters considered are the ratio of wall thickness to its width $(0.02 \le D \le 0.3)$ ( 0.02 ≤ D ≤ 0.3 ) , the wall to porous thermal conductivity ratio $(0.1 \le k_\mathrm{r} \le 10.0)$ ( 0.1 ≤ k r ≤ 10.0 ) , the internal heating $(0 \le \gamma \le 5)$ ( 0 ≤ γ ≤ 5 ) , and the local heating exponent parameters $(1 \le \lambda \le 20)$ ( 1 ≤ λ ≤ 20 ) . It is found that the average Nusselt number on the hot and cold interfaces increases with increasing the radiation intensity. Very high non-uniformity heating does not affect the average Nusselt number at very thick walls.     

Mixtures of foam and paste: suspensions of bubbles in yield stress fluids

We study the rheological behavior of mixtures of foams and pastes, which can be described as suspensions of bubbles in yield stress fluids. Model systems are designed by mixing monodisperse aqueous foams and concentrated emulsions. The elastic modulus of the bubble suspensions is found to depend on the elastic capillary number $\textit{Ca}_{_G}$ , defined as the ratio of the paste elastic modulus to the bubble capillary pressure. For values of $\textit{Ca}_{_G}$ larger than $\simeq 0.5$ , the dimensionless elastic modulus of the aerated material decreases as the bubble volume fraction $\phi $ increases, suggesting that bubbles behave as soft elastic inclusions. Consistently, this decrease is all the sharper as $\textit{Ca}_{_G}$ is high, which accounts for the softening of the bubbles as compared to the paste. By contrast, we find that the yield stress of most studied materials is not modified by the presence of bubbles. This suggests that their plastic behavior is governed by the plastic capillary number $\textit{Ca}_{\tau_y}$ , defined as the ratio of the paste yield stress to the bubble capillary pressure. At low $\textit{Ca}_{\tau_y}$ values, bubbles behave as nondeformable inclusions, and we predict that the suspension dimensionless yield stress should remain close to unity, in agreement with our data up to $\textit{Ca}_{\tau_y}=0.2$ . When preparing systems with a larger target value of $\textit{Ca}_{\tau_y}$ , we observe bubble breakup during mixing, which means that they have been deformed by shear. It then seems that a critical value $\textit{Ca}_{\tau_y}\simeq 0.2$ is never exceeded in the final material. These observations might imply that, in bubble suspensions prepared by mixing a foam and a paste, the suspension yield stress is always close to that of the paste surrounding the bubbles. Finally, at the highest $\phi $ investigated, the yield stress is shown to increase abruptly with $\phi $ : this is interpreted as a “foamy yield stress fluid” regime, which takes place when the paste mesoscopic constitutive elements (here, the oil droplets) are strongly confined in the films between the bubbles.     

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.     

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.     

Invariant Parallels, Invariant Meridians and Limit Cycles of Polynomial Vector Fields on Some 2-Dimensional Algebraic Tori in \mathbb R ^3

We consider the polynomial vector fields of arbitrary degree in $\mathbb R ^3$ R 3 having the 2-dimensional algebraic torus $$\begin{aligned} \mathbb T ^2(l,m,n)=\{(x,y,z)\in \mathbb R ^3 : (x^{2l}+y^{2m}-r^2)^2+z^{2n}-1=0\}, \end{aligned}$$ T 2 ( l , m , n ) = { ( x , y , z ) ∈ R 3 : ( x 2 l + y 2 m - r 2 ) 2 + z 2 n - 1 = 0 } , where $l,m$ l , m , and $n$ n positive integers, and $r\in (1,\infty )$ r ∈ ( 1 , ∞ ) , invariant by their flow. We study the possible configurations of invariant meridians and parallels that these vector fields can exhibit on $\mathbb T ^2(l,m,n)$ T 2 ( l , m , n ) . Furthermore, we analyze when these invariant meridians or parallels are limit cycles.