Front Tracking Approach and Application to Multi-Fluid Flow

Author(s):  
Manasa Ranjan Behera ◽  
K. Murali

Multiphase flows simulations using a robust interface-tracking method, are presented. The method is based on writing one set of governing equations for the whole computational domain and treating the different phases as single fluid domain with variable material properties. Interfacial terms are accounted for by adding the appropriate sources as δ functions at the boundary separating the phases. The unsteady Navier-Stokes equations are solved by finite volume method on a fixed, structured grid and the interface, or front, is tracked explicitly by a lower dimensional grid. Interfacial source terms are computed on the front and transferred to the fixed grid. Advection of fluid properties such as density and viscosity is done by following the motion of the front. The method has been implemented for interfacial flow problems, depicting the interface and topology change capturing capability. The representation of the moving interface and its dynamic restructuring, as well as the transfer of information between the moving front and the fixed grid, is discussed. Extensions of the method to density stratified flows, and interfacial movements are then presented.

2005 ◽  
Vol 19 (28n29) ◽  
pp. 1439-1442
Author(s):  
YUZHI SUN ◽  
Z. J. WANG

In this paper, the spectral volume (SV) method is experimented for the Navier-Stokes equations by treating the viscous terms with a mixed formulation named the local Discontinuous Galerkin approach. The formulation of the SV method for the two-dimensional compressible Navier-Stokes equations is described. The solver is used to solve several viscous flow problems to show its potential.


Processes ◽  
2021 ◽  
Vol 9 (5) ◽  
pp. 841
Author(s):  
Yuzhen Jin ◽  
Huang Zhou ◽  
Linhang Zhu ◽  
Zeqing Li

A three-dimensional numerical study of a single droplet splashing vertically on a liquid film is presented. The numerical method is based on the finite volume method (FVM) of Navier–Stokes equations coupled with the volume of fluid (VOF) method, and the adaptive local mesh refinement technology is adopted. It enables the liquid–gas interface to be tracked more accurately, and to be less computationally expensive. The relationship between the diameter of the free rim, the height of the crown with different numbers of collision Weber, and the thickness of the liquid film is explored. The results indicate that the crown height increases as the Weber number increases, and the diameter of the crown rim is inversely proportional to the collision Weber number. It can also be concluded that the dimensionless height of the crown decreases with the increase in the thickness of the dimensionless liquid film, which has little effect on the diameter of the crown rim during its growth.


2010 ◽  
Vol 297-301 ◽  
pp. 924-929
Author(s):  
Inès Bhouri Baouab ◽  
Nejla Mahjoub Said ◽  
Hatem Mhiri ◽  
Georges Le Palec ◽  
Philippe Bournot

The present work consists in a numerical examination of the dispersion of pollutants discharged from a bent chimney and crossing twin similar cubic obstacles placed in the lee side of the source. The resulting flow is assumed to be steady, three-dimensional and turbulent. Its modelling is based upon the resolution of the Navier Stokes equations by means of the finite volume method together with the RSM (Reynolds Stress Model) turbulent model. This examination aims essentially at detailing the wind flow perturbations, the recirculation and turbulence generated by the presence of the twin cubic obstacles placed tandem at different spacing distances (gaps): W = 4 h, W = 2 h and W = 1 h where W is the distance separating both buildings.


1998 ◽  
Vol 14 (3) ◽  
pp. 153-159 ◽  
Author(s):  
Chou-Jiu Tsai ◽  
Ger-Jyh Chen

ABSTRACTIn this study, fluid flow around bluff bodies are studied to examine the vortex shedding phenomenon in conjuction with the geometrical shapes of these vortex shedders. These flow phenomena are numerically simulated. A finite volume method is employed to solve the incompressible two-dimensional Navier-Stokes equations. Thus, quantitative descriptions of the vortex shedding phenomenon in the near wake were made, which lead to a detailed description of the vortex shedding mechanism. Streamline contours, figures of lift coefficent, and figures of drag coefficent in various time, are presented, respectively, for a physical description.


Author(s):  
Elisabeth Longatte

This work is concerned with the modelling of the interaction of a fluid with a rigid or a flexible elastic cylinder in the presence of axial or cross-flow. A partitioned procedure is involved to perform the computation of the fully-coupled fluid solid system. The fluid flow is governed by the incompressible Navier-Stokes equations and modeled by using a fractional step scheme combined with a co-located finite volume method for space discretisation. The motion of the fluid domain is accounted for by a moving mesh strategy through an Arbitrary Lagrangian-Eulerian (ALE) formulation. Solid dyncamics is modeled by a finite element method in the linear elasticity framework and a fixed point method is used for the fluid solid system computation. In the present work two examples are presented to show the method robustness and efficiency.


Author(s):  
Wolfgang Höhn

During the design of the compressor and turbine stages of today’s aeroengines, aerodynamically induced vibrations become increasingly important since higher blade load and better efficiency are desired. In this paper the development of a method based on the unsteady, compressible Navier-Stokes equations in two dimensions is described in order to study the physics of flutter for unsteady viscous flow around cascaded vibrating blades at stall. The governing equations are solved by a finite difference technique in boundary fitted coordinates. The numerical scheme uses the Advection Upstream Splitting Method to discretize the convective terms and central differences discretizing the viscous terms of the fully non-linear Navier-Stokes equations on a moving H-type mesh. The unsteady governing equations are explicitly and implicitly marched in time in a time-accurate way using a four stage Runge-Kutta scheme on a parallel computer or an implicit scheme of the Beam-Warming type on a single processor. Turbulence is modelled using the Baldwin-Lomax turbulence model. The blade flutter phenomenon is simulated by imposing a harmonic motion on the blade, which consists of harmonic body translation in two directions and a rotation, allowing an interblade phase angle between neighboring blades. Non-reflecting boundary conditions are used for the unsteady analysis at inlet and outlet of the computational domain. The computations are performed on multiple blade passages in order to account for nonlinear effects. A subsonic massively stalled unsteady flow case in a compressor cascade is studied. The results, compared with experiments and the predictions of other researchers, show reasonable agreement for inviscid and viscous flow cases for the investigated flow situations with respect to the Steady and unsteady pressure distribution on the blade in separated flow areas as well as the aeroelastic damping. The results show the applicability of the scheme for stalled flow around cascaded blades. As expected the viscous and inviscid computations show different results in regions where viscous effects are important, i.e. in separated flow areas. In particular, different predictions for inviscid and viscous flow for the aerodynamic damping for the investigated flow cases are found.


Sign in / Sign up

Export Citation Format

Share Document