scholarly journals Computational investigation of liquid-solid slurry flow through an expansion in a rectangular duct

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

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.


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].


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.


Author(s):  
Wei Yao ◽  
Christophe Morel

In this paper, a multidimensional two-fluid model with additional turbulence k–ε equations is used to predict the two-phase parameters distribution in freon R12 boiling flow. The 3D module of the CATHARE code is used for numerical calculation. The DEBORA experiment has been chosen to evaluate our models. The radial profiles of the outlet parameters were measured by means of an optical probe. The comparison of the radial profiles of void fraction, liquid temperature, gas velocity and volumetric interfacial area at the end of the heated section shows that the multidimensional two-fluid model with proper constitutive relations can yield reasonably predicted results in boiling conditions. Sensitivity tests show that the turbulent dispersion force, which involves the void fraction gradient, plays an important role in determining the void fraction distribution; and the turbulence eddy viscosity is a significant factor to influence the liquid temperature distribution.


Author(s):  
Youn-Gyu Jung ◽  
Moon-Sun Chung ◽  
Sung-Jae Yi

This study discusses on the implementation of an upwind method for a one-dimensional two-fluid model including the surface tension effect in the momentum equations. This model consists of a complete set of six equations including two-mass, two-momentum, and two-internal energy conservation equations having all real eigenvalues. Based on this equation system with upwind numerical method, the present authors first make a pilot code and then solve some benchmark problems to verify whether this model and numerical method is able to properly solve some fundamental one-dimensional two-phase flow problems or not.


Sign in / Sign up

Export Citation Format

Share Document