Extremal nozzle contours for gas flows with particle lag

The determination of the extremal nozzle contour for gas flow without foreign particles has been carried out in several studies [1–6], based on the calculation of the flow field using the method of characteristics.In [7, 8] the equations are derived for the characteristics and the relations along the streamlines which are required for calculating two-dimensional gas flow with foreign particles. The variational problem for two-phase flow in the two-dimensional formulation may be solved by the method of Guderley and Armitage [9] with the use of equations given in [7] or [8]; however this method is very tedious, even with the use of high-speed computers.In [10, 11] studies are made of two-phase one-dimensional flows by expanding the unknown functions in series in a small parameter, defined by the particle dimensions. In [12] a solution is given for the variational problem (in the one-dimensional formulation) of designing the contour of a nozzle with maximal impulse. However that study does not take account of the static term appearing in the impulse and the solution is obtained in relative cumbersome form. Moreover, the question of account for the losses due to nonparallelism and nonuniformity of the discharge was not considered.The present paper considers in the one-dimensional formulation the flow of a two-phase medium in a Laval nozzle with small particle lags (in velocity and temperature). The variational problem of determining the maximal nozzle impulse is formulated along the nozzle contour for fixed geometric expansion ratio. The impulse losses due to nonparallelism of the discharge are simulated by a function which depends on the ordinates which are variable along the contour and on the slope of the tangent to the contour.The author wishes to thank Yu. D. Shmyglevskii and A. N. Kraiko for helpful discussions and V. K. Starkov for carrying out the calculations on the computer.     

One class of self-similar flows of a radiating gas

This article discusses self-similar statements of the problem of the motion of a completely radiating and absorbing gas. The field of radiation is assumed to be quasi-steady-state, and the contribution of the radiation to the internal energy, as well as the pressure and the viscosity of the medium, are not taken into account. The presence of local thermodynamic equilibrium is assumed. The absorption coefficient is approximated by a power function of the pressure and the density. Scattering of the radiation is not taken into account. Under these assumptions, there exist self-similar statements of the problem for one-dimensional unsteady-state flows (a strong detonation, the problem of plug-flow, motion under the effect of a radiation source, and others) and two-dimensional steady-state flows (flow in a diffuser, supersonic flow around a wedge or a cone). It is shown that there exists a non steady-state spherically symmetrical flow depending on four parameters; this flow is adiabatic in spite of the presence of radiation. This article is made up of seven sections. It is shown in the first section that the presence of radiation leads to the appearance of new dimensional constants, entering into the equations of the problem. The second section is devoted to self-similar nonsteady-state one-dimensional flows. The third section contains a detailed study of one class of such flows. In a partial case, adiabatic flows of a radiating gas are obtained. In the fourth and fifth sections, a detailed analysis is made of the initial and boundary conditions from the point of view of dimensionality. The sixth section describes self-similar two-dimensional steady-state flows of a radiating-absorbing gas. The seventh section consists of remarks with respect to approximations of the transfer equation.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 4, pp. 8–22, July–August, 1970.     

稠密颗粒床中点火传火过程的两维两相流数值模拟

本文结合火炮工程背景,建立伴随化学反应的稠密颗粒群气固两相非定常两维轴对称流和点火器内的非定常一维两相流数学模型,并考虑点火器对主装药床的耦合作用。采用MacCormack两步差分格式,模拟点火传火过程中火焰传播及其物理量沿轴向和径向变化规律。部分计算结果与实验结果相吻合。发现点火开始阶段存在明显的径向效应,当火焰波传到固壁后,轴向流占主要优势。     

Complementarity of CFD, experimentation and reactor models for solving challenging fluidization problems

Experimentalists, numerical modellers and reactor modellers need to work together, not only just for validation of numerical codes, but also to shed fundamental light on each other＇s problems and underlying assumptions. Several examples are given, Experimental gas axial dispersion data provide a means of choosing the most appropriate boundary condition （no slip, partial slip or full slip） for particles at the wall of fluidized beds. CFD simulations help to identify how close ＂two-dimensional＂ experimental columns are to being truly two-dimensional and to representing three-dimensional columns. CFD also can be used to provide a more rational means of establishing assumptions needed in the modelling of two-phase fluidized bed reactors, for example how to deal with cases where there is a change in molar flow （and hence volumetric flow） as a result of chemical reactions.     

Fixed stream-tube method for solving two-phase plane flow problems and its theoretical analysis

The fixed stream-tube method widely adopted in engineering field for giving an approximate solution to the two-dimensional problems of two-phase flow through porous media is summarized and an improvement has been made in this paper. Its core part, i.e., the fluid displacement within a one-dimensional stream tube with variable cross-sectional area under a given pressure difference across the tube is thoroughly studied. The existence and uniqueness of solution are proved, the exact solution, numerical solution and its convergence, stability analyses are given in this paper.     

Visualisation of air–water bubbly column flow using array Ultrasonic Velocity Profiler

In the present work, an experimental study of bubbly two-phase flow in a rectangular bubble column was performed using two ultrasonic array sensors, which can measure the instantaneous velocity of gas bubbles on multiple measurement lines. After the sound pressure distribution of sensors had been evaluated with a needle hydrophone technique, the array sensors were applied to two-phase bubble column. To assess the accuracy of the measurement system with array sensors for one and two-dimensional velocity, a simultaneous measurement was performed with an optical measurement technique called particle image velocimetry(PIV). Experimental results showed that accuracy of the measurement system with array sensors is under 10% for one-dimensional velocity profile measurement compared with PIV technique. The accuracy of the system was estimated to be under 20% along the mean flow direction in the case of two-dimensional vector mapping.     

Numerical simulation of 1-D unsteady two-phase flows with shocks

In the present paper, random-choice method (RCM) and second-order GRP difference method, which are high resolution methods used for pure gas flows with shocks, are extended and employed to study the problem of one-dimensional unsteady two-phase flows. The two-phase shock wave and the flow field behind it in a dusty gas shock tube are calculated and the time-dependent change of the fiow parameters for the gas antiparticle phase are obtained. The numerical results indicate that both the two methods can give the relaxation structure of the two-phase shocks with a sharp discontinuous front and that the GRP method has the advantages of less time-consuming and higher accuracy over the RCM method.     

Bifurcation phenomena in two-phase natural circulation

Multiple steady-state solutions have been found theoretically for a two-phase natural circulation loop. This bifurcation phenomenon is attributed to the nonmonotonic behavior of the buoyancy and friction forces, and the discontinuities due to transitions between two-phase flow patterns. The paper outlines a general and consistent modeling method to describe two-phase natural circulation at steady state. The solution of the one-dimensional governing conservation equations is performed by a combination of analytical and numerical stages. The method is implicit, to treat the highly nonlinear problem by multinested iteration loops. Results are presented and discussed for the multiple solutions of the flow rate, temperature, quality and void distributions and two-phase elevations. Pressure effects are also studied. A comparison with available data shows that the theoretical model simulates the important characteristics of two-phase thermosyphons.     

Solution of the direct problem of a laval nozzle working with a two-phase medium

A study is made of a method of numerical solution of the system of ordinary differential equations describing the flow of a two-phase medium in a Laval nozzle in the one-dimensional approximation. The paper is concerned with the direct problem, in which the law of variation of the area of the transverse section of the nozzle and functional relationships between the entrance parameters are given and the unknown parameters of the two-phase flow are found along the length of the nozzle.     

A method for calculating the two-dimensional temperature field of a reservoir with thermal injection

In [1–3] and other studies devoted to the determination of the temperature field of an oil reservoir when injecting into it a fluid with a temperature different from that of the initial reservoir temperature, the one-dimensional fluid flow (linear or radial) was considered in the case of an injection gallery or a single injection well in the reservoir. The problem was formulated in [4] with an arbitrary distribution of wells, but the solution was obtained only for the integral thermal-loss characteristic. To evaluate the coverage of the reservoir by the thermal effect, we must know the temperature distribution in the multi-well reservoir system at any instant of time. In this paper we propose a method for calculating the reservoir temperature field in the case of two-dimensional fluid flow using the simplifying assumptions which were used earlier by Lauwerier [1] and other authors to describe the thermal phenomena in a reservoir.     

One-, two-, and three-dimensional solutions for counter-current steady-state two-phase subsurface flow

This paper derives analytical solutions for steady-state one-dimensional (1-D), two-dimensional (2-D), and three-dimensional (3-D) two-phase immiscible subsurface flow for a counter-current problem. Since the governing equations are highly nonlinear, 2-D and 3-D derivations are generally difficult to obtain. The primary purpose for the solutions is to test finite difference/volume/element computer programs for accuracy and scalability using architectures ranging from PCs to parallel high performance computers. This derivation is accomplished by first solving for the saturation of water in terms of a function that is a solution to Laplace’s equation to achieve a set of partial differential equations that allows some degree of latitude in the choice of boundary conditions. Separation of variables and Fourier series are used to obtain the final solution. The test problem consists of a rectangular block of soil where specified pressure is applied at the top and bottom of the sample, and no-flow boundary conditions are imposed on the sides. The pressure at the top of the sample is a step function that allows the testing of adaptive meshing or concentration of grid points in action zones.     

Flux difference splitting for open-channel flows

A finite difference scheme based on flux difference splitting is presented for the solution of the one-dimensional shallow-water equations in open channels, together with an extension to two-dimensional flows. A linearized problem, analogous to that of Riemann for gas dynamics, is defined and a scheme, based on numerical characteristic decomposition, is presented for obtaining approximate solutions to the linearized problem. The method of upwind differencing is used for the resulting scalar problems, together with a flux limiter for obtaining a second-order scheme which avoids non-physical, spurious oscillations. The scheme is applied to a one-dimensional dam-break problem, and to a problem of flow in a river whose geometry induces a region of supercritical flow. The scheme is also applied to a two-dimensional dam-break problem. The numerical results are compared with the exact solution, or other numerical results, where available.     

Bingham Viscoplastic as a Limit of Non-Newtonian Fluids

A new formulation is proposed for the equations of the Bingham viscoplastic. Global existence of x--periodic solutions is proved. A uniqueness theorem is established in the two-dimensional case. A relation with the G. Duvaut--J. L. Lions variational inequality is discussed, and a result on equivalence is obtained. The question of interaction between fluid-rigid phases is studied when the initial state is rigid. A free-boundary problem that describes two-phase one-dimensional flows is considered.     

On the various forms of the conservation equations in two-phase flow

The conservation equations for one-dimensional two-phase flow are derived from first principles. The effects of the radial distributions of velocities, enthalpies, and void fraction are taken into account through the use of correlation coefficients. Several simplified separated-flow model formulations that have appeared in the literature are derived from these equations by specializing the values of the correlation coefficients. The equivalence of these formulations under certain assumptions is demonstrated. Finally, new Lagrangian forms of the conservation equations, written in terms of the velocities of the center of mass, momentum, and energy are presented.     

Dynamics of a liquid film sheared by a co-flowing turbulent gas

Consider the dynamics of a thin laminar liquid film flowing over an inclined wall in the presence of a co-flowing turbulent gas. The solution to the full two-phase flow problem poses substantial technical difficulties. However, by making appropriate assumptions, the solution process can be simplified and can provide valuable insights. The assumptions allow us to solve the gas and liquid problems independently. Solving for the gas flow reduces to finding perturbations to pressure and tangential stresses at the interface, influencing the liquid problem through the boundary conditions. We analyze the effect of gas flow on the liquid problem by developing an integral-boundary-layer model, which is valid up to moderate liquid Reynolds numbers. We seek solitary-wave solutions of this model under the influence of gas flow via a pseudo-arclength continuation method. Our computations demonstrate that as a general trend, the wave speed increases with increasing the gas shear and the liquid flow rate. Further insight into the problem is provided via time-dependent computations of the integral-boundary-layer model.     

Stability of a laminar film flow

The problem of the wave motion of a liquid layer was first investigated by Kapitsa [1, 2], who gave an approximate analysis of the free flow and flow in contact with gas stream, and evaluated the influence of the heat transfer processes on the flow. The problem of the stability of such a flow was studied in detail by Benjamin [3] and Yih [4, 5], These authors proposed seeking the solution of the resulting Orr-Sommerfeld equation in the form of a series in a small parameter and developed a corresponding method of successive approximations. As the small parameter [3–5], they made use of the product of the disturbance wave number and the Reynolds number. In these studies, the tangential stress on the free surface was taken equal to zero, and the fluid film was always considered essentially plane. At the same time, there are certain types of problems of considerable interest in which neither of these assumptions is satisfied. A good example might be the problem on the stability of the annular regime of two-phase flow in pipes and capillaries, when the basic stream of one fluid is separated from the pipe walls by an annular layer of another fluid. In this case, the interface has a finite radius of curvature and the tangential stress on the interface may be significantly different from zero.In the present paper, the problem of the flow stability of a fluid layer with respect to small disturbances of the boundary surface is considered with account for both the finite radius of curvature of the boundary surface and the nonzero hydrodynamic friction at the boundary. The film is assumed to be quite thin. This enables us, firstly, to consider the Reynolds number small, to use the general method of [5], and, second ly, to consider the film thickness sufficiently small in comparison with the radius of curvature of the substrate on which the film lies. Furthermore, for evaluating the stability of the laminar flow of the curved film we can use the results obtained for a plane film with account for the terms which depend on the curvature of the substrate.As a rule, previous studies have considered only one-dimensional disturbances of the boundary surface. In the present paper, in the first approximation, the stability is examined in relation to two-dimensional disturbances of this surface, corresponding to three-dimensional flow disturbances.As an example, the results obtained are applied to the investigation of the stability of the free flow of a layer of fluid over an inclined plane under the sole influence of gravity.     

Modeling of retrograde condensation in the process of steady radial flow through a porous medium

The problem of the steady axisymmetric two-phase flow of a multicomponent mixture through a porous medium with phase transitions is considered. It is shown that the system of equations for the two-phase multicomponent flow process, together with the equations of phase equilibrium, reduces to a system of two ordinary differential equations for the pressures in the gas and liquid phases. A family of numerical solutions is found under certain assumptions concerning the pressure dependence of the molar fraction of the liquid phase.Moscow. Translated from Izvestiya Rossiiskoi Akademii Nauk, Mekhanika Zhidkosti i Gaza, No. 6, pp. 92–97, November–December, 1994.     

Pressure drop caused by flow area changes in capillaries under low flow conditions

Single-phase and two-phase flow pressure drops caused by flow area expansion and contraction were measured using air and water. The test section consisted of two capillaries with 0.84 mm and 1.6 mm diameters. For single-phase flow, the Reynolds numbers defined based on the smaller diameter capillary covered the range 160–11,000. For two-phase flow, the all-liquid Reynolds number based on the smaller capillary varied in the 410–1020 range, and the flow quality varied in the 0.018–0.2 range. The single-phase flow loss coefficients for both flow area expansion and contraction were empirically correlated. For two-phase flow, the data indicated the occurrence of significant velocity slip, and the one-dimensional homogeneous flow model utterly disagreed with the data. For flow area expansion the one-dimensional slip flow model along with an Armand-type slip ratio correlation could predict the data well. For flow area contraction, the one-dimensional slip flow model along with the slip ratio expression of Zivi agreed with the data very well, provided that no vena-contracta was considered.     

Theoretical and experimental investigation of condensation in a centered rarefaction wave

The analysis of the process of spontaneous condensation in one-dimensional formulation is dealt with adequately in many papers. However, in reality supersonic flows are not one-dimensional. The most striking effect of two-dimensionality is manifested in two-phase flows, for example in nozzles, inclined sections of jet turbine grills and rarefaction waves. The investigation of these flows, both in the experimental and theoretical aspect, is a complex problem for which a solution has been found only recently. The results are given in this paper of a theoretical and experimental investigation of spontaneous condensation of water vapor in a centered rarefaction wave formed by flow around a protuberant angle by a hypersonic stream.Translated from Prikladnaya Mekhanika i Tekhnicheskii Fiziki, No. 5, pp. 117–122, September–October, 1971.