Viscous incompressible flow in a narrow channel with one free wall and fluid supply through a porous insert

The problem investigated relates the plane unsteady flow of a viscous incompressible fluid in a narrow channel one of whose walls is free and acted upon by a given load, while the other is rigidly fixed. The fluid enters the channel through a porous insert in the stationary wall. A model of the flow of a thin film of viscous incompressible fluid and Darcy's law for flow in a porous medium are used to find the distribution of fluid pressure and velocity in the channel and the porous insert in the two-dimensional formulation for fairly general boundary conditions in the case where the length of the porous insert exceeds the length of the free wall. In the particular case where the length of the porous insert is equal to the length of the free wall an exact stationary solution of the problem is obtained for a given value of the channel height. The stability of the equilibrium position of the free wall supported on a hydrodynamic fluid film is examined.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 1, pp. 16–24, January–February, 1986.     

Wave flow of a liquid film together with a flow of gas

The wave flow of a thin layer of viscous liquid in conjunction with a flow of gas was considered in a linear formulation earlier [1, 2]. In this paper the problem of the wave flow of a liquid film together with a gas flow is solved in a nonlinear setting. On this basis relationships are derived for calculating the parameters of the film and the hydrodynamic quantities.Ivanovo. Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 1, pp. 12–18, January–February, 1972.     

Stability of quasistatic equilibrium state of an inhomogeneous viscous layer on an inclined plane

A study is made of the stability against small perturbations [1] of a slow flow of an incompressible inhomogeneous linearly viscous liquid under the influence of a force of gravity on an unbounded inclined plane. Problems of such kind arise in glaciology when one estimates the stability of snow on mountain slopes or determines the catastrophic movement of a glacier; the results can also be applied to solifluction phenomena [2, 3]. Equations for perturbations of parallel flows of linearly viscous fluids in the case of a continuous variation of the viscosity and density across the flow were derived in [4]. An attempt to solve the hydrodynamic problem with allowance for a perturbation of the viscosity was made in [5]; however, in the equations for the perturbations, simplifications resulted in the neglect of terms that take into account perturbations of the viscosity. In the quasistatic formulation considered here in the case when allowance is made for perturbation of the density and viscosity, the equation for the perturbation amplitudes is an ordinary differential equation with variable coefficients; analytic solution of the equation is very difficult, even for long-wave perturbations. In this connection a study is made of the stability of a laminar model; the viscosity and density are constant within each layer. A similar hydrodynamic problem in the long-wave approximation was considered in [6]. In the present paper an exact solution is constructed in the quasistatic formulation for a two-layer model; the solution shows that the viscosity of the lower layer has an important influence on the stability. Essentially, instability is observed when the lower layer acts as a lubricant.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 6, pp. 20–24, November–December, 1973.     

OVERLAPPING GRIDS FOR THE SIMULATION OF PARTICLE INTERACTION AT THE MICRO SCALE

This paper describes the use of overlapping grids for the calculation of flow around single and multipleparticle configurations at the micro scale. The basic equations for calculation are those for conservation of mass and momentum which are solved using a common Finite-Volume formulation. The hydrodynamic particle-particle and particle-wall interaction can be calculated by using an overlapping or Chimera grid scheme. With the grid structuring procedure it is possible to use simple and structured grids around the particles and the overall main grid geometry. The particle grids are lapped over the main grid such that they can move independently after each time step without remeshing the whole geometry. The paper gives results for the validation of the code developed for general test cases, for a rotating ellipsoid in simple shear flow, the flow around particles attached to a wall, the motion of a particle in the vicinity of a wall and some results for the flow through a packed bed configuration.     

Exact solutions of boundary-value problems for nonlinear flow in porous media

Exact solutions for flow problems in porous media with a limiting gradient in the case when the flow region in the hodograph plane is a half-strip with a longitudinal cut [1] are known only for two models of the resistance law [2–6]. The present study gives a one-parameter family of flow laws, and argues the possibility of effective determination of exact and approximate analytical solutions on the basis of successive reduction to boundary-value problems for the Laplace equation or for the equation studied in detail in [1]. It should be noted that the characteristics of the flow are determined without additional quadratures.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 3, pp. 107–112, May–June, 1985.     

Instability of combined convective flow in a vertical layer

In the present paper, a study is made of the stability of plane-parallel flow induced by a transverse temperature difference between the boundaries of a layer and a longitudinal pressure gradient. This problem was solved earlier by the author [3] in a purely hydrodynamic formulation without allowance for thermal factors; the results then obtained correspond to the limiting case of small Prandtl numbers. In the paper, a numerical solution to the problem with the complete formulation is given.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 3, pp. 3–9, May–June, 1982.I thank G. Z. Gershuni for supervising the work, and also M. A. Gol'dshtik and V. N. Shtern for a helpful discussion.     

Law of transverse sections for the three-dimensional boundary layer on a thin wing in a hypersonic stream

The investigation of three-dimensional flows in boundary layers is important to determine the aerodynamic characteristics of wings such as the heat fluxes and friction drag. However, the circumstance that interaction of the boundary layer and the wake with an inviscid stream can play a governing role for the formation of the flow diagram as a whole is more important. The three-dimensional flow on a thin delta wing in a hypersonic stream is investigated in this paper. An important singularity of hypersonic flow is the low value of the gas density in the boundary layer as compared with the density on its outer boundary. It is shown that in the general case when the pressure in the wing span direction varies mainly by an order, high transverse velocities originate because of the smallness of the density within the boundary layer. This circumstance permits expansion of the solution for smallspan wings in a series in an appropriate small parameter. The equations in each approximation depend on two variables, while the third—longitudinal—variable enters as a parameter. The zero approximation can be considered as the formulation of the law of transverse plane sections for a three-dimensional boundary layer. As a comparison with the exact solutions calculated for delta wings with power-law distributions of the wing thickness has shown, the first approximation yields a very good approximation. Furthermore, flow modes with a different direction of parabolicity on the whole wing, as well as zones in which interaction with the external stream should absolutely be taken into account, are found.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 2, pp. 75–84, March–April, 1976.     

Experimental investigation of turbulence in a supersonic flow

The structure of disturbances carried by the flow into the working section of a supersonic wind tunnel has been investigated by means of a constant-current hot-wire anemometer. In order to generate the disturbances grids consisting of round rods were introduced upstream from the nozzle throat. It was found that in the working section the disturbances consist of non-correlating vortex, entropy and acoustic modes. The latter is generated by the boundary layer on the nozzle walls and the first two by the grids. The spectral compositions of the various modes are compared. Because of the presence of grid turbulence the point of laminar-turbulent transition in the boundary layer on a flat plate varied widely.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 3, pp. 119–124, May–June, 1990.     

Implementation of some higher-order convection schemes on non-uniform grids

A generalized formulation is applied to implement the quadratic upstream interpolation (QUICK) scheme, the second-order upwind (SOU) scheme and the second-order hybrid scheme (SHYBRID) on non-uniform grids. The implementation method is simple. The accuracy and efficiency of these higher-order schemes on non-uniform grids are assessed. Three well-known bench mark convection-diffusion problems and a fluid flow problem are revisited using non-uniform grids. These are: (1) transport of a scalar tracer by a uniform velocity field; (2) heat transport in a recirculating flow; (3) two-dimensional non-linear Burgers equations; and (4) a two-dimensional incompressible Navier-Stokes flow which is similar to the classical lid-driven cavity flow. The known exact solutions of the last three problems make it possible to thoroughly evaluate accuracies of various uniform and non-uniform grids. Higher accuracy is obtained for fewer grid points on non-uniform grids. The order of accuracy of the examined schemes is maintained for some tested problems if the distribution of non-uniform grid points is properly chosen.     

A Comparison Study Between an Adaptive Quadtree Grid and Uniform Grid Upscaling for Reservoir Simulation

An adaptive quadtree grid generation algorithm is developed and applied for tracer and multiphase flow in channelized heterogeneous porous media. Adaptivity was guided using two different approaches. In the first approach, wavelet transformation was used to generate a refinement field based on permeability variations. The second approach uses flow information based on the solution of an initial-time fine-scale problem. The resulting grids were compared with uniform grid upscaling. For uniform upscaling, two commonly applied methods were used: renormalization upscaling and local-global upscaling. The velocities obtained by adaptive grid and uniformly upscaled grids, were downscaled. This procedure allows us to separate the upscaling errors, on adaptive and uniform grids, from the numerical dispersion errors resulting from solving the saturation equation on a coarse grid. The simulation results obtained by solving on flow-based adaptive quadtree grids for the case of a single phase flow show reasonable agreement with more computationally demanding fine-scale models and local-global upscaled models. For the multiphase case, the agreement is less evident, especially in piston-like displacement cases with sharp frontal movement. In such cases a non-iterative transmissibility upscaling procedure for adaptive grid is shown to significantly reduce the errors and make the adaptive grid comparable to iterative local-global upscaling. Furthermore, existence of barriers in a porous medium complicates both upscaling and grid adaptivity. This issue is addressed by adapting the grid using a combination of flow information and a permeability based heuristic criterion.     

Numerical diffusion for flow‐aligned unstructured grids with application to estuarine modeling

The benefits of unstructured grids in hydrodynamic models are well understood but in many cases lead to greater numerical diffusion compared with methods available on structured grids. The flexible nature of unstructured grids, however, allows for the orientation of the grid to align locally with the dominant flow direction and thus decrease numerical diffusion. We investigate the relationship between grid alignment and diffusive errors in the context of scalar transport in a triangular, unstructured, 3‐D hydrodynamic code. Analytical results are presented for the 2‐D anisotropic numerical diffusion tensor and verified against idealized simulations. Results from two physically realistic estuarine simulations, differing only in grid alignment, show significant changes in gradients of salinity. Changes in scalar gradients are reflective of reduced numerical diffusion interacting with the complex 3‐D structure of the transporting flow. We also describe a method for utilizing flow fields from an unaligned grid to generate a flow‐aligned grid with minimal supervision. Copyright © 2013 John Wiley & Sons, Ltd.     

Dynamic model of liquid motion in a porous medium

A theory which takes account of the role of inertial effects in liquid motion in a porous medium is developed. For a compressible liquid, not only the hydrodynamic equations but also the thermodynamic transfer equation are formulated. The initiation, propagation, and dissipation characteristics of eddy motion are considered. Matching conditions at the interface between the media and boundary conditions are obtained. An approximate formulation of the problem is given, isolating in the porous medium a basic flow region in which the classical Darcy law is valid. As an illustration, weakly perturbed liquid flow in a plane channel with an insert of porous material is investigated.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 6, pp. 89–95, November–December, 1978.     

Flow of a stream containing charged inertial particles past a conducting sphere

The problem of the flow of a hydrodynamic stream containing electrically charged particles past a conducting sphere is solved. The influence of the volume density of the electric charge and the potential of the sphere on the capture coefficient is determined for different values of the inertia parameter of the particles and different drag laws for them. It is shown that for an earthed sphere can appreciably exceed unity. In the formulation of the problem, it is assumed that the region of electrogasdynamic flow is bounded by two electrode grids, which simulate the exit of the source of charged particles and the surrounding electric conditions. The velocity field near the sphere is assumed irrotational. The problem is solved numerically.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 4, pp. 21–26, July–August, 1982.I thank A. B. Vatazhin for assistance in the work.     

Some effects of strong blowing of gas into a supersonic flow

The intense evaporation of bodies moving in the atmospheres of planets at high supersonic velocities has been partly simulated both theoretically [1–5] (numerical calculations of strong blowing in the framework of the Navier-Stokes equations were also made at the Scientific-Research Institute of Mechanics at the Moscow State University by É. A. Gershbein and A. F. Kolesnikov [6]) as well as experimentally [7–9]. Below, the results are given of investigations of strong blowing of gas from the flat end of a cylinder into a supersonic flow at Reynolds numbers such that the mixing layer separating the blown and the oncoming gas is fairly thin. In this case, the mixing layer can be regarded as a contact surface, so that the problem of blowing can be solved in the framework of Euler's equations. The results of a numerical solution are compared with the results of experiments on the separation and profile of the shock wave, the thickness of the blowing layer on the axis, and also on the pressure distribution on the end of the cylinder. It was established experimentally, and then confirmed numerically that there is a downwash of the blown gas on the periphery of a porous end. It is shown that for the same blowing parameter K, which is equal to the ratio of the dynamic head of the blown gas to the dynamic head of the oncoming gas, and for a given distribution of K over the surface of the body the contact surface tends to a certain limiting position with increasing Mach number of the oncoming flow, i.e., the profile of the contact surface is stabilized. The influence of the adiabatic exponent on the thickness of the blowing layer is estimated. The present investigations continue earlier experimental studies, the main results of which have been presented in [9].Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 1, pp. 91–98, January–February, 1980.     

Analysis of the hydrodynamic interaction between cascades of thin profiles taking account of vortex wake evolution

The papers [1–5] are devoted to an investigation of aspects of the hydrodynamic interaction of cascades of profiles in a nonlinear formulation: it is shown experimentally in [1] and theoretically in [2] that the free vortex sheet ruptures upon meeting a profile; taking account of the evolution of vortex wakes, the flows around two cascades of solid profiles of infinitesimal [3] and finite [4] density are computed; results of an experimental investigation of the dynamic reactions of the flow on two mutually moving cascades of thin profiles are presented in [5]. The interference between two cascades of thin profiles in an inviscid, incompressible fluid flow is examined in this paper, where a modified method from [6] is used.Translated from Zhurnal Prikladnoi MekhaniM i Tekhnicheskoi Fiziki, No. 4, pp. 61–65, July–August, 1976.The author is grateful to D. H. Gorelov for discussing the research.     

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.     

Decay of turbulence behind a double grid of cylinders with opposite motion of the grids

The results are given of an experimental investigation into the decay of turbulence in a flow behind a double grid of cylinders with the grids moving in opposite directions. It is shown that the distribution of the turbulence parameters behind such a grid is more uniform than for other known methods of generating turbulence.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 3, pp. 17–25, May–June, 1979.     

Impact of a strip of compressible liquid on a barrier

The exact solution of the plane problem of the impact of a finite liquid strip on a rigid barrier is obtained in the linearized formulation. The velocity components, the pressure and other elements of the flow are determined by means of a velocity potential that satisfies a two-dimensional wave equation. The final expressions for them are given in terms of elementary functions that clearly reflect the wave nature of the motion. The exact solution has been thoroughly analyzed in numerous particular cases. It is shown directly that in the limit the solution of the wave problem tends to the solution of the analogous problem of the impact of an incompressible strip obtained in [1]. A logarithmic singularity of the velocity parallel to the barrier in the corner of the strip is identified. A one-dimensional model of the motion, which describes the behavior of the compressible liquid in a thin layer on impact and makes it possible to obtain a simple solution averaging the exact wave solution, is proposed. Inefficient series solutions are refined and certain numerical data on the impact characteristics for a semi-infinite compressible liquid strip, previously considered in [2–4] in connection with the study of the earthquake resistance of a dam retaining water in a semi-infinite basin, are improved. The solution obtained can be used to estimate the forces involved in the collision of solids and liquids. It would appear to be useful for developing correct and reliable numerical methods of solving the nonlinear problems of fluid impact on solids often examined in the literature [5].Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 6, pp. 138–145, November–December, 1990.The results were obtained by the author under the scientific supervision of B. M. Malyshev (deceased).     

Laminar boundary layer on an oscillating wedge

A study is made of the nonstationary laminar boundary layer on a sharp wedge over which a compressible perfect gas flows; the wedge executes slow harmonic oscillations about its front point. It is assumed that the perturbations due to the oscillations are small, and the problem is solved in the linear approximation. It is also assumed that the thickness of the boundary layer is small compared with the thickness of the complete perturbed region. Then in a first approximation the influence of the boundary layer on the exterior inviscid flow can be ignored, and the parameters on the outer boundary of the boundary layer can be taken equal to their values on the body for the case of inviscid flow over the wedge. They are determined from the solution to the inviscid problem that is exact in the framework of the linear formulation. The wall is assumed to be isothermal. The dependence of the viscosity on the temperature is linear. Under these assumptions, the problem of calculating the nonstationary perturbations in the boundary layer on the wedge is a self-similar problem.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 4, pp. 146–151, July–August, 1980.