Effective Parameter Interpretation and Extrapolation of Dispersion Simulations by Means of a Simple Two-Velocity Model

The investigation of dispersion by microscopic simulations yields a lot of detailed information. To identify characteristic behaviours, it is useful to condense this information into a few effective parameters, which describe the transport process in the model geometry on a larger scale. For this purpose, a very simple two-velocity model has been developed, which models the transition from reversible to irreversible spreading of a tracer volume. It is shown that this model is very similar to Taylor–Aris dispersion and that it is quite suitable to approximate the time dependence of dispersion. The model is applied to characterize the effect of dead end pores on dispersion with a single correlation parameter. Up to Péclet numbers of about 500, 'hold-up'-dispersion similar to Taylor–Aris-dispersion is found. The simulations have been performed by the lattice Bhatnagar–Gross–Krook (BGK) method, which is a particular type of cellular automata and therefore allows an easy implementation of complicated geometries. The fully irreversible asymptotic dispersion is reached in an exponential process, the parameters of which can be identified by the two-velocity model after the mixing has noticeably begun. These are used to extrapolate the process which reduces the computational effort by about one order of magnitude.     

Solution of the problem of a flow of rarefied gas of constant density around a semiinfinite plate by the integral diffusion method

Several papers [1–4] have proposed approximate diffusion models which can be used to examine the transport process in a rarefied gas where the mean free path is large and transport is not determined by the local gradient of the particular quantity.In this paper the integral diffusion model [2] is used to solve the problem of determination of the friction stress and velocity of a flow of an incompressible gas around a plane semi-infinite plate in the whole range of Knudsen numbers. The obtained solution is compared with published solutions and experimental data [9].     

Stability of flows of the Hamel type in channels with permeable walls

The stability of nonparallel flows of a viscous incompressible fluid in an expanding channel with permeable walls is studied. The fluid is supplied to the channel through the walls with a constant velocity v₀ and through the entrance cross section, where a Hamel velocity profile is assigned. The resulting flow in the channel depends on the ratio of flow rates of the mixing streams. This flow was studied through the solution of the Navier—Stokes equations by the finite-difference method. It is shown that for strong enough injection of fluid through the permeable walls and at a distance from the initial cross section of the channel the flow approaches the vortical flow of an ideal fluid studied in [1]. The steady-state solutions obtained were studied for stability in a linear approximation using a modified Orr—Sommerfeld equation in which the nonparallel nature of the flow and of the channel walls were taken into account. Such an approach to the study of the stability of nonparallel flows was used in [2] for self-similar Berman flow in a channel and in [3] for non-self-similar flows obtained through a numerical solution of the Navier—Stokes equations. The critical parameters *, R*, and Cr* at the point of loss of stability are presented as functions of the Reynolds number R₀, characterizing the injection of fluid through the walls, and the parameter , characterizing the type of Hamel flow. A comparison is made with the results of [4] on the stability of Hamel flows with R₀ = 0.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 6, pp. 125–129, November–December, 1977.The author thanks G.I. Petrov for a discussion of the results of the work at a seminar at the Institute of Mechanics of Moscow State University.     

Convective and absolute instabilities in counter-rotating spiral Poiseuille flow

We present results of an experimental study on the stability of Taylor–Couette flow in case of counter-rotating cylinders and an imposed axial through flow. We are able to confirm results form recent numerical investigations done by Pinter et al. [24] by measuring the absolute and convective stability boundaries of both propagating Taylor vortices (PTV) and spiral vortices (SPI). Thus our work shows that these theoretical concepts from hydrodynamic stability in open flows apply to experimental counter-rotating Taylor–Couette systems with an imposed axial through flow. PACS 47.20.-k, 05.45.-a, 47.15.fe     

Experimental investigation of flow in a trench

An experimental investigation was made of the flow of a viscous incompressible liquid in a trench of square transverse cross section, using a laser Doppler velocimeter. The investigation was made with two values of the Reynolds number Re, corresponding to laminar and turbulent flow conditions in the channel. The experimental data show that a core with a constant vorticity is formed in the trench, that a jet propagates near the walls of the trench, and that there are secondary eddies in the corners of the trench. The motion of a viscous liquid in a trench of rectangular cross section is part of a broad class of breakaway flows. Experimental data on the investigation of flow in trenches are extremely few. A majority of the existing information is limited to visual observations [1–4]. In [2, 5, 6] the question of the unstable character of flow in trenches was discussed. Quantitative measurements of stable eddy flows in trenches were made in [7–9] using a thermoanemometer, and in [7] measurements were made of the pressure at the bottom and walls of trenches; there are data on the distribution of the velocity in the middle sections of trenches. In [8] the mean velocity, the intensity of the turbulence, and the stress of the turbulent flow were obtained in several sections parallel to the side walls of the trench, In [9] a measurement was made of the velocities also in two cross sections of a trench in which one component of the velocity prevails. A brief analysis of the existing experimental results shows that these data are insufficient to form a detailed representation of the character of flow in a trench.Translated from Zhurnal Prikladnoi Mekhaniki i Tekhnicheskoi Fiziki, No. 2, pp. 76–86, March–April, 1976.     

Longitudinal dispersion associated with the motion of an adsorbed mixture in a column filled with porous granules

The effect of the radial nonuniformity of the linear flow velocity on the radial nonuniformity of the adsorbate concentration distribution over the cross section of a column and on the longitudinal dispersion is estimated for the diffusion model, and the legitimacy (for relatively small adsorber radii) of using a one-dimensional diffusion model for describing the interphase mass transfer in cylindrical adsorbers is thereby determined.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 4, pp. 93–101, July–August, 1987.     

Sorbing and Non-Sorbing Solute Migration in Rough Fractures with a Multi-Species LGA Model: Dispersion Dependence on Retardation and Roughness

A multi-species Lattice-Gas cellular automaton model was applied to the study of the migration of sorbing and non-sorbing tracers in 2D smooth fractures and in a series of increasing roughness fractures. A tenfold increase of the dispersion of the non-sorbing tracer was calculated in the highest roughness case compared to the smooth fracture. Up to a threefold increase of the dispersion of the sorbing tracers was calculated compared to the non-sorbing tracer. These enhanced dispersions were found to be of a Fickian form and were interpreted in terms of the classical Taylor–Aris dispersion. The effects of roughness and retardation over the increase of dispersion were identified and quantified through a semi-empirical relation. These effects were found additive and independent.     

Direct numerical simulation of turbulent Taylor–Couette flow

The direct numerical simulation (DNS) of the Taylor–Couette flow in the fully turbulent regime is described. The numerical method extends the work by Quadrio and Luchini [M. Quadrio, P. Luchini, Eur. J. Mech. B/Fluids 21 (2002) 413–427], and is based on a parallel computer code which uses mixed spatial discretization (spectral schemes in the homogeneous directions, and fourth-order, compact explicit finite-difference schemes in the radial direction). A DNS is carried out to simulate for the first time the turbulent Taylor–Couette flow in the turbulent regime. Statistical quantities are computed to complement the existing experimental information, with a view to compare it to planar, pressure-driven turbulent flow at the same value of the Reynolds number. The main source for differences in flow statistics between plane and curved-wall flows is attributed to the presence of large-scale rotating structures generated by curvature effects.     

Calculation by the Krylov-Bogolyubov method of the velocity profile and fluxes in the nonisothermal flow of a rarefied gas in a cylindrical capillary

The Krylov-Bogolyubov numerical method is used to solve integral transfer equations obtained from the kinetic equation with the BHC (Bhatnagar—Cross—Krook) model of the collision operator. The velocity profiles and the thermal-creep flows and Poiseuille flow are calculated in different modes of flow under conditions of incomplete accommodation of the tangential momentum of molecules at the wall.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 6, pp. 143–150, November–December, 1978.     

Qualitative mathematical analysis of the Richards equation

The Richards equation is widely used as a model for the flow of water in unsaturated soils. For modelling one-dimensional flow in a homogeneous soil, this equation can be cast in the form of a specific nonlinear partial differential equation with a time derivative and one spatial derivative. This paper is a survey of recent progress in the pure mathematical analysis of this last equation. The emphasis is on the interpretation of the results of the analysis. These are explained in terms of the qualitative behaviour of the flow of water in an unsaturated soil which is described by the Richards equation.Nomenclature a coefficient in second-order diffusion term of equation - b coefficient in first-order advection term of equation - D soil-moisture diffusivity [L²T^-1] - h pressure head [L] - H quarter-plane domain for Cauchy-Dirichlet problem [L] x [T] - K hydraulic conductivity scalar [LT^–1] - K hydraulic conductivity tensor [LT^–1] - q soil-moisture flux scalar [LT^–1] - q soil-moisture flux vector [LT^–1] - r dummy variable - R rectangle [L] x [T] - s dummy variable - s* representative value of dummy variable - S half-plane domain for Cauchy problem [L] x [T] - t time [T] - u unknown solution of partial differential equation - u₀ initial-value function - v soil-moisture velocity scalar [LT^–1] - v soil-moisture velocity vector [LT^–1]     

Theory of vortical helical ideal gas flows in laval nozzles

An asymptotic solution is found for the direct problem of the motion of an arbitrarily vortical helical ideal gas flow in a nozzle. The solution is constructed in the form of double series in powers of parameters characterizing the curvature of the nozzle wall at the critical section and the intensity of stream vorticity. The solution obtained is compared with available theoretical results of other authors. In particular, it is shown that it permits extension of the known Hall result for the untwisted flow in the transonic domain [1]. The behavior of the sonic line as a function of the vorticity distribution and the radius of curvature of the nozzle wall is analyzed. Spiral flows in nozzles have been investigated by analytic methods in [2–5] in a one-dimensional formulation and under the assumption of weak vorticity. Such flows have been studied by numerical methods in a quasi-one-dimensional approximation in [6, 7]. An analogous problem has recently been solved in an exact formulation by the relaxation method [8, 9]. A number of important nonuniform effects for practice have consequently been clarified and the boundedness of the analytical approach used in [2–7] is shown.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 2, pp. 126–137, March–April, 1978.The authors are grateful to A. N. Kraiko for discussing the research and for valuable remarks.     

Quasihomogeneous model of cavitation flows in diffuser channels

Starting with the experiments carried out by Reynolds in 1894, the flow in Venturi tubes has traditionally been used to study and demonstrate various forms of cavitation. Numerous authors have carried out experimental research on the various flow regimes in diffuser channels [1–7] or have investigated theoretical models of such flows [6, 8]. The occurrence and development of cavitation is closely associated with the phenomenon of turbulent separation complicated by the presence of two-phase flow in the dissipation zone. For a long time these effects were considered separately, until Gogish and Stepanov [9] proposed a single model of cavitation and separation based on the theory of intense interaction of an incompressible potential flow and a turbulent cavitation layer of variable density and embracing the various stages of cavitation. The object of this study is to demonstrate the possibilities of this model with reference to the simple example of flows accompanied by cavitation and separation in plane and axisymmetric diffuser channels of the Venturi tube type with straight and curved walls. The dissipative flow near the walls is described by a quasihomogeneous model of turbulent two-phase flow, in which the presence of two phases is taken into account only by varying the mean density. The potential core of the flow is considered in the one-dimensional formulation. The displacement thickness serves as the flow interaction parameter. The conditions of ocurrence and development of circulatory flows are determined. Examples of symmetrical and nonsymmetrical flows are presented.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 5, pp. 47–54, September–October, 1986.     

Quasi-one-dimensional approximation in two-dimensional problems of gas dynamics

A useful means of constructing approximate flow models is the hydraulic (for two-dimensional problems quasi-one-dimensional) approach, based on averaging the initial nonuniform flows over some direction or cross section [1]. In this case, at the expense of a rougher model it is possible to reduce the dimensionality of the problem. Here, this approach is extended to unsteady two-dimensional gas-dynamic processes; certain problems (flow around a cone or a blunt body, jet flows) are considered in the framework of the quasi-one-dimensional model obtained, and results are compared with the solutions of the corresponding two-dimensional problems.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 2, pp. 136–143, March–April, 1989.     

Detonation Wave Propagation in Rotational Gas Flows

This paper studies the propagation of detonation and shock waves in vortex gas flows, in which the initial pressure, density, and velocity are generally functions of the coordinate — the distance from the symmetry axis. Rotational axisymmetric flow having a transverse velocity component in addition to a nonuniform longitudinal velocity is considered. The possibility of propagation of Chapman–Jouguet detonation waves in rotating flows is analyzed. A necessary conditions for the existence of a Chapman–Jouguet wave is obtained.     

Nonequilibrium stretching dynamics of dilute and entangled linear polymers in extensional flow

We propose an extension of the FENE-CR model for dilute polymer solutions [M.D. Chilcott, J.M. Rallison, Creeping flow of dilute polymer solutions past cylinders and spheres, J. Non-Newtonian Fluid Mech. 29 (1988) 382–432] and the Rouse-CCR tube model for linear entangled polymers [A.E. Likhtman, R.S. Graham, Simple constitutive equation for linear polymer melts derived from molecular theory: Rolie–Poly equation, J. Non-Newtonian Fluid Mech. 114 (2003) 1–12], to describe the nonequilibrium stretching dynamics of polymer chains in strong extensional flows. The resulting models, designed to capture the progressive changes in the average internal structure (kinked state) of the polymer chain, include an ‘effective’ maximum contour length that depends on local flow dynamics. The rheological behavior of the modified models is compared with various results already published in the literature for entangled polystyrene solutions, and for the Kramers chain model (dilute polymer solutions). It is shown that the FENE-CR model with an ‘effective’ maximum contour length is able to describe correctly the hysteretic behavior in stress versus birefringence in start-up of uniaxial extensional flow and subsequent relaxation also observed and computed by Doyle et al. [P.S. Doyle, E.S.G. Shaqfeh, G.H. McKinley, S.H. Spiegelberg, Relaxation of dilute polymer solutions following extensional flow, J. Non-Newtonian Fluid Mech. 76 (1998) 79–110] and Li and Larson [L. Li, R.G. Larson, Excluded volume effects on the birefringence and stress of dilute polymer solutions in extensional flow, Rheol. Acta 39 (2000) 419–427] using Brownian dynamics simulations of bead–spring model. The Rolie–Poly model with an ‘effective’ maximum contour length exhibits a less pronounced hysteretic behavior in stress versus birefringence in start-up of uniaxial extensional flow and subsequent relaxation.     

Simulation of ordered structures in open turbulent flows

Extensive experimental material [1–4] indicates that ordered (coherent) structures play an important part in determining the nature of the flow, the generation of Reynolds stresses and turbulence energy, and the transport of heat, momentum, and passive admixtures in a turbulent flow. In the present paper, a model is constructed for describing coherent structures in which, given the profile of the mean velocity, one can determine the characteristic sizes, the propagation velocities, and also the frequency and amplitude characteristics of these ordered motions. The model is based on the analogy between the ordered formations and secondary flows in a subsidiary laminar flow whose velocity profile is the same as the turbulent profile of the mean velocity. The influence of small-scale pulsations is described by the introduction of the coefficient of turbulent viscosity. In the framework of the model, numerical calculations are made for two-dimensional turbulent flows in a mixing layer, a jet, and a wake behind a cylinder. The results of the calculations are compared with experimental data.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 4, pp. 45–52, July–August, 1981.     

Computation of base flow characteristics with uniform blowing

Blowing at bluff body base was considered under different conditions and for small amount of blowing this problem was solved using dividing streamline model [1]. The effect of supersonic blowing on the flow characteristics of the external supersonic stream was studied in [2–4]. The procedure and results of the solution to the problem of subsonic blowing of a homogeneous fluid at the base of a body in supersonic flow are discussed in this paper. Analysis of experimental results (see, e.g., [5]) shows that within a certain range of blowing rate the pressure distribution along the viscous region differs very little from the pressure in the free stream ahead of the base section. In this range the flow in the blown subsonic jet and in the mixing zones can be described approximately by slender channel flow. This approximation is used in the computation of nozzle flows with smooth wall inclination [6, 7]. On the other hand, boundary layer equations are used to compute separated stationary flows with developed recirculation regions [8] in order to describe the flow at the throat of the wake. The presence of blowing has significant effect on the flow structure in the base region. An increasing blowing rate reduces the size of the recirculation region [9] and increases base pressure. This leads to a widening of the flow region at the throat, usually described by boundary-layer approximations. At a certain blowing rate the recirculation region completely disappears which makes it possible to use boundary-layer equations to describe the flow in the entire viscous region in the immediate neighborhood of the base section.Translated from Zhurnal Prikladnoi Mekhaniki i Tekhnicheskoi Fiziki, No. 1, pp. 76–81, January–February, 1984.     

Subsonic flow of a compressible gas around a semiinfinite flat plate in a pipe

The subsonic flow of an ideal compressible gas around the rear end of a semiinfinite flat plate in a pipe is considered. The flow pattern is similar to that assumed by Efros [1, 2] with a return stream for the cavitational flows. Fal'kovich's method [3] is used to solve the problem and this makes it possible to obtain the solution to the problems of the gas streams at several typical velocities. The method is a generalization of that of Chaplygin [4] for flow problems at one typical velocity.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 4, pp. 101–108, July–August, 1970.     

Transition flow of a fibrous suspension in a pipe as the flow of an anisotropic fluid

The transition flow is considered of a fibrous suspension in a pipe. The flow region consists of two subregions: at the center of the flow a plug formed by interwoven fibers and fluid moves as a rigid body; between the solid wall and the plug is a boundary layer in which the suspension is a mixture of the liquid phase and fibers separated from the plug [1–3]. In the boundary region the suspension is simulated as an anisotropic Ericksen—Leslie fluid [4, 5] which satisfies certain additional conditions. Equations are obtained for the velocity profile and drag coefficient of the pipe, which are both qualitatively and quantitatively in good agreement with the experimental results [6–8]. Within the framework of the model, a mechanism is found for reducing the drag in the flow of a fibrous suspension as compared to the drag of its liquid phase.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 5, pp. 91–98, September–October, 1985.     

Heat transfer enhancement with deformation of the velocity profile

The problem of enhancing the heat transfer in channels and boundary layers by the appropriate deformation of the fluid velocity profile is considered. The resulting additional hydraulic losses, the price of heat transfer enhancement, are determined. The possibilities of controlling heat transfer by redistributing the fluid velocity in channels are demonstrated with reference to flows at low Prandtl numbers. Laminar and turbulent liquid and gas flows with heat transfer in channels and boundary layers are numerically modeled on the basis of modern models of turbulence (flow development in channels with different initial velocity profiles, flows with wall roughness and boundary layer flows with forces acting on the flow to cause deformation of the velocity profile). In all cases it is found that the heat transfer can be enhanced only at the expense of a considerable increase in the hydaulic losses. A class of self-similar thermal problems for flows in plane diffusers is formulated. The eigenfunctions — temperature modes — for various velocity profiles are determined with allowance for the nonuniqueness of the solution of the classical dynamical problem for a plane diffuser and the corresponding heat transfer coefficients are found.Translated from Izvestiya Rossiiskoi Akademii Nauk, Mekhanika Zhidkosti i Gaza, No.4, pp. 94–105, May–June, 1993.The authors are grateful to A. Yu. Klimenko for useful discussions.