A Eulerian/Lagrangian model to calculate the evolution of a water droplet spray

We introduce a Eulerian/Lagrangian model to compute the evolution of a spray of water droplets inside a complex geometry. To take into account the complex geometry we define a rectangular mesh and we relate each mesh node to a node function which depends on the location of the node. The time-dependent incompressible and turbulent Navier-Stokes equations are solved using a projection method. The droplets are regarded as individual entities and we use a Lagrangian approach to compute the evolution of the spray. We establish the exchange laws related to mass and heat transfer for a droplet by introducing a mass transfer coefficient and a heat transfer coefficient. The numerical results from our model are compared with those from the literature in the case of a falling droplet in the atmosphere and from experimental investigation in a wind tunnel in the case of a polydisperse spray. The comparison is fairly good. We present the computation of a water droplet spray inside a complex and realistic geometry and determine the characteristics of the spray in the vicinity of obstacles.     

Motion of a solidcontaininga spherical damper

For the purpose of modeling the motion of a solid with a cavity filled with a viscous fluid, M. A. Lavrent'ev [1] has proposed a model in the form of a solid with a spherical cavity in which another solid, spherical in shape, is enclosed. The sphere is separated from the cavity walls by a small, clearance in which viscous forces act (a lubricating film). This simple model with a finite number of degrees of freedom possesses certain mechanical properties of a solid with a cavity containing a viscous fluid. Study of this model is therefore of interest.The present paper examines certain properties of the model, which will be termed a solid with a damper. It is shown that for a highviscosity lubricant the motion of a solid with a damper can be described by the same equations which pertain to the motion of a solid with a spherical cavity filled with a high-viscosity fluid. Expressions relating the parameters of the systems are obtained. If these relations are fulfilled, the systems become mechanically equivalent.The steady motions of a free solid with a damper and their stability conditions are determined.These motions and stability conditions hold for a body with a cavity filled with a viscous fluid [2].     

Motion of a variable-volume sphere in an ideal fluid near a plane surface

The motion of a spherical cavity in a fluid is investigated. The radius of the sphere varies under the action of a constant pressure at infinity. The problems of the collapse of a cavity moving in an unbounded fluid and of the collapse of a cavity near a plane are solved in the exact formulation. The occurrence of an initial translational velocity or the presence of a solid surface, by contrast with the collapse of a sphere at rest in an unbounded fluid [1], yields a limiting radius at which the process of collapse ceases. A sphere initially at rest near a plane always comes into contact with the plane as a result of collapse. The radius and velocities at which the sphere arrives the plane are calculated for various initial distances from the latter. The possible mechanism of the action of a cavitation bubble on a solid surface is discussed.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 5, pp. 94–103, September–October, 1971.     

Measurement of the time-dependent principal-strain directions in birefringent materials

In this paper a simple method is presented to analyze the time variation of the principal-strain directions. The experimental arrangement to be used consists of a dark-field linear polariscope, a beam of monochromatic parallel light with a small diameter penetrating the measuring point with normal incidence. The variation of the light intensity behind the analyzer was measured by a photocell and a CRO. According to the method introduced here, two variations of light intensity were measured by choosing two orientations of the crossed system of polarizer and analyzer. The difference between the orientations is 45 deg. This method was applied to the analysis of the resulting principal directions in a region in which a superposition of a compression and a shear wave occurred. These waves were produced in a photoelastic foil, which was cemented on a brass rod impacted longitudinally. The compressional main pulse generated in the rod has a duration of about 25 μs. This experimental arrangement will be used later to investigate the behavior of the photoelastic material under short-time loading.     

Effect of boundaries of various shapes on vortex dynamics

The motion of a vortex near a boundary of arbitrary shape is considered within the framework of a two-dimensional problem. Integrable differential equations of motion are obtained. Two forms of the algebraic equation of the vortex trajectories are derived. Examples of vortex motion near a straight-line boundary, in a channel, in an angular domain, in the neighborhood of a flat edge, in a round basin, and near a parabolic boundary.     

Bifurcations of rotating structures in a parabolic functional differential equation

We investigate the Andronov-Hopf bifurcation of the birth of a periodic solution from a space-homogeneous stationary solution of the Neumann problem on a disk for a parabolic equation with a transformation of space variables in the case where this transformation is the composition of a rotation by a constant angle and a radial contraction. Under general assumptions, we prove a theorem on the existence of a rotating structure, deduce conditions for its orbital stability, and construct its asymptotic form. __________ Translated from Neliniini Kolyvannya, Vol. 9, No. 2, pp. 155–169, April–June, 2006.     

The extension and oscillation of a non-Hooke's law spring

We derive a wave equation for small-amplitude, undamped, extensional oscillation of a spring-mass system consisting of a mass suspended on a spring governed by a quadratic force-extension relationship. We justify this quadratic model using a Taylor series expansion of the general elasticity equations for a helical spring. Transformation of the equation of motion of the spring leads to a separable wave equation with the spacial component being a transformation of Bessel's equation. The model is successful in predicting static extension and period of oscillation of a helical wire spring for which the wave equation based on Hooke's law is inadequate.     

Lattice of plates in an unsteady supersonic flow

The supersonic unsteady flow of a gas past a lattice of thin, slightly curved profiles has been investigated in several studies. The paper [1] is devoted to an evaluation of the effect of wind tunnel walls on the unsteady aerodynamic characteristics of a profile, and [2] investigates the effects of the boundaries of a free jet. These cases are equivalent respectively to the anti-phase and in-phase oscillations of the profiles of an unstaggered grid. In [3] consideration is given to a more general case of gas flow past a profile in a wind tunnel with perforated walls. Flow past a lattice of profiles with stagger is studied in [4], where the magnitude of the stagger angle is limited by the condition that the lattice leading edge is located in the undisturbed stream.In the present paper we present a method of calculation of the unsteady supersonic flow of a gas past a lattice of profiles with arbitrary stagger. As an example the results are presented of the calculation of the aerodynamic forces and moments acting on an oscillating profile in a wind tunnel with solid walls and in a free jet.     

Open boundary condition symposium benchmark solution: Stratified flow over a backward-facing step

This paper presents a detailed numerical solution to a simplified version of two-dimensional stratified flow over a backward-facing step with a Froude number of 16/9, a Reynolds number of 800 and a Prandtl number of 1—one of the Open Boundary Condition Symposium test problems. The steady state solution was derived by integrating the time-dependent Boussinesq equations forward in time using a semi-implicit finite-element-based model on a 38400-element mesh. In addition to presenting the results derived on this grid, the paper also presents the results of a Richardson extrapolation calculation for a set of ‘key’ parameters. It is hoped that this solution can be used as a baseline to compare the performance of the various techniques discussed at the Open Boundary Condition Symposium.     

Three-dimensional isothermal lava flows over a non-axisymmetric conical surface

Within the thin-layer approximation for a highly-viscous heavy incompressible fluid, a hydrodynamicmodel of a 3D isothermal lava flow over a non-axisymmetric conical surface is constructed. Using analytical methods, a self-similar solution for the law of leading-edge propagation is obtained in the case of a flow from a non-axisymmetric source located at the apex of a conical surface with smoothly varying properties. In the case of a flow over a substantially non-axisymmetric surface, it is shown that there exists a self-similar solution for the free-surface shape and the law of leading-edge motion. This solution is studied numerically for particular examples of the substrate surface and the source. In the general case of a non-self-similar flow over a substantially non-axisymmetric conical surface, a local analytical solution is obtained for the free-surface shape and the velocity field near the leading flow front.     

Interaction of a heavy fluid flow in a channel with a transverse jet barrier

An experimental and numerical analysis of the interaction between a plane horizontal water flow in a rectangular channel (free water current) and a plane thin water jet (water jet curtain) is presented; the jet flows out vertically from either a slot nozzle in the bottom of the channel or the crest of a rigid spillway at a velocity appreciably (several times) greater than the water velocity in the channel. Numerical calculations were carried out using the STAR-CD software package preliminarily tested against the experimental data obtained. The dependence of the water level in the channel at a certain distance ahead of the jet barrier on the main jet parameters and the water flow rate in the horizontal channel is studied. It is found that in the region of the interface between the flows both steady and unsteady (self-oscillatory) flow patterns can be realized. Steady stream/jet interaction patterns of the “ejection” and “ejection-spillway” types are distinguished and a criterion separating these regimes is obtained. The notion of a rigid spillway equivalent to a jet curtain is introduced and an approximate dependence of its height on the relevant parameters of the problem is derived. The possibility of effectively controlling the water level ahead of a rigid spillway with a sharp edge by means of a plane water jet flowing from its crest is investigated. The boundary of transition to self-oscillation interaction patterns in the region of the flow interface is determined. The structure of these flows and a possible mechanism of their generation are described. Within the framework of the inviscid incompressible fluid model in the approximate formulation for a “thin” jet, an analytical dependence of the greatest possible depth of a reservoir filled with a heavy fluid at rest and screened by a vertical jet barrier on the jet parameters is obtained.     

Dynamical regimes of fiber spinning

An analytical solution is obtained that describes fiber spinning with a given force on the receiving bobbin. As an example, a calculation is made of the response of the final fiber section to a periodically varying draw force; a solution is constructed that describes the propagation along the fiber of a finite perturbation associated with a change in the conditions at the spinneret for a fixed draw force. The problem of the small perturbations of a fiber spun at a given rate onto a bobbin is reduced to a linear integrodifferential equation with retardation whose characteristic equation determines the region of the “draw resonance” instability. The reasons for the occurrence for the instability are elucidated.     

Flow of a gas behind a Chapman-Jouguet detonation wave in the vicinity of a line of discontinuity of a surface wave curvature

This paper discusses questions of constructing a solution of the gasdynamic equations near a line of curvature discontinuity at the surface of a detonation wave, propagating under Chapman—Jouguet conditions. It describes the construction of the solution in two cases: in a flow arising with the initiation of a detonation along a half-plane in a quiescent homogeneous combustible gas and in a flow arising with the initiation of a detonation along a half-line under these same conditions.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 1, pp. 120–126, January–February, 1978.     

Draft of a nozzle with the circulating outflow of a vortical flow of gas

The presence of circulation in an outflowing gas leads to a change in the working parameters of a nozzle. The question of the mass flow rate and the draft of a nozzle without a diffusor (a point) for twisted flows has been studied theoretically and experimentally [1–6]. The use of nozzles with a supersonic part introduces a considerable degree of complication into the method for the analytical calculation of the draft characteristics and the program for their experimental investigation. In [2, 7], a theory of a nozzle is formulated for a model of a potential circulating flow of gas; in [5, 8], an electronic computer was used to solve the complete system of the equations of gasdynamics for the motion of a rotating flow along a nozzle; in [7, 9], an investigation was made of a variational problem of the shaping of a diffusor for a circulation flow. The calculation of the draft, carried out in the above-mentioned communications (with the exception of [2], in which a study was made of a partial model of an eddyless rotational motion), is bound up with labor-consuming computer calculations. In the present article, in a development of [3, 6], a quasi-one-dimensional theory of a supersonic nozzle for a vortical flow of gas is formulated and verified experimentally.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 5, pp. 145–149, September–October, 1975.     

On configurational balance in slender bodies

We propose a new derivation of the evolution equation of a sharp, coherent interface in a two-phase body having elongated shape, a body which we regard as a one-dimensional micropolar continuum. To this aim, we introduce a system of forces acting at the interface, and we apply the method of virtual powers to derive a balance law involving these forces. By exploiting the dissipation inequality, we manage to write this balance law in terms of a scalar field whose form is reminiscent of a well-known expression for the configurational stress in three dimensional micropolar continua.     

Transient responses of beams and plates subject to travelling load. Miscellaneous results

The steady-state response of a free and infinite Timoshenko beam is specified analytically in terms of non-dimensional displacements and stresses. The beam is supposed loaded by a travelling concentrated force or a moving step load. By a validated explicit numerical calculation, it is shown how a load travelling on a beam at constant velocity, from defined time and abscissa, generates a response which evolves towards the steady-state solution for a part, and towards a quantified transient solution for another part. Asymptotic values are given for the transient displacements and stresses according to the time and the speed of the loading. The solution is also found for a plate subject to a pressure, which spreads respecting the cylindrical symmetry. It is possible to identify in the response a part which follows the pressure front, and which is comparable with the steady-state response of a beam, and another transient part, which generates displacements and stresses with a much less oscillating character. An asymptotic solution is also presented for the plate.The whole series of the results makes it possible to better understand qualitatively the beginning of the transient response of a beam or of a plate to a moving load, and also makes it possible to estimate the stresses and displacements without needing specialised numerical codes.     

Function of region

The purpose of this paper is to extend points function and interval functions theoretics to an arbitrary region. For this, the new theory, the contraction of a region, and the retraction of a region; the extension of a region, and the kernel-preserving extension of a region are established by the author. Starting from these concepts, the new definitions of a region function is given. And a kernel (i.e. fixed point) of a region function is connected with a stable centre of defining region of such a region function. Thereby, the region theoretics and algorithms are established.In applications, to find a stable centre of a region, the author has utilized the measure theoretics of matrice defined by Hartfiel⁽⁷⁾ and other authors. The measure problems of coefficient matrice of system of equations of linear algebra associated with some region are discussed.     

On the problem of vortex interaction with a plane

We consider the well-known problem of the interaction of a vortex filament with a perpendicular plane in a viscous incompressible fluid. In this study, the vortex filament is represented by a semi-infinite rotating needle. Different models are considered: a zero-radius needle and fixed and movable in the axial direction needles of a finite radius. The ranges of the existence of the solution are found, and the correspondence of the flow around a finite-radius needle to that around a zero-radius needle, as the needle radius decreases, is studied.     

Displacement of yield-stress fluids in a fracture

We consider a displacement of several yield-stress fluids in a Hele-Shaw cell. The topic is relevant to the development of a model for the flow of multiple phases inside a narrow fracture with application to hydraulically fracturing a hydrocarbon-bearing underground formation. Existing models for fracturing flows include only pure power-law models without yield stress, and the present work is aimed at filling this gap. The fluids are assumed to be immiscible and incompressible. We consider fluid advection in a plane channel in the presence of density gradients. Gravity is taken into account, so that there can be slumping and gravitational convection. We use the lubrication approximation so that governing equations are reduced to a 2D width-averaged system formed by the quasi-linear elliptic equation for pressure and transport equations for volume concentrations of fluids. The numerical solution is obtained using a finite-difference method. The pressure equation is solved using an iterative algorithm and the Multigrid method, while the transport equations are solved using a second-order TVD flux-limiting scheme with the superbee limiter. This numerical model is validated against three different sets of experiments: (i) gravitational slumping of fluids in a closed Hele-Shaw cell, (ii) viscous fingering of fluids with a high viscosity contrast due to the Saffman–Taylor (S–T) instability in a Hele-Shaw cell at microgravity conditions, (iii) displacement of Bingham fluids in a Hele-Shaw cell with the development of fingers due to the S–T instability. Good agreement is observed between simulations and laboratory data. The model is then used to investigate the joint effect of fingering and slumping. Numerical simulations show that the slumping rate of yield-stress fluid is significantly less pronounced than that of a Newtonian fluid with the same density and viscosity. If a low-viscosity Newtonian fluid is injected after a yield-stress one, the S–T instability at the interface leads to the development of fingers. As a result, fingers penetrating into a fluid with a finite yield stress locally decrease the pressure gradient and unyielded zones develop as a consequence.