Effect of the Nature of the Interaction between a Gas and a Surface on the Relations for the Knudsen Layer in the Case of Strong Evaporation

The effect of the temperature accommodation coefficient _T on the relations at the Knudsen layer edge is investigated for strong evaporation using the moment method. An explicit expression for the dimensionless density as a function of the temperature and the Mach number M is obtained for 0 < _T < 1. For _T = 0 the entire solution is obtained in explicit form. It is shown that for = 0 and a condensation coefficient << 1 the temperature outside the Knudsen layer changes sharply as M varies from 0 to a certain value much less than unity after which the temperature ceases to depend on . For the model of specular reflection of the molecules from the surface the density and the temperature outside the Knudsen layer are found in explicit form as functions of the Mach number.     

Regular reflection of plane detonation waves from an elastic half-space

An effective method for the approximate solution of the Eq. [1] for the intensity of a reflected shock wave in the case of oblique incidence of a detonation wave on an elastic half-space is described; the elastic half-space is described by a certain specific form of the equation of state. Formulas relating the front and particle velocities behind the transmitted wave front to physical parameters are derived. Values of the wave intensity and other quantities determined with the aid of a Ural-2 computer are cited.The author of [1, 2] investigated the regular reflection of shock waves from the boundary between two bodies. In the present paper we solve the analogous problem in the case of oblique incidence of a detonation wave on an elastic half-space. The detonation wave deforms the elastic half-space, which assumes the position OK₁ (Fig. 1) forming the angle to the initial direction KO of the halfspace boundary. We assume that the acoustic stiffness of the halfspace is larger than the acoustic stiffness of the explosive. In this case, both reflected wave 2 and transmitted wave 3 are shock waves [3]. Let us denote the velocities of propagation of the detonation, reflected, and transmitted waves by U_i(i=1, 2, 3), respectively; let the pressure be p_i and let the density bep _i(i=0, 1, 2, 3, 4). The quantities U₁, ₁, ₀, and ₄ are given. We determine the intensities of waves 2 and 3, their velocities of propagation, and the angles ₂, ₃ and . The parameters are constant within each of the domains a, b, c, d, and e. In domains a and e the medium is stationary, i.e., u₀=u₄ =0. The basic equations of the problem express the conditions at the wave fronts and the dynamic and kinematic relationships.     

Stability of convective motion caused by inhomogeneously distributed internal heat sources

The stability of steady convective plane-parallel flow in a vertical layer of viscous incompressible liquid of thickness h is investigated. The motion is caused by heat sources distributed in the liquid with volume density Q = Q₀exp (^–x) (the x axis is taken perpendicular to the boundary layer). The region of instability is determined for various values of the Prandtl number and the parameter N = h characterizing the inhomogeneity of the internal sources. It is shown that with increase in N there is qualitative rearrangement of the stability limit for perturbations of hydrodynamic type and incremental thermal waves.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 3, pp. 140–144, May–June, 1977.     

Electrical fluctuations and entry flow instabilities during capillary flow of aqueous poly (ethylene oxide) solutions

Capillary flow of poly (ethylene oxide) solutions generates voltage fluctuations (noise) between electrodes placed on both sides of the capillary. The noise has a 1/f type spectral distribution, the value of increasing with (shear rate) to a limiting level. Within certain ranges, two sets of harmonic peaks appear in the spectra. It is demonstrated that these peaks are related to the frequency components of the instabilities in the entry flow region. The lower frequency set of harmonics corresponds to axial oscillations (pulsations) of the flow, while the higher frequency peaks are associated with the transverse oscillation of the stream lines in the vicinity of the entry. The corresponding frequencies were measured by visual counting and by spectral analysis of laser light transmitted through the entry region during injection of a coloured solution.The noise measurements were carried out using both platinum and reversible Ag/AgCl electrodes. In the latter case the streaming potential was also measured; its variation due to pressure fluctuations in the instable flow region appears to provide a plausible explanation of flow-induced noise phenomena as observed in elastic solutions.Dedicated to Prof. Dr. J. Schurz on the occasion of his sixtieth birthday.     

Elastic Waves in Swelling Porous Media

Time harmonic waves in a swelling porous elastic medium of infinite extent and consisting of solid, liquid and gas phases have been studied. Employing Eringen’s theory of swelling porous media, it has been shown that there exist three dilatational and two shear waves propagating with distinct velocities. The velocities of these waves are found to be frequency dependent and complex valued, showing that the waves are attenuating in nature. Here, the appearance of an additional shear wave is new and arises due to swelling phenomena of the medium, which disappears in the absence of swelling. The reflection phenomenon of an incident dilatational wave from a stress-free plane boundary of a porous elastic half-space has been investigated for two types of boundary surfaces: (i) surface having open pores and (ii) surface having sealed pores. Using appropriate boundary conditions for these boundary surfaces, the equations giving the reflection coefficients corresponding to various reflected waves are presented. Numerical computations are performed for a specific model consisting of sandstone, water and carbon dioxide as solid, liquid and gas phases, respectively, of the porous medium. The variations of phase speeds and their corresponding attenuation coefficients are depicted against frequency parameter for all the existing waves. The variations of reflection coefficients and corresponding energy ratios against the angle of incidence are also computed and depicted graphically. It has been shown that in a limiting case, Eringen’s theory of swelling porous media reduces to Tuncay and Corapcioglu theory of porous media containing two immiscible fluids. The various numerical results under these two theories have been compared graphically.     

Stability of stationary traveling waves on the surface of a vertical film of viscous fluid

The stability of stationary traveling waves of the first and second families with respect to infinitesimal perturbations of arbitrary wavelength is subjected to a detailed numerical investigation. The existence of a unique region of stability of the first family is established for wave numbers (₁, ₁) corresponding to the optimal wave regime. There are several regions of stability of the second family ( _k , _k),k=2,3,..., lying close to the local flow rate maxima. In the regions of instability the growth rates of perturbations of the first family are several times greater than for the second family. This difference increases with increase in the Reynolds number. The calculations make it possible to explain a number of experimental observations.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 3, pp. 33–41, May–June, 1989.The authors are grateful to V. Ya. Shkadov for his constant interest, and to A. G. Kulikovskii, A. A. Barmin and their seminar participants for useful discussions and suggestions.     

New concepts in relative permeabilities at high capillary numbers for surfactant flooding

Flooding oil reservoirs with surfactant solutions can increase the amount of oil that can be recovered. Macroscopic modelling of the process requires relative permeabilities to be functions of saturation and capillary number. With only limited experimental data, relative permeabilities have usually been assumed to be linear functions of saturation at high capillary numbers. The experimental data is reviewed, some of which suggest that this assumption is not necessarily correct. The basis for the assumption is therefore reviewed and it is concluded that the linear model corresponds to microscopically segregated flow in the porous medium. Based on new but equally plausible complementary assumptions about the flow pattern, a mixed flow model is derived. These models are then shown to be limiting cases of a droplet model which represents the mixing scale within the porous medium and gives a physical basis for interpolating between the models. The models are based on physical concepts of flow in a porous medium and so the approach described here represents a significant improvement in the understanding of high capillary number flow. This is shown by the fact that fewer parameters are needed to describe experimental data.Notation A total cross-sectional area assigned to capillary bundle - A ⁽ⁱ⁾ physical cross-sectional area of tube i - c ⁽ⁱ⁾ ordered configurational label for droplets in tube i - c configuration label for tube i (order not considered) - D defined by Equation (26) - E(...) expectation value with respect to the trinomial distribution - S _r ⁽⁾ fractional flow of phase - k absolute permeability - k _r relative permeability of phase - k _r ⁰ endpoint relative permeability of phase - L capillary tube length in bundle model - m ⁽ⁱ⁾ number of droplets of phase a occupying tube i - n exponent for phase a in Equation (2) - N number of droplets in bundle model - N _c capillary number - p pressure - p(c') probability of configuration c - Q ⁽ⁱ⁾ total volume flow rate in tube i - S saturation of phase - S flowing saturation of phase - S _r residual saturation of phase - S _r ⁽⁾ saturations when fractional flow of phase is 1 in the case of varying residual saturations for three-phase flow ( ) - t _c residence time for droplet configuration c - v ⁽ⁱ⁾ total fluid velocity in bundle tube i - , phase label - p pressure differential across capillary bundle - ⁽ⁱ⁾ tube conductivity defined by Equation (7) - viscosity of phase - interfacial tension - gradient operator - ... average over tube droplet configurations     

Wave propagation in transversely isotropic porous piezoelectric materials

Wave propagation in porous piezoelectric material (PPM), having crystal symmetry 6 mm, is studied analytically. Christoffel equation is derived for the propagation of plane harmonic waves in such a medium. The roots of this equation give four complex wave velocities which can propagate in such materials. The phase velocities of propagation and the attenuation quality factors of all these waves are described in terms of complex wave velocities. Phase velocities and attenuation of the waves in PPM depend on the phase direction. Numerical results are computed for the PPM BaTiO₃. The variation of phase velocity and attenuation quality factor with phase direction, porosity and the wave frequency is studied. The effects of anisotropy and piezoelectric coupling are also studied. The phase velocities of two quasi dilatational waves and one quasi shear waves get affected due to piezoelectric coupling while that of type 2 quasi shear wave remain unaffected. The phase velocities of all the four waves show non-dispersive behavior after certain critical high frequency. The phase velocity of all waves decreases with porosity while attenuation of respective waves increases with porosity of the medium. The characteristic curves, including slowness curves, velocity curves, and the attenuation curves, are also studied in this paper.     

Absolute determination of the isochromatic parameter by load-stepping photoelasticity

A new approach to full-field automated photoelasticity is presented in which a circular polariscope is used to enable the isochromatic phase value () to be determined unambiguously and without input of a known isochromatic value obtained using an auxiliary technique. Values of cos are obtained from light-field and dark-field images for three loads of small incremental steps. Using a relatively straight-forward procedure, ramped phase maps for are produced which can be unwrapped using conventional techniques. The resulting distribution of is then found absolutely using information provided by which is the incremental change in the isochromatic phase value between the load steps. The results obtained for disk-in-compression tests presented here in comparison with theoretical solutions demonstrate that the technique is both simple to use and very accurate. A similar approach may be adopted using three wavelengths instead of three load steps.     

Reservoir transport equations by compositional approach

The objective of this paper is to present an overview of the fundamental equations governing transport phenomena in compressible reservoirs. A general mathematical model is presented for important thermo-mechanical processes operative in a reservoir. Such a formulation includes equations governing multiphase fluid (gas-water-hydrocarbon) flow, energy transport, and reservoir skeleton deformation. The model allows phase changes due to gas solubility. Furthermore, Terzaghi's concept of effective stress and stress-strain relations are incorporated into the general model. The functional relations among various model parameters which cause the nonlinearity of the system of equations are explained within the context of reservoir engineering principles. Simplified equations and appropriate boundary conditions have also been presented for various cases. It has been demonstrated that various well-known equations such as Jacob, Terzaghi, Buckley-Leverett, Richards, solute transport, black-oil, and Biot equations are simplifications of the compositional model.Notation List B reservoir thickness - B formation volume factor of phase - Cⁱ mass of component i dissolved per total volume of solution - C ⁱ mass fraction of component i in phase - C heat capacity of phase at constant volume - C_p heat capacity of phase at constant pressure - D ⁱ hydrodynamic dispersion coefficient of component i in phase - D_MTf thermal liquid diffusivity for fluid f - F = F(x, y, z, t) defines the boundary surface - f_p fractional flow of phase - g gravitational acceleration - H_p enthalpy per unit mass of phase - J_p volumetric flux of phase - k_rf relative permeability to fluid f - k₀ absolute permeability of the medium - M_p ⁱ mass of component i in phase - n porosity - N rate of accretion - P_f pressure in fluid f - p_ca capillary pressure between phases and =p-p - Rⁱ rate of mass transfer of component i from phase to phase - Rⁱ _source source rate of component i within phase - S saturation of phase - s gas solubility - T temperature - t time - U displacement vector - u velocity in the x-direction - v velocity in the y-direction - V volume of phase - V_s velocity of soil solids - W_i body force in coordinate direction i - x horizontal coordinate - z vertical coordinate Greek Letters _p volumetric coefficient of compressibility - _T volumetric coefficient of thermal expansion - _ij Kronecker delta - volumetric strain - _m thermal conductivity of the whole matrix - internal energy per unit mass of phase - gf suction head - density of phase - _ij tensor of total stresses - _ij tensor of effective stresses - volumetric content of phase - f viscosity of fluid f     

On the Non-Isentropic Sound Waves Propagation in a Cylindrical Tube Filled with Saturated Porous Media

A rigid frame, cylindrical capillary theory of sound propagation in porous media that includes the nonlinear effects of the Forchheimer type is laid out by using variational solutions. It is shown that the five main parameters governing the propagation of sound waves in a fluid contained in rigid cylindrical tubes filled with a saturated porous media are: the shear wave number, , the reduced frequency parameter, , the porosity, ε, Darcy number, , and Forchheimer number, . The manner in which the flow influences the attenuation and the phase velocities of the forward and backward propagating non-isentropic acoustic waves is deduced. It is found that the inclusion of the solid matrix increases wave’s attenuations and phase velocities for both forward and backward sound waves, while increasing the porosity and the reduced frequency number decreased attenuation and increased phase velocities. The effect of the steady flow is found to decrease the attenuation and phase velocities for forward sound waves, and enhance them for the backward sound waves. This work is done during a sabbatical leave year granted form the University of Jordan to Dr. Hamzeh Duwairi for the academic year 2007/2008 at the German Jordanian University.     

On the viscosity-temperature behavior of polymer melts

Based on the free volume concept and the equation by Doolittle, an empirical equation is offered for the flow activation energy, E ^*, for polymer melts for the range of over 150°C above glass transition temperature, T _g. This E ^* represents the temperature coefficient of viscosity for the Newtonian region which is also equal to the value measured at constant shear stress for non-Newtonian flow. Data show that the E ^* of linear polymers approaches a constant value for a temperature range above T _g+150°C. Data on 17 polymers are correlated. The proposed equation for this region predicts the E ^* of polymer melts from the volume expansion coefficient, _l, above T _g and also from the T _g.Correlations have also been developed between E ^* and _l and between E ^* and T _g by simplifying the equation by use of the Simha-Boyer expression. A polymer having a lower _l or higher T _g generally has a higher E ^*. However, more satisfactory results are obtained by calculating E ^* from both _l and T _g. The E ^* calculated is found to agree with measurements within the experimental precision of about ±1 Kcal/mole.The effects of polymer composition, molecular weight, branching and microstructure on E ^* are also discussed. These factors influence E ^* in the way in which they effect _l and T _g.     

Mixed convection flow in narrow vertical ducts

The mixed convection flow in a vertical duct is analysed under the assumption that , the ratio of the duct width to the length over which the wall is heated, is small. It is assumed that a fully developed Poiseuille flow has already been set up in the duct before heat from the wall causes this to be changed by the action of the buoyancy forces, as measured by a buoyancy parameter . An analytical solution is derived for the case when the Reynolds numberRe, based on the duct width, is of 0 (1). This is extended to the case whenRe is 0 (^–1) by numerical integrations of the governing equations for a range of values of representing both aiding and opposing flows. The limiting cases, || 1 andR=Re of 0 (1), andR and both large, with of 0 (R ^1/3) are considered further. Finally, the free convection limit, large with R of 0 (1), is discussed.

---

Mischkonvektion in engen senkrechten Rohren

Zusammenfassung Mischkonvektion in einem senkrechten Rohr wird unter der Voraussetzung untersucht, daß das Verhältnis der Rohrbreite zur Länge, über welche die Wand beheizt wird, klein ist. Es wird angenommen, daß sich bereits eine voll entwickelte Poiseuille-Strömung in dem Rohr eingestellt hat, bevor Antriebskräfte, gemessen mit dem Auftriebsparameter , aufgrund der Wandbeheizung die Strömung verändern. Es wird eine analytische Lösung für den Fall erhalten, daß die mit der Rohrbreite als charakteristische Länge gebildete Reynolds-ZahlRe konstant ist. Dies wird mittels einer numerischen Integration der wichtigsten Gleichungen auf den FallRe =f (^–1) sowohl für Gleich- als auch für Gegenstrom ausgedehnt. Weiterhin werden die beiden Grenzfälle betrachtet, wenn || 1 undR=Re konstant ist, sowieR und beide groß mit proportionalR ^1/3. Schließlich wird der Grenzfall der freien Konvektion, großes mit konstantem R, diskutiert.

---

Nomenclature g acceleration due to gravity - Gr Grashof number - G modified Grashof number - h duct width - l length of the heated section of the duct wall - p pressure - Pr Prandtl number - Q flow rate through the duct - Q ₀ heat transfer on the wally=0 - Q ₁ heat transfer on the wally=1 - Re Reynolds number - R modified Reynolds number - T temperature of the fluid - T ₀ ambient temperature - T applied temperature difference - u, velocity component in thex-direction - v, velocity component in they-direction - x, co-ordinate measuring distance along the duct - y, co-ordinate measuring distance across the duct - buoyancy parameter - ₀ modified buoyancy parameter, ₀=R ^–1/3 - coefficient of thermal expansion - ratio of duct width to heated length, =h/l - (non-dimensional) temperature - _w applied temperature on the wally=0 - kinematic viscosity - density of the fluid - ₀ shear stress on the wally=0 - ₁ shear stress on the wally=1 - stream function     

Hypersonic flow over a narrow delta plate at large angles of attack

The problem of hypersonic flow over a flat delta plate with a high sweepback anglex at angles of attack close to /2 is solved using a numerical algorithm based on transition to the conical solution. The existence of conical flow at /2 with the velocity vector directed towards the apex of the plate is established. Values ofC _p/sin² and the thickness of the shock layer in the plane of symmetry of the plate are given as functions of the hypersonic similarity parameterk=tan tanx. A comparison of the calculated and experimental data shows that they are in good agreement.Translated from Izvestiya Rossiiskoi Akademii Nauk, Mekhanika Zhidkosti i Gaza, No.5, pp. 183–185, September–October, 1992.     

Numerical analysis of the branching of couette flow between concentric cylinders for various wave numbers

The character of stability loss of the circular Couette flow, when the Reynolds number R passes through the critical value R₀, is investigated within a broad range of variation of the wave numbers. The Lyapunov-Schmidt method is used [1, 2]; the boundary-value problems for ordinary differential equations arising in the case of its realization are solved numerically on a computer. It is shown that the branching character substantially depends on the wave number . For all a, excluding a certain interval (₁, ₂), the usual postcritical branching takes place: at a small supercriticality the circular flow loses stability and is softly excited into a secondary stationary flow — stable Taylor vortices. For wave numbers from the interval (₁,₂) a hard excitation of Taylor vortices takes place: at a small subcriticality R=R₀–² the secondary mode is unstable and merges with the Couette flow for 0; however, for a small supercriticality in the neighborhood of a circular flow there exist no stationary modes which are different.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 3, pp. 132–135, May–June, 1976.     

Darcy-Forchheimer natural,forced and mixed convection heat transfer in non-Newtonian power-law fluid-saturated porous media

The governing equation for Darcy-Forchheimer flow of non-Newtonian inelastic power-law fluid through porous media has been derived from first principles. Using this equation, the problem of Darcy-Forchheimer natural, forced, and mixed convection within the porous media saturated with a power-law fluid has been solved using the approximate integral method. It is observed that a similarity solution exists specifically for only the case of an isothermal vertical flat plate embedded in the porous media. The results based on the approximate method, when compared with existing exact solutions show an agreement of within a maximum error bound of 2.5%.Nomenclature A cross-sectional area - b _i coefficient in the chosen temperature profile - B ₁ coefficient in the profile for the dimensionless boundary layer thickness - C coefficient in the modified Forchheimer term for power-law fluids - C ₁ coefficient in the Oseen approximation which depends essentially on pore geometry - C _i coefficient depending essentially on pore geometry - C _D drag coefficient - C _t coefficient in the expression forK ^* - d particle diameter (for irregular shaped particles, it is characteristic length for average-size particle) - f _p resistance or drag on a single particle - F _R total resistance to flow offered byN particles in the porous media - g acceleration due to gravity - g _x component of the acceleration due to gravity in thex-direction - Grashof number based on permeability for power-law fluids - K intrinsic permeability of the porous media - K ^* modified permeability of the porous media for flow of power-law fluids - l _c characteristic length - m exponent in the gravity field - n power-law index of the inelastic non-Newtonian fluid - N total number of particles - Nux,_D,F local Nusselt number for Darcy forced convection flow - Nux,_D-F,F local Nusselt number for Darcy-Forchheimer forced convection flow - Nux,_D,M local Nusselt number for Darcy mixed convection flow - Nux,_D-F,M local Nusselt number for Darcy-Forchheimer mixed convection flow - Nux,_D,N local Nusselt number for Darcy natural convection flow - Nux,_D-F,N local Nusselt number for Darcy-Forchheimer natural convection flow - pressure - p exponent in the wall temperature variation - _Pe c characteristic Péclet number - _Pe x local Péclet number for forced convection flow - _Pe x modified local Péclet number for mixed convection flow - _Ra c characteristic Rayleigh number - _Ra x local Rayleigh number for Darcy natural convection flow - _Ra x local Rayleigh number for Darcy-Forchheimer natural convection flow - Re convectional Reynolds number for power-law fluids - Reynolds number based on permeability for power-law fluids - T temperature - T _e ambient constant temperature - T w,_ref constant reference wall surface temperature - T _w(X) variable wall surface temperature - T _w temperature difference equal toT w,_ref–T _e - T ₁ term in the Darcy-Forchheimer natural convection regime for Newtonian fluids - T ₂ term in the Darcy-Forchheimer natural convection regime for non-Newtonian fluids (first approximation) - T _N term in the Darcy/Forchheimer natural convection regime for non-Newtonian fluids (second approximation) - u Darcian or superficial velocity - u ₁ dimensionless velocity profile - u _e external forced convection flow velocity - u _s seepage velocity (local average velocity of flow around the particle) - u _w wall slip velocity - U c _M characteristic velocity for mixed convection - U c _N characteristic velocity for natural convection - x, y boundary-layer coordinates - x ₁,y ₁ dimensionless boundary layer coordinates - X coefficient which is a function of flow behaviour indexn for power-law fluids - effective thermal diffusivity of the porous medium - shape factor which takes a value of/4 for spheres - shape factor which takes a value of/6 for spheres - ₀ expansion coefficient of the fluid - _T boundary-layer thickness - T ₁ dimensionless boundary layer thickness - porosity of the medium - similarity variable - dimensionless temperature difference - coefficient which is a function of the geometry of the porous media (it takes a value of 3 for a single sphere in an infinite fluid) - ₀ viscosity of Newtonian fluid - ^* fluid consistency of the inelastic non-Newtonian power-law fluid - constant equal toX(2 ^2–n )/ - density of the fluid - dimensionless wall temperature difference     

Peristaltic pumping of a second-order fluid in a planar channel

The peristaltic motion of a non-Newtonian fluid represented by the constitutive equation for a second-order fluid was studied for the case of a planar channel with harmonically undulating extensible walls. A perturbation series for the parameter ( half-width of channel/wave length) obtained explicit terms of 0(²), 0(²Re²) and 0(₁Re²) respectively representing curvature, inertia and the non-Newtonian character of the fluid. Numerical computations were performed and compared to the perturbation analysis in order to determine the range of validity of the terms.Presented at the second conference Recent Developments in Structured Continua, May 23–25, 1990, in Sherbrooke, Québec, Canada     

Linear shear flow past a porous particle

Linear shear flow past a porous spherical particle is studied using a generalized boundary condition proposed by Jones. The torque on a porous sphere rotating in a quiescent fluid is calculated. Streamlines patterns are illustrated for the case of a particle freely suspended in a simple shear flow. These patterns are shown to differ significantly from those associated with an impermeable rigid sphere. Finally, an expression for the effective viscosity of a dilute suspension of porous spherical particles is obtained.Nomenclature A, B dimensionless flow parameter - a radius of the porous sphere - C, E, F constants of integration - d shear strength - d constant rate of deformation of ambient field - e rate of strain tensor - f, g functions of distance - k permeability of the porous medium - n unit normal vector - p pressure - p unit vector - Q coefficient of spherical harmonic - q filter velocity within the porous medium - r polar spherical coordinate - S _p surface of porous particle - S, T, T* coefficients of spherical harmonics - T torque exerted on the particle - u fluid velocity vector - x cartesian coordinates - dimensionless constant - , polar spherical coordinates - dimensionless flow parameter - viscosity of the fluid - stress tensor - rotational velocity of the particle - rotational velocity of the ambient field.     

The Sensitivity of Tracer Concentration to Nonuniform Permeability and Porosity

Concentration of tracer at a producing well is a function of time and of the flow field that transports the tracer from one well to another. Small localized changes in the permeability and porosity fields result in small changes in the observed concentration history. Calculation of the effect of changes in the petrophysical parameter fields on observed concentration is an integral part of most approaches to solving the inverse problem, i.e. calculation of permeability and porosity fields that are consistent with observed tracer concentrations. This paper does not address the inverse problem directly, only the forward problem, but the results are presented as semi-analytic formulas which are intended to provide physical insight into the limitations of the inverse procedure. In particular, for a simple example of dipole flow between two wells in a nearly homogeneous infinite 2D porous medium, the width of the region of significant influence is shown to scale as d where is the dispersivity and d is half the interwell distance. Also, for fixed injection and production rates the influence on concentration of variation in porosity and in log-permeability are shown to be similar in magnitude.     

Thermodynamically complete equations of state for α-, ε-, and γ-iron

The numerical model of phase transition in iron in stress waves described in [1] contains equations of state with a limited range of applicability. They do not consider thermal excitation of conduction electrons and the presence of and — -triple point on the phase equilibrium curve, the effect of which should appear in shock loading of porous or preheated specimens. The present study will offer thermodynamically complete equations of state for the -, -, -phases of iron, free of these shortcomings.Translated from Zhurnal Prikladnoi Mekhaniki i Tekhnicheskoi Fiziki, No. 3, pp. 112–114, May–June, 1986.