Non-stationary compressible flow in ducts with varying cross-section

Author(s):  
O Igra ◽  
L Wang ◽  
J Falcovitz

A comparative study of the interaction of shock wave or rarefaction wave with a converging duct separating long constant cross-section segments is presented. Quasi-one-dimensional computations are compared with fully two-dimensional computations. It is observed that in some cases the two-dimensional results approach the respective one-dimensional approximations over long times, while in other cases the two-dimensional computed flow is genuinely two-dimensional and cannot be reduced to a one-dimensional equivalent. In the latter cases, significant errors are incurred by analysing the flow using the quasi-one- dimensional duct flow approximation.

2010 ◽  
Vol 7 ◽  
pp. 90-97
Author(s):  
M.N. Galimzianov ◽  
I.A. Chiglintsev ◽  
U.O. Agisheva ◽  
V.A. Buzina

Formation of gas hydrates under shock wave impact on bubble media (two-dimensional case) The dynamics of plane one-dimensional shock waves applied to the available experimental data for the water–freon media is studied on the base of the theoretical model of the bubble liquid improved with taking into account possible hydrate formation. The scheme of accounting of the bubble crushing in a shock wave that is one of the main factors in the hydrate formation intensification with increasing shock wave amplitude is proposed.


1999 ◽  
Vol 67 (2) ◽  
pp. 274-281 ◽  
Author(s):  
D. Das ◽  
J. H. Arakeri

In this paper we give a procedure to obtain analytical solutions for unsteady laminar flow in an infinitely long pipe with circular cross section, and in an infinitely long two-dimensional channel, created by an arbitrary but given volume flow rate with time. In the literature, solutions have been reported when the pressure gradient variation with time is prescribed but not when the volume flow rate variation is. We present some examples: (a) the flow rate has a trapezoidal variation with time, (b) impulsively started flow, (c) fully developed flow in a pipe is impulsively blocked, and (d) starting from rest the volume flow rate oscillates sinusoidally. [S0021-8936(00)01702-5]


1975 ◽  
Vol 15 (04) ◽  
pp. 269-276 ◽  
Author(s):  
J.R. Kyte ◽  
D.W. Berry

Abstract This paper presents an improved procedure for calculating dynamic pseudo junctions that may be used in two-dimensional, areal reservoir simulations to approximate three-dimensional reservoir behavior. Comparison of one-dimensional areal and two-dimensional vertical cross-sectional results for two example problems shows that the new pseudos accurately transfer problems shows that the new pseudos accurately transfer the effects of vertical variations in reservoir properties, fluid pressures, and saturations from the properties, fluid pressures, and saturations from the cross-sectional model to the areal model. The procedure for calculating dynamic pseudo-relative permeability accounts for differences in computing block lengths between the areal and cross-sectional models. Dynamic pseudo-capillary pressure transfers the effects of pseudo-capillary pressure transfers the effects of different pressure gradients in different layers of the cross-sectional model to the areal model. Introduction Jacks et al. have published procedures for calculating dynamic pseudo-relative permeabilities fro m vertical cross-section model runs. Their procedures for calculating pseudo functions are procedures for calculating pseudo functions are more widely applicable than other published approaches. They demonstrated that, in some cases, the derived pseudo functions could be used to simulate three-dimensional reservoir behavior using two-dimensional areal simulators. For our purposes, an areal simulator is characterized by purposes, an areal simulator is characterized by having only one computing block in the vertical dimension. The objectives of this paper are to present an improved procedure for calculating dynamic pseudo functions, including a dynamic pseudo-capillary pressure, and to demonstrate that the new procedure pressure, and to demonstrate that the new procedure generally is more applicable than any of the previously published approaches. The new pseudos previously published approaches. The new pseudos are similar to those derived by jacks et al. in that they are calculated from two-dimensional, vertical cross-section runs. They differ because (1) they account for differences in computing block lengths between the cross-sectional and areal models, and (2) they transfer the effects of different flow potentials in different layers of the cross-sectional potentials in different layers of the cross-sectional model to the areal model. Differences between cross-sectional and areal model block lengths are sometimes desirable to reduce data handling and computing costs for two-dimensional, areal model runs. For very large reservoirs, even when vertical calculations are eliminated by using pseudo functions, as many as 50,000 computing blocks might be required in the two-dimensional areal model to minimize important errors caused by numerical dispersion. The new pseudos, of course, cannot control numerical pseudos, of course, cannot control numerical dispersion in the cross-sectional runs. This is done by using a sufficiently large number of computing blocks along die length of the cross-section. The new pseudos then insure that no additional dispersion will occur in the areal model, regardless of the areal computing block lengths. Using this approach, the number of computing blocks in the two-dimensional areal model is reduced by a factor equal to the square of the ratio of the block lengths for the cross-sectional and areal models. The new pseudos do not prevent some loss in areal flow-pattern definition when the number of computing blocks in the two-dimensional areal model is reduced. A study of this problem and associated errors is beyond the scope of this paper. Our experience suggests that, for very large reservoirs with flank water injection, 1,000 or 2,000 blocks provide satisfactory definition. Many more blocks provide satisfactory definition. Many more blocks might be required for large reservoirs with much more intricate areal flow patterns. The next section presents comparative results for cross-sectional and one-dimensional areal models. These results demonstrate the reliability of the new pseudo functions and illustrate their advantages pseudo functions and illustrate their advantages over previously derived pseudos for certain situations. The relationship between two-dimensional, vertical cross-sectional and one-dimensional areal reservoir simulators has been published previously and will not be repeated here in any detail. Ideally, the pseudo functions should reproduce two-dimensional, vertical cross-sectional results when they are used in the corresponding one-dimensional areal model. SPEJ P. 269


1970 ◽  
Vol 10 (04) ◽  
pp. 405-417 ◽  
Author(s):  
N.D. Shutler

Abstract This paper describes a numerical mathematical model that is a significant extension of a previously published one-dimensional model of the steamflood published one-dimensional model of the steamflood process. process. The model describes the simultaneous flow of the three phases - oil, water and gas - in two dimensions. Interphase mass transfer between water and gas phases is allowed, but the oil is assumed nonvolatile and the hydrocarbon gas insoluble in the liquid phases. The model allows two-dimensional heat convection within the reservoir and two-dimensional heat conduction in a vertical cross-section spanning the oil sand and adjacent strata. Example calculations are presented which, on comparison with experimental results, tend to validate the model. Steam overriding due to gravity effects is shown to significantly reduce oil recovery efficiency in a thick system while jailing to do so in a thinner system. A study of the effect of capillary pressure indicates that failure to scale capillary forces in laboratory models of thick sands may lead to optimistic recovery predictions, while properly scaled capillary forces may be sufficiently low as to play no important role in oil recovery. Calculations made with and without vertical permeability show that failure to account for vertical fluid flow can lead to predictions of pessimistic oil recovery efficiency. pessimistic oil recovery efficiency Introduction Mathematical tools of varying complexity have been used in studying the steamflood process. A "simplified" class of mathematical models has served primarily as aids in engineering design. A more comprehensive class of models has improved understanding of the nature of the process. The model described in this report is of the latter class, but it is more comprehensive than any previously published model. published model. All previously available calculations of the steamflood process are confined to one space dimension in their treatments of fluid flow. Thus all previous models necessary ignore all effects of gravity reservoir heterogeneity, and nonuniform initial fluid-phase distributions on fluid flow in a second dimension. This model, an extension of a previously published model accounts for heat and previously published model accounts for heat and fluid transfer in two space dimensions and, hence, can evaluate these effects on simultaneous horizontal and vertical flow. While the model can describe the areal performance of a steamflood (in which case the heat transfer is described in three dimensions), this aspect will not be considered in this paper. Rather, this paper will describe the model in its application to a vertical cross-section through the reservoir and will consider some preliminary investigations to demonstrate the importance of being able to simultaneously account for horizontal and vertical fluid flow. Mathematical details are given in appendices. MATHEMATICAL DESCRIPTION OF STEAMFLOODING Darcy's law provides expressions for the velocities of the three phases (oil, water and gas), which, when combined with oil, water and gas mass balances give the partial differential equations governing Now of the three phases within a reservoir sand: OIL PHASE ..(1) WATER PHASE ..(2) SPEJ P. 405


1963 ◽  
Vol 3 (01) ◽  
pp. 19-27 ◽  
Author(s):  
P.M. Blair ◽  
D.W. Peaceman

Abstract The shape and position of the gas-oil transition zone during downdip displacement of oil by gas has been calculated using flow equations which include the effects of gravity, relative permeability, capillary pressure and compressibility of the fluids. The calculations treat the problem in two space dimensions, and results are compared with data from a laboratory model tilted at 30 degrees and 60 degrees from the horizontal on displacements near and above the maximum rate at which gravity segregation prevents channeling of the gas along the top of the stratum. The good agreement between calculated and experimental results demonstrates the validity of the technique as well as that of the flow equations. Introduction Knowledge of the fluid distribution and movement in and oil reservoirs important in producing operations and estimation of reserves. The history of the oil industry has included steady progress in improving the accuracy of calculations which provide the required knowledge. The earliest method of calculating reservoir performance consisted of material-balance equations based on the assumption that all properties were uniform throughout a reservoir. For many reservoirs such a simple formulation is still the most useful. However, when large pressure and saturation gradients exist in a reservoir, the assumption of uniform values throughout may lead to significant error. To reduce these errors, Buckley and Leverett introduced a displacement equation which considers pressure and saturation gradients. Methods available at that time permitted solutions to the Buckley-Leverett equation in one space dimension; these solutions have been very useful in solving many problems related to the production of oil. However, the one-dimensional methods are not adequate for systems in which saturations vary in directions other than the direction of flow. An example of such a system is the case of gas displacing oil down a dipping stratum in which the gas-oil contact becomes significantly tilted. Of course, the Buckley-Leverett displacement method cannot predict the tilt of the gas-oil contact. Recent improvements of the one-dimensional Buckley-Leverett method achieve some success in predicting the tilt of the gas-oil contact at sufficiently low flow rates. However, at rates high enough that the viscous pressure gradient nearly equals or exceeds the gravity gradient, even these improved one-dimensional methods incorrectly predict the shape and velocity of the contact. Further progress in estimating such fluid movements in a reservoir appears to require consideration of the problem in more than one space dimension. The recent two-dimensional method of Douglas, Peaceman and Rachford appears adaptable to calculate changes with time of the saturation distribution in a vertical cross-section of a reservoir. The movement of saturation contours should represent the moving fluid contacts and include the effects of crossflow due to gravity, as well as variations in the rock and fluid properties. The nonlinear nature of the equations used in the method has prevented proof of the validity of the solutions. Douglas, Peaceman and Rachford made some comparisons with experiment but did not include cases in which gravity was important nor cases involving displacement by the nonwetting phase. Forthesereasons, atestof the two-dimensional method for a case in which these factors are included would be very desirable. The test selected was a comparison of calculated results with those from a carefully controlled laboratory experiment on a model with measured physical properties. The model selected was one in which gas displaced oil down a tilted, rectangular sand pack. The model can be thought of as representing a vertical cross-section taken parallel to the dip of a reservoir. The displacement thus simulates gas displacing oil downdip that might result from gas-cap expansion or gas injection. SPEJ P. 19^


2021 ◽  
Vol 1036 ◽  
pp. 432-441
Author(s):  
Ming Qi Li ◽  
Hui Su Chen

The ITZ (interfacial transition zone) in concrete has very high porosity and permeability, which affects concrete’s macroscopic mechanical properties and transport properties. Two-dimensional (2D) areal analysis and one-dimensional (1D) linear analysis are usually used to study ITZ’s microstructure. However, 3D microstructure is difficult to be characterized by 1D and 2D information. For example, 2D cross-section planes do not always intercept both the ITZ and the corresponding aggregate, which causes some ITZ regions are ignored by researchers. Therefore, ITZ’s volume and thickness will be misestimated, and leads to the misestimation of the diffusivity. In this paper, the effect of aggregate’s shape on the misestimation of ITZ thickness t is studied. The results reveal that the misestimation increases with the increasing sphericity s of aggregates.


1974 ◽  
Vol 96 (1) ◽  
pp. 21-24 ◽  
Author(s):  
A. T. Olson

Using Walz’s approximation method for boundary layer calculation, along with a one-dimensional treatment of the compressible inviscid core flow, discharge coefficients for small nozzle to pipe diameter ratios have been calculated. Discharge co-efficients calculated for the ASME long radius nozzle agree with those recommended by the ASME Power Test Code. In addition, experimental confirmation of an indicated Mach number effect has been achieved in a nozzle modified to minimize two-dimensional effects.


1995 ◽  
Vol 53 (3) ◽  
pp. 267-276 ◽  
Author(s):  
J. Ashkenazy ◽  
D. Zoler

A quasi-one-dimensional model for the steady-state flow of a plasma in an ablative discharge capillary is presented for capillaries with non-constant cross- section. It is demonstrated that small modifications of the capillary geometry can lead to significant changes in the plasma exit parameters. In this respect, the possibility of obtaining an extended range of plasma parameters makes this type of ablative capillary a useful source of plasma for a variety of applications. Numerical solution of the equations of the model for the critical-flow case allows evaluation of the main fluid-dynamic and thermodynamic parameters of the plasma inside the capillary and at its exit.


2003 ◽  
Vol 125 (3) ◽  
pp. 529-537 ◽  
Author(s):  
J. Lo´pez ◽  
F. Faura ◽  
J. Herna´ndez ◽  
P. Go´mez

During the initial slow stage of the injection process in high-pressure die casting machines with horizontal cold chamber, a plunger pushes the molten metal which partially fills the injection chamber, causing the formation of a gravity wave. The evolution of the wave surface profile, which depends on the plunger acceleration law, may trap air in the molten metal, causing porosity when the metal solidifies. In this work, a one-dimensional shallow-water model, which is solved numerically using the method of characteristics, and a three-dimensional numerical model, based on a finite element formulation and the volume of fluid (VOF) method for treating the free surface, are used to analyze the flow of molten metal in an injection chamber of circular cross section. The results for the evolution of the free surface obtained from both models for different plunger motion laws and initial filling fractions of the injection chamber were in good agreement for broad ranges of operating conditions. The existence of a critical plunger speed, above which the reflection of the wave of molten metal against the chamber ceiling might appreciably increase air entrapment effects, is investigated. The results for the wave profiles in chambers of circular cross section are compared with those obtained in an equivalent two-dimensional configuration of the injection chamber, for which the shallow-water model is solved analytically. It is shown how the results obtained by applying the one-dimensional model to a two-dimensional chamber configuration can be used to reproduce, with an acceptable degree of accuracy, the salient characteristics of the flow of molten metal in a real injection chamber of circular cross section.


Sign in / Sign up

Export Citation Format

Share Document