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.     

Asymptotic models of solidification in cooling thin-layer flows of a highly viscous fluid

Asymptotic models are constructed for the solidification process in a highly viscous film flow on the surface of a cone with a given mass supply at the cone apex. In the thin-layer approximation, the problem is reduced to two parabolic equations for the temperatures of the liquid and the solid coupled with an ordinary differential equation for the solidification front. For large Péclet numbers, an analytical steady-state solution for the solidification front is found. A nondimensional parameter which makes it possible to distinguish flows (i) without a solid crust, (ii) with a steady-state solid crust, and (iii) with complete solidification is determined. For finite Péclet numbers and large Stefan numbers, an analytical transient solution is found and the time of complete flow solidification is determined. In the general case, when all the governing parameters are of the order of unity, the original system of equations is studied numerically. The solutions obtained are qualitatively compared with the data of field observations for lava flows produced by extrusive volcanic eruptions.     

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.     

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.     

Flow of a Nonlinearly Viscous Fluid over the Surface of a Rotating Flat Disk

The flow of a nonlinearly viscous (power-law) fluid over the surface of a rotating flat disk is investigated. A solution form which makes it possible to reduce the complete system of partial differential equations to a system of ordinary differential equations is found. This system is integrated using the Runge-Kutta method and reduction to a Cauchy problem on the basis of Newton's method. The velocity and pressure fields in a power-law fluid film flowing over the surface of a rotating flat disk are found numerically.     

Unsteady boundary layer flow due to a stretching surface in a rotating fluid

The induced unsteady flow due to a stretching surface in a rotating fluid, where the unsteadiness is caused by the suddenly stretched surface is studied in this paper. After a similarity transformation, the unsteady Navier–Stokes equations have been solved numerically using the Keller-box method. Also, the perturbation solution for small times as well as the asymptotic solution for large times, when the flow becomes steady, has been obtained. It is found that there is a smooth transition from the small time solution to the large time or steady state solution.     

Simulation of interface and free surface flows in a viscous fluid using adapting quadtree grids

A new adaptive quadtree method for simulating laminar viscous fluid problems with free surfaces and interfaces is presented in this paper. The Navier–Stokes equations are solved with a SIMPLE‐type scheme coupled with the Compressive Interface Capturing Scheme for Arbitrary Meshes (CICSAM) (Numerical prediction of two fluid systems with sharp interfaces, Ph.D. Thesis, Imperial College of Science, Technology and Medicine, London, 1997) volume of fluid (VoF) method and PLIC reconstruction of the volume fraction field during refinement and derefinement processes. The method is demonstrated for interface advection cases in translating and shearing flow fields and found to provide high interface resolution at low computational cost. The new method is also applied to simulation of the collapse of a water column and the results are in excellent agreement with other published data. The quadtree grids adapt to follow the movement of the free surface, whilst maintaining a band of the smallest cells surrounding the surface. The calculation is made on uniform and adapting quadtree grids and the accuracy of the quadtree calculation is shown to be the same as that made on the equivalent uniform grid. Copyright © 2004 John Wiley & Sons, Ltd.     

Submerged shocks in conical gas flows in the presence of a mach shock-wave configuration

The results of a numerical study of a new type of singularities in the Mach shock-wave structure realized in supersonic nonsymmetric conical flows over V-wings with a bow shock attached to the leading edges are presented. Within the framework of the ideal gas model we study the changes in the shock system on transition, with increase in the sweep angle, from the region of nonsymmetric Mach interaction of the shocks attached to the leading edges of the wing to the region of special flow patterns, where on the windward cantilever surface a rarefaction flow is realized rather than a flow with an internal shock. It is shown, in particular, that in the region with special wing flow patterns a Mach system of shocks with a submerged shock proceeding from the branch point above the windward cantilever may exist.     

Formation of the free surface of a viscous fluid volume inside a rotating horizontal cylinder

Two-dimensional viscous flow with a free surface in a horizontal cylinder rotating at a constant speed is investigated numerically using the boundary element method. It is shown that in the initial stage of rotation of the cylinder four different variants of the behavior of the free surface can be realized in the stage of transition from horizontal to steady-state form.     

Three-dimensional flow of Oldroyd-B fluid over surface with convective boundary conditions

The present study addresses the three-dimensional flow of an Oldroyd-B fluid over a stretching surface with convective boundary conditions. The problem formulation is presented using the conservation laws of mass, momentum, and energy. The solutions to the dimensionless problems are computed. The convergence of series solutions by the homotopy analysis method (HAM) is discussed graphically and numerically. The graphs are plotted for various parameters of the temperature profile. The series solutions are verified by providing a comparison in a limiting case. The numerical values of the local Nusselt number are analyzed.     

Slip effects and heat transfer analysis in a viscous fluid over an oscillatory stretching surface

In this study, we investigate the heat transfer problem in a viscous fluid over an oscillatory infinite sheet with slip condition. The sheet is moved back and forth in its own plane. The derived problem involves a dimensionless parameter indicating the relative magnitude of frequency to sheet stretching rate. A system of non‐linear partial differential equations is solved numerically using the finite‐difference scheme, in which a coordinate transformation is employed to transform the semi‐infinite physical space to a bounded computational domain. The physical features of interesting parameters on the velocity and temperature distributions are shown graphically and discussed. The values of the skin‐friction coefficient and the local Nusselt number are given in tabular form. Copyright © 2008 John Wiley & Sons, Ltd.     

Three-dimensional motion of a viscous incompressible fluid in a narrow tube

An approximate solution of the problem of unsteady motion of a viscous incompressible fluid in a long narrow deformable tube at low Reynolds numbers is obtained. Pressure oscillations and tube deformation are shown to be related by an integrodifferential equation. The solution obtained extends the Poiseuille solution in elliptic tubes to the case of comparatively arbitrary small deformations in terms of the tube length and angle. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 50, No. 4, pp. 28–32, July–August, 2009.     

The development of a Cartesian cut cell method for incompressible viscous flows

This paper describes the extension of the Cartesian cut cell method to applications involving unsteady incompressible viscous fluid flow. The underlying scheme is based on the solution of the full Navier–Stokes equations for a variable density fluid system using the artificial compressibility technique together with a Jameson‐type dual time iteration. The computational domain encompasses two fluid regions and the interface between them is treated as a contact discontinuity in the density field, thereby eliminating the need for special free surface tracking procedures. The Cartesian cut cell technique is used for fitting the complex geometry of solid boundaries across a stationary background Cartesian grid which is located inside the computational domain. A time accurate solution is achieved by using an implicit dual‐time iteration technique based on a slope‐limited, high‐order, Godunov‐type scheme for the inviscid fluxes, while the viscous fluxes are estimated using central differencing. Validation of the new technique is by modelling the unsteady Couette flow and the Rayleigh–Taylor instability problems. Finally, a test case for wave run‐up and overtopping over an impermeable sea dike is performed. Copyright © 2006 John Wiley & Sons, Ltd.     

Stability and nonlinear wavy regimes in downward film flows on a corrugated surface

The linear and nonlinear stability of downward viscous film flows on a corrugated surface to freesurface perturbations is analyzed theoretically. The study is performed with the use of an integral approach in ranges of parameters where the calculated results and the corresponding solutions of Navier-Stokes equations (downward wavy flow on a smooth wall and waveless flow along a corrugated surface) are in good agreement. It is demonstrated that, for moderate Reynolds numbers, there is a range of corrugation parameters (amplitude and period) where all linear perturbations of the free surface decay. For high Reynolds numbers, the waveless downward flow is unstable. Various nonlinear wavy regimes induced by varying the corrugation amplitude are determined. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 48, No. 1, pp. 110–120, January–February, 2007.     

Two-dimensional models of magma flows in a volcanic conduit taking the magma compressibility and thermal effects into account

Two steady-state models of magma flow in a conduit are considered, with and without allowance for magma compressibility. As distinct from studies [{xc1}–{xc6}], in which either simplified equations were solved or unrealistic values of the parameters were used, in the present study the complete systems of equations are solved and the values of the parameters correspond to magma flow in a volcanic conduit. The secondary flows obtained in [{xc5}] for model conditions are not formed when the magma is simulated by an incompressible fluid and all the terms of the equations are taken into account. When the magma compressibility is taken into account, in the isothermal case and for constant magma viscosity the entire flow is adequately described by the one-dimensional isothermalmodel, although this approach is not formally applicable.     

Mixed convection unsteady stagnation point flow over a stretching sheet with heat transfer in the presence of variable free stream

An analysis is performed for the unsteady mixed convection flow of an incompressible viscous fluid about a stagnation point on a stretching sheet in the presence of a variable free stream. The equations of motion and energy are transformed into the ordinary differential equations by using similarity transformations. Homotopy analysis method is used for the solution of the governing problem. The results have been discussed by plots. The present values of the function are shown very close to the previous limiting solutions. Copyright © 2011 John Wiley & Sons, Ltd.     

Free surface Stokes flows obstructed by multiple obstacles

Gravity‐driven Stokes flow down an inclined plane over and around multiple obstacles is considered. The flow problem is formulated in terms of a boundary integral equation and solved using the boundary element method. A Hermitian radial basis function (RBF) is used for the interpolation of the free surface, generation of the unit normal and curvature, and to prescribe the far‐field conditions. For flow over an obstacle, hemispheres are taken. For flow around an obstacle, circular cylinders are modelled and the contact angle condition on the obstacle/free surface intersection specified using the RBF formulation. Explicit profiles are produced for flow over and around two obstacles placed in various locations relative to one another. Interaction due to two obstacles is given by comparisons made with the profiles for flow over and around individual obstacles. In general, when the obstacles are separated by a sufficiently large distance the flow profiles are identical to a single obstacle analysis. For flow over and around two obstacles in‐line with the incident flow, effects of the governing parameters are examined, with variations in plane inclination angle, Bond number, obstacle size, and in the case of obstacles intersecting the free surface, static contact angle is considered. Finally flows over and around three obstacles are modelled. Copyright © 2009 John Wiley & Sons, Ltd.     

Three‐dimensional transient free‐surface flow of viscous fluids inside cavities of arbitrary shape

The three‐dimensional transient free‐surface flow inside cavities of arbitrary shape is examined in this study. An adaptive (Lagrangian) boundary‐element approach is proposed for the general three‐dimensional simulation of confined free‐surface flow of viscous incompressible fluids. The method is stable as it includes remeshing capabilities of the deforming free‐surface, and thus can handle large deformations. A simple algorithm is developed for mesh refinement of the deforming free‐surface mesh. Smooth transition between large and small elements is achieved without significant degradation of the aspect ratio of the elements in the mesh. The method is used to determine the flow field and free‐surface evolution inside cubic, rectangular and cylindrical containers. These problems illustrate the transient nature of the flow during the mixing process. Surface tension effects are also explored. Copyright © 2003 John Wiley & Sons, Ltd.