Towards a numerical simulation of a high temperature flow

Two mathematical models of unsteady plasma flows based on 3D finite volume methods for Navier‐Stokes equations are presented. In the first, cold jet model, the nondimensional NS equations are solved by the artificial compressibility method, and the results include time‐varying velocity fields for different Reynolds numbers. The second formulation is a compressible Stegger‐Warming flux vector splitting scheme that takes into account the non‐perfect character of the plasma fluid. Temperature and composition dependent gas properties are to be evaluated from real data.     

On the Crenation of Coated Droplets

The dynamic response of a small, very viscous, coated droplet is examined as additional fluid is deposited into a compressible surface layer. At early times, the surface material behaves like an incompressible fluid and a number of crenations form. These protuberances decrease in size and in number as the compressibility of the coating takes effect, until finally, the droplet again becomes spherical. This sequence of events and its dependence on the compressibility and rheological properties of the fluids are studied. Similar ruffling effects are observed in cell biology, and some of the possible implications of the fluid dynamical model are discussed.     

Numerical simulation of cavitation dynamics using a cavitation-induced-momentum-defect (CIMD) correction approach

A new unsteady cavitation event tracking model is developed for predicting vapor dynamics occurring in multi-dimensional incompressible flows. The procedure solves incompressible Navier–Stokes equations for the liquid phase supplemented with an additional vapor transport equation for the vapor phase. The novel cavitation-induced-momentum-defect (CIMD) correction methodology developed in this study accounts for cavitation inception and collapse events as relevant momentum-source terms in the liquid phase momentum equations. The model tracks cavitation zones and applies compressibility effects, employing homogeneous equilibrium model (HEM) assumptions, in constructing the source term of the vapor transport model. Effects of vapor phase accumulation and diffusion are incorporated by detailed relaxation models. A modified RNG k–ε model, including the effects of compressibility in the vapor regions, is employed for modeling turbulence effects. Numerical simulations are carried out using a finite volume methodology available within the framework of commercial CFD software code Fluent v.6.2. Simulation results are in good qualitative agreement with experiments for unsteady cloud cavitation behavior in planar nozzle flows. Multitude of mechanisms such as formation of vortex cavities, vapor cluster shedding and coalescence, cavity pinch off are sharply captured by the CIMD approach. Our results indicate the profound influence of re-entrant jet motion and adverse pressure gradients on the cavitation dynamics.     

Unsteady flows of a viscoelastic fluid with the fractional Burgers’ model

The aim of this work is to discuss some unidirectional flows of a viscoelastic fluid between two parallel plates with fractional Burgers’ fluid model. The exact analytical solutions for Plane Poiseuille and Plane Couette flows are obtained by using the finite Fourier sine transform and the Laplace transform. Moreover, the graphs are plotted to show the effects of different parameters on the velocity field.     

Spherical shock wave generated by a moving piston in mixture of a non-ideal gas and small solid particles under a gravitational field

Self-similar solutions are obtained for one-dimensional isothermal and adiabatic unsteady flows behind a strong spherical shock wave propagating in a dusty gas. The shock is assumed to be driven out by a moving piston and the dusty gas to be a mixture of a non-ideal (or perfect) gas and small solid particles, in which solid particles are continuously distributed. It is assumed that the equilibrium flow-conditions are maintained and variable energy input is continuously supplied by the piston. The medium is under the influence of the gravitational field due to a heavy nucleus at the origin (Roche model). The effects of an increase in the mass concentration of solid particles, the ratio of the density of the solid particles to the initial density of the gas, the gravitational parameter and the parameter of non-idealness of the gas in the mixture, are investigated. It is shown that due to presence of gravitational field the compressibility of the medium at any point in the flow-field behind the shock decreases and all other flow-variables and the shock strength increase. A comparison has also been made between the isothermal and adiabatic flows. It is investigated that the singularity in the density and compressibility distributions near the piston in the case of adiabatic flow are removed when the flow is isothermal.     

Constraints on equation of state for cavitating flows with thermodynamic effects

Thermodynamic effects play an important role in the cavitation dynamics of cryogenics fluids. Such flows are characterized by strong variations in fluid properties with the temperature. A compressible, multiphase, one-fluid solver was developed to study and to predict thermodynamic effects in cavitating flows. To close the system, a cavitation model is proposed to capture metastable behaviours of fluids and non isothermal thermodynamic path. The thermodynamical consistency based on entropy conditions and the evolution of the mixture speed of sound are investigated. These constraints are applied to other models. The considered working fluid is the refrigerant R-114.     

Modified k – ω Model for Simulation of Cavitating Flows

Although cavitating flows are generally turbulent, the interaction of turbulence and cavitation dynamics is still not well understood. In general, two‐equation models are employed, which were originally developed for single‐phase flows. Therefore they fail by handling cavitation based flow phenomena with very high density variations (dependent on operating condition up to 40000:1). This sudden change of the density causes strong pressure gradients, secondary flows and local compressibility. The aim of this study is to enhance the Wilcox's k‐? model with empirical correlations in order to simulate turbulent cavitating flows more precisely and effciently.     

Numerical Solution of Some Cardiovascular Problem

This paper deals with problem of numerical solution of flows through vessel with bypass. One could use model of Navier‐Stokes equations and find solution by using multistage Runge‐Kutta method together with time dependent artificial compressibility method. Some results of numerical solution of cardiovascular problems are presented: stationary and unstationary 2D flows in a vessel and a bypass. (© 2004 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Application of an exactly incompressible lattice Bhatnagar-Gross-Krook method to internal pressure driven and depth-averaged flows

Modifications to the lattice Bhatnagar-Gross-Krook (BGK) methodeliminate compressibility errors from the simulation of macroscopicfluids governed by the steady-state Navier-Stokes equation.We describe the means by which this modified scheme makes itpossible to compute both pressure drops and flow fields, bringinginternal, pressure driven, isothermal flows within the scopeof the method. The scheme is applied to the simulation of flowof a viscous incompressible fluid past a sudden expansion inan infinite aspect ratio duct (back-facing step geometry). Further,we devise the means by which an exactly incompressible (El)LBGK scheme may be extended and calibrated to calculate depth-averaged flow fields in ducts of uniform depth. By analyticalcalculation and by simulation we demonstrate the validity ofour approach.     

Finite-element stream function solutions of transonic rotational internal and external flows

The finite-element method is applied to the stream function formulation of transonic flows. Numerical dissipation, necessary for the calculation of mixed flows with shocks, is introduced via the artificial compressibility method. The classical problem of double-valuedness of the mass flux versus Mach number is resolved by direct integration of the vorticity equation. Solutions are obtained for isolated airfoils, the blade-to-blade cascade equation, as well as the radial equilibrium equation governing the hub-to-shroud through flow in turbomachinery.     

Numerical Solution of Newtonian and Non-Newtonian Flows

This paper deals with the numerical solution of Newtonian and non-Newtonian flows. The flows are supposed to be laminar, viscous, incompressible and steady. The model used for non-Newtonian fluids is some variant of power-law. Governing equations in this model are incompressible Navier-Stokes equations. For numerical solution one could use artificial compressibility method with three stage Runge-Kutta finite volume method in cell centered formulation for discretization of space derivatives. Following cases of flows are solwed: flow through a bypass connected to main channel in 2D and 3D and non-Newtonian flow through branching channels in 2D. These results are presented for 2D and 3D case. This problem could have an application in the area of biomedicine. (© 2006 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Nano boundary layers over stretching surfaces

In this paper, we present similarity solutions for the nano boundary layer flows with Navier boundary condition. We consider viscous flows over a two-dimensional stretching surface and an axisymmetric stretching surface. The resulting nonlinear ordinary differential equations are solved analytically by the Homotopy Analysis Method. Numerical solutions are obtained by using a boundary value problem solver, and are shown to agree well with the analytical solutions. The effects of the slip parameter K and the suction parameter s on the fluid velocity and on the tangential stress are investigated and discussed. As expected, we find that for such fluid flows at nano scales, the shear stress at the wall decreases (in an absolute sense) with an increase in the slip parameter K.     

Exact solution for rotating flows of a generalized Burgers’ fluid in a porous space

This article considers the oscillatory flows of a generalized Burgers’ fluid on an infinite insulating plate when the fluid is permeated by a transverse magnetic field. The effects of Hall current are taken into account. Modified Darcy’s law for a generalized Burgers’ fluid has been used to discuss the flows in a porous medium. The governing time dependent equations in a rotating frame are first developed and then solved for the two problems. The influence of various emerging parameters is discussed through various graphs. The solutions for the Newtonian, second grade, Maxwell, Oldroyd-B and Burgers’ fluids can be obtained from our solutions as the limiting cases.     

The effect of seepage on the stress–strain state of rock near a borehole

The effect of the strain-strength and seepage properties of rock and the compressibility of the percolating fluid on the dimensions of the rock fracture zones, which occur in oil and gas boreholes when the bottom hole pressure is reduced, is investigated. The seepage is considered basing on the stationary formulation of the problem, which enables the general case to be investigated. It is shown that in the case of unsteady flow, the stresses on the boundary of the rock fracture zone and, as a consequence, on its dimensions, are independent of the nature of the pressure distribution in the stratum, and are determined solely by the pressure of the percolating fluid on the boundary of this zone. It is established that an increase in the compressibility of the percolating fluid leads to an increase in the dimensions of the rock fracture zone.     

Barotropic instability of shear flows

We consider barotropic instability of shear flows for incompressible fluids with Coriolis effects. For a class of shear flows, we develop a new method to find the sharp stability conditions. We study the flow with Sinus profile in details and obtain the sharp stability boundary in the whole parameter space, which corrects previous results in the fluid literature. Our new results are confirmed by more accurate numerical computation. The addition of the Coriolis force is found to bring fundamental changes to the stability of shear flows. Moreover, we study dynamical behaviors near the shear flows, including the bifurcation of nontrivial traveling wave solutions and the linear inviscid damping. The first ingredient of our proof is a careful classification of the neutral modes. The second one is to write the linearized fluid equation in a Hamiltonian form and then use an instability index theory for general Hamiltonian partial differential equations. The last one is to study the singular and nonresonant neutral modes using Sturm-Liouville theory and hypergeometric functions.     

Initial-boundary value problem of three phase flow in porous media

We study the mathematical model of three phase compressible flows through porous media. Under the condition that the rock, water and oil are incompressible, and the compressibility of gas is small, we present a finite element scheme to the initial-boundary value problem of the nonlinear system of equations, then by the convergence of the scheme we prove that the problem admits a weak solution.     

Visualization of vortical flows in computational fluid dynamics

The concepts and methods of the visual representation of fluid dynamics computations of vortical flows are studied. Approaches to the visualization of vortical flows based on the use of various definitions of a vortex and various tests for its identification are discussed. Examples of the visual representation of solutions to some fluid dynamics problems related to the computation of vortical flows in jets, channels, and cavities and of the computation of separated flows occurring in flows around bodies of various shapes are discussed.     

Optimization methods for pipeline transportation of natural gas with variable specific gravity and compressibility

In this paper, the problem of flow maximization in pipeline systems for transmission of natural gas is addressed. We extend previously suggested models by incorporating the variation in pipeline flow capacities with gas specific gravity and compressibility. Flow capacities are modeled as functions of pressure, compressibility and specific gravity by the commonly-used Weymouth equation, and the California Natural Gas Association method is used to model compressibility as a function of specific gravity and pressure. The sources feeding the transmission network do not necessarily supply gas with equal specific gravity. In our model, it is assumed that when different flow streams enter a junction point, the specific gravity of the resulting flow is a weighted average of the specific gravities of entering flows. We also assume the temperature to be constant, and the system to be in steady state. Since the proposed model is non-convex, and global optimization hence can be time consuming, we also propose a heuristic method based on an iterative scheme in which a simpler NLP model is solved in each iteration. Computational experiments are conducted in order to assess the computability of the model by applying a global optimizer, and to evaluate the performance of the heuristic approach. When applied to a wide set of test instances, the heuristic method provides solutions with deviations less than 10% from optimality, and in many instances turns out to be exact. We also report several experiments demonstrating that letting the compressibility and the specific gravity be global constants can lead to significant errors in the estimates of the total network capacity.     

On the pollution effect of quasi‐compressibility methods in magneto‐hydrodynamics and reactive flows

The study of the asymptotics of quasi‐compressibility methods is an essential tool to construct and improve numerical schemes in hydrodynamics, like e.g., stabilized finite element methods or time‐splitting projection methods. The goal of this work is to illustrate different asymptotical solution behaviour for equations in magneto‐hydrodynamics and those that describe reactive flows for standard quasi‐compressibility methods that influences the design of numerical algorithms. Copyright © 1999 John Wiley & Sons, Ltd.