Slow gas motions and the negative drag of a strongly heated spherical particle

In the slow flows of a strongly and nonuniformly heated gas, in the continuum regime (Kn → 0) thermal stresses may be present. The theory of slow nonisothermal continuum gas flows with account for thermal stresses was developed in 1969–1974. The action of the thermal stresses on the gas results in certain paradoxical effects, including the reversal of the direction of the force exerted on a spherical particle in Stokes flow. The propulsion force effect is manifested at large but finite temperature differences between the particle and the gas. This study is devoted to the thermal-stress effect on the drag of a strongly heated spherical particle traveling slowly in a gas for small Knudsen numbers (M ～ Kn → 0), small but finite Reynolds numbers (Re ≤ 1), a linear temperature dependence of the transport coefficients µ ∝ T, and large but finite temperature differences ((T _w ? T _∞ )/T _M8 ～ 1). Two different systems of equations are solved numerically: the simplified Navier-Stokes equations and the modified Navier-Stokes equations with account for the thermal stresses.     

认识稀薄气体动力学

以通俗易懂的方式介绍了空气动力学当气体间断分子效应显著时发展起来的特殊分文——稀薄气体动力学、讨论了非平衡现象与稀薄气体动力学的关系．通过与8速度气体模型的间断Boltzmann方程的对比，解释了Boltzmann方程碰撞项的物理意义和数学困难，简要综述了其一般解法、讨论了分子在物体表面的反射和问题的边界条件，着重介绍了直接模拟Monte Carlo(DSMC)方法和为克服低速稀薄流动(如MEMS中流动)中模拟困难的信息保存(IP)方法。     

Region of applicability and the main features of the generalized Chapman-Enskog method

The present paper is made necessary by the publication of the foregoing paper in this issue by Kolesnichenko [1]. It considers the basic propositions of the generalized Chapman-Enskog method and analyzes the arguments put forward by Kolesnichenko [1] and the validity of the method. The position of the results obtained by Kolesnichenko [14–17] is indicated. Nonequilibrium flows of multiatomic gases in which there occur processes of exchange of internal energy of the molecules in collisions between them and chemical reactions (such processes are called inelastic) are encountered frequently in nature and technology. It is therefore naturally of interest to derive gas-dynamic equations for such flows. The methods of the kinetic theory of gases were first used to obtain equations describing the limiting cases of very fast inelastic processes that take place in times of the order of the molecule-molecule collision times (equilibrium case) and very slow inelastic processes that take place over times of the order of the characteristic flow time (relaxation case). In [2–5], an algorithm was proposed for deriving gas-dynamic equations valid for arbitrary ratios of the rates of the elastic and inelastic processes and reducing to the well-known equations for the limiting cases already mentioned. The algorithm is called the generalized Chapman-Enskog method (abbreviated to the generalized method). The development, modification, and analysis of its properties can be found in [4, 6–13]. In [1], Kolesnichenko has questioned the validity of this algorithm.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 3, pp. 126–136, May–June, 1984.We thank V. A. Rykov for helpful and constructive discussions of the work.     

On the accuracy and efficiency of CFD methods in real gas hypersonics

A study of viscous and inviscid hypersonic flows using generalized upwind methods is presented. A new family of hybrid flux-splitting methods is examined for hypersonic flows. The hybrid method is constructed by the superposition of the flux-vector-splitting (FVS) method and second-order artificial dissipation in the regions of strong shock waves. The conservative variables on the cell faces are calculated by an upwind extrapolation scheme to third-order accuracy. A second-order-accurate scheme is used for the discretization of the viscous terms. The solution of the system of equations is achieved by an implicit unfactored method. In order to reduce the computational time, a local adaptive mesh solution (LAMS) method is proposed. The LAMS method combines the mesh-sequencing technique and local solution of the equations. The local solution of either the Euler or the NAVIER-STOKES equations is applied for the region of the flow field where numerical disturbances die out slowly. Validation of the Euler and NAVIER-STOKES codes is obtained for hypersonic flows around blunt bodies. Real gas effects are introduced via a generalized equation of state.     

Multicomponent Mixture Flow in an Axisymmetric Capillary in the Presence of a Mobile Surface Phase

Within the framework of the density-functional method, the hydrodynamics of a multicomponent mixture in the presence of a mobile surface phase is investigated. A system of interrelated three- and two-dimensional hydrodynamic equations for isothermal volume and surface flows is formulated. For flow in a thin axisymmetric capillary, the principal term of the asymptotic expansion of the solution in powers of the characteristic capillary radius to length ratio is obtained. For slow motion, the solution is found in quadratures.     

Radiative heat transfer in a hypersonic gas flow generated a spherical source past a plane

An asymptotic solution of radiative gas-dynamic equations for stationary interaction of two hypersonic gas flows emanating from two identical spherical sources is obtained. Under the assumption that the gas in the shock layer is in local thermodynamic equilibrium and volume emission (energy loss for radiation) occurs there, analytical expressions for the distributions of gas-dynamic functions and temperature are derived. The shock wave shape and the radiant flux on the contact plane are examined as functions of problem parameters.     

Three-dimensional non-isothermal extrusion flows

A three-dimensional (3-D) non-isothermal study of viscous free-surface flows with exponential dependence of viscosity on temperature is presented. The effects of non-isothermal conditions and/or geometry on the extrudate shape are investigated with a fully three-dimensional finite element/Galerkin formulation. Apart from the well known thermally induced extrudate swelling phenomenon, bending and distortion of the extrudate may occur because of temperature differences and/or geometric asymmetries. A temperature difference across the die can be imposed by heating or cooling the die walls, but can also arise because of asymmetric viscous heat generation due to the die geometry. Temperature differences affect velocity profiles because of the temperature dependence of viscosity and lead to extrudate bending, an effect known as kneeing in the fiber spinning industry. It is also shown numerically and confirmed experimentally that the die geometry induces extrudate bending even in the case of isothermal Newtonian flows.     

On basic equations of turbulent swirling gas-solid flows and their application in cyclones

The basic equations of turbulent gas-solid flows are derived by using the pseudo-fluid model of particle phase with a refined two-phase turbulence model. These equations are then applied to swirling gas-particle flows for analyzing the collection efficiency in cyclone separators.     

Investigation of slow monatomic gas flows past a circular cylinder

Slow low-Knudsen-number monatomic-gas flow past a circular cylinder is numerically investigated on the basis of a model kinetic equation. The gas flow is described by a new kinetic equation, from which the continuum equations for slow nonisothermal gas flows containing temperature stresses follow rigorously. It is shown that a closed convective-flow region arises near a nonuniformly heated cylinder in a slow gas flow if the flow impinges on the hot side of its surface. Using a new model of the Boltzmann equation makes it possible to study gas flows both in continuum and rarefied flow regimes.     

On Dissipative Fluxes in a Mixture of Polyatomic Gases

From the system of linearized kinetic equations, the diffusion and heat fluxes and the stress tensor for slow steady-state flows of a mixture of polyatomic gases are obtained and the transfer coefficients are found.     

Regular submodels of types (1,2) and (1,1) of the equations of gas dynamics

Partially invariant solutions of types (1, 2) and (1, 1) for gas-dynamic equations are regularly divided into two classes: for the first class, the invariant independent variable is the time, i.e., this class contains barochronic solutions, and for the second class, the invariant variable necessarily depends on spatial coordinates. The barochronic submodel of gas-dynamic equations, as well as a passive subsystem for solutions of the second class, is integrated in finite form. In the latter case, the invariant subsystem is reduced to an ordinary differential equation and quadratures. Integration of the submodels is illustrated by a number of examples. The following common properties of barochronic gas flows are described: rectilinear trajectories of gas particles, the possibility of collapse of density on a manifold, and stratification of the space of events. Lavrent’ev Institute of Hydrodynamics, Siberian Division, Russian Academy of Sciences, Novosibirsk 630090. Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 40, No. 2, pp. 40–49, March–April, 1999.     

Macroscopic thermodynamics of flowing polymeric liquids

The thermodynamics and mechanics of non-isothermal polymeric fluids are examined within the auspices of a new methodology wherein the laws of physics and principles of mechanics which are applicable to these thermodynamic systems are imbedded in a definite mathematical structure of a general, abstract equation. Such a concept allows new insight to be obtained concerning some aspects of non-isothermal flows of polymeric fluids, and permits a consistent expression and interpretation of other thermodynamic theories for these systems which have been developed over the past forty years. A major portion of this article is devoted to demonstrating the above statements, and in so doing some common misconceptions occurring in a significant fraction of the literature regarding this subject are exposed. The definite mathematical structure of the new methodology permits the thermodynamically consistent generalization of isothermal, incompressible models of polymeric fluids to non-isothermal, compressible conditions. Doing thus reproduces, corrects, and extends non-isothermal models which have been developed over the years, and also allows for simpler (but equivalent) representations of these models in terms of alternate variables with a clearer connection to the microstructure of the material than the stress tensor and heat flux vector fields. Furthermore, a generalization of the GENERIC structure is proposed that accommodates interactions between phenomena of differing parities, which impose antisymmetry upon the corresponding elements of the dissipative operator matrix. Received: 28 July 1998 Accepted: 6 November 1998     

Steady-State MHD Flows in Nozzles with an External Longitudinal Magnetic Field

Steady-state MHD flows in channels of the nozzle type in the presence of an external longitudinal magnetic field can be divided into two significantly different classes. Subcritical flows, in which the Alfvén velocity calculated from the longitudinal magnetic field is less than the plasma velocity, have mainly the same properties as flows in a transverse magnetic self-field and their quantitative characteristics depend only slightly on the longitudinal magnetic field strength. Supercritical flows with the opposite inequality for the velocities correspond to strong longitudinal magnetic field. The main difference is the transitions between different forms of energy (kinetic, thermal, and electromagnetic). The present study contains a classification of possible flows, namely, sub- and supercritical and sub-, super-, and transonic flows with respect to the fast and slow magnetosonic and Alfvén velocities. Examples of these flows are given. The effect of the problem parameters on the flow properties is investigated.     

Advances in the study of separated flows

This paper consists of three parts. The first deals with the separation conditions for three dimensional steady viscous separated flows, in which the behaviour of separated flow described by Navier-Stokes equations and boundary layer equations is studied. The second part involves an application of differential topology to qualitative analysis of flow fields. Here the distribution rule of singular points on the separation line is studied. The last part discusses the numerical method solving Navier-Stokes equations for separated flows. The obtained computational results are analysed by the above mentioned theories and methods.     

三角翼上分离及涡流的数值模拟

综述了三角翼(包括双三角翼,边条-三角翼,近距耦合鸭翼)上分离及涡流问题的数值模拟进展,介绍了用Euler方程、N-S方程的数值方法模拟不同三角翼,在不同攻角、不同来流M数等多种条件下的复杂涡流形态和流动结构,研究了控制方程、网格、湍流模型和计算方法等对计算结果的影响。      

Buoyancy-induced flow of non-Newtonian fluids over a non-isothermal body of arbitrary shape in a fluid-saturated porous medium

The buoyancy-induced flows of non-Newtonian fluids over non-isothermal bodies of arbitrary shape within saturated porous media have been treated using the boundary layer approximations and the power-law model to characterize the non-Newtonian fluid behavior. Upon introducing a general similarity transformation which considers both the geometrical effect and the wall temperature effect on the development of the boundary layer length scale, the governing equations for a non-isothermal body of arbitrary shape have been reduced to those for a vertical flat plate. The transformed equations reveal that a plane or axisymmetric body of arbitrary shape possesses its corresponding family of the wall temperature distributions which permit similarity solutions. Numerical integrations were carried out using the Runge-Kutta-Gill method, and the results of the heat transfer function were presented once for all plane and axisymmetric bodies. As illustrations, local wall heat flux distributions were discussed for wedges, cones, spheres, circular cylinders and other geometries. Furthermore, an approximate formula based on the Karman-Pohlhausen integral relation has been presented for speedy and sufficiently accurate estimation of heat transfer rates.     

A practical implementation of turbulence models for the computation of three-dimensional separated flows

An upwind MUSCL-type implicit scheme for the three-dimensional Navier-Stokes equations is presented and details on the implementation for three-dimensional flows of a ‘diagonal’ upwind implicit operator are developed. Turbulence models for separated flows are also described with an emphasis on the numerical specificities of the Johnson-King non-equilibrium model. Good predictions of separated two- and three-dimensional flows are demonstrated.     

Possibilities and Limitations of Computer Simulations of Industrial Turbulent Dispersed Multiphase Flows

We give an overview on the usage of computer simulations in industrial turbulent dispersed multiphase flows. We present a few examples of industrial flows: bubble columns and bubbly pipe flows, stirred tanks, cyclones, and a fluid catalytic cracking unit. The fluid catalytic cracking unit is used to illustrate the complexity of the physical phenomena involved, and the possibilities and limitations of the different approaches used: Eulerian–Lagrangian (particle-tracking) and Eulerian–Eulerian (two-fluid). In the first approach, the continuous phase is solved using either RANS simulations (Reynolds-Averaged Navier–Stokes simulations) or DNS/LES (Direct Numerical Simulations/Large-Eddy Simulations), and the individual particles are tracked. In the second approach, the dispersed phase is averaged, leading to two sets equations, which are quite similar to the RANS equations of single-phase flows. The Eulerian–Eulerian approach is the most commonly used in industrial applications, however, it requires a significant amount of modelling. Eulerian–Lagrangian RANS can be simpler to use; in particular in situations involving complex boundary conditions, polydisperse flows and agglomeration/breakup. The key issue for the success of the simulations is to have good models for the complex physics involved. A major weakness is the lack of good models for: the turbulence modification promoted by the particles, the inter-particle interactions, and the near-wall effects. Eulerian–Lagrangian DNS/LES can play an important role as a research tool, in order to get a better physical understanding, and to improve the models used in the RANS simulations (either Eulerian–Eulerian or Eulerian–Lagrangian).     

Numerical simulation of the vertical structure of discontinuous flows

A numerical method to solve the Reynolds‐averaged Navier–Stokes equations with the presence of discontinuities is outlined and discussed. The pressure is decomposed into the sum of a hydrostatic component and a hydrodynamic component. The numerical technique is based upon the classical staggered grids and semi‐implicit finite difference methods applied for quasi‐ and non‐hydrostatic flows. The advection terms in the momentum equations are approximated in order to conserve mass and momentum following the principles recently developed for the numerical simulation of shallow water flows with large gradients. Conservation of these properties is the most important aspect to represent near local discontinuities in the solution, following from sharp bottom gradients or hydraulic jumps. The model is applied to reproduce the flow over a step where a hydraulic jump forms downstream. The hydrostatic pressure assumption fails to represent this type of flow mainly because of the pressure deviation from the hydrostatic values downstream the step. Fairly accurate results are obtained from the numerical model compared with experimental data. Deviation from the data is found to be inherent to the standard k–ε model implemented. Copyright © 2001 John Wiley & Sons, Ltd.