Assessment of Components of the Subcooled Boiling Model for CFD Simulations

Author(s):  
Deoras Prabhudharwadkar ◽  
Martin A. Lopez de Bertodano ◽  
John Buchanan ◽  
Avinash Vaidheeswaran

This paper describes the details of validation of heat and mass transfer models used for subcooled boiling simulation with a CFD two-fluid model. This research was focused on assessment of the wall heat flux partitioning model using the state-of-the-art multidimensional experimental data available in the literature. Various constitutive relations used to close the vapor generation rate at the heated wall were studied and the best suited combination of these was obtained. The current study was restricted to vertical flows through pipe and annulus geometries. Three data sets from the literature were considered: first with R12 at about 26 bar pressure, second with water at atmospheric pressure and third with R113 at 2.69 bar pressure. In these data sets, the bubble diameter distribution across the ducts was measured. Bubble diameter estimation brings in the largest uncertainty in the two-fluid model predictions and hence using the data with known bubble sizes allowed to focus on assessment of other parameters which use constitutive relations to model vapor generation rate, e.g. bubble nucleation site density and bubble departure frequency at the wall. The simulations were carried out using the CFD code CFX-12. The R12 data used here corresponds to fluid-vapor density ratio which is equivalent to that of water-steam at 150 bar. Therefore the density ratio varies over two orders of magnitude. The surface tension also varies over a wide range from 0.0017 to 0.057 N/m. The ratio of the flow channel hydraulic diameter to the bubble diameter in these simulations varied between 4 and 40. The two-fluid model was modified, for cases involving bubble sizes too large to be represented using the continuum assumption, in order to obtain satisfactory results. Hence, the model that has been developed for this study is applicable for a wide range of physical conditions and bubble sizes.

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):  
Xiaodong Sun ◽  
Yang Liu ◽  
Basar Ozar ◽  
Mamoru Ishii ◽  
Joseph M. Kelly

To apply the two-fluid model to a wide range of flow regimes in gas-liquid two-phase flows, the gas phase is categorized into two groups: small spherical/distorted bubbles as Group 1 and large cap/slug/churn-turbulent bubbles as Group 2 in the modeling of interfacial area transport. The interfacial transfer terms of momentum and energy for the gas phase are then divided into two groups accordingly in the implementation of the two-group interfacial area transport equation to the two-fluid model. Thus, the drag coefficients and the interfacial heat transfer for each group bubbles need to be developed. An approach has been sought for evaluating the drag coefficients of each bubble group based on a comprehensive experimental data base obtained in air-water upward flows in various size round pipes. Comparisons have been made with the theory of the drag coefficients and it was found that the agreement is not very satisfactory although the general trends can be predicted by the current approach.


Energies ◽  
2021 ◽  
Vol 14 (18) ◽  
pp. 5776
Author(s):  
Yongzhong Zeng ◽  
Weilin Xu

Bubble flow can be simulated by the two-fluid model and the multi-fluid model based on the Eulerian method. In this paper, the gas phase was further divided into several groups of dispersed phases according to the diameter by using the Eulerian-Eulerian (E-E) multi-fluid model. The diameters of bubbles in each group were considered to be the same, and their distributions were reorganized according to a specific probability density function. The experimental data of two kinds of bubble flow with different characteristics were used to verify the model. With the help of the open-source CFD software, OpenFOAM-7.x (OpenFOAM-7.0, produced by OpenFOAM foundation, Reading, England), the influences of the group number, the probability distribution function, and the parameters of different bubble diameters on the calculation results were studied. Meanwhile, the numerical simulation results were compared with the two-fluid model and the experimental data. The results show that for the bubble flow with the unimodal distribution, both the multi-fluid model and the two-fluid model can obtain the distribution of gas volume fraction along the pipe radius. The calculation results of the multi-fluid model agree with the experimental data, while those of the two-fluid model differ greatly from the experimental data, which verifies the advantage of the multi-fluid model in calculating the distribution of gas volume fraction in the polydisperse bubble flow. Meanwhile, the multi-fluid model can be used to accurately predict the distribution of the parameters of each phase of the bubble flow if the reasonable bubble diameter distribution is provided and the appropriate interphase force calculation model is determined.


2020 ◽  
Vol 501 (2) ◽  
pp. 1940-1950
Author(s):  
A Alharbi ◽  
I Ballai ◽  
V Fedun ◽  
G Verth

ABSTRACT The plasma dynamics at frequencies comparable with collisional frequency between various species has to be described in multifluid framework, where collisional interaction between particles is an important ingredient. In our study, we will assume that charged particles are strongly coupled, meaning that they form a single fluid that interacts with neutrals, therefore we will employ a two-fluid model. Here, we aim to investigate the evolutionary equation of slow sausage waves propagating in a gravitationally stratified flux tube in the two-fluid solar atmosphere in a strongly ionized limit using an initial value analysis. Due to the collisional interaction between massive particles (ions and neutrals), the governing equations are coupled. Solutions are sought in the strongly ionized limit and the density ratio between neutrals and charged particles is a small parameter. This limit is relevant to the upper part of the chromosphere. Our results show that slow sausage waves associated with charged particles propagate such that their possible frequency is affected by a cut-off due to the gravitational stratification. In contrast, for neutral acoustic waves the cut-off value applies on their wavelength and only small wavelength waves are able to propagate. Slow modes associated with neutrals are driven by the collisional coupling with ions.


2003 ◽  
Vol 125 (4) ◽  
pp. 660-669 ◽  
Author(s):  
Venkatraman Iyer ◽  
John Abraham

A two fluid Eulerian-liquid Eulerian-gas (ELEG) model for diesel sprays is developed. It is employed to carry out computations for diesel sprays under a wide range of ambient and injection conditions. Computed and measured results are compared to assess the accuracy of the model in the far field, i.e., at axial distances greater than 300 orifice diameters, and in the near field, i.e., at axial distances less than 100 orifice diameters. In the far field, the comparisons are of drop mean velocities and drop fluctuation velocities and in the near field they are of entrainment velocities and entrainment constants. Adequate agreement is obtained quantitatively, within 30 percent, and qualitatively as parameters are changed. Unlike in traditional Lagrangian-drop Eulerian-fluid (LDEF) approaches that are employed for diesel spray computations, adequate resolution can be employed in the near field to achieve numerical grid independence when the two-fluid model is employed. A major source of uncertainty in the near field is in the modeling of liquid jet breakup and atomization.


Author(s):  
Ik Kyu Park ◽  
Seung Jun Lee ◽  
Soo Hyoung Kim ◽  
Hyungrae Kim ◽  
Jae Ryong Lee ◽  
...  

This paper presents a set of numerical procedure to innovate CFD code into a PWR component analysis code. A porous media approach is adapted to two-fluid model and conductor model, and a pack of constitutive relations close the numerical model into a PWR component analysis code. The separate verification calculations on conductor model and porous media approach, and the validation calculation for the integrated component-scale code are introduced.


Sign in / Sign up

Export Citation Format

Share Document