A Self-Similar Solution to the Problem of Lava Dome Growth on an Arbitrary Conical Surface

Within the framework of the asymptotic thin-film equations for a highly viscous heavy Newtonian fluid, a hydrodynamic model of non-axisymmetric lava dome growth on a conical surface is constructed. A new class of self-similar solutions describing the flow on a conical surface with finite inclination to the horizontal and a point mass supply at the apex is found analytically for power-law or exponential growth of the liquid volume with time. For a conical surface with a small inclination to the horizontal, the free-surface shape is found numerically. The asymptotics of this solution are compared with the solutions describing the flow on a horizontal plane and a conical surface with finite inclination to the horizontal.     

Depth-sensing indentation of a transversely isotropic elastic layer: Second-order asymptotic models for canonical indenters

Simple analytical approximations to the frictionless indentation problem for a transversely isotropic layer are obtained for spherical, conical, and pyramidal indenters as well as for axisymmetric indenters of power-low profile and self-similar non-axisymmetric indenters. These approximations are asymptotically exact in the small-contact limit. The results obtained are validated in the case of an isotropic layer for spherical and conical indenters.     

Free-surface deformation due to spiral flow owing to a source/sink and a vortex in Stokes flow

The free-surface shape and cusp formation are analyzed by considering a viscous flow arising from the superposition of a source/sink and vortex below the free surface where the strength of the source and vortex are arbitrary. In the analysis, Stokes’ approximation is used and surface tension effects are included, but gravity is neglected. The solution is obtained analytically by using conformal mapping and complex function theory. From the solution, shapes of the free surface are obtained, and the formation of a cusp on the free surface is discussed. Above some critical capillary number with a sink, the free-surface shape becomes singular and an apparent cusp should form on the free surface below a real fluid. On the other hand, no cusp would occur for sources of zero or positive strength. Typical streamline patterns are also shown for some capillary numbers. As the capillary number vanishes, the solution is reduced to a linearized potential flow solution.     

Non-Isothermal Lava Flows over a Conical Surface

Using the equations of a non-isothermal thin layer of viscous fluid with an exponential dependence of the viscosity on temperature, a family of hydrodynamic models of a cooling lava flow over a conical surface in the presence of mass supply is constructed. These models correspond to asymptotically different rates of heat exchange with the ambient medium. The evolution of the free-surface shape and the temperature fields is investigated numerically for a stationary mass supply. Using the matched asymptotic expansions method, solutions valid both near and very far from the mass supply region are constructed. The solutions obtained are compared with known analytical solutions for isothermal flow.__________Translated from Izvestiya Rossiiskoi Academii Nauk, Mekhanika Zhidkosti i Gaza, No. 2, 2005, pp. 62–75.Original Russian Text Copyright © 2005 by Osiptsov.     

Steady Film Flow of a Highly Viscous Heavy Fluid with Mass Supply

Asymptotic models of a thin layer of highly viscous heavy incompressible Newtonian fluid are constructed for steady axisymmetric (plane) flow on a curved rigid surface with distributed or point mass supply on a surface section near the axis (plane) of symmetry. Examples of analytical and numerical investigations of the free-surface shape and hydrodynamic-parameter fields are given. The models constructed are generalized for the case of a viscoplastic fluid and solutions which can be used for describing extrusive volcanic eruptions are obtained.     

On the theory of supersonic inviscid flow separation in gasdynamic problems

A general schematic flow representation that explains the mechanism of inviscid gas separation in time-dependent and three-dimensional gas flows is presented. The scenario of gas flow separation from a body surface or a mixing layer is described as a vortex which induces in the flowfield a velocity opposing to that of the main flow, thus decelerating it. Within the framework of this scenario the analytical conditions of separation are obtained for conical and self-similar gas flows which coincide with the results of experimental and numerical simulations.     

Fresh water injection into a geothermal reservoir which contains superheated vapor and solid phase salt

A mathematical model is proposed for the process with two unknown phase transition interfaces. At the leading interface, evaporation occurs and a salt precipitates. The interface that moves at a lower velocity is a dissolution surface. In the isothermal approximation the self-similar solution is obtained. It is shown that fresh water injection may lead to solid precipitate transfer from one region to another, substantially increasing the precipitate content behind the evaporation front. This significantly decreases permeability and, hence, the flow velocity through permeable rocks. With approaching the critical parameter region the branches of the obtained two-valued self-similar solution also approach and at critical values merge, which corresponds to the disappearance of the solution. The non-existence of the self-similar solution can be interpreted as the filling of the pores with salt precipitate and the flow arrest in the geothermal reservoir.     

Two-dimensional spreading of a cloud of conducting gas in a magnetic field self-similar solution

A self-similar solution of the problem on the spreading in a magnetic field of a cloud of conducting gas, having the shape of a cylinder of noncircular cross section, is constructed. The cylindrical surface of the gas is restrained by a nonconducting sheath that spreads according to a prescribed law. The shape of the transverse cross section of the cylindrical cloud is determined from the solution. Cross sections obtained for a concrete case are represented in graphic form.Translated from Zhurnal Prikladnoi Mekhaniki i Tekhnicheskoi Fiziki, No. 2, pp. 29–36, March–April, 1973.In conclusion the author thanks V. I. Khonichev for assistance with the numerical calculations.     

Transcritical Flow over a Ramp in an Open Channel

Steady free-surface flows over an elongated obstacle located on the channel bottom are studied theoretically and experimentally. For determining the free-surface shape and the main flow parameters, the first and second shallow-water approximations are used. In the second approximation, a solution describing a smooth transition from the subcritical to superctitical flow is found. A mathematical model of the hydraulic jump behind the ramp is constructed. The results of the mathematical modeling are compared with experimental data.     

Long-lasting effect of initial configuration in gravitational spreading of material fronts

We present the results from laboratory experiments and fully resolved simulations pertaining to finite-release turbulent density flows with a non-axisymmetric initial shape. First, we demonstrate that the effects of the initial shape influence the current’s evolution well into the long-time phase which would corresponds to the inertial self-similar phase in the case of planar or axisymmetric configurations. Then, we identify the physical mechanisms responsible for this dependence and propose a new model capable of capturing the dynamics of such releases. Finally, we show that this dependence on the initial configuration is robust for various types of gravity currents over a wide range of parameters such as Reynolds number, density ratio, and aspect ratio.     

Complementary mild-slope equations in a two-layer fluid

Mild-slope (MS) type equations are depth-integrated models, which predict under appropriate conditions refraction and diffraction of linear time-harmonic water waves. By using a streamfunction formulation instead of a velocity potential one, the complementary mild-slope equation (CMSE) was shown to give better agreement with exact linear theory compared to other MS-type equations. The main goal of this work is to extend the CMSE model for solving two-layer flow with a free-surface. In order to allow for an exact reference, an analytical solution for a two-layer fluid over a sloping plane beach is derived. This analytical solution is used for validating the results of the approximated MS-type models. It is found that the two-layer CMSE model performs better than the potential based one. In addition, the new model is used for investigating the scattering of linear surface water waves and interfacial ones over variable bathymetry.     

Numerical simulations of granular free-surface flows using smoothed particle hydrodynamics

This paper presents a two-dimensional SPH model designed to simulate free-surface flows of dense granular materials. Smoothed particle hydrodynamics (SPH) is a mesh-free numerical method based on a Lagrangian discretization of the continuum mass and momentum conservation equations. The rheology of dense granular materials is modelled using a new local constitutive law recently proposed by Jop et al. (Nature, 2006). Of the viscoplastic class, this law is characterized by an apparent viscosity depending both on the local strain rate and local pressure. Validation test cases of the model in steady and unsteady configurations are presented. For steady cases (vertical chute flow and uniform free-surface layers on inclines), excellent agreement with analytical predictions is obtained. In the unsteady case, the simulations satisfactorily capture the dynamics of gravity-driven surges observed in experiments, including behaviours that are very specific to granular materials. Among the various parameters involved in the computations, the influence of SPH particle configuration within the flow and of the threshold viscosity used in the regularization of the constitutive yield criterion are particularly discussed.     

Interaction of viscous wakes with a free surface

The interaction of laminar wakes with.free-surface waves generated by a moving body beneath the surface of an incompressible viscous fluid of infinite depth was investigated analytically. The analysis was based on the steady Oseen equations for disturbed flows.The kinematic and dynamic boundary conditions were linearized for the small-amplitude free-surface waves. The effect of the moving body was mathematically modeled as an Oseenlet.The disturbed flow was regarded as the sum of an unbounded singular Oseen flow which represents the effect of the viscous wake and a bounded regular Oseen flow which represents the influence of the free surface. The exact solution for the free-surface waves was obtained by the method of integral transforms. The asymptotic representation with additive corrections for the free-surface waves was derived by means of Lighthill‘s two-stage scheme. The symmetric solution obtained shows that the amplitudes of the free-surface waves are exponentially damped by the presences of viscosity and submergence depth.     

Viscous Spreading of Non-Newtonian Gravity Currents on a Plane

A gravity current originated by a power-law viscous fluid propagating on a horizontal rigid plane below a fluid of lower density is examined. The intruding fluid is considered to have a pure Ostwald power-law constitutive equation. The set of equations governing the flow is presented, under the assumption of buoyancy-viscous balance and negligible inertial forces. The conditions under which the above assumptions are valid are examined and a self-similar solution in terms of a nonlinear ordinary differential equation is derived. For the release of a time-variable volume of fluid, the shape of the gravity current is determined numerically using an approximate analytical solution derived close to the current front as a starting condition. A closed-form analytical expression is derived for the special case of the release of a fixed volume of fluid. The space-time development of the gravity current is discussed for different flow behavior indexes.     

Control of unsteady vortical lift on an airfoil by leading-edge blowing-suction

The effects of leading-edge blowing-suction on the vortex flow past an airfoil at high incidence are investigated numerically by solving the Navier-Stokes equations. The results indicate that the frequency of the flowfield excited by the periodic blowing-suction locks into the forcing frequency, which is half of the dominant frequency for the flow past a fixed airfoil without injection. In that case, a well-developed primary leading-edge vortex occupies the upper surface of the airfoil and the largest lift augmentation is obtained. The project supported by the National Defence Research Fund of China     

An analytical theory of heated duct flows in supersonic combustors

One-dimensional analytical theory is developed for supersonic duct flow with variation of cross section, wall friction, heat addition, and relations between the inlet and outlet flow parameters are obtained. By introducing a selfsimilar parameter, effects of heat releasing, wall friction, and change in cross section area on the flow can be normalized and a self-similar solution of the flow equations can be found. Based on the result of self-similar solution, the sufficient and necessary condition for the occurrence of thermal choking is derived. A relation of the maximum heat addition leading to thermal choking of the duct flow is derived as functions of area ratio, wall friction, and mass addition, which is an extension of the classic Rayleigh flow theory, where the effects of wall friction and mass addition are not considered. The present work is expected to provide fundamentals for developing an integral analytical theory for ramjets and scramjets.     

A Computational Method for Free-Surface Viscous Flows Through Porous Media

In this paper we have investigated free-surface unsteady viscous flow through porous media using a finite difference technique via primitive formulation. A transformation has been used to convert the irregular flow domain into a rectangular domain. Two illustrations have been presented in case of a standing wave problem and a cavity flow problem with a free surface. The method is found satisfactory for numerical solution for Reynolds number upto 200.     

Thermo-electromechanical behavior of functionally graded piezoelectric hollow cylinder under non-axisymmetric loads

This paper presents an analytical solution of a thick walled cylinder com- posed of a functionally graded piezoelectric material (FGPM) and subjected to a uniform electric field and non-axisymmetric thermo-mechanical loads. All material properties, except Poisson's ratio that is assumed to be constant, obey the same power law. An exact solution for the resulting Navier equations is developed by the separation of variables and complex Fourier series. Stress and strain distributions and a displacement field through the cylinder are obtained by this technique. To examine the analytical approach, different examples are solved by this method, and the results are discussed.     

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.     

Free surface deformation due to a source or a sink in Stokes flow

Free surface shape and cusp formation are analyzed by considering two-dimensional viscous flow due to a line source or a line sink below the free surface where the strength of source/sink is given arbitrarily. In the analysis, the Stokes' approximation is used and surface tension effects are included, but gravity is neglected. The solution is obtained analytically by using conformal mapping and complex function theory. From the solution, shapes of the free surface are shown and the formation of a cusp on the free surface is discussed. As the capillary number decreases in negative, the free surface shape becomes singular and in a real fluid a cusp should form on the free surface below some negative critical capillary number. Typically, streamline patterns for some capillary numbers are also shown. As the small capillary number vanishes, the solution is reduced to a linearized potential flow solution.