高温真实气体底部流动的NS方程数值计算

本文数值模拟了高超音速飞行时钝锥的底部流动。采用轴对称NS方程并考虑真实气体效应。湍流模型采用修正的Baldwin-Lomax涡粘性代数模型,数值方法空间离散对流项采用显式NND格式,粘性项采用中心差。时间离散采用三阶的龙格-库塔法。真实气体模型采用考虑七种组分四种反应的汉森模型。给出了底部流场的压力和温度分布及各组分的浓度分布。可以看出在近底部区域高速流-绕过拐角就产生一回流旋涡区。由于温度变化很大,气体的热力学特性受气体离解、复合和振动能激发的影响。所以整个流动过程变得十分复杂。     

不同发射深度下导弹水下点火气水流体动力计算

从流体动力角度研究了不同发射深度下，导弹水下点火这一非定常非线性过程。整个系统分为外部水流场、喷管流场和燃气泡流场三个区域加以考虑。水流场采用不可压势流模型，用边界元方法求解；喷管内流场采用非定常一元流动模型，用特征线差分法求解，并设置了激波检测功能；燃气泡采用基于质量和能量守恒的零维计算模型。在时间域中用步进方法实现了三个流场的耦合求解。给出了四种发射深度下的数值计算结果，展示了导弹水下点火的一     

Burnett simulations of gas flow in microchannels

The Burnett equations with slip boundary conditions are used to model the gas flow in microchannels in transition flow regime. As the Navier-Stokes equations are not appropriate to model the gas flow in this regime, the higher-order Burnett equations are adopted in the present study. In earlier studies, convergent solutions of the Burnett equations of microPoiseuille flow could only be obtained when Knudsen number is less than 0.2. By using a relaxation method on the boundary values, convergent solutions of the Burnett equations can be obtained even when Knudsen number reaches 0.4. The solutions of Burnett equations agree very well with experimental data and direct simulation Monte Carlo (DSMC) results. The pressure distributions and velocity profiles are then discussed in detail.     

Microscopic Mechanisms of Oil Recovery By Near-Miscible Gas Injection

As gas flooding becomes a more viable means of enhanced oil recovery, it is important to identify and understand the pore-scale flow mechanisms, both for the development of improved gas flooding applications and for the predicting phase mobilisation under secondary and tertiary gas flooding. The purpose of this study was to visually investigate the pore-level mechanisms of oil recovery by near-miscible secondary and tertiary gas floods. High-pressure glass micromodels and model fluids representing a near-miscible fluid system were used for the flow experiments. A new pore-scale recovery mechanism was identified which significantly contributed to oil recovery through enhanced flow and cross-flow between the bypassed pores and the injected gas. This mechanism is strongly related to a very low gas/oil interfacial tension (IFT), perfect wetting conditions and simultaneous flow of gas and oil in the same pore, all of which occur as the gas/oil critical point is approached. The results of this study helps us to better understand the pore-scale mechanisms of oil recovery in very low-IFT (near-miscible) systems. In particular we show that in near-miscible gas floods, behind the main gas front, the recovery of the oil continues by cross-flow from the bypassed pores into the main flow stream and as a result almost all of the oil, which has been contacted by the gas, could be recovered. Our observations in high-pressure micromodel experiments have demonstrated that this mechanism can only occur in near-miscible processes (as opposed to immiscible and completely miscible processes), which makes oil displacement by near-miscible gas floods a very effective process.     

Simulations for gas flows in microgeometries using the direct simulation Monte Carlo method

Micro gas flows are often encountered in MEMS devices and classical CFD could not accurately predict the flow and thermal behavior due to the high Knudsen number. Therefore, the gas flow in microgeometries was investigated using the direct simulation Monte Carlo (DSMC) method. New treatments for boundary conditions are verified by simulations of micro-Poiseuille flow, compared with the previous boundary treatments, and slip analytical solutions of the continuum theory. The orifice flow and the corner flow in microchannels are simulated using the modified DSMC codes. The predictions were compared with existing experimental phenomena as well as predictions using continuum theory. The results showed that the channel geometry significantly affects the microgas flow. In the orifice flow, the flow separation occurred at very small Reynolds numbers. In the corner flow, no flow separation occurred even with a high driving pressure. The DSMC results agreed well with existing experimental information.     

A Numerical Study of Microscale Flow Behavior in Tight Gas and Shale Gas Reservoir Systems

Various attempts have been made to model flow in shale gas systems. However, there is currently little consensus regarding the impact of molecular and Knudsen diffusion on flow behavior over time in such systems. Direct measurement or model-based estimation of matrix permeability for these “ultra-tight” reservoirs has proven unreliable. The composition of gas produced from tight gas and shale gas reservoirs varies with time for a variety of reasons. The cause of flowing gas compositional change typically cited is selective desorption of gases from the surface of the kerogen in the case of shale. However, other drivers for gas fractionation are important when pore throat dimensions are small enough. Pore throat diameters on the order of molecular mean free path lengths will create non-Darcy flow conditions, where permeability becomes a strong function of pressure. At the low permeabilities found in shale gas systems, the dusty-gas model for flow should be used, which couples diffusion to advective flow. In this study we implement the dusty-gas model into a fluid flow modeling tool based on the TOUGH+ family of codes. We examine the effects of Knudsen diffusion on gas composition in ultra-tight rock. We show that for very small average pore throat diameters, lighter gases are preferentially produced at concentrations significantly higher than in situ conditions. Furthermore, we illustrate a methodology which uses measurements of gas composition to more uniquely determine the permeability of tight reservoirs. We also describe how gas composition measurement could be used to identify flow boundaries in these reservoir systems. We discuss how new measurement techniques and data collection practices should be implemented in order to take advantage of this method. Our contributions include a new, fit-for-purpose numerical model based on the TOUGH+ code capable of characterizing transport effects including permeability adjustment and diffusion in micro- and nano-scale porous media.     

Analytical Modeling for Gravity Segregation in Gas Improved Oil Recovery of Tilted Reservoirs

Stone’s model for gravity segregation in gas improved oil recovery (IOR) indicates the distance that injected gas and water travel together before the segregation being completed (length of complete segregation). This model is very useful for co-injection of water and gas into horizontal depleted reservoirs. A proof by Rossen and van Duijn showed that Stone’s model applies to steady-state gas–liquid flow, and also foam flow, in horizontal reservoirs as long as the standard assumptions of fractional flow theory (incompressible flow, Newtonian mobilities, local equilibrium) are applied. However, until now, there has been no analytical study on the length of segregation when co-injection of water and gas occurs in tilted reservoirs. In this article, in order to extent the validity of Stone’s model to tilted reservoirs, governing equations of fluids displacement based on fractional flow theory are solved by the method of characteristics, MOC. The results are then compared to Stone’s model and to the results of a three-dimensional finite-difference compositional reservoir simulator. This study shows that Stone’s model should be corrected for tilted reservoirs and that the presented math proof can model gravity segregation in gas IOR of tilted reservoirs, appropriately. The effect of co-injecting of water and gas into tilted reservoirs on recovery efficiency is also examined.     

The direct simulation Monte Carlo method using unstructured adaptive mesh and its application

The implementation of an adaptive mesh‐embedding (h‐refinement) scheme using unstructured grid in two‐dimensional direct simulation Monte Carlo (DSMC) method is reported. In this technique, local isotropic refinement is used to introduce new mesh where the local cell Knudsen number is less than some preset value. This simple scheme, however, has several severe consequences affecting the performance of the DSMC method. Thus, we have applied a technique to remove the hanging node, by introducing the an‐isotropic refinement in the interfacial cells between refined and non‐refined cells. Not only does this remedy increase a negligible amount of work, but it also removes all the difficulties presented in the originals scheme. We have tested the proposed scheme for argon gas in a high‐speed driven cavity flow. The results show an improved flow resolution as compared with that of un‐adaptive mesh. Finally, we have used triangular adaptive mesh to compute a near‐continuum gas flow, a hypersonic flow over a cylinder. The results show fairly good agreement with previous studies. In summary, the proposed simple mesh adaptation is very useful in computing rarefied gas flows, which involve both complicated geometry and highly non‐uniform density variations throughout the flow field. Copyright © 2002 John Wiley & Sons, Ltd.     

Effect of compressibility and nonisothermicity on the efficiency of film cooling in a turbulent boundary layer

Most studies of film cooling deal with the analysis of the thermal shielding efficiency in an incompressible gas flow with constant physical properties. In practice, however, film cooling is used for machine elements exposed to high-speed, high-temperature gas flows. These conditions have received relatively little attention [1–3]. In this paper the influence of these factors on thermal shielding efficiency is analyzed by the method proposed in [4]. It is shown that the effect of compressibility and nonisothermicity on the thermal shielding efficiency is not very great.     

Driver gas contaminationin a detonation-driven reflected-shock tunnel

A computational study has been carried out to examine the effects of driver gas contamination in the NASA HYPULSE facility at GASL when operating with a detonation driver in reflected-shock tunnel mode. Unlike high-enthalpy shock tunnels which use helium as a driver gas, the driver gas in a detonation driver consists of a mixture of water vapour and argon, which has very different chemical and thermodynamic properties than those of helium. The purpose of the present work is to quantitatively evaluate the effects of driver gas contamination on the flow properties in the test section. Two computational analyses have been performed. The first analysis examined the nozzle flow under the influence of a prescribed level of driver gas contamination. In the second analysis, the transient development of the driver gas leakage in the reflected-shock region in the shock tube is studied. The unique flow features brought about by the detonation-driver gas and the method for detecting the contamination are discussed.Received: 6 August 2002, Revised: 22 September 2003, Published online: 10 February 2004Correspondence to: R.S.M. ChueAn abridged version of this paper was presented at the 23rd International Symposium on Shock Waves at Fort Worth, Texas, July 23 to 27, 2001     

The elimination of severe slugging-experiments and modeling

Severe slugging can occur in a pipeline-riser system operating at low liquid and gas rates. The flow of gas into the riser can be blocked by liquid accumulation at the base of the riser. This can cause formation of liquid slugs of a length equal to or longer than the height of the riser. A cyclic process results in which a period of no liquid production into the separator occurs, followed by a period of very high liquid production. This study is an experimental and theoretical investigation of two methods for eliminating this undesirable phenomenon, using choking and gas lift. Choking was found to effectively eliminate or reduce the severity of the slugging. However, the system pressure might increase to some extent. Gas lift can also eliminate severe slugging. While choking reduces the velocities in the riser, gas lift increases the velocities, approaching annular flow. It was found that a relatively large amount of gas was needed before gas injection would completely stabilize the flow through the riser. However, gas injection reduces the slug length and cycle time, causing a more continuous production and a lower system pressure. Theoretical models for the elimination of severe slugging by gas lift and choking have been developed. The models enable the prediction of the flow behavior in the riser. One model is capable of predicting the unstable flow conditions for severe slugging based on a static force balance. The second method is a simplified transient model based on the assumption of a quasi-equilibrium force balance. This model can be used to estimate the characteristics of the flow, such as slug length and cycle time. The models were tested against new severe slugging data acquired in this study. An excellent agreement between the experimental data and the theoretical models was found.     

Analytical Solutions to Gas Flow Problems in Low Permeability Porous Media

In most of conventional porous media the flow of gas is basically controlled by the permeability and the contribution of gas flow due to gas diffusion is ignored. The diffusion effect may have significant impact on gas flow behavior, especially in low permeability porous media. In this study, a dual mechanism based on Darcy flow as well as diffusion is presented for the gas flow in homogeneous porous media. Then, a novel form of pseudo pressure function was defined. This study presents a set of novel analytical solutions developed for analyzing steady-state and transient gas flow through porous media including effective diffusion. The analytical solutions are obtained using the real gas pseudo pressure function that incorporates the effective diffusion. Furthermore, the conventional assumption was used for linearizing the gas flow equation. As application examples, the new analytical solutions have been used to design new laboratory and field testing method to determine the porous media parameters. The proposed laboratory analysis method is also used to analyze data from steady-state flow tests of three core plugs. Then, permeability (k) and effective diffusion coefficient (D _e) was determined; however, the new method allows one to analyze data from both transient and steady-state tests in various flow geometries.     

On the accuracy of gas flow rate measurements in gas–liquid pipe flows by cross-correlating dual wire-mesh sensor signals

This article presents an assessment of the accuracy of gas flow rate measurement in gas–liquid pipe flows by cross-correlating dual wire-mesh sensor signals. The differences between the estimated and the actual gas superficial velocities in different flow regimes were analyzed. It was demonstrated that this gas flow rate measurement method is susceptible to significant systematic errors, some of which are inherent to the use of cross-correlation and others which are specific to wire-mesh sensors. It was concluded that this method would be accurate only for flow conditions within narrow ranges.     

The Solution by the Wave Curve Method of Three-Phase Flow in Virgin Reservoirs

There are two goals of this study. The first is to provide an introduction to the wave curve method for finding the analytic solution of a porous medium injection problem. Similar to fractional and chromatographic flow theory, the wave curve method is based on the method of characteristics, but it is applicable to an expanded range of physical processes in porous medium flow. The second goal is to solve injection problems for immiscible three-phase flow, as described by Corey’s model, in which a mixture of gas and water is injected into a porous medium containing oil and irreducible water. In particular we determine, for any choice of the phase viscosities, the proportion of the injected fluids that maximizes recovery around breakthrough time. Numerical simulations are performed to compare our solutions for Corey’s model with those of other models. For the injection problems we consider, solutions for Corey’s model are very similar to those for Stone’s model, despite the presence of an elliptic region in the latter; and they are very different from those for the Juanes-Patzek model, which preserves strict hyperbolicity. A nice feature of our analytical method is that it facilitates explaining both differences and similarities among the solutions for the three models considered.     

Mesoscopic Simulation of Rarefied Flow in Narrow Channels and Porous Media

Rarefied gas flow in channels and computer-aided reconstructions of porous media is simulated using the direct simulation Monte Carlo (DSMC) method and a modified lattice-Boltzmann (LB) method that can account for rarefaction effects. An increase of two orders of magnitude was noted for the gas permeability as the Knudsen number increased from 0.1 to 10. It was found that incorporation of a Bosanquet expression for the viscosity in the dusty gas flux equations leads to the recovery of the well known Klinkenberg expression for the gas permeability, revealing an explicit relation of the, thus far empirical, permeability correction factor to the fluid and structure properties. The expression for the effective gas viscosity in the transition flow regime is also incorporated in the LB method, which is then validated against the DSMC method by comparing predictions for the velocity profiles in straight channels over the entire transition flow regime. The new LB method offers the additional advantages of simplicity in the code implementation and great savings in computational time and memory compared to the DSMC method. It is shown that the rough adjustment of a single parameter suffices to make the LB method suitable for the reliable prediction of the gas permeability in porous media over the entire transition flow regime.     

Thermohydrodynamic characteristics of combined double-diffusive radiation convection heat transfer in a cavity

This work deals with the numerical analysis of a radiating gas flow caused by both temperature and buoyancy concentration gradients in a square cavity; in this regard, the set of governing equations, including conservation of mass, momentum, species, and energy are solved by a numerical technique. In terms of radiation, since the fluid is considered as a semitransparent medium, the radiative term in the energy equation appears and is calculated by numerical solving of the radiative transfer equation (RTE). Furthermore, all of the surrounding cavity walls are considered to be opaque, gray, and diffuse with constant emissivity. All of the flow equations are solved by the finite difference method (FDM) and the RTE by the discrete ordinate one (DOM). In the present study, an attempt is made to verify the optical thickness effects on flow, thermal behavior, and mass transform in a cavity flow, such that reciprocating trends were seen in this manner. Our numerical results show that the thermal field in double-diffusive convection flow reaches very fast its steady-state situation in comparison to the concentration distribution. Besides, it is found that the thermohydrodynamic characteristics of a double-diffusive convection flow of a radiating gas are much affected by optical thickness.     

Non-newtonian liquid-air stratified flow through horizontal tubes—ii

Non-Newtonian liquid gas stratified flow data were obtained using 0.052 and 0.025 m dia horizontal circular ducts. Unless the liquid velocity was very low, the flow pattern generally observed was non-uniform stratified flow having an interfacial level gradient between the two phases. The Heywood-Charles model is valid for predicting the pressure drop and liquid holdup in pseudoplastic (shear thinning) non-Newtonian liquid-gas uniform stratified flow. Two-phase drag reduction, which is predicted by the Heywood-Charles model did not occur because there was a transition to semi-slug flow before the model criteria were reached. Interfacial liquid and gas shear stresses were compared.     

Investigation of the influence of the conditions of swirling on the structure of a two-phase flow in an expanding channel

Both experimental and theoretical investigations show that the main feature in the structure of sufficiently strongly swirling gas flows is the presence in then of reverse circulation regions whose configuration depends very strongly both on the law of swirling of the flow and the conditions at the entrance as well as on the channel geometry [1–6]. In expanding channels, the occurrence of such regions is most probable in the axial region [7, 8]. In short annular channels for which the characteristic transverse and longitudinal dimensions are of the same order, reverse flow arises in the exit part of the channel along its inner wall [6, 9]. Hitherto, the investigations have been made for single-phase gas flows. The present paper reports a numerical investigation of the influence of particles of a condensed phase on the intensity of the reverse flow and the structure of the gas flow in an annular expanding channel under conditions of thermal, mass, and mechanical interaction of the phases. The method of stabilization was used to solve the boundary-value problem. The system of equations describing the axisymmetric unsteady flow of the two-phase medium was integrated by means of Godunov's difference scheme [10, 11]. The calculations were made for different conditions of injection of the particles of the condensed phase into the channel.     

Effect of gas flow swirling on coating deposition by the cold gas-dynamic spray method

The effect of gas flow swirling on the process of coating deposition onto a target by the cold gas-dynamic spray method is studied experimentally and numerically. Flow swirling is found to change the gas flow field and to reduce the gas flow rate under typical conditions of cold gas-dynamic spray. In a non-swirled flow, the shape of the deposited spot is similar to a sharp cone. In contrast, the deposited spot in a swirled flow is shaped as a crater without particles at the center of this crater. It is found that this effect is caused by centrifugal forces acting on particles in a swirled gas flow.     

Mathematical Model of Coalbed Gas Flow with Klinkenberg Effects in Multi-Physical Fields and its Analytic Solution

The deep-mining coal seam impacted by high in situ stress, where Klinkenberg effects for gas flow were very obvious due to low gas permeability, could be regarded as a porous and tight gas-bearing media. Moreover, the Klinkenberg effects had a significant effect on gas flow behavior of deep-mining coal seam. Based on the gas flow properties of deep-mining coal seams affected by in situ stress field, geothermal temperature field and geo-electric field, a new mathematical model of coalbed gas flow, which reflected the impact of Klinkenberg effects on coalbed gas flow properties in multi-physical fields, was developed by establishing the flow equation, state equation, and continuity equation and content equation of coalbed gas. The analytic solution was derived for the model of one-dimensional steady coalbed gas flow with Klinkenberg effects affected by in situ stress field and geothermal temperature field, and a sensitivity analysis of its physical parameters was carried out by comparing available analytic solutions and the measured values. The results show that the analytic solutions of this model of coalbed gas flow with Klinkenberg effects are closer to the measured values compared to those without Klinkenberg effects, and this model can reflect more accurately gas flow of deep-mining coal seams. Moreover, the analytic solution of this model is more sensitive to the change of Klinkenberg factor b and temperature grad G than depth h.