Multi-Scale Near-Wall Averaging and Two-Phase Logarithmic Law for a CFD Two-Fluid Model

Author(s):  
Avinash Vaidheeswaran ◽  
John R. Buchanan ◽  
Paul Guilbert ◽  
Martin Lopez de Bertodano

A considerable amount of work has been done in the past to improve the solution methodology using the two-fluid model in the near-wall region. This includes the works of Larrateguy et al. [1], and Moraga et al. [2], based on a multi-scale bubble-center averaging technique. However one shortcoming is that the primitive variables must be recovered from the bubble-center averaged variables. This makes it difficult to implement it in a commercial CFD code. The current research focuses on an engineering approach to overcome this issue. A multi-scale near-wall averaging technique is proposed which separates the effects of bubble dynamics from its geometry in this region. In addition, the averaged volume fraction profile makes the CFD approach consistent with the modified logarithmic law of Marie et al. [3]. A step function volume fraction distribution was assumed in the near-wall region while developing the theory. However, the volume fraction prediction obtained from CFD calculations is not uniform in this region. The proposed near-wall averaging technique resolves this issue and makes the CFD implementation of the modified wall function approach consistent with the theory of Marie et al. [3].

2016 ◽  
Vol 139 (1) ◽  
Author(s):  
Avinash Vaidheeswaran ◽  
Deoras Prabhudharwadkar ◽  
Paul Guilbert ◽  
John R. Buchanan ◽  
Martin Lopez de Bertodano

A new two-fluid model averaging in the near-wall region is proposed to ensure consistent matching of the two-phase k–ε turbulence model with the two-phase logarithmic law of the wall (Marie J. L., Moursali, E., and Tran-Cong, S., 1997, “Similarity Law and Turbulence Intensity Profiles in a Bubbly Boundary Layer,” Int. J. Multiphase Flow, 23(2), pp. 227–247). The void fraction distribution obtained with the averaging procedure is seen to conform to the two-phase wall function approach which is based on a double step function void fraction distribution. In particular, the proposed averaging technique is shown to achieve grid convergence in the near-wall region, which could not be obtained otherwise. Computational fluid dynamics (CFD) results with the proposed technique are in good agreement with experiments on upward bubbly flows over a flat plate, and upward and downward flows in pipes. An additional advantage of the proposed technique is that it replaces the wall force model, which has a significant degree of uncertainty in turbulent flow modeling, with a simpler geometric constraint.


2014 ◽  
Vol 62 (3) ◽  
pp. 234-240 ◽  
Author(s):  
Gianandrea Vittorio Messa ◽  
Stefano Malavasi

Abstract The flow of a mixture of liquid and solid particles at medium and high volume fraction through an expansion in a rectangular duct is considered. In order to improve the modelling of the phenomenon with respect to a previous investigation (Messa and Malavasi, 2013), use is made of a two-fluid model specifically derived for dense flows that we developed and implemented in the PHOENICS code via user-defined subroutines. Due to the lack of experimental data, the two-fluid model was validated in the horizontal pipe case, reporting good agreement with measurements from different authors for fully-suspended flows. A 3D system is simulated in order to account for the effect of side walls. A wider range of the parameters characterizing the mixture (particle size, particle density, and delivered solid volume fraction) is considered. A parametric analysis is performed to investigate the role played by the key physical mechanisms on the development of the two-phase flow for different compositions of the mixture. The main focuses are the distribution of the particles in the system and the pressure recovery


2020 ◽  
Vol 142 (6) ◽  
Author(s):  
Jamel Chahed ◽  
Lucien Masbernat

Abstract A two-fluid model with second-order turbulence closure is used for the simulation of a turbulent bubbly boundary layer. The turbulence model is based on the decomposition of the Reynolds stress tensor in the liquid phase into two parts: a turbulent part and a pseudo-turbulent part. The reduction in second-order turbulence closure in the near-wall region is interpreted according to a modified wall logarithmic law. Numerical simulations of bubbly boundary layer developing on a vertical flat plate were performed in order to analyze the bubbles effect on the liquid turbulence structure and to evaluate the respective roles of turbulence and of interfacial forces in the near-wall distribution of the void fraction. The two-fluid model with the second-order turbulence closure succeeds in reproducing the diminution of the turbulent intensity observed in the near-wall region of bubbly boundary layer and the increase in turbulence outside the boundary layer. The analysis of the interfacial force in the near-wall zone has led to the development of relatively simple formulation of the lift-wall force in the logarithmic zone that depends on dimensionless distances to the wall. After appropriate adjustment, this formulation makes it possible to reproduce the shape of the near-wall void fraction peaking observed in bubbly boundary layer experiments.


Author(s):  
Raphael V. N. de Freitas ◽  
Carina N. Sondermann ◽  
Rodrigo A. C. Patricio ◽  
Aline B. Figueiredo ◽  
Gustavo C. R. Bodstein ◽  
...  

Numerical simulation is a very useful tool for the prediction of physical quantities in two-phase flows. One important application is the study of oil-gas flows in pipelines, which is necessary for the proper selection of the equipment connected to the line during the pipeline design stage and also during the pipeline operation stage. The understanding of the phenomena present in this type of flow is more crucial under the occurrence of undesired effects in the duct, such as hydrate formation, fluid leakage, PIG passage, and valve shutdown. An efficient manner to model two-phase flows in long pipelines regarding a compromise between numerical accuracy and cost is the use of a one-dimensional two-fluid model, discretized with an appropriate numerical method. A two-fluid model consists of a system of non-linear partial differential equations that represent the mass, momentum and energy conservation principles, written for each phase. Depending on the two-fluid model employed, the system of equations may lose hyperbolicity and render the initial-boundary-value problem illposed. This paper uses an unconditionally hyperbolic two-fluid model for solving two-phase flows in pipelines in order to guarantee that the solution presents physical consistency. The mathematical model here referred to as the 5E2P (five equations and two pressures) comprises two equations of continuity and two momentum conservation equations, one for each phase, and one equation for the transport of the volume fraction. A priori this model considers two distinct pressures, one for each phase, and correlates them through a pressure relaxation procedure. This paper presents simulation cases for stratified two-phase flows in horizontal pipelines solved with the 5E2P coupled with the flux corrected transport method. The objective is to evaluate the numerical model capacity to adequately describe the velocities, pressures and volume fraction distributions along the duct.


2021 ◽  
Vol 33 (3) ◽  
pp. 033324
Author(s):  
Alejandro Clausse ◽  
Martín López de Bertodano

2013 ◽  
Author(s):  
Sung Chan Cho ◽  
Yun Wang

In this paper, two-phase flow dynamics in a micro channel with various wall conditions are both experimentally and theoretically investigated. Annulus, wavy and slug flow patterns are observed and location of liquid phase on different wall condition is visualized. The impact of flow structure on two-phase pressure drop is explained. Two-phase pressure drop is compared to a two-fluid model with relative permeability correlation. Optimization of correlation is conducted for each experimental case and theoretical solution for the flows in a circular channel is developed for annulus flow pattern showing a good match with experimental data in homogeneous channel case.


Author(s):  
David Heinze ◽  
Thomas Schulenberg ◽  
Lars Behnke

A simulation model for the direct contact condensation of steam in subcooled water is presented that allows determination of major parameters of the process, such as the jet penetration length. Entrainment of water by the steam jet is modeled based on the Kelvin–Helmholtz and Rayleigh–Taylor instability theories. Primary atomization due to acceleration of interfacial waves and secondary atomization due to aerodynamic forces account for the initial size of entrained droplets. The resulting steam-water two-phase flow is simulated based on a one-dimensional two-fluid model. An interfacial area transport equation is used to track changes of the interfacial area density due to droplet entrainment and steam condensation. Interfacial heat and mass transfer rates during condensation are calculated using the two-resistance model. The resulting two-phase flow equations constitute a system of ordinary differential equations, which is solved by means of the explicit Runge–Kutta–Fehlberg algorithm. The simulation results are in good qualitative agreement with published experimental data over a wide range of pool temperatures and mass flow rates.


Author(s):  
Hiroyuki Yoshida ◽  
Takeharu Misawa ◽  
Kazuyuki Takase

Two-fluid model can simulate two phase flow less computational cost than inter-face tracking method and particle interaction method. Therefore, two-fluid model is useful for thermal hydraulic analysis in large-scale domain such as a rod bundle. Japan Atomic Energy Agency (JAEA) develops three dimensional two-fluid model analysis code ACE-3D, which adopts boundary fitted coordinate system in order to simulate complex shape channel flow. In this paper, boiling two-phase flow analysis in a tight lattice rod bundle is performed by ACE-3D code. The parallel computation using 126CPUs is applied to this analysis. In the results, the void fraction, which distributes in outermost region of rod bundle, is lower than that in center region of rod bundle. At height z = 0.5 m, void fraction in the gap region is higher in comparison with that in center region of the subchannel. However, at height of z = 1.1m, higher void fraction distribution exists in center region of the subchannel in comparison with the gap region. The tendency of void fraction to concentrate in the gap region at vicinity of boiling starting point, and to move into subchannel as water goes through rod bundle, is qualitatively agreement with the measurement results by neutron radiography. To evaluate effects of two-phase flow model used in ACE-3D code, numerical simulation of boiling two-phase in tight lattice rod bundle with no lift force model (neglecting lift force acting on bubbles) is also performed. From the comparison of numerical results, it is concluded that the effects of lift force model are not so large on overall void fraction distribution in tight lattice rod bundle. However, higher void fraction distribution in center region of the subchannel was not observed in this simulation. It is concluded that the lift force model is important for local void fraction distribution in rod bundles.


Sign in / Sign up

Export Citation Format

Share Document