Nonlinear Pressure Diffusion in Flow of Compressible Liquids Through Porous Media

In the flow of liquids through porous media, nonlinear effects arise from the dependence of the fluid density, porosity, and permeability on pore pressure, which are commonly approximated by simple exponential functions. The resulting flow equation contains a squared gradient term and an exponential dependence of the hydraulic diffusivity on pressure. In the limiting case where the porosity and permeability moduli are comparable, the diffusivity is constant, and the squared gradient term can be removed by introducing a new variable y, depending exponentially on pressure. The published transformations that have been used for this purpose are shown to be special cases of the Cole–Hopf transformation, differing in the choice of integration constants. Application of Laplace transformation to the linear diffusion equation satisfied by y is considered, with particular reference to the effects of the transformation on the boundary conditions. The minimum fluid compressibilities at which nonlinear effects become significant are determined for steady flow between parallel planes and cylinders at constant pressure. Calculations show that the liquid densities obtained from the simple compressibility equation of state agree to within 1% with those obtained from the highly accurate Wagner-Pru?  equation of state at pressures to 20 MPa and temperatures approaching 600 K, suggesting possible applications to some geothermal systems.     

Stabilization and evolution of parameters of a symmetric laminar flow in a plane channel with sudden expansion

A problem is formulated for computing the fields of parameters of a stationary laminar symmetric flow. A two-dimensional flow in a channel with a sudden change in the cross-sectional area is computed. The evolution of a three-dimensional perturbation inserted into the channel at the initial stage of computations is analyzed. It is demonstrated that the parameters of a two-dimensional flow in the channel at a Reynolds number Re = 50 become stabilized at a dimensionless time t > 20, whereas the steady state is reached under the same conditions at t ≈ 100. At a distance of approximately 10h (h is the channel width at the entrance), the flow becomes one-dimensional, but the streamwise component of the velocity vector remains a function of the streamwise coordinate owing to flow compressibility. __________ Translated from PrikladnayaMekhanika i Tekhnicheskaya Fizika, Vol. 48, No. 1, pp. 35–42, January–February, 2007.     

Microstructural effects on the permeability of periodic fibrous porous media

An analytical–numerical approach is presented for computing the macroscopic permeability of fibrous porous media taking into account their microstructure. A finite element (FE) based model for viscous, incompressible flow through a regular array of cylinders/fibers is employed for predicting the permeability associated with this type of media. High resolution data, obtained from our simulations, are utilized for validating the commonly used semi-analytical models of drag relations from which the permeability is often derived. The effect of porosity, or volume fraction, on the macroscopic permeability is studied. Also microstructure parameters like particle shape, orientation and unit cell staggered angle are varied. The results are compared with the Carman–Kozeny (CK) equation and the Kozeny factor (often assumed to be constant) dependence on the microstructural parameters is reported and used as an attempt to predict a closed form relation for the permeability in a variety of structures, shapes and wide range of porosities.     

Numerical simulation of turbulent free-surface flow in curved channel

This paper presents first results of numerical simulation of turbulent free-surface flow. Simple implementation of surface capturing method is based on the variable density approach. The flow is treated as if there is only one fluid, but with variable material properties (density, viscosity). The switch in these values is done by a function resulting from the mass conservation principle. This approach simplifies the implementation of turbulence model. In this case the SST k−ω model was chosen in modification given by Hellsten.Numerical solution was carried out by finite-volume method with explicit Runge-Kutta time-integration. The artificial compressibility method was used for time-marching search for steady state solution. The whole model was tested on horizontally placed square-sectioned 90^∘ bend, which was partially filled by the water. The main goal of this study was to demonstrate the applicability of this model and solution method for capturing the water-air interface as well as for predicting the turbulent effects in both fluids.     

Assessment of morphology in transient and steady state regimes resulting from phase separation in polystyrene/poly(vinyl methyl ether) blend

Morphology development after phase separation in polystyrene (PS)/poly(vinyl methyl ether) (PVME) blend was assessed both in transient and in steady state regimes. Phase segregation was evaluated by various techniques including optical microscopy, light transmission, dynamic scanning calorimetry and rheological analyses. The steady state particle size resulted from phase separation was determined experimentally and then compared to the predictions of both the emulsion models that assume zero-thickness interfacial boundaries and to the asymptotic value of the Cahn–Hilliard theory that assumes rather a diffuse interphase.     

The Effect of Void Structure on the Permeability of Fibrous Networks

A Kozeny–Carman-based model of permeability for fibrous networks is proposed: the original model is extended by incorporating information about the local structure of the void space. Furthermore, it is demonstrated how in practice this added structural information can be retrieved from a three-dimensional digital image of a fibrous material. The proposed model is then validated for both foam- and water-deposited laboratory sheets of bleached kraft pulp (Scots pine) and chemi-thermo-mechanical pulp (CTMP, Norway spruce). The validation is carried out by comparing the model predictions against computationally determined permeability values. The related fluid-flow simulations are executed using the lattice-Boltzmann method together with high-resolution X-ray microtomography images. For both pulp materials, the sample sets had nearly equal porosities, but deviated substantially in their permeabilities. The proposed model was shown to improve prediction of permeability for the fibrous materials considered: the deviation between the predicted and computationally determined values was no more than 8%.     

A turbulent plane offset jet with small offset ratio

 Mean velocities and turbulence characteristics of a turbulent plane offset jet with a small offset ratio of 2.125 have been studied using laser Doppler anemometry (LDA). Static pressure measurements highlight the importance of side plates in enhancing two-dimensionality of the jet. The spatial distributions of turbulence intensities and Reynolds shear stress show a high turbulence recirculating flow region close to the nozzle plate between the jet and the offset plate. The LDA results have been used to examine the capability of three different turbulence models (i.e. k–ɛ, RNG and Reynolds stress) in predicting the velocity field of this jet. While all three models are able to predict qualitatively the recirculation, converging and reattachment regions observed experimentally, the standard k–ɛ turbulence model predicts a reattachment length that best agrees with the experimentally determined value. Received: 11 September 1996/Accepted: 30 May 1997     

Marangoni Mixed Convection Boundary Layer Flow

The paper deals with a steady coupled dissipative layer, called Marangoni mixed convection boundary layer, which can be formed along the interface of two immiscible fluids, in surface driven flows. The mixed convection boundary layer is generated when besides the Marangoni effects there are also buoyancy effects due to gravity and external pressure gradient effects. We shall use a model proposed by Golia and Viviani (L’ Aerotecnica missili e Spazio 64 (1985) 29–35, Meccanica 21 (1986) 200–204) wherein the Marangoni coupling condition has been included into the boundary conditions at the interface. The similarity equations are first determined, and the pertinent equations are solved numerically for some values of the governing parameters and the features of the flow and temperature fields as well as the interface velocity and heat transfer at the interface are analysed and discussed.     

Structure and rheology of heat-set gels of globular proteins

The structure and the rheology of systems resulting from heating at 80°C isoelectric solutions of bovine serum albumin (BSA) in the concentration range 10–200mg/ml were studied.  Small-angle neutron scattering measurements view the systems as being formed of large aggregates of micrometric size with a close packed arrangement of denatured protein molecules. No indication of a fractal structure stands out.  The viscoelastic behaviour is linear up to about 5% strain, except in the BSA concentration range 30–90mg/ml where the linearity limit is below 1% strain. The viscoelastic response was analysed in the linear domain, or as close as possible to it, by combining the results of dynamic and creep recovery measurements. The dependence on concentration of the steady state viscosity, of the steady state compliance, and of the average retardation time shows a marked change around a concentration C ₀∼50mg/ml, corresponding probably to a percolation threshold. Received: 10 March 1998 Accepted: 1 May 1998     

Modeling of Breakthrough Curves for Adsorption of Propane, <Emphasis Type="Italic">n</Emphasis>-Butane,and Iso-Butane Mixture on 5A Molecular Sieve Zeolite

Breakthrough curves for the adsorption of propane, n-butane, and iso-butane mixture on 5A molecular sieve zeolite were obtained experimentally and theoretically at a constant temperature of 301 K. The equilibrium model and linear driving force model were used to predict the experimental breakthrough curves for this multicomponent mixture. The equilibrium model gave a satisfactory fit for experimental data. The model equations were solved by a numerical method based on backward finite difference with a fixed griding technique. The effects of feed flow rate (0.552–3.496 l/min), feed concentration (60.72–141.68 mmol/l), and adsorbates composition (58.75–75.32%) on these breakthrough curves were examined.     

Flow near the permeable boundary of a porous medium: An experimental investigation using LDA

 Fluid flow at the interface of a porous medium and an open channel is the governing phenomenon in a number of processes of industrial importance. Traditionally, this has been modeled by applying the Brinkman’s modification of Darcy’s law to obtain the velocity profile in terms of an additional parameter known as the “apparent viscosity” or the “slip coefficient”. To test this ad hoc approach, a detailed experimental investigation of the flow was conducted using Laser Doppler Anemometry (LDA) in the close vicinity of the permeable boundary of a porous medium. The porous medium used in the experiments consisted of a network of continuous glass strands woven together in a random fashion. A Hele–Shaw cell was partially filled with a fibrous preform such that an open channel flow is coupled with the Darcy flow inside the preform through the permeable interface of the preform. The open channel portion of the Hele–Shaw cell also acts as an ideal porous medium of known in-plane permeability which is much higher than the permeability of the fibrous porous medium. A viscous fluid is injected at a constant flow rate through the above arrangement and a saturated and steady flow is established through the cell. Using LDA, steady state velocity profiles are accurately measured by traversing across the cell in the direction perpendicular to the flow. A series of experiments were conducted in which fluid viscosity, flow rate, solid volume fraction of the porous medium and depth of the Hele–Shaw cell were varied. For each and every case in which the conditions for Hele–Shaw approximation were valid, the depth of the boundary layer zone or the screening length inside the fibrous preform was found to be of the order of the channel depth. This is much larger as compared to the Brinkman’s prediction of the screening length which is of the order of √K, where K is the permeability of the fibrous porous medium. Based on this finding, we modified the boundary condition in the Brinkman’s solution and found that the velocity profile results compared well with the experimental data for the planar geometry and the fibrous preforms for volume fractions of 7%, 14% and 21% for Hele–Shaw cell depths of 1.6 and 3.175 mm. For a cell depth of 4.8 cm, in which the Hele–Shaw approximation was not valid, the boundary layer thickness or the screening length was found to be less than the mold or channel depth but was still much larger than the Brinkman’s prediction. Received: 10 May 1996 / Accepted: 26 August 1996     

Pore-Scale Modeling of Electrokinetics in Geomaterials

Pore-scale finite-volume continuum models of electrokinetic processes are used to predict the Debye lengths, velocity, and potential profiles for two-dimensional arrays of circles, ellipses and squares with different orientations. The pore-scale continuum model solves the coupled Navier–Stokes, Poisson, and Nernst–Planck equations to characterize the electro-osmotic pressure and streaming potentials developed on the application of an external voltage and pressure difference, respectively. This model is used to predict the macroscale permeabilities of geomaterials via the widely used Carmen–Kozeny equation and through the electrokinetic coupling coefficients. The permeability results for a two-dimensional X-ray tomography-derived sand microstructure are within the same order of magnitude as the experimentally calculated values. The effect of the particle aspect ratio and orientation on the electrokinetic coupling coefficients and subsequently the electrical and hydraulic tortuosity of the porous media has been determined. These calculations suggest a highly tortuous geomaterial can be efficient for applications like decontamination and desalination.

     

Dynamic method of determining aerodynamic characteristics of models from on experimental results obtained in short-duration wind tunnels

A technique is developed to recover the loads acting in model testing in short-duration wind tunnels. The model together with the load-measuring device is assumed to be a linear dynamic system. Normal responses of the system to unit loads are determined experimentally. The solution of a system of integral Volterra equations of the first kind is found in the class of piecewise-constant functions. Examples are given on solving a test problem with a known exact solution and a problem of determining the aerodynamic characteristics of the HB-2 reference model on the basis of loads obtained in wind-tunnel measurements. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 47, No. 5, pp. 47–55, September–October, 2006.     

Effects of Viscous Dissipation on Unsteady Free Convection in a Fluid Past a Vertical Plate Immersed in a Porous Medium

The effects of viscous dissipation on unsteady free convection from an isothermal vertical flat plate in a fluid saturated porous medium are examined numerically. The Darcy–Brinkman–Forchheimer model is employed to describe the flow field. A new model of viscous dissipation is used for the Darcy–Brinkman–Forchheimer model of porous media. The simultaneous development of the momentum and thermal boundary layers are obtained by using a finite difference method. Boundary layer and Boussinesq approximation have been incorporated. Numerical calculations are carried out for various parameters entering into the problem. Velocity and temperature profiles as well as local friction factor and local Nusselt number are shown graphically. It is found that as time approaches infinity, the values of friction factor and heat transfer coefficient approach steady state.     

Rheological study of segmented polyurethane and polycaprolactone blends

The dynamic rheological properties of segmented polyurethane based on polycaprolactone diol (PU-PCL diol) and poly (ɛ-caprolactone) (PCL) blends were experimentally examined and theoretically analyzed using Palierne model. PU–PCL diol was melt blended with PCL in three different ratios of 20, 30, and 40%. Initial modeling attempts failed to fit the experimental data of these blends, as the model overpredicted their rheological data. This failure is believed to be due to partial dissolution of PCL in PU. According to our obtained results from differential scanning calorimetry and dynamic thermal mechanical analysis, pure PU–PCL diol is diluted by its homopolymer. By calculating the complex modulus of this diluted matrix using a simple mixing rule and its incorporation into the model, good fitting of the Palierne equation to the experimental data was obtained.     

Influence of elastic material compressibility on parameters of an expanding spherical stress wave

We investigated the influence of elastic material compressibility on parameters of an expanding spherical stress wave. The material compressibility is represented by Poisson’s ratio, ν, in this paper. The stress wave is generated by a pressure produced inside a spherical cavity surrounded by the isotropic elastic material. The analytical closed form formulae determining the dynamic state of the mechanical parameters (displacement, particle velocity, strains, stresses, and material density) in the material have been derived. These formulae were obtained for surge pressure p(t) = p ₀ = const inside the cavity. From analysis of these formulae, it is shown that the Poisson’s ratio substantially influences the course of material parameters in space and time. All parameters intensively decrease in space together with an increase of the Lagrangian coordinate, r. On the contrary, these parameters oscillate versus time around their static values. These oscillations decay in the course of time. We can mark out two ranges of parameter ν values in which vibrations of the parameters are “damped” at a different rate. Thus, Poisson’s ratio in the range below about 0.4 causes intense decay of parameter oscillations. On the other hand in the range 0.4 < ν < 0.5, i.e. in quasi-incompressible materials, the “damping” of parameter vibrations is very low. In the limiting case when ν = 0.5, i.e. in the incompressible material, “damping” vanishes, and the parameters harmonically oscillate around their static values. The abnormal behaviour of the material occurs in the range 0.4 < ν < 0.5. In this case, an insignificant increase of Poisson’s ratio causes a considerable increase of the parameter vibration amplitude and decrease of vibration “damping”.      

Mechanism of compressive stress formation during weak shock waves impact with granular materials

 Certain aspects of wave propagation and the dynamic reaction of a granular material when subjected to a long-duration impulse load are studied. In the majority of studies published on this subject the unsteady pressure behavior at the end-wall covered by a layer of granular material was observed and documented. However, up to now little attention was given to explaining the physical mechanism of this process. Experimental results, obtained in the course of this study, regarding the pressure fields inside granular layers of different materials, clearly show that the compaction effect strongly depends on the characteristics of the medium. This phenomenon manifests itself by changing the gas-particle interaction in the course of the gas filtration, and by variation in the contribution of the different forces and effective stress, σ, to the energy exchange between the gas, the particles and the shock-tube wall. The material permeability,  f, the relative density, ν, and the particle response time, τ _p , are the most important parameters affecting the stress formation at the end-wall covered by the granular layer. In addition to the effect of the material parameters, the effective stress, σ, was found to strongly depend on the granular layer height, h. Based on detailed pressure measurements a qualitative analysis regarding the role of the particle rearrangement in the formation of the unsteady peak at the end-wall was performed. The phenomenology of the particle–particle interaction includes rotation and consolidation of the granules and movement or sliding of particle planes within the layer over each other. Most of these processes are frictional in their nature. They are related to the energy losses and affect the profile and magnitude of the compressive stress as measured at the shock-tube end-wall covered by the granular layer. Received: 10 June 1996/Accepted: 15 October 1996     

An empirical constitutive law for concentrated colloidal suspensions in the approach of the glass transition

Concentrated, non-crystallizing colloidal suspensions in their approach of the glass state exhibit distinct dynamics patterns. These patterns suggest a powerlaw rheological constitutive model for near-glass viscoelasticity, as presented here. The rheological parameters used for this model originate in the mode-coupling theory. The proposed constitutive model provides explicit expressions for the steady shear viscosity, the steady normal stress coefficient, the modulus-compliance relation, and the α peak of G″. The relaxation pattern distinctly differs from gelation.     

Transient analysis of a tennis racket using PC-based finite elements and experimental techniques

A method has been developed to facilitate the fully stressed design optimization of a tennis racket. The method consists of a PC-based finite-element model with experimental verification, and transient analysis using experimentally determined dynamic loading data. Results are obtained and discussed. Paper was presented at the 1988 SEM Spring Conference held in Portland, OR on June 5–10.     

Insights into the Relationships Among Capillary Pressure,Saturation, Interfacial Area and Relative Permeability Using Pore-Network Modeling

To gain insight in relationships among capillary pressure, interfacial area, saturation, and relative permeability in two-phase flow in porous media, we have developed two types of pore-network models. The first one, called tube model, has only one element type, namely pore throats. The second one is a sphere-and-tube model with both pore bodies and pore throats. We have shown that the two models produce distinctly different curves for capillary pressure and relative permeability. In particular, we find that the tube model cannot reproduce hysteresis. We have investigated some basic issues such as effect of network size, network dimension, and different trapping assumptions in the two networks. We have also obtained curves of fluid–fluid interfacial area versus saturation. We show that the trend of relationship between interfacial area and saturation is largely influenced by trapping assumptions. Through simulating primary and scanning drainage and imbibition cycles, we have generated two surfaces fitted to capillary pressure, saturation, and interfacial area (P _c –S _w –a _nw ) points as well as to relative permeability, saturation, and interfacial area (k _r –S _w –a _nw ) points. The two fitted three-dimensional surfaces show very good correlation with the data points. We have fitted two different surfaces to P _c –S _w –a _nw points for drainage and imbibition separately. The two surfaces do not completely coincide. But, their mean absolute difference decreases with increasing overlap in the statistical distributions of pore bodies and pore throats. We have shown that interfacial area can be considered as an essential variable for diminishing or eliminating the hysteresis observed in capillary pressure–saturation (P _c –S _w ) and the relative permeability–saturation (k _r –S _w ) curves.