Vertical immersion of a cylindrical floating solid

This paper is devoted to a study of the vortex-free (irrotational) motion of an ideal incompressible liquid during the vertical immersion of a cylindrical solid. In contrast to problems of impact [1] and the entry of a solid into water [2], the case here treated deals with continuous immersion involving a change in the shape of the free surface but with a constant width of the wetted surface of the solid. The coefficients of the time-dependent power series for the velocity potential, the equation of the free surface, and the pressure on the solid are determined, allowing for all the terms in the Cauchy-Lagrange equation. The results of calculations relating to the immersion of a bottom with an elliptical shape of the submerged part are presented.Translated from Izvestiya Akademii Nauk SSSR. Mekhanika Zhidkosti i Gaza, No. 2, pp. 27–31, March–April, 1973.     

Nonlinear waves on the surface of a charged liquid. Instability,bifurcation, and nonequilibrium shapes of the charged surface

Plane nonlinear waves in shallow water are described by the Kortewegde Vries equation [1–3]. The present paper contains theoretical investigations of nonlinear waves and nonlinear equilibrium shapes on the surface of a charged liquid. The influence of the field on the velocity and shape of a hydrodynamic soliton is considered. The bifurcation of the equilibrium shapes is investigated. Problems of the equilibrium shapes of a charged liquid are solved in the nonlinear formulation of the dynamics of nonlinear solitary forms (lunes, trenches) on the surface.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 3, pp. 94–102, May–June, 1984.     

Stability of the equilibrium of a liquid with multiply connected free surface in closed systems

A study is made of the equilibrium of a capillary liquid that partly fills a container and is in contact with several isolated gas cavities. The conditions of its stability are obtained with allowance for the change in the internal energy of the gas when the free surface is perturbed. The analysis of these conditions is accompanied by the solution of examples.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 3, pp. 4–9, May–June, 1984.     

Free surface asymptotics for thermocapillary flow at high marangoni numbers

Formal asymptotic expansions of the solution of the steady-state problem of incompressible flow in an unbounded region under the influence of a given temperature gradient along the free boundary are constructed for high Marangoni numbers. In the boundary layer near the free surface the flow satisfies a system of nonlinear equations for which in the neighborhood of the critical point self-similar solutions are found. Outside the boundary layer the slow flow approximately satisfies the equations of an inviscid fluid. A free surface equation, which when the temperature gradient vanishes determines the equilibrium free surface of the capillary fluid, is obtained. The surface of a gas bubble contiguous with a rigid wall and the shape of the capillary meniscus in the presence of nonuniform heating of the free boundary are calculated.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 3, pp. 61–67, May–June, 1989.     

A modified finite element method for determining equilibrium capillary surfaces of liquids with specified volumes

This paper describes a modified finite element method (MFEM) for determining the static equilibrium shape of the capillary surface of a liquid with a prescribed volume constrained by rigid boundaries with arbitrary shapes. It is assumed that the liquid is in static equilibrium under the influence of surface tension, adhesion, and gravity forces. This problem can be solved by employing the conventional FEM; however, a major difficulty arises due to the presence of the volume (integral) constraint and usually requires the use of the Lagrange multiplier method, the sequential unconstrained minimization technique, or the augmented Lagrange multiplier method. With the MFEM, the space variables defining the equilibrium surfaces (or curves) are expanded in terms of parametric interpolation functions, which are designed such that the boundary conditions and the integral constraint equation are automatically satisfied during each iteration of a direct numerical search process. Hence, there is no need to include Lagrange multipliers and/or penalty factors and the problem can be treated more simply as one involving unconstrained optimization. This investigation indicates that the MFEM is more efficient and reliable than the other methods. Results are presented for several case study problems involving liquid solder drops. Copyright © 2000 John Wiley & Sons, Ltd.     

Branching of rotationally symmetric solutions describing flows of a viscous liquid with a free surface

Equilibrium shapes of a liquid, situated on the outer or inner surface of a rigid cylinder and rotating together with it as a solid body, are studied. We determine the principal part of the solution of the equilibrium equation for small deviations of the determining parameter from the critical value. The bifurcation of rotationally symmetric motions with a free boundary in a body force field is also investigated.Translated from Zhurnal Prikladnoi Mekhaniki i Tekhnicheskoi Fiziki, No. 2, pp. 127–134, March–April, 1973.The authors thank V. Kh. Izakson for his discussion of the work.     

Transient boundary layer effect in underground oil and gas reservoir depletion processes

Modern methods of exploiting underground deposits of oil and gas are characterized by conditions that vary slowly with time, which makes it possible to treat them over an extended period as near-equilibrium processes and use for their analysis the effective methods of perturbation theory. At the same time, during a brief initial period these systems display essentially nonequilibrium behavior leading to a transient boundary layer effect. For closed reservoir depletion problems a measure of the degree of nonequilibrium of the reservoir system is introduced and for real deposits shown to be small, the existence of a boundary layer is established, and the exterior and interior problems are formulated, together with the matching conditions. The general form of the exterior asymptotic behavior, in which the space variables and time are separated, is established and the initial parabolic system is reduced to a linear Poisson equation. Examples of problem solving are given.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 6, pp. 72–77, November–December, 1985.     

Thermocapillary instability of the equilibrium of a flat layer containing internal heat sources

In the absence of body forces, a factor which has a strong influence on the equilibrium stability of a nonuniformly heated liquid is the dependence of the coefficient of surface tension on the temperature and the thermocapillary effect generated by it. If the equilibrium temperature gradient is sufficiently great, then the presence of the thermocapillary forces on the free surface can lead to the occurrence of convective motion. The monotonie instability of the equilibrium of a flat layer was investigated in [1–3]. Analysis of nonmonotonic disturbances [4] showed that in the case of an undeformable free surface there is no oscillatory instability. In [5] it was found that oscillatory instability is possible if there is a nonlinear dependence of the coefficient of surface tension on the temperature. The present paper is devoted to numerical investigation of the equilibrium stability of a flat layer with respect to arbitrary disturbances. It is shown that for a deformable free boundary there appears an additional neutral curve, which corresponds to monotonie capillary instability. In addition, when the capillary convection mechanism is taken into account, there appears an oscillatory instability, which becomes the most dangerous in the region of small Prandtl and wave numbers.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 2, pp. 27–31, March–April, 1991.I thank V. K. Andreev for a helpful discussion of the work.     

不定边界问题的边界元和有限元数学模型

对于不可压缩有势流动,有两类典型的不定边界或可动边界问题,即定常型不定边界如常水位闸门出流、过水坝或过水闸、水堰绕流及射流等,非定常型的不定边界如装液容器内的流体晃动问题.它们的共同特点是自由面位置与形状事先是的,需在计算过程中调整网格作自由面拟合.且由于自由面条件呈强烈非线性,给不确定数值计算带来困难.本文综述了两类不定边界问题的有限元和边界元模式,简述了笔者的一些计算经验.      

Dynamic problem for a two-layer medium with a vibrational source at the boundary interface

The present study is concerned with an analysis of gravitational and acoustic waves which are excited by a vibrational source deeply placed in a liquid covered by ice. An analysis of the rigidity characteristics of ice modeled by an elastic layer or by a Kirchhoff plate is done by factorization of the solution to the integral equation equivalent to an initially combined boundary value problem. The uncombined boundary condition is used to solve problems for unrestricted ice fields in [1–3], whereas combined conditions with vibrational sources positioned at the boundary of the medium are used in [4].Translated from Zhurnal Prikladnoi Mekhaniki, No. 3, pp. 125–129, May–June, 1986.     

Influence of surface tension on the stability of heterogeneous equilibrium of barotropic liquid phases

On the basis of the results of earlier work of the author [1] a study is made of the equilibrium and stability of a two-phase single-component heterogeneous liquid system with respect to perturbations of arbitrary shape. Allowance is made for the influence of surface tension, which plays a critical part in the formation of nucleating centers of a new phase [2]. Conditions of equilibrium are derived, and also a criterion of radial stability of a nucleating center of a new phase bounded by a closed spherical boundary. It is shown that radial stability of spherical nucleating centers also guarantees stability with respect to perturbations of arbitrary shape. The part played by the finite size of the system and the boundary conditions is elucidated. For this, two different cases are studied: a) a system under a constant external pressure, b) a system with fixed volume. In the first case, all equilibrium states are unstable. In the second, there are both unstable and stable configurations (depending on the corresponding values of two dimensionless parameters). The equation of the hyperbola of neutral stability is derived. The limits of a very small coefficient of surface tension and a very large size of the container are considered. The first situation corresponds to stable configurations, the second to unstable. For simplicity, the considered systems are assumed to be isothermal, and the equilibrium and stability are analyzed on the basis of the mechanical analog of Gibbs's principle, namely, the principle of a minimum of the mechanical potential energy of the barotropic heterogeneous liquid system. The case of nonisothermal perturbations leads to similar results, but the expressions for the corresponding dimensionless parameters are more cumbersome and less physically perspicuous.     

Nonlinear problems of the thermal conductivity equation

The asymptotic behavior of solutions of parabolic equations at infinite times has been investigated for various cases [1–6]. Two initial boundary-value problems are considered in this paper. The solution of the thermal conductivity equation with a nonlinear right-hand side is found, including also nonlinear boundary conditions. It is shown that the solution of the corresponding problem tends either to a stable, steady-state solution, or to a periodic solution, depending on the initial values of the functions and constants appearing in the conditions of the problem. Other papers [7, 8] are devoted to finding the periodic solutions of these two problems encountered in hydrodynamics (diffusion, underground hydrodynamics), and to studying the asymptotic behavior of the corresponding initial boundary problems.Translated from Zhurnal Prikladnoi Mekhaniki i Tekhnicheskoi Fiziki, No. 3, pp. 123–128, May–June, 1972.     

Reduced numerical modeling of flows involving liquid-crystalline polymers

Kinetic theory models involving the Fokker–Planck equation can be accurately discretized using a mesh support (Finite Elements, Finite Differences, Finite Volumes, Spectral Techniques, …). However, these techniques involve a high number of approximation functions. In the finite element framework, widely used in complex flow simulations, each approximation function has only local support and is related to a node that defines the associated degree of freedom. In the technique proposed here, a reduced approximation basis is constructed. The new shape functions have extended support and are defined in the whole domain in an appropriate manner (the most characteristic functions related to the model solution). Thus, the number of degrees of freedom involved in the solution of the Fokker–Planck equation is very significantly reduced. The construction of those new approximation functions is done with an ‘a priori’ approach, which combines a basis reduction (using the Karhunen–Loève decomposition) with a basis enrichment based on the use of some Krylov subspaces. This paper analyzes the application of model reduction to the simulation of non-linear kinetic theory models involving complex behaviors, such as those coming from stability analysis, complex geometries and coupled models. We apply our model reduction approach to the Doi's classical constitutive equation for viscoelasticity of liquid-crystalline polymer.     

Waves on the surface of a thin layer of viscous liquid

Under the assumption that the boundary layer approximation for the original equations is valid, we show the possibility of the existence of progressive waves on the surface of a vertically flowing film when surface tension is neglected. From the system of equations obtained for a thin layer of viscous liquid flowing down an inclined plane, one equation for perturbations of a thin film follows. Steady solutions of this equation allow periodic discontinuous solutions of the roll-wave type.Translated from Zhurnal Prikladnoi Mekhaniki i Teknicheskoi Fiziki, No. 2, pp. 109–113-March–April, 1973.     

Vibrations of a gas-filled spherical cavity in pseudoplastic and dilatant liquids

Questions of the dynamics of bubbles in a liquid are connected with problems of cavitation [1]. In connection with cavitation phenomena in non-Newtonian media, in particular in polymeric liquids [2, 3], a study is made of the pulsations of a bubble in a polymeric liquid with an exponential rheological law. The equation of the motion of the boundary of the gas cavity is integrated numerically; here, the cases of pseudo-plastic and dilatant liquids are discussed separately. The results obtained can be used in the analysis of acoustical cavitation in aqueous solutions of polymers.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 1, pp. 146–148, January–February, 1975.     

Nonaxisymmetric equilibrium shapes of a drop on a plane

A nonuniform temperature distribution, the presence of surface-active substances and impurities, and also other factors lead to a change in the wetting angle along a plane. A study is made of the influence of a small perturbation of the equilibrium contact angle on the shape of the free surface of the liquid. Two cases are considered: a surface of small slope in a gravity field and a nearly spherical shape under conditions of weightlessness. The equilibrium shapes of a liquid drop on an inclined plane under conditions of hysteresis of the wetting are also obtained.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 4, pp. 164–167, July–August, 1983.I thank I. E. Tarapov and I, I, Ievlev for constant interest in the work and valuable comments.     

Numerical simulation of the dynamics of a liquid in a moving axisymmetric cavity

The problem of the motion of an ideal liquid with a free surface in a cavity within a rigid body has been most fully studied in the linear formulation [1, 2]. In the nonlinear formulation, the problem has been solved by the small-parameter method [3] and numerically [4–7]. However, the limitations inherent in these methods make it impossible to take into account simultaneously the large magnitude and the threedimensional nature of the displacements of the liquid in the moving cavity. In the present paper, a numerical method is proposed for calculating such liquid motions. The results of numerical calculations for spherical and cylindrical cavities are given.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 2, pp. 174–177, March–April, 1984.     

Intense evaporation of a gas from a two-dimensional periodic surface

Evaporation (or condensation) of a gas is said to be intense when the normal component of the velocity of the gas in the Knudsen layer has a value of the order of the thermal velocity of a molecule, c_T=(2kT/m)^1/2. In this case the distribution function of the molecules with respect to their velocities in the Knudsen layer differs from the equilibrium (Maxwellian) value by its own magnitude. As a result of this, over the thickness of the Knudsen layer the macroparameters also vary by their own magnitudes. So in order to obtain the correct boundary conditions for the Euler gas dynamic equations, it is necessary to solve the nonlinear Boltzmann equation in the Knudsen layer. The problem of obtaining such boundary conditions for the case of a plane surface was considered in [1–11]. In the present study this problem is solved for a two-dimensional periodic surface in the case when the dimensions of the inhomogeneities are of the order of the mean free path of the molecules and the inhomogeneities have a rectangular shape. The flow in the Knudsen layer becomes two-dimensional, and this leads to a considerable complication of the solution of the problem.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 2, pp. 132–139, March–April, 1985.In conclusion the author would like to express his gratitude to V. A. Zharov for his valuable advice, and also V. S. Galkin, M. N. Kogan, and N. K. Makashev for discussion of the results obtained.     

Dynamics of a flat cavity in a low-viscosity liquid

The problem of the motion of a cavity in a plane-parallel flow of an ideal liquid, taking account of surface tension, was first discussed in [1], in which an exact equation was obtained describing the equilibrium form of the cavity. In [2] an analysis was made of this equation, and, in a particular case, the existence of an analytical solution was demonstrated. Articles [3, 4] give the results of numerical solutions. In the present article, the cavity is defined by an infinite set of generalized coordinates, and Lagrange equations determining the dynamics of the cavity are given in explicit form. The problem discussed in [1–4] is reduced to the problem of seeking a minimum of a function of an infinite number of variables. The explicit form of this function is found. In distinction from [1–4], on the basis of the Lagrauge equations, a study is also made of the unsteady-state motion of the cavity. The dynamic equations are generalized for the case of a cavity moving in a heavy viscous liquid with surface tension at large Reynolds numbers. Under these circumstances, the steady-state motion of the cavity is determined from an infinite system of algebraic equations written in explicit form. An exact solution of the dynamic equations is obtained for an elliptical cavity in the case of an ideal liquid. An approximation of the cavity by an ellipse is used to find the approximate analytical dependence of the Weber number on the deformation, and a comparison is made with numerical calculations [3, 4]. The problem of the motion of an elliptical cavity is considered in a manner analogous to the problem of an ellipsoidal cavity for an axisymmetric flow [5, 6]. In distinction from [6], the equilibrium form of a flat cavity in a heavy viscous liquid becomes unstable if the ratio of the axes of the cavity is greater than 2.06.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 5, pp. 15–23, September–October, 1973.The author thanks G. Yu. Stepanov for his useful observations.     

One-dimensional propagation of a monochromatic light pulse in absorbing media

A study is made of the equations describing the propagation of a monochromatic pulse of radiation of arbitrary shape in absorbing media — a plasma and an absorbing two-level photodissociable medium. Exact analytic solutions are found for a wide variety of boundary conditions. The discussion is carried through for problems with plane, cylindrical, and spherical symmetry. The formulas obtained can be used directly to compare calculation and experiment.Translated from Zhurnal Prikladnoi Mekhaniki i Tekhnicheskoi Fiziki, No. 5, pp. 177–180, September–October, 1972.In conclusion the authors thank O. N. Krokhina for a discussion of the problems considered in the paper.