Effect of viscous dissipation on nonisothermal high-viscosity magma flow in a volcanic conduit

A two-dimensional nonisothermal model of magma flow in a volcanic conduit is proposed. The model makes it possible to investigate the effect of the processes of viscous dissipation and heat conduction on the magma flow. It is established that the effect of these processes is significant, particularly in the case of high flow rates. It is shown that in this case the conduit resistance calculated from the Poiseuille formula widely used in one-dimensional models is highly overestimated. This is related to the formation of a strongly heated fluid layer with reduced viscosity in the near-wall conduit zone. Within the framework of the proposed model it is possible to describe eruptions with flow rates which are several times higher than the flow rates obtained within the framework of one-dimensional models.Translated from Izvestiya Rossiiskoi Academii Nauk, Mekhanika Zhidkosti i Gaza, No. 6, 2004, pp. 21–32. Original Russian Text Copyright © 2004 by Barmin, Vedeneeva, and Melnik.     

Indentation of a punch into an elastic strip with friction and adhesion

In the contact interaction between elastic bodies with friction taken into account, the contact region splits, as a rule, into adhesion and sliding regions {xc[1]}. Contact with adhesion and sliding was first considered by L. A. Galin {xc[2]} in the problem of indentation of a punch with a rectilinear foundation into an elastic half-plane, who obtained an approximate solution of this problem [{xc2}, {xc3}]. Galin's problem was further studied in [{xc4}–{xc9}].     

Investigation of a Gas Flow with Solid Particles in a Supersonic Nozzle

A two-phase flow with high Reynolds numbers in the subsonic, transonic, and supersonic parts of the nozzle is considered within the framework of the Prandtl model, i.e., the flow is divided into an inviscid core and a thin boundary layer. Mutual influence of the gas and solid particles is taken into account. The Euler equations are solved for the gas in the flow core, and the boundary-layer equations are used in the near-wall region. The particle motion in the inviscid region is described by the Lagrangian approach, and trajectories and temperatures of particle packets are tracked. The behavior of particles in the boundary layer is described by the Euler equations for volume-averaged parameters of particles. The computed particle-velocity distributions are compared with experiments in a plane nozzle. It is noted that particles inserted in the subsonic part of the nozzle are focused at the nozzle centerline, which leads to substantial flow deceleration in the supersonic part of the nozzle. The effect of various boundary conditions for the flow of particles in the inviscid region is considered. For an axisymmetric nozzle, the influence of the contour of the subsonic part of the nozzle, the loading ratio, and the particle diameter on the particle-flow parameters in the inviscid region and in the boundary layer is studied. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 46, No. 6, pp. 65–77, November–December, 2005.     

Laminar axisymmetric flow of a weakly compressible viscoelastic fluid

The combined effects of weak compressibility and viscoelasticity in steady, isothermal, laminar axisymmetric Poiseuille flow are investigated. Viscoelasticity is taken into account by employing the Oldroyd-B constitutive model. The fluid is assumed to be weakly compressible with a density that varies linearly with pressure. The flow problem is solved using a regular perturbation scheme in terms of the dimensionless isothermal compressibility parameter. The sequence of partial differential equations resulting from the perturbation procedure is solved analytically up to second order. The two-dimensional solution reveals the effects of compressibility and the other dimensionless numbers and parameters in the flow. Expressions for the average pressure drop, the volumetric flow rate, the total axial stress, as well as for the skin friction factor are also derived and discussed. The validity of other techniques used to obtain approximate solutions of weakly compressible flows is also discussed in conjunction with the present results.     

Periodic structures developing with account for dispersion in a turbulence model

Wave processes are studied within the framework of a turbulence model that describes the reaction-diffusion processes in physicochemical hydrodynamics [{xc1}–{xc5}]. For certain parameters of the equation, exact analytical traveling-wave solutions in the form of kinks are obtained. In the general case, the wave processes can be analyzed using numerical simulation. It is confirmed that for a zero dispersion coefficient the nonlinear wave processes are disordered. It is established that when the dispersion terms are taken into account, as for the Kuramoto-Sivashinsky equation, periodic structures develop in the system starting from a certain threshold dispersion-coefficient value.     

Unsteady one-dimensional water and steam flows through a porous medium with allowance for phase transitions

Fluid flow through a porous medium is considered with allowance for heat conduction and phase transition processes. The one-dimensional problem of the breakdown of an arbitrary discontinuity is solved with reference to the processes of combined nonisothermal water and steam flow through the porous medium. It is assumed that there are two-phase zones of water and steam flow through the porous medium to the left and right of the initial discontinuity. Six qualitatively different discontinuous solutions with internal single-phase water or steam zones are constructed and domains corresponding to each of the solutions are found in the determining parameter space. For the parameters considered a solution of the breakdown problem exists and is unique when the requirements for the existence of a discontinuity structure are satisfied [{xc1}].     

MHD stagnation point flow towards heated shrinking surface subjected to heat generation/absorption

The magnetohydrodynamic(MHD) stagnation point flow of micropolar fluids towards a heated shrinking surface is analyzed.The effects of viscous dissipation and internal heat generation/absorption are taken into account.Two explicit cases,i.e.,the prescribed surface temperature(PST) and the prescribed heat flux(PHF),are discussed.The boundary layer flow and energy equations are solved by employing the homotopy analysis method.The quantities of physical interest are examined through the presentation of plots/tabulated values.It is noticed that the existence of the solutions for high shrinking parameters is associated closely with the applied magnetic field.     

Fluid Flows past Spheroids at Moderate Reynolds Numbers

Axisymmetric viscous fluid flows past spheroids are considered. The time-independent complete Navier-Stokes equations written in a spherical coordinate system are used for describing the flow. The problem is solved by the stabilization method on the basis of a variable direction scheme. The input coordinate system is transformed in order to construct a regular computational grid. As a result of the numerical investigation, the stream patterns of flow past elongated and oblate spheroids are obtained for various values of the determining parameters. Numerical values of the dimensions of the circulation zone and the drag coefficient are given for various values of the spheroid semi-axis ratio in the domain of moderate Reynolds numbers.     

MHD stagnation-point flow and heat transfer towards stretching sheet with induced magnetic field

The problem of the steady magnetohydrodynamic (MHD) stagnation-point flow of an incompressible viscous fluid over a stretching sheet is studied. The effect of an induced magnetic field is taken into account. The nonlinear partial differential equations are transformed into ordinary differential equations via the similarity transformation. The transformed boundary layer equations are solved numerically using the shooting method. Numerical results are obtained for various magnetic parameters and Prandtl numbers. The effects of the induced magnetic field on the skin friction coefficient, the local Nusselt number, the velocity, and the temperature profiles are presented graphically and discussed in detail.     

配置点谱方法-人工压缩法(SCM-ACM)求解同心圆筒内流体流动

发展了配置点谱方法SCM(Spectral collocation method)和人工压缩法ACM(Artificial compressibility method)相结合的SCM-ACM数值方法,计算了柱坐标系下稳态不可压缩流动N-S方程组。选取典型的同心圆筒间旋转流动Taylor-Couette流作为测试对象,首先,采用人工压缩法获得人工压缩格式的非稳态可压缩流动控制方程;再将控制方程中的空间偏微分项用配置点谱方法进行离散,得到矩阵形式的代数方程;编写了SCM-ACM求解不可压缩流动问题的程序;最后,通过与公开发表的Taylor-Couette流的计算结果对比,验证了求解程序的有效性。结果证明,本文发展的SCM-ACM数值方法能够用于求解圆筒内不可压缩流体流动问题,该方法既保留了谱方法指数收敛的特性,也具有ACM形式简单和易于实施的特点。本文发展的SCM-ACM数值方法为求解柱坐标下不可压缩流体流动问题提供了一种新的选择。     

Application of a thermodynamically compatible two‐phase flow model to the high‐resolution simulations of compressible gas–magma flow

This paper reports on the application and development of a fully hyperbolic and fully conservative two‐phase flow model for the simulation of gas and magma flow within volcanic processes. The model solves a set of mixture conservation equations for the gas and magma two‐phase flow with velocity non‐equilibrium. In this model, the effect of the relative velocity is introduced by a kinetic constitutive equation with other equations for volume and mass fractions of the gas phase. The model is examined numerically by the widely used finite volume Godunov methods of centered‐type. Using the Riemann problem, we numerically simulate wave propagation and the development of shocks and rarefactions in volcanic eruptions. These simulations are of magma fragmentation type where the relative velocity continues to dominate. A series of test cases whose solution contains features relevant to gas–magma mixtures are conducted. In particular, numerical results indicate that the model implementation predicts key features of the relative velocity within volcanic processes without any mathematical or physical simplifications. Simulation results are sharply and accurately provided without any spurious oscillations in all of the flow variables. The numerical methods and results are also compared with other numerical methods available in the literature. It is found that the provided resolutions are more accurate for the considered test cases. Copyright © 2014 John Wiley & Sons, Ltd.     

On one model of cyclic ejections of magma in the process of the explosive volcanic eruption

The dynamics of the structure of a cavitating magma flow behind the decompression wave front is experimentally studied by the method of hydrodynamic shock tubes. It is demonstrated that a discrete system of intensely cavitating zones with alternation of low and high densities of the gas phase can be formed at a certain regime of shock-wave loading of the examined fluid sample. Based on the results of a numerical analysis of the formation of an anomalous zone in the cavitating magma flow with anomalously high flow characteristics exceeding the values of these characteristics outside this zone at least by an order of magnitude, a model of an instantaneous transformation of the cavitating magma in the anomalous zone to a gas-droplet system, its ejections, and formation of a free surface on the interface is proposed. A numerical analysis shows that the characteristic wave structure and the anomalous saturation zone are fairly rapidly reconstructed in the vicinity of this free surface of the flow part remaining in the conduit after the ejections, and the above-mentioned jumps of flow characteristics are again formed in the anomalous zone.     

A Study of Magma Flows in the Case of Non-Equilibrium Diffusion of Water in the Melt

All volcanic eruptions are accompanied by the degasification of the magma, which results in the growth of gas bubbles in the silicate melt. The number and growth rate of the bubbles determine the character of the eruption. When the free-gas concentration is low, the eruption is weak and takes a short time. At high concentrations, the melt is fragmented and a gas-particle mixture is formed. This results in the catastrophic intensification of the eruption. In this study, for describing the magma flow in a volcano conduit with account for the mass transfer between the bubbles and the melt, we construct models in which the conduit magma flow is considered simultaneously with the dynamics of gas-bubble growth in the melt. The models describe the magma eruptions with equilibrium or weak growth of the bubbles and also with a moderate rate of bubble growth, which is more typical of volcanic systems. Using the models constructed, an intense steady-state eruption and the problem of the evolution of the eruption to the steady-state regime after the rupture of a plug near the top of the conduit are considered. The effect of gas diffusion on the intensity and duration of the eruption is indicated.     

Gasified magma flow in a volcanic conduit

A closed model for describing the magma flow in a volcanic conduit in the eruption regime characterized by quasisteady venting of a gas suspension jet is proposed.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 5, pp. 35–43, September–October, 1990.     

How to predict viscous fingering in three component flow

In a WAG process (Water Alternate Gas), water and a miscible solvent (gas) are injected into a reservoir containing water and oil. The solvent will finger through the oil, leading to early breakthrough and poor recovery. Compared with a miscible flood, when only solvent is injected, fingering is supressed by the simultaneous injection of water, since this reduces the apparent mobility contrast between the injected and displaced fluids. The fingering in a miscible flood, with only hydrocarbon flowing, can be modelled successfully using a Todd and Longstaff fractional flow. In this paper, we demonstrate how to modify the effective Todd and Longstaff mobility ratio self-consistently to account for fingering in three component systems. The resultant empirical equations of flow are solved exactly in one dimension and are in excellent agreement with the averaged saturation and concentration profiles computed using two dimensional high resolution simulation, for a variety of injected water saturations, in both secondary and tertiary displacements.     

Magnetohydrodynamic effects on blood flow through an irregular stenosis

This paper investigates the magnetohydrodynamic (MHD) effects on the blood flow. Rheological properties of blood have been taken into account through the constitutive equations of a micropolar fluid. Unsteady nonlinear differential system is solved numerically by employing finite difference method. Explicit results of axial velocity, flow rate and wall shear stress are obtained and analyzed. It is found that an applied magnetic field reduces the blood flow rate. Copyright © 2010 John Wiley & Sons, Ltd.     

Scaling group transformation for MHD boundary layer flow over permeable stretching sheet in presence of slip flow with Newtonian heating effects

Taking into account the slip flow effects, Newtonian heating, and thermal radiation, two-dimensional magnetohydrodynamic （MHD） flows and heat transfer past a permeable stretching sheet are investigated numerically. We use one parameter group transformation to develop similarity transformation. By using the similarity transformation, we transform the governing boundary layer equations along with the boundary conditions into ordinary differential equations with relevant boundary conditions. The obtained ordinary differential equations are solved with the fourth-fifth order Runge-Kutta- Fehlberg method using MAPLE 13. The present paper is compared with a published one. Good agreement is obtained. Numerical results for dimensionless velocity, temperature distributions, skin friction factor, and heat transfer rates are discussed for various values of controlling parameters.     

A Derivation of the Volume-Averaged Boussinesq Equations for Flow in Porous Media with Viscous Dissipation

The volume-averaged equations are derived for convective flow in porous media. In the thermal energy equation viscous dissipation is taken into account, and a suitable form is obtained which is valid when Brinkman effects are significant.