Fractal geometry model for effective thermal conductivity of three-phase porous media

An approximate model, a fractal geometry model, for the effective thermal conductivity of three-phase/unsaturated porous media is proposed based on the thermal-electrical analogy technique and on statistical self-similarity of porous media. The proposed thermal conductivity model is expressed as a function of porosity (related to stage n of Sierpinski carpet), ratio of areas, ratio of component thermal conductivities, and saturation. The recursive algorithm for the thermal conductivity by the proposed model is presented and found to be quite simple. The model predictions are compared with the existing measurements. Good agreement is found between the present model predictions and the existing experimental data. This verifies the validity of the proposed model. (C) 2004 American Institute of Physics.

Fractal dimensions for unsaturated porous media

The analytical expressions of the fractal dimensions for wetting and non-wetting phases for unsaturated porous media are derived and are found to be a function of porosity, maximum and minimum pore sizes as well as saturation. There is no empirical constant in the proposed fractal dimensions. It is also found that the fractal dimensions increase with porosity of a medium and are meaningful only in a certain range of saturation S-w, i.e. S-w > S-min for wetting phase and S-w < S-max for non-wetting phase at a given porosity, based on real porous media for requirements from both fractal theory and experimental observations. The present analysis of the fractal dimensions is verified to be consistent with the existing experimental observations and it makes possible to analyze the transport properties such as permeability, thermal dispersion in unsaturated porous media by fractal theory and technique.

Wave propagation and the frequency domain Green's functions in viscoelastic Biot/squirt (BISQ) media

In this paper, we examine the characteristics of elastic wave propagation in viscoelastic porous media, which contain simultaneously both the Biot-flow and the squirt-flow mechanisms (BISQ). The frequency-domain Green's functions for viscoelastic BISQ media are then derived based on the classic potential function methods. Our numerical results show that S-waves are only affected by viscoelasticity, but not by squirt-flows. However, the phase velocity and attenuation of fast P-waves are seriously influenced by both viscoelasticity and squirt-flows; and there exist two peaks in the attenuation-frequency variations of fast P-waves. In the low-frequency range, the squirt-flow characteristic length, not viscoelasticity, affects the phase velocity of slow P-waves, whereas it is opposite in the high-frequency range. As to the contribution of potential functions of two types of compressional waves to the Green's function, the squirt-flow length has a small effect, and the effects of viscoelastic parameter are mainly in the higher frequency range. Crown Copyright (C) 2006 Published by Elsevier Ltd. All rights reserved.

Stress-distribution around a rigid line in dissimilar media

The elastic plane problem of a rigid line inclusion between two dissimilar media was considered. By solving the Riemann-Hilbert problem, the closed-form solution was obtained and the stress distribution around the rigid line was investigated. It was found that the modulus of the singular behavior of the stress remains proportional to the inverse square root of the distance from the rigid line end, but the stresses possess a pronounced oscillatory character as in the case of an interfacial crack tip.

Effects of contact resistance on thermal-conductivity of composite media with a periodic structure

The thermal conductivity of periodic composite media with spherical or cylindrical inclusions embedded in a homogeneous matrix is discussed. Using Green functions, we show that the Rayleigh identity can be generalized to deal with thermal properties ot these systems. A new calculating method for effective conductivity of composite media is proposed. Useful formulae for effective thermal conductivity are derived, and meanings of contact resistance in engineering problems are explained.

An asymptotic analysis for one dimensional stress wave propagation in damaged viscoelastic media

A regular perturbation technique is suggested to deal with the problem of one dimensional stress wave propagation in viscoelastic media with damage. Based upon the first order asymptotic solution obtained, the characteristics of wave attenuation are studied. In fact, there exist three different time-dependent phenomena featuring the dynamic response of the materials, the first expressing the characteristics of wave propagation, the second indicating the innate effect of visco-elastic matrix and the third coming from the time dependent damage. The comparision of first order asymptotic solution with the numerical results calculated by a finite difference procedure shows that the perturbation expansion technique may offer a useful approach to the problem concerned.

Thermal-conductivity of periodic composite media with spherical inclusions

The thermal conductivity of periodic composite media with spherical inclusions embedded in a homogeneous matrix is discussed. Using Green's function, we show that the Rayleigh identity can be generalized to deal with the thermal properties of these systems. A technique for calculating effective thermal conductivities is proposed. Systems with cubic symmetries (including simple cubic, body centered cubic and face centered cubic symmetry) are investigated in detail, and useful formulae for evaluating effective thermal conductivities are derived.

A general theory on media with randomly distributed inclusions .1. The average field behaviors

In this paper, a theory is developed to calculate the average strain field in the materials with randomly distributed inclusions. Many previous researches investigating the average field behaviors were based upon Mori and Tanaka's idea. Since they were restricted to studying those materials with uniform distributions of inclusions they did not need detailed statistical information of random microstructures, and could use the volume average to replace the ensemble average. To study more general materials with randomly distributed inclusions, the number density function is introduced in formulating the average field equation in this research. Both uniform and nonuniform distributions of inclusions are taken into account in detail.

The asymptotic near-tip solution for mode-iii crack in steady growth in power hardening media

Two local solutions, one perpendicular and one parallel to the direction of solar gravitational field, are discussed. The influence of gravity on the gas-dynamical process driven by the piston is discussed in terms of characteristic theory, and the flow field is given quantitatively. For a typical piston trajectory similar to the one for an eruptive prominence, the velocity of the shock front which locates ahead the transient front is nearly constant or slightly accelerated, and the width of the compressed flow region may be kept nearly constant or increased linearly, depending on the velocity distribution of the piston. Based on these results, the major features of the transient may be explained. Some of the fine structure of the transient is also shown, which may be compared in detail with observations.

Transient Flows in Low Permeability Porous Media with Composite Start-Up Pressure Gradient

A new mathematical model for the transient flow in the composite low permeability is established. It is solved by FEM with different boundary conditions such as infinite, circular closed and constant pressure boundary conditions. The typical curves for transient wellbore pressure have been presented. It is shown that the pressure and pressure derivative curves with composite start-up pressure gradients have different slopes which are depended on the start-up pressure gradients and the mobility radios in different regions. The boundary effects are the same as the normal reservoirs without start-up pressure gradients. The study provides a new tool to analyze the transient pressure test data in the low permeability reservoir.

Propagation of Explosion Wave in Hard-Soft-Hard Layered Media

The paper presents a reasonable analysis for dynamic response and failure process of a plane multi-layered media, which are subjected to a blast loading. This blast loading is induced by a cylindric explosive put on the center of top surface of the layered media. With the help of numerical simulation technique provided by LS-DYNA software, the whole process of explosion wave propagation and attenuation can be revealed. The feature of local failure around the blasting site is also discussed in some detail. Our focus will be on the explosion wave attenuation for the hard-soft-hard sandwich layers. As seen in the paper, the computational results are delivered in a feasible way by comparing with experimental data.

Numerical Sensitivity Analysis and Partial Similarity of Porous Media Flows

A numerical optimisation approach to identify dominant dimensionless variables in porous media flows by sensitivity analysis is proposed. We have validated the approach at first by examining a simple oil reservoir theoretically and numerically as well. A more complex water-flooding reservoir is examined based on sensitivity analysis of oil recovery to the similarity parameters, thus demonstrating the feasibility of the proposed approach to identify dominant similarity parameters for water-oil two-phase flows.

Catastrophic Rupture Induced Damage Coalescence in Heterogeneous Brittle Media

In heterogeneous brittle media, the evolution of damage is strongly influenced by the multiscale coupling effect. To better understand this effect, we perform a detailed investigation of the damage evolution, with particular attention focused on the catastrophe transition. We use an adaptive multiscale finite-element model (MFEM) to simulate the damage evolution and the catastrophic failure of heterogeneous brittle media. Both plane stress and plane strain cases are investigated for a heterogeneous medium whose initial shear strength follows the Weibull distribution. Damage is induced through the application of the Coulomb failure criterion to each element, and the element mesh is refined where the failure criterion is met. We found that as damage accumulates, there is a stronger and stronger nonlinear increase in stress and the stress redistribution distance. The coupling of the dynamic stress redistribution and the heterogeneity at different scales result in an inverse cascade of damage cluster size, which represents rapid coalescence of damage at the catastrophe transition.

Comparisons of static, quasi-static and dynamic 3D porous media scale network models for two-phase immiscible flow in porous media

A dynamic 3D pore-scale network model is formulated for investigating the effect of interfacial tension and oil-water viscosity during chemical flooding. The model takes into account both viscous and capillary forces in analyzing the impact of chemical properties on flow behavior or displacement configuration, while the static model with conventional invasion percolation algorithm incorporates the capillary pressure only. From comparisons of simulation results from these models. it indicates that the static pore scale network model can be used successfully when the capillary number is low. With the capillary increases due to the enhancement of water viscosity or decrease of interfacial tension, only the quasi-static and dynamic model can give insight into the displacement mechanisms.