Experimental exploration of underexpanded supersonic jets

Two underexpanded free jets at fully expanded Mach numbers $M_\mathrm{j}$ = 1.15 and 1.50 are studied. Schlieren visualizations as well as measurements of static pressure, Pitot pressure and velocity are performed. All these experimental techniques are associated to obtain an accurate picture of the jet flow development. In particular, expansion, compression and neutral zones have been identified in each shock cell. Particle lag is considered by integrating the equation of motion for particles in a fluid flow and it is found that the laser Doppler velocimetry is suitable for investigating shock-containing jets. Even downstream of the normal shock arising in the $M_\mathrm{j}$ = 1.50 jet, the measured gradual velocity decrease is shown to be relevant.     

Measurements of Transitional and Turbulent Flow in a Randomly Packed Bed of Spheres with Particle Image Velocimetry

Particle image velocimetry (PIV) has been used to investigate transitional and turbulent flow in a randomly packed bed of mono-sized transparent spheres at particle Reynolds number, \(20<{{ Re}}_{\mathrm{p}}< 3220\). The refractive index of the liquid is matched with the spheres to provide optical access to the flow within the bed without distortions. Integrated pressure drop data yield that Darcy law is valid at \({{ Re}}_{\mathrm{p}} \approx 80\). The PIV measurements show that the velocity fluctuations increase and that the time-averaged velocity distribution start to change at lower \({{ Re}}_{\mathrm{p}}\). The probability for relatively low and high velocities decreases with \({{ Re}}_{\mathrm{p}}\) and recirculation zones that appear in inertia dominated flows are suppressed by the turbulent flow at higher \({{ Re}}_{\mathrm{p}}\). Hence there is a maximum of recirculation at about \({{ Re}}_{\mathrm{p}} \approx 400\). Finally, statistical analysis of the spatial distribution of time-averaged velocities shows that the velocity distribution is clearly and weakly self-similar with respect to \({{ Re}}_{\mathrm{p}}\) for turbulent and laminar flow, respectively.     

An experimental investigation of a wall jet in a torus

Results of an experimental investigation of the velocity profile of a turbulent gas injected in a toroidal configuration are presented. The measurements surprisingly show that it is possible to describe the radial distribution of the azimuthal velocity in terms of a plane wall jet discharging in an external stream. The growth of the inner boundary layer, the width of the jet, and the velocity profile are in accordance with the known experimental data on this subject. A fundamentally different relation has been deduced for the decay of the maximum velocity. Up to now Sigalla's formula \(U_m /U_j \propto \sqrt {a/x}\) is generally accepted. Our data based on an essentially extended range of x/a, correlate with the exponential relation $$U_m /U_j = exp\left[ { - 154(Re_x )^{ - 0.777} \frac{x}{a}} \right].$$      

Simultaneous measurement of velocity and pressure in a plane jet

A new technique was developed for the simultaneous measurement of velocity and pressure in turbulent flows. To accomplish this objective, a new probe (hereafter called the combined probe) that consists of an X-type hot-wire probe and a newly devised pressure probe was developed. The pressure probe was miniaturized by the MEMS fabrication process and by using a 0.1-in. microphone as a pressure sensor for improving the spatial resolution. This pressure probe was placed between two hot-wire sensors of which the X-type hot-wire probe was composed. The pressure probe was given a hemispherical tip, like that of a pitot tube, because an earlier pressure probe with a conical tip suffered from a reduction in spatial resolution. The spatial arrangement of the pressure probe and the hot-wire probe for the combined probe was carefully determined, because there was a risk that the measurement accuracy of one probe will be influenced by disturbances caused by the other probe when the two probes were placed very close to each other. Therefore, the combined probe was arranged to engender no noticeable interference between the velocity data and the pressure data measured by their respective probes. As one application of this combined probe, simultaneous measurements of pressure and two components of instantaneous velocity were performed in a plane jet. The turbulent energy budget and the cross-correlation coefficient of velocity and pressure in the intermittent region of the plane jet were estimated. The results show that the mean streamwise velocity, velocity fluctuation, and pressure fluctuation profiles were consistent with those measured individually using the X-type hot-wire probe or pressure probe. Moreover, it was shown that the integral value of the diffusion term (which should theoretically be equal to zero) in the turbulent energy transport equation was closer to zero than previous reports (Bradbury in J Fluid Mech 23(Part 1):31–64, 1965). In addition, the time variation of the cross-correlation coefficient in the intermittent region supports the vortex structure model predicted in previous studies (Browne et?al. in J Fluid Mech 149:355–373, 1984; Tanaka et?al. JSME Int J Ser B 49(4):899–905, 2006; Sakai et?al. J Fluid Sci Technol 2(3):611–622, 2007).     

Super-equilibrium increase of chemical reaction rate in the detonation front and other effects in the detonation wave initiated by a shock wave

In the present work the problem of detonation wave formation in a shock tube was considered in one-dimensional formulation. The Monte Carlo non-stationary method of statistical simulation (MCNMSS), also known as DSMC, was used for simulation. The method automatically takes into account all details of mass and heat transfer. At an initial moment, the low-pressure channel (LPC) of the shock tube was filled with gas A while the high-pressure chamber (HPC) was filled with gas C. The cross-sections of the HPC and LPC, as well as the temperatures of gases A and C were equal to each other. At the beginning of the simulation the ratio of pressures in the HPC and LPC was equal to 100. It was assumed that chemical reactions $\mathrm{{A}}+\mathrm{{M}} \rightarrow \mathrm{{B}}+\mathrm{{M}}$ ( $\mathrm{{M}}=\mathrm{{A}},\, \mathrm{{B}}$ and $\mathrm{{C}}$ ) took place. The ratio of molecular masses of gases $\mathrm{{A}},\, \mathrm{{B}}$ and $\mathrm{{C}}$ was taken as 20:20:1. Different reaction thresholds were considered. For the case of a low reaction threshold, the velocity of the resulting detonation wave was found to be higher than the Chapman–Jouguet velocity. A region with constant values of flow parameters inside product was observed. An increase of the reaction threshold led to disappearance of this region and gave rise to something similar to an expansion wave, with peaks of flow parameters at the leading part of the detonation wave. The values of these peaks were found to be constant in time. The velocity of the detonation wave became appreciably lower than the Chapman–Jouguet velocity. Further increase of the reaction threshold led to disappearance of detonation. The reactions $\mathrm{{A}}+\mathrm{{B}} \rightarrow \mathrm{{B}}+\mathrm{{B}}$ and $\mathrm{{A}}+\mathrm{{C}}\rightarrow \mathrm{{B}}+\mathrm{{C}}$ turned out to be very important for initiation of detonation.     

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.     

Mean Flow and Turbulence Measurements in a Turbulent Free Cruciform Jet

The results of measurements of all three components of the mean velocity vector, the Reynolds normal and primary shear stresses and the mean static pressure in a turbulent free jet, issuing from a sharp-edged cruciform orifice, are presented in this paper. The measurements were made with an x-array hot-wire probe and a pitot-static tube in the near flow field of the jet. The Reynolds number, based upon the equivalent diameter of the orifice, was 1.70 × 10⁵. In addition to the quantities measured directly, the mean streamwise centreline velocity decay, the jet half-velocity widths, the jet spreading rate, the mean streamwise vorticity, the mass entrainment rate, the integral momentum flux and the one-dimensional energy spectra have been derived from the measured data. The results show that the mean streamwise centreline velocity decay rate of the cruciform jet is higher than that of a round jet issuing from an orifice with the same exit area as that of the cruciform orifice. The mean streamwise velocity field changed shape continuously from a cruciform close to the orifice exit plane to circular at 12 and half equivalent diameters downstream. The mean streamwise vorticity field, up to about three equivalent diameters downstream of the orifice exit plane, consists of four pairs of counter-rotating cells, which are aligned with the four edges in the centre of the cruciform orifice.     

Simultaneous measurement of fluctuating velocity and pressure in the near wake of a circular cylinder

Simultaneous measurement of fluctuating velocity and pressure by a static-pressure probe and a hot-wire probe was performed in the near wake of a circular cylinder, in order to strengthen reliability of the measurement technique. Effect of geometry of the static-pressure probe was systematically investigated, and validity of the measurement results was addressed by quantitative comparison with reference data by a large-eddy simulation. Interference between the probes was found to mainly depend on the diameter of the pressure probe and only weakly on the length. A certain time lag between the velocity and pressure signals was detected in the experiment, and the measurement results of velocity–pressure correlation $\overline{up}$ and $\overline{vp}$ obtained with the correction of the time lag were in good agreement with the computational results. It was also found that the measurement of $\overline{vp}$ is extremely sensitive to a small time lag between the velocity and pressure signals, while that of $\overline{up}$ is not.     

Unsteady MHD double-diffusive convection boundary-layer flow past a radiate hot vertical surface in porous media in the presence of chemical reaction and heat sink

This paper presents an analytical study of the unsteady MHD free convective heat and mass transfer flow of a viscous, incompressible, gray, absorbing-emitting but non-scattering, optically-thick and electrically conducting fluid occupying a semi-infinite porous regime adjacent to an infinite moving hot vertical plate with constant velocity. We employ a Darcian viscous flow model for the porous medium the Rosseland diffusion approximation is used to describe the radiative heat flux in the energy equation. The homogeneous chemical reaction of first order is accounted in mass diffusion equation. The governing equations are solved in closed form by Laplace-transform technique. A parametric study of all involved parameters is conducted and representative set of numerical results for the velocity, temperature, concentration, shear stress function $\frac{\partial u}{\partial y} \vert_{y=0}$ , temperature gradient $\frac{\partial \theta }{ \partial y}\vert_{y=0}$ , and concentration gradient $\frac{ \partial \phi }{\partial y}\vert_{y=0}$ is illustrated graphically and physical aspects of the problem are discussed.     

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.     

Direct Numerical Simulation of Head-On Quenching of Statistically Planar Turbulent Premixed Methane-Air Flames Using a Detailed Chemical Mechanism

A three-dimensional compressible Direct Numerical Simulation (DNS) analysis has been carried out for head-on quenching of a statistically planar stoichiometric methane-air flame by an isothermal inert wall. A multi-step chemical mechanism for methane-air combustion is used for the purpose of detailed chemistry DNS. For head-on quenching of stoichiometric methane-air flames, the mass fractions of major reactant species such as methane and oxygen tend to vanish at the wall during flame quenching. The absence of \(\text {OH}\) at the wall gives rise to accumulation of carbon monoxide during flame quenching because \(\text {CO}\) cannot be oxidised anymore. Furthermore, it has been found that low-temperature reactions give rise to accumulation of \(\text {HO}_{2}\) and \(\mathrm {H}_{2}\mathrm {O}_{2}\) at the wall during flame quenching. Moreover, these low temperature reactions are responsible for non-zero heat release rate at the wall during flame-wall interaction. In order to perform an in-depth comparison between simple and detailed chemistry DNS results, a corresponding simulation has been carried out for the same turbulence parameters for a representative single-step Arrhenius type irreversible chemical mechanism. In the corresponding simple chemistry simulation, heat release rate vanishes once the flame reaches a threshold distance from the wall. The distributions of reaction progress variable c and non-dimensional temperature T are found to be identical to each other away from the wall for the simple chemistry simulation but this equality does not hold during head-on quenching. The inequality between c (defined based on \(\text {CH}_{4}\) mass fraction) and T holds both away from and close to the wall for the detailed chemistry simulation but it becomes particularly prominent in the near-wall region. The temporal evolutions of wall heat flux and wall Peclet number (i.e. normalised wall-normal distance of \(T = 0.9\) isosurface) for both simple and detailed chemistry laminar and turbulent cases have been found to be qualitatively similar. However, small differences have been observed in the numerical values of the maximum normalised wall heat flux magnitude \(\left ({\Phi }_{\max } \right )_{\mathrm {L}}\) and the minimum Peclet number \((Pe_{\min })_{\mathrm {L}}\) obtained from simple and detailed chemistry based laminar head-on quenching calculations. Detailed explanations have been provided for the observed differences in behaviours of \(\left ({\Phi }_{\max }\right )_{\mathrm {L}}\) and \((Pe_{\min })_{\mathrm {L}}\). The usual Flame Surface Density (FSD) and scalar dissipation rate (SDR) based reaction rate closures do not adequately predict the mean reaction rate of reaction progress variable in the near-wall region for both simple and detailed chemistry simulations. It has been found that recently proposed FSD and SDR based reaction rate closures based on a-priori DNS analysis of simple chemistry data perform satisfactorily also for the detailed chemistry case both away from and close to the wall without any adjustment to the model parameters.     

Statistical behavior of post-shock overpressure past grid turbulence

When a shock wave ejected from the exit of a 5.4-mm inner diameter, stainless steel tube propagated through grid turbulence across a distance of 215 mm, which is 5–15 times larger than its integral length scale \(L_{u}\) , and was normally incident onto a flat surface; the peak value of post-shock overpressure, \(\Delta P_{\mathrm{peak}}\) , at a shock Mach number of 1.0009 on the flat surface experienced a standard deviation of up to about 9 % of its ensemble average. This value was more than 40 times larger than the dynamic pressure fluctuation corresponding to the maximum value of the root-mean-square velocity fluctuation, \(u^{\prime }= 1.2~\hbox {m}/\hbox {s}\) . By varying \(u^{\prime }\) and \(L_{u}\) , the statistical behavior of \(\Delta P_{\mathrm{peak}}\) was obtained after at least 500 runs were performed for each condition. The standard deviation of \(\Delta P_{\mathrm{peak}}\) due to the turbulence was almost proportional to \(u^{{\prime }}\) . Although the overpressure modulations at two points 200 mm apart were independent of each other, we observed a weak positive correlation between the peak overpressure difference and the relative arrival time difference.     

An experimental study on turbulence modification in the near-wall boundary layer of a dilute gas-particle channel flow

Turbulence modifications of a dilute gas-particle flow are experimentally investigated in the lower boundary layer of a horizontal channel by means of a simultaneous two-phase PIV measurement technique. The measurements are conducted in the near-wall region with y ⁺?<?250 at Re _τ (based on the wall friction velocity u _τ and half channel height h)?=?430. High spatial resolution and small interrogation window are used to minimize the PIV measurement uncertainty due to the velocity gradient near the wall. Polythene beads with the diameter of 60?μm (d _p ⁺ ?=?1.71, normalized by the fluid kinematic viscosity ν and u _τ) are used as dispersed phase, and three low mass loading ratios (Φ _m ) ranging from 10^?4 to 10^?3 are tested. It is found that the addition of the particles noticeably modifies the mean velocity and turbulent intensities of the gas-phase, as well as the turbulence coherent structures, even at Φ _m ?=?0.025?%. Particle inertia changes the viscous sublayer of the gas turbulence with a smaller thickness and a larger streamwise velocity gradient, which increases the peak value of the streamwise fluctuation velocity ( $ u_{\text{rms}}^{ + } $ ) of the gas-phase with its location shifting to the wall. Particle sedimentation increases the roughness of the bottom wall, which significantly increases the wall-normal fluctuation velocity ( $ v_{\text{rms}}^{ + } $ ) and Reynolds shear stress ( $ - \langle u^{ \prime } v^{\prime } \rangle^{ + } $ ) of the gas-phase in the inner region of the boundary layer (y ⁺?<?10). Under effect of particle–wall collision, the Q2 events (ejections) of the gas-phase are slightly increased by particles, while the Q4 events (sweeps) are obviously decreased. The spatial scale of the coherent structures near the wall shrinks remarkably with the presence of the particles, which may be attributed to the intensified crossing-trajectory effects due to particle saltation near the bottom wall. Meanwhile, the $ v_{\text{rms}}^{ + } $ and $ - \langle u^{ \prime } v^{\prime } \rangle^{ + } $ of the gas-phase are significantly reduced in the outer region of the boundary layer (y ⁺?>?20).     

Multi-component Reynolds-averaged Navier–Stokes simulations of Richtmyer–Meshkov instability and mixing induced by reshock at different times

Turbulent mixing generated by shock-driven acceleration of a perturbed interface is simulated using a new multi-component Reynolds-averaged Navier–Stokes (RANS) model closed with a two-equation $K$ – $\epsilon $ model. The model is implemented in a hydrodynamics code using a third-order weighted essentially non-oscillatory finite-difference method for the advection terms and a second-order central difference method for the gradients in the source and diffusion terms. In the present reshocked Richtmyer–Meshkov instability and mixing study, an incident shock with Mach number $M\!a_{\mathrm{s}}=1.20$ is generated in air and progresses into a sulfur hexafluoride test section. The time evolution of the predicted mixing layer widths corresponding to six shock tube test section lengths are compared with experimental measurements and three-dimensional multi-mode numerical simulations. The mixing layer widths are also compared with the analytical self-similar power-law solution of the simplified model equations prior to reshock. A set of model coefficients and initial conditions specific to these six experiments is established, for which the widths before and after reshock agree very well with experimental and numerical simulation data. A second set of general coefficients that accommodates a broader range of incident shock Mach numbers, Atwood numbers, and test section lengths is also established by incorporating additional experimental data for $M\!a_{\mathrm{s}}=1.24$ , $1.50$ , and $1.98$ with $At=0.67$ and $M\!a_{\mathrm{s}}=1.45$ with $At=-0.67$ and previous RANS modeling. Terms in the budgets of the turbulent kinetic energy and dissipation rate equations are examined to evaluate the relative importance of turbulence production, dissipation and diffusion mechanisms during mixing. Convergence results for the mixing layer widths, mean fields, and turbulent fields under grid refinement are presented for each of the $M\!a_{\mathrm{s}}=1.20$ cases.     

Measurements of Turbulent Jet Mixing in a Turbulent Co-Flow Including the Influence of Periodic Forcing and Heating

In this work, the turbulent mixing of a confined coaxial jet in air is investigated by means of simultaneous particle image velocimetry and planar laser induced fluorescence of the acetone seeded flow injection. The jet is injected into a turbulent duct flow at atmospheric pressure through a 90 ^° pipe bend. Measurements are conducted in a small scale windtunnel at constant mass flow rates and three modes of operation: isothermal steady jet injection at a Dean number of 20000 (R e _d =32000), pulsed isothermal injection at a Womersley number of 65 and steady injection at elevated jet temperatures of ΔT=50 K and ΔT=100 K. The experiment is aimed at providing statistically converged quantities of velocity, mass fraction, turbulent fluctuations and turbulent mass flux at several downstream locations. Stochastic error convergence over the number of samples is assessed within the outer turbulent shear layer. From 3000 samples the statistical error of time-averaged velocity and mass fraction is below 1 % while the error of Reynolds shear stress and turbulent mass flux components is in the of range 5-6 %. Profiles of axial velocity and turbulence intensity immediately downstream of the bend exit are in good agreement with hot-wire measurements from literature. During pulsed jet injection strong asymmetric growing of shear layer vortices lead to a skewed mass fraction profile in comparison with steady injection. Phase averaging of single shot PLIF-PIV measurements allows to track the asymmetric shear layer vortex evolvement and flow breakdown during a pulsation cycle with a resolution of 10^°. Steady injection with increased jet temperature supports mixing downstream from 6 nozzle diameters onward.     

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

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.     

Fracture Propagation Driven by Fluid Outflow from a Low-Permeability Aquifer

Deep saline aquifers are promising geological reservoirs for $\mathrm {CO}_2$ CO 2 sequestration if they do not leak. The absence of leakage is provided by the caprock integrity. However, $\mathrm {CO}_2$ CO 2 injection operations may change the geomechanical stresses and cause fracturing of the caprock. We present a model for the propagation of a fracture in the caprock driven by the outflow of fluid from a low-permeability aquifer. We show that to describe the fracture propagation, it is necessary to solve the pressure diffusion problem in the aquifer. We solve the problem numerically for the two-dimensional domain and show that, after a relatively short time, the solution is close to that of one-dimensional problem, which can be solved analytically. We use the relations derived in the hydraulic fracture literature to relate the width of the fracture to its length and the flux into it, which allows us to obtain an analytical expression for the fracture length as a function of time. Using these results we predict the propagation of a hypothetical fracture at the In Salah $\mathrm {CO}_2$ CO 2 injection site to be as fast as a typical hydraulic fracture. We also show that the hydrostatic and geostatic effects cause the increase of the driving force for the fracture propagation and, therefore, our solution serves as an estimate from below. Numerical estimates show that if a fracture appears, it is likely that it will become a pathway for $\mathrm {CO}_2$ CO 2 leakage.     

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.