The use of drift-flux techniques for the analysis of horizontal two-phase flows

New data is presented for horizontal air/water two-phase flow having various flow regimes. It is shown that drift-flux models are able to correlate these data and that the drift velocity, V_gj, is normally finite.     

The Riemann problem for a hyperbolic model of two-phase flow in conservative form

This article is to continue the present author's work (International Journal of Computational Fluid Dynamics (2009) 23 (9), 623–641) on studying the structure of solutions of the Riemann problem for a system of three conservation laws governing two-phase flows. While existing solutions are limited and found quite recently for the Baer and Nunziato equations, this article presents the first instance of an exact solution of the Riemann problem for two-phase flow in gas–liquid mixture. To demonstrate the structure of the solution, we use a hyperbolic conservative model with mechanical equilibrium and without velocity equilibrium. The Riemann problem solution for the model equations comprises a set of elementary waves, rarefaction and discontinuous waves of various types. In particular, such a solution treats both the wave structure and the intermediate states of the two-phase gas–liquid mixture. The resulting exact Riemann solver is fully non-linear, direct and complete. On this basis then, we use locally the exact Riemann solver for the two-phase flow in gas–liquid mixture within the framework of finite volume upwind Godunov methods. In order to demonstrate the effectiveness and accuracy of the proposed solver, we consider a series of test problems selected from the open literature and compare the exact and numerical results by using upwind Godunov methods, showing excellent oscillation-free results in two-phase fluid flow problems.     

The Riemann problem for non-ideal isentropic compressible two phase flows

We consider the Riemann problem for a five-equation, two-pressure (5E2P) model of non-ideal isentropic compressible gas–liquid two-phase flows. This system is more complex due to the extended thermodynamics model for van der Waals gases, that is, typical real gases for gas phase and Tait׳s equation of state for liquid phase. The overall model is strictly hyperbolic and non-conservative form. We investigate the structure of Riemann problem and construct the solution for it. To construct solution of Riemann problem approximately assuming that all waves corresponding to the genuinely non-linear characteristic fields are rarefaction and then we discuss their properties. Lastly, we discuss numerical examples and study the solution influenced by the van der Waals excluded volume.     

Solution of Riemann problem for dusty gas flow

A direct approach is used to solve the Riemann problem for a quasilinear hyperbolic system of equations governing the one dimensional unsteady planar flow of an isentropic, inviscid compressible fluid in the presence of dust particles. The elementary wave solutions of the Riemann problem, that is, shock waves, rarefaction waves and contact discontinuities are derived and their properties are discussed for a dusty gas. The generalised Riemann invariants are used to find the solution between rarefaction wave and the contact discontinuity and also inside rarefaction fan. Unlike the ordinary gasdynamic case, the solution inside the rarefaction waves in dusty gas cannot be obtained directly and explicitly; indeed, it requires an extra iteration procedure. Although the case of dusty gas is more complex than the ordinary gas dynamics case, all the parallel results for compressive waves remain identical. We also compare/contrast the nature of the solution in an ordinary gasdynamics and the dusty gas flow case.     

An extension of Osher's Riemann solver for chemical and vibrational non-equilibrium gas flows

In this paper we study an extension of Osher's Riemann solver to mixtures of perfect gases whose equation of state is of the form encountered in hypersonic applications. As classically, one needs to compute the Riemann invariants of the system to evaluate Osher's numerical flux. For the case of interest here it is impossible in general to derive simple enough expressions which can lead to an efficient calculation of fluxes. The key point here is the definition of approximate Riemann invariants to alleviate this difficulty. Some of the properties of this new numerical flux are discussed. We give 1D and 2D applications to illustrate the robustness and capability of this new solver. We show by numerical examples that the main properties of Osher's solver are preserved; in particular, no entropy fix is needed even for hypersonic applications.     

On the composite Riemann problem for multi‐material fluid flows

We develop an Eulerian fixed grid numerical method for calculating multi‐material fluid flows. This approach relates to the class of interface capturing methods. The fluid is treated as a heterogeneous mixture of constituent materials, and the material interface is implicitly captured by a region of mixed cells that have arisen owing to numerical diffusion. To suppress this numerical diffusion, we propose a composite Riemann problem (CRP), which describes the decay of an initial discontinuity in the presence of a contact point between two different fluids, which is located off the initial discontinuity point. The solution to the CRP serves to calculate multi‐material no mixed numerical flux without introducing any material diffusion. We discuss the CRP solution and its implementation in the multi‐material fluid Godunov method. Numerical results show that a simple framework of the CRP greatly improves capturing material interfaces in the Godunov method and reproduces many of the advantages of more complicated interface tracking multi‐material treatments. Copyright © 2014 John Wiley & Sons, Ltd.     

Modified ghost fluid method for three-dimensional compressible multimaterial flows with interfaces exhibiting large curvature and topological change

In order to capture the material interface dynamics, especially under the impact of strong shocks, the key feature of the modified ghost fluid method (MGFM) is to construct a multimaterial Riemann problem normal to the interface and use its solution to define interface conditions. However, such process sometimes may not be easily or accurately implemented when the multidimensional interfaces come into contact or undergo significant deformations. In this article, a three-dimensional interface treating procedure is developed for a wide range of compressible multimaterial flows. It utilizes the MGFM with an explicit approximate Riemann problem solver to define interface conditions. More importantly, a weighted average technique is designed to optimize the treatment for interfaces exhibiting large curvature and topological change. This remedies two defects of the traditional approach in these extreme cases. One is that the normal directions of interfacial ghost nodes may not be easily calculated. The other is that the interface conditions may not be accurately defined. The numerical methodology is validated through several typical problems, including gas-liquid Riemann problem and shock-bubble/droplet interaction. These results indicate that the developed method is capable of dealing with interfacial evolutions in three dimensions, especially when interfaces undergo merger, fragmentation, and other complex changes.     

An approximate solution of the Riemann problem for a realisable second-moment turbulent closure

Abstract. An approximate solution of the Riemann problem associated with a realisable and objective turbulent second-moment closure, which is valid for compressible flows, is examined. The main features of the continuous model are first recalled. An entropy inequality is exhibited, and the structure of waves associated with the non-conservative hyperbolic convective system is briefly described. Using a linear path to connect states through shocks, approximate jump conditions are derived, and the existence and uniqueness of the one-dimensional Riemann problem solution is then proven. This result enables to construct exact or approximate Riemann-type solvers. An approximate Riemann solver, which is based on Gallou?t's recent proposal is eventually presented. Some computations of shock tube problems are then discussed. Received 2 March 1999 / Accepted 24 August 2000     

On the Riemann problem for liquid or gas-liquid media

Self-similar solutions to the Riemann problem for water with the modified Tait equation of state are presented. The methods of Smoller for gas dynamics are employed to reduce the problem to the solution of a single non-linear equation. The same methods are used for solving the Riemann problem at a gas-water interface. In both cases the method of interval bisections affords a solution technique free of problems with convergence.     

Generalized Riemann problem for gas dynamic combustion

The generalized Riemann problem for gas dynamic combustion in a neighborhood of the origin t > 0 in the (x, t) plane is considered. Under the modified entropy conditions, the unique solutions are constructed, which are the limits of the selfsimilar Zeldovich-von Neumann-Dring (ZND) combustion model. The results show that, for some cases, there are intrinsical differences between the structures of the perturbed Riemann solutions and the corresponding Riemann solutions. Especially, a strong detonation in the...     

A two-scale second-order moment two-phase turbulence model for simulating dense gas-particle flows

A two-scale second-order moment two-phase turbulence model accounting for inter-particle collision is developed, based on the concepts of particle large-scale fluctuation due to turbulence and particle small-scale fluctuation due to collision and through a unified treatment of these two kinds of fluctuations. The proposed model is used to simulate gas-particle flows in a channel and in a downer. Simulation results are in agreement with the experimental results reported in references and are near the results obtained using the single-scale second-order moment two-phase turbulence model superposed with a particle collision model (USM-θ model) in most regions. The project supported by the Special Funds for Major State Basic Research, China (G-1999-0222-08), and the Postdoctoral Science Foundation (2004036239) The English text was polished by Keren Wang     

High-resolution algorithms of sharp interface treatment for compressible two-phase flows

In this paper, a kind of arbitrary high order derivatives (ADER) scheme based on the generalised Riemann problem is proposed to simulate multi-material flows by a coupling ghost fluid method. The states at cell interfaces are reconstructed by interpolating polynomials which are piece-wise smooth functions. The states are treated as the equivalent of the left and right states of the Riemann problem. The contact solvers are extrapolated in the vicinity of contact points to facilitate ghost fluids. The numerical method is applied to compressible flows with sharp discontinuities, such as the collision of two fluids of different physical states and gas–liquid two-phase flows. The numerical results demonstrate that unexpected physical oscillations through the contact discontinuities can be prevented effectively and the sharp interface can be captured efficiently.     

A contribution to the great Riemann solver debate

The aims of this paper are threefold: to increase the level of awareness within the shock-capturing community of the fact that many Godunov-type methods contain subtle flaws that can cause spurious solutions to be computed; to identify one mechanism that might thwart attempts to produce very-high-resolution simulations; and to proffer a simple strategy for overcoming the specific failings of individual Riemann solvers.     

The Riemann problem for nonlinear degenerate wave equations

This paper studies the Riemann problem for a system of nonlinear degenerate wave equations in elasticity. Since the stress function is neither convex nor concave, the shock condition is degenerate. By introducing a degenerate shock under the generalized shock condition, the global solutions are constructively obtained case by case.     

Derivation of a well-posed and multidimensional drift-flux model for boiling flows

In this Note, we derive a multidimensional drift-flux model for boiling flows. Within this framework, the distribution parameter is no longer a scalar but a tensor that might account for the medium anisotropy and the flow regime. A new model for the drift-velocity vector is also derived. It intrinsically takes into account the effect of the friction pressure loss on the buoyancy force. On the other hand, we show that most drift-flux models might exhibit a singularity for large void fraction. In order to avoid this singularity, a remedy based on a simplified three field approach is proposed. To cite this article: O. Grégoire, M. Martin, C. R. Mecanique 333 (2005).     

Explicit solution to the exact Riemann problem and application in nonlinear shallow‐water equations

The Riemann solver is the fundamental building block in the Godunov‐type formulation of many nonlinear fluid‐flow problems involving discontinuities. While existing solvers are obtained either iteratively or through approximations of the Riemann problem, this paper reports an explicit analytical solution to the exact Riemann problem. The present approach uses the homotopy analysis method to solve the nonlinear algebraic equations resulting from the Riemann problem. A deformation equation defines a continuous variation from an initial approximation to the exact solution through an embedding parameter. A Taylor series expansion of the exact solution about the embedding parameter provides a series solution in recursive form with the initial approximation as the zeroth‐order term. For the nonlinear shallow‐water equations, a sensitivity analysis shows fast convergence of the series solution and the first three terms provide highly accurate results. The proposed Riemann solver is implemented in an existing finite‐volume model with a Godunov‐type scheme. The model correctly describes the formation of shocks and rarefaction fans for both one and two‐dimensional dam‐break problems, thereby verifying the proposed Riemann solver for general implementation. Copyright © 2007 John Wiley & Sons, Ltd.     

Direct Riemann solvers for the time-dependent Euler equations

Approaches for finding direct, approximate solutions to the Riemann problem are presented. These result in three approximate Riemann solvers. Here we discuss the time-dependent Euler equations but the ideas are applicable to other systems. The approximate solvers are (i) assessed on local Riemann problems with exact solutions and (ii) used in conjunction with the Weighted Average Flux (WAF) method to solve the two-dimensional Euler equations numerically. The resulting numerical technique is assessed on a shock reflection problem. Comparison with experimental observation is carried out.     

A numerical method for one‐dimensional compressible multiphase flows on moving meshes

This paper is devoted to the numerical approximation of a hyperbolic non‐equilibrium multiphase flow model with different velocities on moving meshes. Such goal poses several difficulties. The presence of different flow velocities in conjunction with cell velocities poses difficulties for upwinding fluxes. Another issue is related to the presence of non‐conservative terms. To solve these difficulties, the discrete equations method (J. Comput. Phys. 2003; 186 (2):361–396; J. Fluid. Mech. 2003; 495 :283–321; J. Comput. Phys. 2004; 196 :490–538; J. Comput. Phys. 2005; 205 :567–610) is employed and generalized to the context of moving cells. The complementary conservation laws, available for the mixture, are used to determine the velocities of the cells boundaries. With these extensions, an accurate and robust multiphase flow method on moving meshes is obtained and validated over several test problems with exact or experimental solutions. Copyright © 2007 John Wiley & Sons, Ltd.     

Explicit algebraic Reynolds stress model for shock-dominated flows

Shock waves drastically alter the nature of Reynolds stresses in a turbulent flow, and conventional turbulence models cannot reproduce this effect. In the present study, we employ explicit algebraic Reynolds stress model (EARSM) to predict the Reynolds stress anisotropy generated by a shockwave. The model by Wallin and Johansson (2000) is used as the baseline model. It is found to over-predict the post-shock Reynolds stresses in canonical shock turbulence interaction. The budget of the transport equation of Reynolds stresses computed using linear interaction analysis shows that the unsteady shock distortion mechanism and the pressure–velocity correlations are important. We propose improvement to the baseline model using linear interaction analysis results and redistribute the turbulent kinetic energy between the principle Reynolds stresses. The new model matches DNS data for the amplification of Reynolds stresses across the shock and their post-shock evolution, for a range of Mach numbers. It is applied to oblique shock/boundary-layer interaction at Mach 5. Significant improvements are observed in predicting surface pressure and skin friction coefficient, with respect to experimental measurements.     

低浓度固液两相流的颗粒相动理学模型

用广义Fokker-Planck扩散模型描述液相湍动对颗粒的挟带作用，用修正的BGK模型描述粒间碰撞效应，建立了封闭的颗粒相PDF输运方程．运用Chapman-Enskog迭代法求得方程的二阶近似解，获得颗粒相脉动速度二阶矩和三阶矩闭合关系．模型与颗粒流模型相容，与液相湍流闭合模型是否相容依赖于扩散模型的具体形式，并据此比较了不同的涡一颗粒作用模型．模型与二维明渠流轻质沙和天然沙试验资料符合很好．表明细小粒径颗粒能够充分跟随水流运动；大粒径颗粒的相间平均速度差和壁面滑移速度明显，近壁区内的颗粒沿流向和垂向脉动强度都可能大于水流，并存在一定程度的颗粒碰撞效应．