Numerical Simulation of the Physical–Chemical–Thermal Processes During Hydration Reaction of the Calcium Oxide/Calcium Hydroxide System in an Indirect Reactor

Transport in Porous Media - Solutal convection in a horizontal layer filled with porous media with horizontal seepage of a mixture is investigated considering the solute immobilization and...     

Thermal Dispersion–Radiation Effects on Non-Darcy Natural Convection in a Fluid Saturated Porous Medium

The effects of thermal dispersion and thermal radiation on the non-Darcy natural convection over a vertical flat plate in a fluid saturated porous medium are studied. Forchheimer extension is considered in the flow equations. The coefficient of thermal diffusivity has been assumed to be the sum of molecular diffusivity and the dispersion thermal diffusivity due to mechanical dispersion. Rosseland approximation is used to describe the radiative heat flux in the energy equation. Similarity solution for the transformed governing equations is obtained. Numerical results for the details of the velocity and temperature profiles which are shown on graphs have been presented. The combined effect of thermal dispersion and thermal radiation, for the two cases Darcy and non-Darcy porous medium, on the heat transfer rate which are entered in tables is discussed.     

Non-linear Evolution of P-waves in Viscous–Elastic Granular Saturated Media

The Frenkel–Biot P-wave of the first type is a seismic longitudinal wave observed in rocks fully saturated with oil, water or high-pressure gas. The P-wave of the second type is observed in unsaturated soils and other porous media saturated with gas of low pressure. Their models include properties of the skeleton, that is, its elastic modules and its own viscosity. If the non-linear terms are accounted for, the asymptotic analysis, usual for weak non-linear waves, might be applied to get the wave spectrum evolution. The wetness of grains contacts in soils and such components of oil as tars or bitumen, which attached to the skeleton, can be described by generalized viscous–elastic stress–strain connections. The latter are nominated in such a way that creates the narrow frequency interval of wave of negative dissipation where the non-linear terms begin to play the main role besides the neutral stability for waves of zero wave number. The corresponding case, relevant to single continuum model, was analyzed in the literature. Here it is shown that the interpenetrating continua with interaction of the Darcy type provide the dissipation sink in the wave evolution equation. This generalization, (Tribelsky, M.I.: Phys. Rev. Lett. (2007, submitted)), can stabilize the asymptotic solution of the evolution equation, where the dispersion terms are omitted. The asymptotic solution of the equation is invariant to initial conditions and it means a transformation of initial wave spectra to unique one while wave is spreading in the viscous–elastic medium under consideration. This explains the phenomenon, observed in wave tests at marine beach, when any dynamics action (impact, explosion, and ultrasound action) created at some distance a wave of a single frequency (~25 Hz).     

Thermal–mechanical coupled analysis of a brake disk rotor

The main purpose of this study is to analyze the thermomechanical behavior of the dry contact between the brake disk and pads during the braking phase. The simulation strategy is based on computer code ANSYS11. The modeling of transient temperature in the disk is actually used to identify the factor of geometric design of the disk to install the ventilation system in vehicles The thermal–structural analysis is then used with coupling to determine the deformation and the Von Mises stress established in the disk, the contact pressure distribution in pads. The results are satisfactory when compared to those of the specialized literature.     

The Onset of Darcy–Brinkman Electroconvection in a Dielectric Fluid Saturated Porous Layer

The combined effect of a vertical AC electric field and the boundaries on the onset of Darcy–Brinkman convection in a dielectric fluid saturated porous layer heated either from below or above is investigated using linear stability theory. The isothermal bounding surfaces of the porous layer are considered to be either rigid or free. It is established that the principle of exchange of stability is valid irrespective of the nature of velocity boundary conditions. The eigenvalue problem is solved exactly for free–free (F/F) boundaries and numerically using the Galerkin technique for rigid–rigid (R/R) and lower-rigid and upper-free (F/R) boundaries. It is observed that all the boundaries exhibit qualitatively similar results. The presence of electric field is emphasized on the stability of the system and it is shown that increasing the AC electric Rayleigh number R _ea is to facilitate the transfer of heat more effectively and to hasten the onset of Darcy–Brinkman convection. Whereas, increase in the ratio of viscosities Λ and the inverse Darcy number Da ⁻¹ is to delay the onset of Darcy–Brinkman electroconvection. Besides, increasing R _ea and Da ⁻¹ as well as decreasing Λ are to reduce the size of convection cells.     

Effect of Conduction in Bottom Wall on B��nard Convection in a Porous Enclosure with Localized Heating and Lateral Cooling

Darcy-Bénard convection in a square porous enclosure with a localized heating from below and lateral cooling is studied numerically in the present paper. A finite-thickness bottom wall is locally heated, the top wall is kept at a lower temperature than the bottom wall temperature, and the lateral walls are cooled. The finite difference method has been used to solve the dimensionless governing equations. The analysis in the undergoing numerical investigation is performed in the following ranges of the associated dimensionless groups: the heat source length?? ${0.2\leq H \leq 0.9}$ , the wall thickness?? ${0.05\leq D \leq 0.4}$ , the thermal conductivity ratio?? ${0.8\leq K_{\rm r} \leq 9.8}$ , and the Biot number?? ${0.1\leq Bi \leq 1.1}$ . It is observed that the heat transfer rate could increase with increasing heat source lengths, thermal conductivity ratio, and cooling intensity. There exists a critical wall thickness for a high wall conductivity below which the increasing wall thickness increases the heat transfer rate and above which the increasing wall thickness decreases the heat transfer rate.     

Role of pH in Electro-Osmosis: Experimental Study on NaCl–Water Saturated Kaolinite

Electro-osmosis experiments were conducted on rigid cylindrical samples containing 0.01 M NaCl–water saturated Speswhite kaolinite. It is experimentally found that the electro-osmotic permeability is pH-dependent. It is also experimentally found that time and spatial variations of the sample pH and of the pore water pressure correlate. This is qualitatively confirmed by a simple analysis that couples the electro-osmotic and hydraulic flows through the pH-dependent electro-osmotic permeability. However quantitative agreement between the experimental and numerical values of the pore water pressure is not obtained throughout the whole sample. This suggests that the hydraulic permeability may also depend on the pH.     

Multi–Front Phase Transitions During Nonisothermal Filtration of Live Paraffin–Base Crude

By numerical modeling of nonisothermal filtration of multicomponent oil with allowance for the Joule—Thomson effect, adiabatic effect, and heat of phase transitions that occur during oil degassing and paraffin crystallization, the formation of profiles of phase–saturation, concentrations of oil components, and temperature in oil beds is studied. It is shown that consideration of many components results in occurrence of phase–transition fronts during degassing of oil components and paraffin crystallization. In turn, paraffin crystallization gives rise to temperature oscillations. Depending on the initial paraffin concentration and on the ratio of phase–transition heats for oil degassing and paraffin crystallization in an oil bed, either decaying or steady–state temperature oscillations are observed.     

An Upscaling–Static-Downscaling Scheme for Simulation of Enhanced Oil Recovery Processes

We investigate whether upscaling errors for EOR simulation can be reduced by an upscaling–static-downscaling method where the scales of simulation for the pressure and saturation/concentration switch between coarse simulation model and fine geological model. We apply a static downscaling that has been previously shown to be reliable for water flooding. We use the same algorithm of static downscaling for EOR processes that have been used for water flooding. Different EOR processes are considered: polymer, surfactant and thermal. This range of flooding processes ensures that we are examining more physically complicated systems than water flooding. For these processes, one major difference from water flooding is existence of a secondary front. The effective capturing of this front is a criterion of accuracy for upscaling because, for this front, the coupling of dispersion with the fractional flow creates excessive smearing. A scheme for numerical dispersion control is implemented to both upscaled and downscaled models to determine and reduce the sensitivity to dispersion errors.     

Instability of Hadley–Prats Flow with Viscous Heating in a Horizontal Porous Layer

The onset of convective rolls instability in a horizontal porous layer subject to a basic temperature gradient inclined with respect to gravity is investigated. The basic velocity has a linear profile with a non-vanishing mass flow rate, i.e., it is the superposition of a Hadley-type flow and a uniform flow. The influence of the viscous heating contribution on the critical conditions for the onset of the instability is assessed. There are four governing parameters: a horizontal Rayleigh number and a vertical Rayleigh number defining the intensity of the inclined temperature gradient, a Péclet number associated with the basic horizontal flow rate, and a Gebhart number associated with the viscous dissipation effect. The critical wave number and the critical vertical Rayleigh number are evaluated for assigned values of the horizontal Rayleigh number, of the Péclet number, and of the Gebhart number. The linear stability analysis is performed with reference either to transverse or to longitudinal roll disturbances. It is shown that generally the longitudinal rolls represent the preferred mode of instability.     

A Benchmark Solution for Infiltration and Adsorption of Polluted Water Into Unsaturated–Saturated Porous Media

We discuss the numerical modeling of the infiltration of contaminated water into unsaturated porous media. A system with contaminant transport, dispersion, and adsorption is considered. The mathematical model for unsaturated flow is based on Richards nonlinear and degenerate equation. Nonlinear adsorption is represented by adsorption isotherms and kinetic rates. An accurate numerical method is constructed in 1D which can be a good candidate for the solution of inverse problems to determine model parameters in the adsorption part of the model. Our numerical solution is based on the method of lines (MOL method) where space discretization leads to the corresponding system of ODEs. We substantially use the numerical modeling of interfaces, separating fully saturated, partially saturated, and dry zones in the underground. Finally, in a series of numerical experiments and in comparisons with HYDRUS (?imunek et al., The HYDRUS-1D software package for simulating the one-dimensional movement of water, heat, and multiple solutes in variably/saturated media, version 2.0, Rep. IGWMC-TPS-70, 202 pp., Int. Groundwater Model. Cent., Colo. Sch of Mines, Golden, Colo), we demonstrate the effectiveness of our method.     

Impacts of Limestone Vertical Permeability Heterogeneity on Fluid–Rock Interaction During CCS

Transport in Porous Media - Within the context of subsurface CO2 storage, “live brine” (formed by the mixing of CO2 and brine) in limestone rock environments can cause mineral...     

Coordinated Variation of Contact Angles During Mobilization of Double Liquid–Gas Interfaces in a Microcapillary

Transport in Porous Media - Effectively mobilizing displacement and predicting mobilization pressure in a porous-type reservoir filled with bubbles or blobs require the knowledge of variation of...     

Kinetic Approach to Prediction of the Life of Aluminum Alloys under Various Thermal–Temporal Loading Conditions

Results of investigation of the life of D16 T, AK41 T1, and 1201 T1 aluminum alloys are generalized on the basis of the kinetic concept of failure. The life is studied under creep at constant loads and loads increasing with different rates and at different temperatures. The temperature is varied within the range of 473–77 K, and the duration of tests ranges from fractions of a second to ten thousand hours. Information on the effect of internalstress relaxation on the life of alloys is obtained. A method for predicting the life with allowance for relaxation processes in solids is verified experimentally.     

Simulation of Structure Change in Porous Media During Gas–Solid Reactions Using Cellular Automata Model

Transport in Porous Media - In some gas–solid reactions, a new solid substance is produced. The product acts as a shield and prevents the collision between gas and solid reactants which...     

Heating and Melting of Asphalt–Paraffin Plugs in Oil‐Well Equipment Using an Electromagnetic Radiation Source Operating in a Periodic Mode

This paper considers the elimination of asphalt–paraffin plugs in wellbore equipment using a highfrequency radiation source which is energized and deenergized periodically. The dynamics of melting of plugs is analyzed numerically. The time of removal of plugs is determined with variation in the offduty ratio of the operating cycle of the highfrequency oscillator and the time of its continuous operation in a cycle.     

A Highly Conductive Porous Medium for Solid–Gas Reactions: Effect of the Dispersed Phase on the Thermal Tortuosity

The heat transfer in a highly conductive material constituted by a graphite matrix in which a granular phase is dispersed is studied. The effective thermal conductivity of this anisotropic porous composite medium used in solid–gas reactors can vary largely with the component fractions. The effect of the dispersed grains on the deformable structure of the matrix is considered. A model developed on the basis of thermal tortuosity by analogy with mass transfer is adequately correlated with experimental results.