Application of the Density-Functional Method for Numerical Simulation of Flows of Multispecies Multiphase Mixtures

Numerical examples of application of the density functional used to describe isothermal flows of two-phase two-species mixtures are given. The following flows are calculated in a two-dimensional formulation: impact of a drop on a liquid layer, breakdown of a drop in the velocity field of the Couette flow, formation of the wetting angle of a drop on a solid surface, and development of the Rayleigh–Taylor and Kelvin–Helmholtz instabilities at the gas–liquid interface.     

Empirical slip and viscosity model performance for microscale gas flow

For the simple geometries of Couette and Poiseuille flows, the velocity profile maintains a similar shape from continuum to free molecular flow. Therefore, modifications to the fluid viscosity and slip boundary conditions can improve the continuum based Navier–Stokes solution in the non‐continuum non‐equilibrium regime. In this investigation, the optimal modifications are found by a linear least‐squares fit of the Navier–Stokes solution to the non‐equilibrium solution obtained using the direct simulation Monte Carlo (DSMC) method. Models are then constructed for the Knudsen number dependence of the viscosity correction and the slip model from a database of DSMC solutions for Couette and Poiseuille flows of argon and nitrogen gas, with Knudsen numbers ranging from 0.01 to 10. Finally, the accuracy of the models is measured for non‐equilibrium cases both in and outside the DSMC database. Flows outside the database include: combined Couette and Poiseuille flow, partial wall accommodation, helium gas, and non‐zero convective acceleration. The models reproduce the velocity profiles in the DSMC database within an L₂ error norm of 3% for Couette flows and 7% for Poiseuille flows. However, the errors in the model predictions outside the database are up to five times larger. Copyright © 2005 John Wiley & Sons, Ltd.     

Two-phase couette flow

A study is made of plane laminar Couette flow, in which foreign particles are injected through the upper boundary. The effect of the particles on friction and heat transfer is analyzed on the basis of the equations of two-fluid theory. A two-phase boundary layer on a plate has been considered in [1, 2] with the effect of the particles on the gas flow field neglected. A solution has been obtained in [3] for a laminar boundary layer on a plate with allowance for the dynamic and thermal effects of the particles on the gas parameters. There are also solutions for the case of the impulsive motion of a plate in a two-phase medium [4–6], and local rotation of the particles is taken into account in [5, 6]. The simplest model accounting for the effect of the particles on friction and heat transfer for the general case, when the particles are not in equilibrium with the gas at the outer edge of the boundary layer, is Couette flow. This type of flow with particle injection and a fixed surface has been considered in [7] under the assumptions of constant gas viscosity and the simplest drag and heat-transfer law. A solution for an accelerated Couette flow without particle injection and with a wall has been obtained in [6]. In the present paper fairly general assumptions are used to obtain a numerical solution of the problem of two-phase Couette flow with particle injection, and simple formulas useful for estimating the effect of the particles on friction and heat transfer are also obtained.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 3, pp. 42–46, May–June, 1976.     

Solution of the temperature jump (Knudsen layer flow) and linear heat-transfer problems for two parallel plates in a rarefied gas

The Monte Carlo method [1, 2] is used to solve the linearized Boltzmann equation for the problem of heat transfer between parallel plates with a wall temperature jump (Knudsen layer flow). The linear Couette problem can be separated into two problems: the problem of pure shear and the problem of heat transfer between two parallel plates. The Knudsen layer problem is also linear [3] and, like the Couette problem, can be separated into the velocity slip and temperature jump problems. The problems of pure shear and velocity slip have been examined in [2].The temperature jump problem was examined in [4] for a model Boltzmann equation. For the linearized Boltzmann equation the problems noted above have been solved either by expanding the distribution function in orthogonal polynomials [5–7], which yields satisfactory results for small Knudsen numbers, or by the method of moments, with an approximation for the distribution function selected from physical considerations in the form of polynomials [8–10]. The solution presented below does not require any assumptions on the form of the distribution function.The concrete calculations were made for a molecular model that we call the Maxwell sphere model. It is assumed that the molecules collide like hard elastic spheres whose sections are inversely proportional to the relative velocity of the colliding molecules. A gas of these molecules is close to Maxwellian or to a gas consisting of pseudo-Maxwell molecules [3].     

Couette flow of a binary mixture of rigid-sphere gases described by the linearized Boltzmann equation

A concise and accurate solution to the problem of plane Couette flow for a binary mixture of rigid-sphere gases described by the linearized Boltzmann equation and general (specular-diffuse) Maxwell boundary conditions for each of the two species of gas particles is developed. An analytical version of the discrete-ordinates method is used to establish the velocity, heat-flow, and shear-stress profiles for both types of particles, as well as the particle-flow and heat-flow rates associated with each of the two species. Accurate numerical results are given for the case of a mixture of helium and argon confined between molybdenum and tantalum plates.     

Nonisothermal Couette flow of a non-Newtonian fluid under the action of a pressure gradient

Nonisothermal Couette flow has been studied in a number of papers [1–11] for various laws of the temperature dependence of viscosity. In [1] the viscosity of the medium was assumed constant; in [2–5] a hyperbolic law of variation of viscosity with temperature was used; in [6–8] the Reynolds relation was assumed; in [9] the investigation was performed for an arbitrary temperature dependence of viscosity. Flows of media with an exponential temperature dependence of viscosity are characterized by large temperature gradients in the flow. This permits the treatment of the temperature variation in the flow of the fluid as a hydrodynamic thermal explosion [8, 10, 11]. The conditions of the formulation of the problem of the articles mentioned were limited by the possibility of obtaining an analytic solution. In the present article we consider nonisothermal Couette flows of a non-Newtonian fluid under the action of a pressure gradient along the plates. The equations for this case do not have an analytic solution. Methods developed in [12–14] for the qualitative study of differential equations in three-dimensional phase spaces were used in the analysis. The calculations were performed by computer.Translated from Zhurnal Prikladnoi Mekhaniki i Tekhnicheskoi Fiziki, No. 3, pp. 26–30, May–June, 1981.     

Three-dimensional laminar boundary layer on a permeable surface in the neighborhood of a symmetry plane

Analytical and numerical methods are used to investigate a three-dimensional laminar boundary layer near symmetry planes of blunt bodies in supersonic gas flows. In the first approximation of an integral method of successive approximation an analytic solution to the problem is obtained that is valid for an impermeable surface, for small values of the blowing parameter, and arbitrary values of the suction parameter. An asymptotic solution is obtained for large values of the blowing or suction parameters in the case when the velocity vector of the blown gas makes an acute angle with the velocity vector of the external flow on the surface of the body. Some results are given of the numerical solution of the problem for bodies of different shapes and a wide range of angles of attack and blowing and suction parameters. The analytic and numerical solutions are compared and the region of applicability of the analytic expressions is estimated. On the basis of the solutions obtained in the present work and that of other authors, a formula is proposed for calculating the heat fluxes to a perfectly catalytic surface at a symmetry plane of blunt bodies in a supersonic flow of dissociated and ionized air at different angles of attack. Flow near symmetry planes on an impermeable surface or for weak blowing was considered earlier in the framework of the theory of a laminar boundary layer in [1–4]. An asymptotic solution to the equations of a three-dimensional boundary layer in the case of strong normal blowing or suction is given in [5, 6].Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 5, pp. 37–48, September–October, 1980.     

Circulating rarefied gas flow in a rotating cylinder with a stationary upper end face

In order to analyze the intensity of the circulating flow in the rarefaction zone near the axis of a rapidly rotating cylinder, the problem of gas flow in a cylinder with a stationary upper end face is solved for intermediate Knudsen numbers by the direct statistical modeling method. The effect of the rarefaction of the gas on the intensity of the secondary flow is investigated. The contribution of the self-excited thermal circulation of the gas to the total circulating flow is analyzed.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 3, pp. 186–188, May–June, 1990.     

Self-similar solution to a problem of axisymmetric motion of gas through a porous medium in the case of a quadratic law of resistance

The results of solution of the self-similar problem of planar flow of gas through a porous medium in the case of a quadratic law of resistance [1] are generalized to the case of axisymmetric motion. The equation in similarity variables for the velocity of isothermal gas flow is reduced to an equation having cylindrical functions as solution. Analytic dependences of the pressure and the gas velocity on the coordinate and time are obtained for a given flow rate of the gas at the coordinate origin and for zero Initial gas pressure in the porous medium.Translated from Izvestlya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza., No. 4, pp. 168–171, July–August, 1982.     

Heat exchange of sphere with blowing for slow flows

An analytic solution is given for the problem of convective heat- and mass-exchange of a sphere with transverse flow of matter along the surface for values of Peclet numbers smaller than one and for blowing velocity smaller than that of the incoming gas flow. The solution for velocity field obtained by the authors in a previously published publication is employed on flow past a sphere with blowing; the method of asymptotic expansions of Acrivos and Taylor is also used. Expressions to the second approximation are determined for temperature field and for the values of local and averaged Nusselt numbers. It is shown that blowing reduces the temperature gradient or the concentrations at the surface.Translated from Zhurnal Prikladnoi Mekhaniki i Tekhnicheskoi Fiziki, No. 4, pp. 89–94, July–August, 1972.     

Prediction of amount of entrained droplets in vertical annular two-phase flow

Prediction of amount of entrained droplets or entrainment fraction in annular two-phase flow is essential for the estimation of dryout condition and analysis of post dryout heat transfer in light water nuclear reactors and steam boilers. In this study, air–water and organic fluid (Freon-113) annular flow entrainment experiments have been carried out in 9.4 and 10.2 mm diameter test sections, respectively. Both the experiments covered three distinct pressure conditions and wide range of liquid and gas flow conditions. The organic fluid experiments simulated high pressure steam–water annular flow conditions. In each experiment, measurements of entrainment fraction, droplet entrainment rate and droplet deposition rate have been performed by using the liquid film extraction method. A simple, explicit and non-dimensional correlation developed by Sawant [Sawant, P.H., Ishii, M., Mori, M., 2008. Droplet entrainment correlation in vertical upward co-current annular two-phase flow. Nucl. Eng. Des. 238 (6), 1342–1352] for the prediction of entrainment fraction is further improved in this study in order to account for the existence of critical gas and liquid flow rates below which no entrainment is possible.Additionally, a new correlation is proposed for the estimation of minimum liquid film flow rate at the maximum entrainment fraction condition. The improved correlation successfully predicted the newly collected air–water and Freon-113 entrainment fraction data. Furthermore, the correlations satisfactorily compared with the air–water, helium–water and air–genklene experimental data measured by Willetts [Willetts, I.P., 1987. Non-aqueous annular two-phase flow. D.Phil. Thesis, University of Oxford]. However, comparison of the correlations with the steam–water data available in literature showed significant discrepancies. It is proposed that these discrepancies might have been caused due to the inadequacy of the liquid film extraction method used to measure the entrainment fraction or due to the change in mechanism of entrainment under high liquid flow conditions.     

Calculation of a supersonic gas jet exhausting into vacuum from a nozzle with an inclined exit

A study is made of the problem of supersonic exhausting of ideal gas into vacuum from a conical nozzle with inclined exit. The solution is sought in a region including part of the flow where the projection of the velocity vector of the gas onto the nozzle axis can be less than the local velocity of sound and take negative values.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 3, pp. 185–186, May–June, 1982.     

Effect of an inhomogeneous magnetic field on the structure of an ionized gas flow

Results of an experimental study of the flow of an ionized gas produced by a shock wave through an inhomogeneous magnetic field are presented. Braking of the gas flow produced by the end currents is determined at two fixed sections of the magnetogasdynamic channel from the value of the isolated shock wave formed in the vicinity of the hemispherical model over which the flow passes. Maximum recorded reduction in Mach number was 30%, and with a magnetogasdynamic interaction parameter greater than 1.5, a transition of supersonic flow to infrasonic at the exit of the magnetic zone was observed. Experimental results were compared with a solution of a model problem which assumed one-dimensional flow in the flow core. The gas used was argon, with a maximum magnetic field induction of 1.5 T.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 5, pp. 174–178, September–October, 1976.     

Rarefied gas shear flow between two movable segments of parallel plates

A study is made of the two-dimensional steady-state rarefied gas flow observed between two parallel plane surfaces of finite and different length when one of the surfaces is fixed and the other moves parallel to itself at a constant velocity, while remaining within the bounds of a given segment with fixed ends (the motion is similar to that of a conveyer belt). This flow can be regarded as a twodimensional counterpart of the classical one-dimensional Couette flow. The corresponding problem is formulated in a rectangular domain for the nonlinear kinetic equation with a model collision operator and is solved by a finite-difference method for various boundary conditions. For simplicity's sake, the flow was studied under conditions such that it can be considered near-isothermal. The gas pressures on each side of the gap formed by the plates may be the same or different. If the pressures on both sides of the gap are equal, then a near-zero-gradient flow develops between the plates. In this case, the greater the plate length, the nearer the flow in the middle of the gap to one-dimensional Couette flow. The end effects are examined, together with the conditions in which the flow in the middle of the domain can be assumed to be practically one-dimensional. In the zero-gradient regime, the system operates, in general, as a pump transferring gas from one side of the gap the other. The ability to pump gas also remains if a small counterpressure exists.Translated from Izvestiya Rossiiskoi Akademii Nauk, Mekhanika Zhidkosti i Gaza, No. 3, pp. 150–155, May–June, 1995.The work was financially supported by the Russian Foundation for Fundamental Research (project No. 93-013-17928).     

Resistance of rough plates in a rotating fluid

Generalizing Navier’s partial slip condition, the flow due to a rough or striated plate moving in a rotating fluid is studied. It is found that the motion of the plate, the fluid surface velocity, and the shear stress are in general not in the same direction. The solution is extended to the case of finite depth, or Couette slip flow in a rotating system. In this case an optimum depth for minimum drag is found. The solutions are also closed form exact solutions of the Navier–Stokes equations. The results are fundamental to flows with Coriolis effects.     

Low speed MHD Couette flow of a slightly rarefied and electrically conducting gas

The effects of an applied magnetic field on the steady, laminar, low speed plane Couette flow of a slightly rarefied and electrically conducting gas are studied. Consideration is given to the slip-flow regime, wherein the gas rarefaction begins to play its important role. The generally accepted method of analysis for slip flows is utilized, i.e. the continuum magnetohydrodynamic equations of motion are used throughout the gas, together with the first and the second order slip velocity and temperature jump boundary conditions. Considerations are further given to (1) the case of zero electric field and (2) the case of a nonconducting channel in which the net current across the channel is zero.     

Development of pulsation regimes in one-dimensional unsteady MHD flows with switching off of the electrical conductivity

This paper investigates the gas flow in an electromagnetic field when the conductivity, being a function of the thermodynamic gas parameters, vanishes during the flow (switching off of the conductivity). In the case of steady supersonic flows in an expanding nozzle it was first shown analytically [1] and then confirmed by numerical experiment [2] that stable steady flow is not possible for all the problem parameters (for example, the values of the magnetic field at the exit). Instead of a steady flow a periodic regime is realized when narrow regions of conducting gas with currents flowing through them detach from the conducting region and propagate down the channel. In these papers the conductivity was assumed to be a function of only the temperature, such that for T T* (T) = 0. In [3, 4] the flows of conducting gas in the channels were calculated both with the given dependence of the gas conductivity on the temperature and on the basis of a three-component model by means of the Saha equation. At the same time, the development of periodic regimes in the flow in the nozzle was observed in both cases, but the mechanism of the origin of the current layers was not explained. The self-similar problem of the withdrawal of a nonconducting piston from a half-space occupied by a conducting gas with a magnetic field was investigated in [5] in a linear formulation. At the same time, regions of the problem parameters (the velocity of the piston and the magnetic field on it) were found when, in spite of the self-similar formulation of the problem, there is no self-similar solution. At the same time, regions exist where several solutions are possible. The possibility of the formation of isothermal rarefaction zones with low electrical conductivity when the Joule heating is balanced by the cooling of the gas on expansion (Butler waves) [6] was not taken into account in this paper, since they are unstable with respect to superheating. However, in the case of flow in a nozzle it was shown [2] that precisely the development of instabilities in these zones leads to the formation of the periodic regime. In the present paper the solution of the self-similar problem is constructed in a nonlinear formulation. The reason for the occurrence of regions in which the solution is multiply valued, which is associated with the process of arrival at self-similar boundary conditions, is explained. It is shown that a quasiperiodic regime can arise in the solution, occurring, in particular, in the regions of the problem parameters where there is no self-similar solution.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 4, pp. 115–122, July–August, 1986.     

Oscillatory rotationally symmetric loss of stability of nonisothermal couette flow

A study is made of the stability of nonisothermal Couette flow — steady flow of a viscous heat conducting fluid between two rotating concentric cylinders heated to different temperatures. The methods of perturbation theory are used to establish conditions sufficient for bifurcation of a neutral curve of oscillatory instability from the neutral curve of monotonic instability. Computer calculations show that for certain values of the parameters of the problem these conditions are realized and there is an oscillatory loss of stability of the nonisothermal Couette flow.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 1, pp. 76–80, January–February, 1984.I thank V. I. Yudovich for constant interest in the work.     

Blast wave propagation in air from a gaseous charge

In the point explosion problem it is assumed that an instantaneous release of finite energy causing shock wave propagation in the ambient gas occurs at a space point. The results of the solution of the problem of such blasts are contained in [1–4]. This point model is applied for the determination of shock wave parameters when the initial pressure in a sphere of finite radius exceeds the ambient air pressure by 2–3 orders of magnitude. The possibility of such a flow simulation at a certain distance from the charge is shown in papers [4, 5] as applied to the blast of a charge of condensed explosive and in [6, 8] as applied to the expansion of a finite volume of strongly compressed hot gas. In certain practical problems the initial pressure in a volume of finite dimensions exceeds atmospheric pressure by a factor 10–15 only. Such cases arise, for example, in the detonation of gaseous fuel-air mixtures. The present paper considers the problem of shock wave propagation in air, caused by explosion of gaseous charge of spherical or cylindrical shape. A numerical solution is obtained in a range of values of the specific energy of the charge characteristic for fuel-air detonation mixtures by means of the method of characteristics without secondary shock wave separation. The influence of the initial conditions of the gas charge explosion (specific energy, nature of initiation, and others) is investigated and compared with the point case with respect to the pressure difference across the shock wave and the positive overpressure pulse.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 3, pp. 110–118, May–June, 1986.     

Interaction of rarefaction waves with wet foam

The interaction of rarefaction waves of different shapes with wet water foams is studied experimentally. It is found that the observed values of the pressure are greater, while the surface velocity is lower than the corresponding values predicted by the pseudogas model. The foam breakdown starts as the pressure decreases by 0.3 atm relative to the initial pressure. During downstream propagation of the rarefaction-wave leading edge the propagation velocity decreases.Using of water-based foams as effective screens for damping blast waves in different technological processes has caused considerable interest in studying wave propagation in such systems. The pressure wave dynamics in a foam have been investigated in much detail, both experimentally and theoretically [1–3]. However, the interaction of rarefaction waves with foam has practically never been studied, although it was mentioned in [4] that the unloading phase following the compression wave phase is one of the factors defining the damaging action of blast waves. Besides blast-wave damping, rarefaction wave propagation takes place if such waves are used to breakup foam in oil-producing wells [5].Below, the interaction of rarefaction waves of different shapes with wet water foams is studied. The vertical shock tube described in detail in [3] was used in these experiments.Brest. Translated from Izvestiya Rossiiskoi Akademii Nauk, Mekhanika Zhidkosti i Gaza, No. 2, pp. 76–82, March–April, 1995.