ANALYSIS AND APPLICATIONS OF A TWO-FLUID MULTI-FIELD HYDRODYNAMIC MODEL FOR CHURN-TURBULENT FLOWS
|
|
- Philip Stewart
- 6 years ago
- Views:
Transcription
1 Proceedings of the st International Conference on Nuclear Engineering ICONE21 July 29 - August 2, 2013, Chengdu, China ICONE ANALYSIS AND APPLICATIONS OF A TWO-FLUID MULTI-FIELD HYDRODYNAMIC MODEL FOR CHURN-TURBULENT FLOWS Gustavo Montoya Helmholtz-Zentrum Dresden-Rossendorf P.O. Box Dresden Germany g.montoya-zabala@hzdr.de Yixiang Liao Helmholtz-Zentrum Dresden-Rossendorf P.O. Box Dresden Germany Dirk Lucas Helmholtz-Zentrum Dresden-Rossendorf P.O. Box Dresden Germany Eckhard Krepper Helmholtz-Zentrum Dresden-Rossendorf P.O. Box Dresden Germany Keywords: CFD modeling, CFD validation, churn-turbulent flow, multiphase flow, two-fluid model ABSTRACT Today Computational Fluid Dynamic (CFD) codes are widely used for industrial applications, mostly in the case of single phase flows in automotive or aircraft engineering, but multiphase flow modeling had gain an increasing importance in the last years. Safety analyses on nuclear power plants require reliable prediction on steam-water flows in case of different accident scenarios. This is particularly true for passive safety systems such as the GEKO component of the KERENA reactor. Here flashing may occur in the riser (Leyer and Wich, 2012). In such case, high gas volume fractions and the churn-turbulent flow regime may ensue. In the past, the codes for the prediction of churn-regime have not shown a very promising behavior. In this paper, a two-fluid multi-field hydrodynamic model has been developed based in the Euler-Euler framework. The main emphasis of this work has been on the modeling and applicability of various interfacial forces between dispersed gaseous phases and the continuous liquid, as well as bubblebubble interactions, and the evolution of different bubble sizes in an adiabatic vertical pipe inside the churn-turbulent flow regime. All the expected mechanistic models that intervene in this flow pattern have been taken into account including drag force, wall force, lift force, turbulent dispersion, and bubble induced turbulence. Bubble breakup and coalescence has been defined (Liao et al., 2011), and in order to design a polydispersed model related to reality, the inhomogeneous MUSIG approach (Krepper et al., 2008) has been used to defined an adequate number of bubble size fractions which are arranged into different groups with their own velocity field. Based on these models, a series of simulations were made on the framework of ANSYS CFX 14.0, and all of the calculations were further validated with experimental data extracted from the TOPFLOW facility at the Helmholtz-Zentrum Dresden- Rossendorf. Different water and gas flow rates were used inside the churn-turbulent flow regime, as well as for the transition from bubbly to churn flow. The calculated cross-section averaged bubble size distributions, gas velocities, and time averaged radial profile for the gas fraction have shown a promising agreement with the experimental data. Nevertheless there are also clear deviations which indicate shortcomings of the present modelling. In order to further improve the modeling of this flow regime, a discussion based on the results will be used to shown a series of limitations of the actual modeling and possible solutions to be implemented in future works. 1. INTRODUCTION The applicability of CFD codes for two-phase flows has always been limited to special cases due to the very complex structure of its interface. Furthermore, the additional increase of computational power will not allow reflecting the details of the interface of gas-liquid with practical relevance. Instead, averaging procedures are commonly used in these cases, which at the same time cause the loss of information on the interfacial structure. In order to apply the wide-spread two-fluid model, closure laws are needed to reflect the interfacial mass, momentum and heat transfer (Lucas et al., 2007). Such closure models strongly depend on the flow pattern. When considering vertical pipe flow, it is known that for low gas volume flow rates, bubbly flow occurs (Qi et al., 2012). With increasing gas volume flow rates more and more large 1 Copyright 2013 by ASME
2 bubbles are generated by bubble coalescence, which leads to a transition to slug or/and churn-turbulent flow. Further increase will finally results in an annular flow pattern. Considering, for example, a heated tube producing steam by evaporation, such as steam generator tubes, all these flow patterns including transitions between them are expected to occur in the system. Up to now, no CFD simulation of such case is possible. This becomes clear, if the complexity of the flow situations is considered. Even in the case of bubble flow, the momentum exchange between liquid and the bubbles strongly depends on the bubble size. An example of this would be the sign change of the lateral lift force action on a bubble at a critical size, which leads to a de-mixing of small and large bubbles. The purpose of this paper is the development, test, and experimental validation of a multifield model of adiabatic gasliquid flows at high gas volume fractions (churn-turbulent flow regime) for which a multiple-size bubble approach has been implemented by dividing each of the gas structures into an specified number of groups, each of which represents a prescribed range of sizes. The simulations were performed using the computational fluid dynamic code from ANSYS, CFX A complete threefield and four-field model, including the continuous liquid field and two to three gas fields representing bubbles of different sizes, were first tested for numerical convergence and then validated against TOPFLOW experimental data (Lucas et al., 2010; Beyer et al., 2008). The simulations were made with the objective of predicting the evolution of bubble concentration from the inlet of the channel to fully-developed conditions eight meters downstream along an adiabatic vertical pipe. 2. NOMENCLATURE Notation C D Drag coefficient. C L Lift coefficient. C cc Interphase drag force coefficient for the FAD turbulence dispersion model. C TT Turbulent dispersion coefficient. C VV Virtual mass coefficient. C W Wall force coefficient. C μ Coefficient for viscocity. d b Bubble diameter. D h,k Maximum horizontal dimension of the bubble. EE k Modified Eötvos number. g Acceleration due to gravity. m k Volumetric mass transfer term into field-k from other fields representing the same phase. i M kk Interfacial momentum transfer per unit time between the fields k and j. p Absolute pressure. ε S c Source term for the turbulent eddy dissipation. r d Drag radius [m]. RR k Reynolds number of the k-phase. S Source term. t Time [s]. v k Phase-k averaged velocity. v m Mass averaged mixture velocity. w rrr Relative velocity. Greek letters α Void fraction. ε Turbulent eddy dissipation term. κ Turbulent kinetic energy. μ Viscocity. ρ k Density of the phase-k. σ tc Turbulent Schmidt number for the continuous phase volume fraction. τ k Total shear stress term. μ τ k Total shear stress term due to viscocity. RR τ k Total shear stress term due to the Reynolds number. ω Turbulent frequency term. Superscripts i Interfacial. D Drag. NN Non-drag. VV Virtual mass. TT Turbulent dispersion. L Lift. W Wall. Subscripts c Continuous phase. d Dispersed phase. j Index. k k-phase. 3. COMPUTATIONAL FLUID DYNAMICS (CFD) MODEL FORMULATION FOR CHURN-TURBULENT FLOWS A three to four-field multiphase CFD modeling was created to simulate two-phase churn-turbulent flow regime in a vertical pipe. Such a multifield modeling concept of multiphase flows is based on ensemble averaging the respective governing equations for each component fluid. When using such modeling along adequate closure laws, this approach should be capable of representing flow regimes from bubbly through churn-turbulent to annular flow. With the present available CFD-models, it is not possible to describe transitions between flow regimes. In order to model the adiabatic flow, the resultant conservation equations shown in (1) and (2) are commonly used, representing the mass and momentum conservation equations respectively. (α k ρ k ) d + (ρ k α k v k ) = m k (1) 2 Copyright 2013 by ASME
3 (α k ρ k v k ) + (ρ t k α k v k v k ) = α k p k + α k τ k + ρ k α k g i + M kk j i + m m,k v m m (2) The drag force is known for playing a major role in the axial direction (Drew, 1992). The drag contribution model use in this research work is given as M D d,k ρ c = 3 α 4 d,k C d D,k v c v d,k v c v d,k (5) b Where m k = m m m,k In the aforementioned equations, α k represents the volume fraction of the field-k, while m m,k is the volumetric mass transfer term from other fields which represents the same phase into field-k. The term τ k = τ μ k + τ Rr k = μ k ( v k + v T k ) 2 μ 3 k v k I is the total shear stress term. i The variable M kk represents the interfacial momentum transfer per unit time between two different fields. 3.1 Description of the Mechanistic Models for Interfacial Forces One of the most important factors for the correct representation of the hydrodynamics on multiphase flows is the adequate designation of mechanistic closure laws which allow the modeling of the principle mass and momentum interactions at the fluid interfaces. Due to the averaging of the conservation equations, some information is lost, but has to be reintroduced by the use of these closure relations. The closure laws objective is to account for the mass and momentum transfer between the different fields and phases while providing the functional form expected for the interfacial forces. Still, the presented models are limited by the need of local condition dependent coefficients derived from the fact that the closure laws have been developed for ideal bubbly flow and are now being applied to churn-turbulent flow conditions. The system is assumed to be under adiabatic conditions which translate into the consideration of momentum exchange only at the bubble interface. The different forces that interact at the interface of the flow can be divided into drag and non-drag components with their respective mechanistic models. M k i = M k D + M k NN (3) The interfacial non-drag forces include the lift force, the wall force, the turbulent dispersion force, and the virtual mass force. For reasons that will be explained later in this chapter, the virtual mass force was not taken into account for the present model. Other contributions such as the Basset force may be presented but neglected due to its small effect (Guillen et al, 2009). The total force acting at the interface of the flow can be expressed as the superposition of several component forces. M k i = M k D + M k VV + M k TT + M k L + M k W (4) For the drag coefficient calculation, the expression proposed by Ishii and Zuber (1979) has been used. The next equations represent the correlations for the drag coefficient for spherical cap, ellipse, and distorted bubbles respectively. C D,k (ccc) = 8 3 (1 α k) 2 C D,k (eeeeeee) = 2 3 EE1 2 (6) C D,k (ddddddddd) = min (C D,k (eeeeeee), C D,k (ccc) ) It is important to note that other researchers have reported some significant overestimation of this drag coefficient so that one has to artificially increase the bubble size for the simulation to match the observed flow data when applied to churnturbulent flow conditions (Chen et al., 1995; Tsuchiya et al., 1997). The non-drag interfacial forces are dominant and control the gas profile in the radial direction. The turbulent volume fraction dispersion force is included with the objective of take into account the increased in mixing due to turbulence in the multiphase flow. The turbulent dispersion force transport the gaseous phase in a direction opposing the void fraction gradient. The model that was used in this research is the one known as Favre average drag model (FAD) proposed by Burns (2004). M c TT = M c TT = C TT C cc v tc α d σ tc α d α c α c (7) In the prior equation, the C cc represents the interface drag force coefficient. Therefore, this model depends on details of the drag force correlation that is being used. The term σ tt refers to the turbulent Schmidt number for the continuous phase volume fraction, which is commonly taken to be 0.9 (CFX- Theory User Manual, ANSYS V14.0). The coefficient for turbulent dispersion (C TT ) is commonly take as. In order to account for the momentum exchange due to lift between the continuous and dispersed gaseous phases, the interfacial lift force closure law is used. L M d,k = C L ρ c α d,k (v d,k v c ) ( v c ) (8) The lift coefficient (C L ) is calculated as proposed by Tomiyama et al. (2002). 3 Copyright 2013 by ASME
4 C L,k = min ttth(0.121rr k ), f(ee k ) fff EE k < 4 C L,k = f(ee k ) fff EE k < 4 (9) C L,k = 0.27 fff EE k > 10 In this case, f(ee k ) is defined as shown in equation (10). The critical bubble diameter (Lucas et al., 2006) at which the lift force change sign can be found by solving the roots of this correlation which gives a value of 5.8 mm under atmospheric pressure and room temperature for air bubbles in water (Tomiyama, 1998). f(ee k ) = 0105EE k 3 159EE k 2 204EE k (10) The modified Eötvös number or Eo k is calculated as EE k = g(ρ c ρ d )D h,k 2 σ (11) From an empirical correlation by Wellek et al. (1966) for aspect ratio of an elliptical shape bubble, the maximum horizontal dimension of the bubble, D h,k can be calculated as D h,k = D k EE k (12) In order to account for the hydrodynamic force produce on a bubble which moves in close proximity to a wall, the interfacial wall force is used. The model chose for this force is the one proposed by Tomiyama (1998). M W d d,k = C b α 1 W 2 y 2 ρ cw 2 rrr n r (13) The wall coefficient (C W ) used in this investigation was the one proposed by Hosokawa et al. (2002), who studied the influence of the Morton number, and proposed a correlation which is valid for both high and low Morton numbers provided that the bubble does not collide with the wall (Rzehak et al., 2012). 7 C W = max RR 1.9, 217EE k (14) k The effect of the wall is generally observed as a sharp spike in the bubble void fraction near the wall, that affects the small bubbles and in a weaker manner also the large bubbles. Another of the interfacial non-drag forces named in equation (4) is the virtual mass or added mass force, which is in charge of accelerating the flow around a bubble as it moves through the liquid. As it has been said prior in this section of the present work, this force has been neglected. The reason for doing this is due to that in churn-turbulent flow regime, one of the distinguished characteristics of the flow is the existence of large-fast rising bubbles. These large bubbles do not have a close wake and so the concept of added mass is not applicable. On the other hand, the small bubbles do have a close wake, but in this regime they suffer strong recirculation, moving downwards near the wall region. For this reason, the inclusion of the added mass contributions to the small bubbles would lead to severe convergence difficulties (Krishna et al., 2001). Furthermore, it has been reported (Guillen et al., 2009) that for churn-turbulent regime this force can be neglected since its effect is very small, and does not affect the overall gas profile. For the simulations presented in this work, it has been assumed that between seventeen to nineteen groups of bubbles of different diameters enter the flow channel. These bubbles will eventually start to coalesce forming a large spectrum of new and larger than the original bubbles. The opposite process will also occur making the large bubbles breakup into smaller ones. A mechanistic and accurate mass transfer model, which should be also computationally effective, is needed in order to be used for modeling the interactions between the different sizes of bubbles (Kumbaro and Podowski, 2006; Monahan et al., 2005). Liao et al. (2011) proposed an approach for bubble breakup and coalescence for the inhomogeneous multiple bubble size group (MUSIG) model. Their main objective was to consider a major number of important mechanisms presented in bubble breakup and coalescence in turbulent gas-liquid mixtures. The interactions that were considered for the coalescence mechanism are coalescence due to turbulence, laminar shear, wake-entrainment, and eddy capture. For bubble breakup they consider the effect due to turbulent fluctuation, laminar shear, and interfacial slip velocity. This approach was implemented in the framework of ANSYS CFX via user fortran routines, and applied to the case of turbulent air-water mixtures in a large vertical pipe in order to be validated for bubbly flow conditions in the aforementioned paper. The implementation used in the actual research work is similar to the one presented there with only slight changes in the calculation of the wake entrainment and turbulence contributions for the coalescence model. The turbulence approach used for the liquid phase is the k- ω based SST model, which accounts for the transport of the turbulent shear stress and produces highly accurate predictions of the onset and the amount of flow separation under adverse pressure gradients. This model behaves better than the baseline (BSL) approach, which even by combining the advantages of the Wilcox and the k-ε models, fails to properly predict the onset and the amount of flow separation from smooth surfaces (CFX-Theory User Manual, ANSYS V.14.0). In general, the k- ω model counts with various advantages over the k-ε approach, such as the treatment when in low Reynolds number close to the wall. This model does not need the complex non-linear damping functions required for the k-ε, and for that matter can be considered as more accurate and robust. Finally, the bubble induced turbulence effect was taken into account by the implementation of the Rzehak and Krepper (2012) approach, which is based in a two equation turbulence model with additional source term description for the particle 4 Copyright 2013 by ASME
5 induced turbulence. The source term describing the bubble effects for the k-equation is shown in the next correlation, and it is formulated as an assumption that all energy lost by the bubble due to drag is converted to turbulent kinetic energy. S c k = F c D. v c v c (15) The source term for the ω-equation is formulated as S c ω = 1 C μ k c S c ε ω c k c S c k (16) In the last equation, C μ has a standard value of 9. This model has been compared with others approaches, such as Sato enhanced turbulence, showing a better prediction on the actual model. 3.2 Description of the Inhomogeneous Multiple Size Group model (MUSIG) Krepper et al. (2005) proposed an association of different dispersed phases and the algebraic multiple size group model in order to combine both, the adequate number of bubble size classes for the simulation of coalescence and breakup and a limited number of dispersed gaseous phases, with the purpose of limiting the computational effort. The gaseous dispersed phase is divided into a number of N velocity groups (or phases), where each velocity group is characterized by its own velocity field. Then, the overall bubble size distribution is further represented by dividing the bubble diameter range within each of the velocity groups j into a number, M j (j = 1 N), bubble sub-size fractions (see Figure 1). The population balance model considers bubble coalescence and breakup that are applied to the sub-size groups, which means that the mass exchange between the sub-size groups can exceed the size ranges assigned to the velocity clusters (Krepper et al., 2009). Figure 1: Improvement of the polydispersed approach: the size fractions M j are assigned to the velocity field V j (Krepper et al., 2008). In order to control the lower and upper boundaries of the intervals for the bubble size fractions, equal distribution of bubble diameter, equal bubble mass, or user definition of the bubble diameter ranges for each distinct bubble diameter fraction could be used. The subdivision should be based on the physics of bubble motion for gas structures of different sizes, for example, dissimilar behavior of distinctively sized bubbles with respect to lift force or turbulent dispersion. Extensive validation made by Krepper et al. (2007 & 2008), have shown that in most cases N = 2 or 3 velocity groups should be sufficient in order to capture the main phenomena, at least in bubbly regimes. 4. DESCRIPTION OF THE EXPERIMENTAL FACILITY The simulations presented in this research work have been validated using an experimental database for air/water twophase flows in vertical stainless steel pipes with nominal diameter of 50 mm and 200 mm, that has been extracted from the thermal-dynamic test facility TOPFLOW at the Helmholtz- Zentrum Dresden-Rossendorf (see figure 2). This database provides information about the evolution of two-phase flow along the pipe height. In order to achieve this, several measurements with varying distances between the gas injection and the measurement plane were made for different combination of gas and water flow rates with superficial gas velocities (j G ) varying from 025 to m/s, and liquid superficial velocities (j L ) from 405 to m/s. Desalinated water was used in the experiments. The pressure at the gas injection is kept constant at 0.25 MPa. The advantage of such boundary condition lies in that the measure data represents the evolution along the pipe as if the gas injection was in a fixed position and the measurement plane varies. The temperature in the system was also kept constant at 30 C. The test section DN200 is equipped with six gas inlet locations which allowed the injection of air or steam via orifices in the pipe wall. To have such injection via wall orifices allow the two-phase flow to rise smoothly to the measurement plane without having any influence from the feeder within the tube at any other location along the flow (Beyer et al., 2008). The liquid phase was supplied from the bottom of the test section. The measurement plane is kept constant at the top of the system, and consists in two wire-mesh sensors. The lower sensor is used for obtaining data on the gas volume fraction profiles and bubble size distributions, while the second one is to determine gas velocities by cross-correlation of both sensors. The results of these experiments are the time averaged radial profiles of the gas fraction, the gas velocity, the time and cross-section averaged bubble sized distributions, and the gas fraction data regarding the resolved bubble sized and spatial distribution. Furthermore, there is also an extensive analysis of the encountered flow patterns and a study of the effect of the orifice s diameter in the gas injection (Beyer et al., 2008). In the present studies, the 1 mm orifices for the injection of the gas were used. As expected from the literature (Ohnuki et al., 2000; Prasser et al., 2007), no slug pattern in the DN200 pipe was observed. It was reported a criterion for the analyzed of the observed flow pattern in the experiments based on the 5 Copyright 2013 by ASME
6 maximum observable diameter, and the visualizing of the void fraction. For bubbly flow this diameter must be less than 50 mm and present a mono-modal bubble sized distribution, while in the churn-turbulent flow there should be a bimodal bubble sized distribution with bubbles of diameter similar to those encountered in the bubbly flow and also bigger gas structures of diameter less than 400 mm. Figure 2: Vertical test section of the TOPFLOW facility with variable gas injection system. 5. COMPUTATIONAL GRID AND BOUNDARY CONDITIONS gas phase have been taken from the experimental data given by TOPFLOW at L/D = 1.1, including gas velocity, mass flow rate and size fractions. In principle, the TOPFLOW data regarding the gas phase was arrange into two velocity groups. The first velocity group contains six bubble size fractions up to 5.8 mm, while the second group is organized into twelve to fourteen bubble size fractions larger than the limit size of the first velocity group. The defined number of size fractions is dependent of the test case. This limitation in the size in order to establish the difference between one gas field and the other is given by the change on the sign of the Tomiyama lift force coefficient. A logarithmic distribution is used for the calculation of the different size groups. The smallest size bubble that can be measured by the wire mesh sensor of the TOPFLOW facility is about 2 mm. Some first calculations will be given when considering an extra field for the separation of very large cap bubbles. In this case, a third velocity group is defined based on the maximum cap bubble limit diameter correlation proposed by Ishii and Zuber (1979). This diameter was calculated to be around 10.8 mm, reducing the second gas field to mere three gas fractions, and leaving the third gas field with eleven size fractions. A more in-depth explanation about the reason for these considerations will be given in the results and discussion section of this paper. The simulations are based to replicate the experimental conditions established in TOPFLOW air-water flow in a vertical test section experiments. An axisymmetric 40 x 800 grid with refined cell size of mm near the wall, and mm near the inlet, were used to represent the 9765m x 8m system. There are 40 nodes along the radial direction and 800 axial nodes. Taking advantage of the symmetry, only a segment of the pipe was modeled, establishing a so-called 2D axisymmetric approximation with its respective symmetry planes (see figure 3). Three cases from the TOPFLOW experimental studies were taken for the validation of the presented model as it can be seen in table 1. Cases 127 and 129 are both under churnturbulent flow regime conditions, while case 118 is considered to be in transition from bubbly to churn-turbulent flow. Table 1: Initial conditions used in the simulations Run number Superficial liquid velocity (m/s) Superficial gas velocity (m/s) The data corresponding to the inlet conditions at the continuous phase in the multiphase modeling (turbulent kinetic energy, and turbulent dissipation) were taken from fully developed single phase simulations made on the same computational grid, while the velocity profile is calculated from the TOPFLOW experiments. The boundary conditions for the Figure 3: 3D pipe geometry, 2D axisymmetric approximation, and 2D axisymmetric mesh with unequal node distribution (segmented in the axial direction). At the wall boundary, a no slip condition was used. A symmetry boundary condition was established along the pipe 6 Copyright 2013 by ASME
7 centerline. The meshing software utilized was ANSYS ICEM CFD. The mesh unequal distribution, with a finer spacing at the inlet and wall of the pipe was established with the intention of capturing the developing flow at the inlet as well as the effect of high gradients of velocity and bubble concentration near the wall of the tube. 5. RESULTS AND DISCUSSION Extensive testing has been performed on the proposed modeling approach using the experimental conditions corresponding to various data series from TOPFLOW experiments with the objective of validating a number of models previously proposed for bubbly flows, and now implemented on churn-turbulent flow regime. The radial concentrations of the different defined gas fields, as well as the gas and liquid velocity magnitude profiles, can be seen in figure 4 to 9 for the complete channel. In order to plot the whole channel, an unequal ratio height to length has been used in the presented images. The left side in each of the domains represents the wall boundary. For all three cases, it can be seen that the small bubbles field or gas 1, remains near the wall along the channel, while the largest bubbles field or gas 2, are distributed in a way in which the maximum value of volumetric concentration is higher at the center of the pipe. It can also be observed in the images how the highest velocity values appear at the center of the pipe for both, liquid and gas phases. The experiment 127, which presents a lower superficial liquid velocity, showed a better distribution for the second gas field across the whole domain. Figure 5: Distributions of velocity for each dispersed and continuous field in experiment 118. Figure 6: Distributions of void fraction for each dispersed field in experiment 127. Figure 4: Distributions of void fraction for each dispersed field in experiment 118. Figure 7: Distributions of velocity for each dispersed and continuous field in experiment Copyright 2013 by ASME
8 Figure 8: Distributions of void fraction for each dispersed field in experiment 129. the TOPFLOW facility experiments are only accurate for sizes larger than 2 mm (Prasser et al., 2007). When the superficial gas velocity is increased in experiment 129 (see figure 12), the behavior of the bubble size distribution remains similar as in 118. There is only a slight increment of the probability of appearance of larger bubbles, which slightly decrease the overall appearance of the smaller ones. Still, the underprediction observed for the smaller bubbles distribution peak is highly reduced as to the observed one in the experiment 118. Once the superficial liquid velocity decreased in experiment 127 as in comparison with 129 (see figure 11), there is the appearance of even larger bubbles in the experimental data. In turn, an overprediction of smaller bubble sizes close to both the inlet boundary (L/D = 1.7) and the outlet (L/D = 39.9) is observable on these case, in contrast to the typical underprediction found in the prior two simulations. This last deviation results from the effect of reducing the liquid superficial velocity to less than half, which in turn also increase the turbulence parameters and for that matter increase the breakup on the whole domain. Bubble Size Distribution (% / mm) Figure 9: Distributions of velocity for each dispersed and continuous field in experiment 129. A quantitative analysis of experiment 118 through the bubble size distribution study (see figure 10) showed an overall good prediction of the experimental trend with a constant underprediction of bubbles with smaller sizes, which is followed by an overprediction of bubbles ranging from about 10 to 35 mm. This last deviation is more visible close to the inlet boundary (L/D = 1.7), and gradually decreases until the highest measurement level (L/D = 39.9). Bubbles which sizes are larger than 40 mm are constantly underpredicted for all measurement heights. The behavior observed for the values smaller than 35 mm can be explained in part due to too high coalescence at the beginning of the pipe, which in turn creates an underprediction of the smaller bubbles. This could be improved by diminishing the coefficient for wake entrainment in the coalescence model. It can also be observed that for the bubble size distribution, the smallest bubble size in the simulations differs from the experiments, but there is no possibility of real comparison since Bubble Size (mm) Figure 10: Validation of bubble size distribution against the experiment 118. Bubble Size Distribution (% / mm) Bubble Size (mm) Figure 11: Validation of bubble size distribution against the experiment Copyright 2013 by ASME
9 Bubble Size Distribution (% / mm) Bubble Size (mm) Figure 12: Validation of bubble size distribution against the experiment 129. In terms of the total radial void fraction, there is a very good agreement between the experiment 127 and the corresponding simulation (see figure 14). Once the liquid superficial velocity is increase in the experiment 129, a good fit with the simulations can be found for the lower levels (L/D = 1.7), but a peak appears between the wall and the center of the pipe for the next shown higher measurement plane (see figure 15). For posterior higher measurement positions these deviation disappear obtaining a similar trend as in experiment 127. While maintaining the same liquid superficial velocity, and decreasing the gas velocity, an analogous behavior of the total radial void fraction is encountered (see figure 13). Total Void Fraction ( ) Figure 13: Validation of the radial total void fraction against the experiment 118. In principle, it appears as if this deviation at L/D = 7.9 could be cause by a too large lift force in comparison with the turbulent dispersion. Furthermore, parametric sensitive studies made by Tselishcheva et al. (2010) have shown that both a negative lift coefficient and a too low drag force coefficient could produce similar behaviors in the total void fraction as the one observed in the highest level of the studied system (L/D = 39.9). As it was said prior to this section, the correlation proposed by Ishii and Zuber (1979) is known for overpredict the drag coefficient at churn-conditions. Total Void Fraction ( ) Figure 14: Validation of the radial total void fraction against the experiment 127. Total Void Fraction ( ) Figure 15: Validation of the radial total void fraction against the experiment 129. A study of the averaged gas velocities for the three experimental cases was also prepared. The observed behavior of the trend in the gas velocities for the three cases is analogous to the overall radial void fraction in terms of the expected deviations. For experiment 127, the velocities for all the evolution of the flow through the pipe qualitatively fit the experimental data very well (see figure 17). For both, experiment 118 (see figure 16) and 129 (see figure 18), the simulation fit almost perfectly the experiments for the height which is closer to the inlet boundary as expected, but for the measurement level L/D = 7.9, the velocity presents a similar peak as the one seem in figure 13 and 15, which attenuates eventually before achieving higher positions (L/D = 39.9). From the bubble size distribution charts, it is already possible to observed that many of the deviations appears to have their base in the misconception than the largest bubbles behave in the same way as it is expected from bubbly flow conditions. For that matter, an approach has been model in 9 Copyright 2013 by ASME
10 which a third dispersed velocity field has been used for further subdivision of the overall gas phase. Average Gas Velocity (m / s) Figure 16: Validation of the gas averaged velocity against the experiment 118. Average Gas Velocity (m / s) Figure 17: Validation of the gas averaged velocity against the experiment 127. In the first set of simulations, the new gas field shares the same closure models as the other gas phases, while in the other shown approach, the modeling has been made by neglecting the lift force in the newly added dispersed field. In this research paper, preliminary results are shown to validate this approach by comparison against the prior mentioned experiment 129. It can be observed that the addition of the third gas phase greatly improves the bubble size distribution in terms of the underprediction encountered with the two-gas field approach, especially for the higher positions of the measurement planes (see figure 19). On the other hand, there was not a visible effect on the simulation by the elimination of the lift force on the third gas phase. Bubble Size Distribution (% / mm) Bubble Size (mm) Figure 19: Validation of bubble size distribution against the experiment 129 using a three gas field model. In terms of the radial total void fraction, it can be observed that the peak encountered at L/D = 7.9 in figure 15, is no longer visible with the three-field simulation (see figure 20). This comes at the cost of the accuracy for the radial total void fraction at the lower measurement levels. The lack of lift force in this case, only shifted slightly the total void fraction to lower values for the higher measurement planes. Average Gas Velocity (m / s) Figure 18: Validation of the gas averaged velocity against the experiment 129. Total Void Fraction ( ) Figure 20: Validation of the radial total void fraction against the experiment 129 using a three gas field model. 10 Copyright 2013 by ASME
11 Finally, the addition of a third gas phase in terms of gas averaged velocity also in this case makes the peak observable in the two-field approach not visible anymore (see figure 21). A constant overprediction of averaged gas velocities for the values at the center of the pipe, and a constant underprediction close to the wall can be seen. The elimination of the lift force for the third gas phase did not produce an important effect in terms of the averaged gas velocity. Average Gas Velocity (m / s) Figure 21: Validation of the gas averaged velocity against the experiment 129 using a three gas field model. These deviations obtained for the higher measurement levels could be the result of a too low drag force coefficient. As prior stated, it has been shown that both the underprediction of the drag coefficient and a negative lift force could create such trends in the prediction of total void fraction and gas velocity. Since the effect of neglecting the lift force appear to have little effect in the third gas-field simulations (see figure 22), studies have to be made to consider the effect of other forces that could be producing the higher concentration of gas-3 at the center of the pipe, independently of the lift force, such as the wall lubrication force. Further studies should be made on the subject. Total Void Fraction ( ) Gas 1 Gas 2 Gas 3 Total Gas Gas 1 w/o Lift Force Gas 2 w/o Lift Force Gas 3 w/o Lift Force Total Gas w/o Lift Force Figure 22: Comparison of the total void fraction for the dispersed phases in the 3-field model with and without lift force at L/D = There are other factors which will be studied in future research works such as the effect of the turbulence modeling. In most of the simulations and modeling approaches defined in the last few years for churn-turbulent flow, including this work, the Reynolds-averaged Navier-Stokes (RANS) turbulence approach is the easiest and favorite choice in terms of description for the turbulence of the continuous phase. It is known (Bestion, 2011) that the RANS approach is not suitable for the overall prediction of the churn-regime due to its incapacity for notfiltering of the large scales produce by the largest bubbles (third gas field), not being able to predict the intermittency produce by the passage of these bubbles. Other factor affecting the third gas field is the effect produced by the confinement of the pipe in the simulation, which has not been taking into account before. A possible future modification in the drag force equation could account for the addition of this effect in the simulations. A different possible approach would be the complete resolution of the interface on the third gas field. A suitable technique for such modeling would be the Generalized Two Phase flow approach or GENTOP (Hänsch et al. 2012) currently being developed at the HZDR. While this is a very promising method, it would also bring many challenges and higher computational time simulations due to that it will require three dimensional transient calculations. 6. CONCLUSIONS The current presented work has shown the achievement of a series of important objectives in the studies of the churnturbulent flow as a significant phenomenon which occur in flashing and GEKO risers for safety research in the nuclear industry. It has been shown that while results using a bimodal gas distribution and bubbly flow closure laws produce promising and relatively good results, they cannot yet predict the expected physics at this flow regime. The modeling flow at churnturbulent flow conditions should be based on at least three gas field approach or trimodal distribution based on small spherical bubbles, larger spherical bubbles, and large deformable cap bubbles. In general, while the acknowledgement and consideration of the existence of a third gas field in comparison of the traditionally expected bimodal distribution in the churnturbulent flow regime (Krishna et al. 2001), predicted the physics behind the experiments in churn-conditions in a much more realistic way than before, there is still not achievement of an optimal modeling approach. This occurs because the expected behavior of the third gas phase should differ of the other two dispersed fields. It is expected that further studies and modifications on the specific closure and turbulence models in order to acknowledge the actual physics behind the churnregime, which are in some specific ways very different to those encountered under bubbly flow conditions, will improve in a greatly manner the current modeling at these flow conditions. 11 Copyright 2013 by ASME
12 Future works will include further studies on the closure laws and its effect on the third gas field, as well as tests on the turbulence modeling for these largest and deformable dispersed structures (churn-turbulent bubbles). The inclusion and study of the confinement effect of the pipe on the third gas field by implementation on the drag force equation will be also a next step on the future simulations. ACKNOWLEDGMENTS This work is carried out in the frame of a current research project funded by the energy company E.ON Kernkraft GmbH. REFERENCES 1. Bestion, D. Applicability of two-phase CFD to nuclear reactor thermalhydraulics and elaboration of Best Practice Guidelines. Nucl. Eng. Des, In Press (2011). 2. Beyer, M.; Lucas, D.; Kussin, J.; Schütz, P. Air-water experiments in a vertical DN200-pipe. FZD-505 TOPFLOW Wissenschaftlich-Technische Berichte (2008). 3. Burns, A.; Frank, T.; Hamill, I.; Shi, J. The Favre Averaged Drag Model for Turbulent Dispersion in Eulerian Multi-phase Flows. 5th International Conference on Multiphase Flow, ICMF04, Japan (2004). 4. Chen, P.; Dudukovic, M.P.; Sanyal, J. Three-Dimensional Simulation of Bubble Column Flows with Bubble Coalescence and Breakup. AIChE J., 51(3): (2005). 5. CFX-Theory User Manual, ANSYS V Drew, D.A. Analytical modeling of multiphase flows in Boiling Heat Transfer. Elsevier, New York (1992). 7. Guillen, D.; Shelley, J.; Antal, S.; Tselishcheva, E.; Podowski, M.; Lucas, D.; Beyer, M. Optimization of a Two-Fluid Hydrodynamic Model of Churn-Turbulent Flows. Proceedings of the 17th International Conference on nuclear Engineering (ICONE17), Paper ICONE (2009). 8. Hänsch, S.; Lucas, D.; Krepper, E.; Höhne, T. A multifield two-fluid concept for transitions between different scales of interfacial structures. Int. J. Mult. Flow 47: (2012). 9. Hosokawa, S.; Tomiyama, A.; Misaki, S.; Hamada, T. Lateral migration of single bubbles due to the presence of wall. Proceedings of ASME Fluids Engineering Division Summer Meeting, Canada (2002). 10. Ishii, M.; Zuber, N. Drag Coefficient and Relative Velocity in Bubbly, Droplet or Particle Flows. AIChE J., 25: (1979). 11. Krepper, E.; Lucas, D.; Prasser, H.-M. On the modelling of bubbly flow in vertical pipes. Nucl. Eng. Des. 235: (2005). 12. Krepper, E., Beyer, M.M Frank, Th.; Lucas, D.; Prasser, H.-M. Application of a population balance approach for polydispersed bubbly flows. 6th Int. Conf. on Multiphase Flow, Leipzig (2007). 13. Krepper, E.; Frank, Th.; Lucas, D.; Prasser, H.-M; Zwart, P.J. The Inhomogeneous MUSIG model for the simulation of polydispersed flow. Nucl. Eng. Des., 238: (2008). 14. Krepper, E.; Ruyer, P.; Beyer, M.; Lucas, D.; Prasser, H.- M.; Seiler, N. CFD simulation of polydispersed bubbly two-phase flow around an obstacle. STNI Article ID (2009). 15. Krishna, R.; van Baten, J.M. Eulerian simulations of bubble columns operating at elevated pressures in the churn turbulent flow regime. Chem. Eng. Sci., 56: (2001). 16. Kumbaro, A.; Podowski, M.Z. The effect of bubble/bubble interactions on local void distribution in two-phase flows. Proceedings 13th International Heat Transfer Conference, Australia (2006). 17. Leyer, S.; Wich, M. The Integral Test Facility Karlstein. Sci. Tech. Nucl. Inst., Article ID: (2012). 18. Liao, Y.; Lucas, D.; Krepper, E.; Schmidtke, M. Development of a generalized coalescence and breakup closure for the inhomogeneous MUSIG model. Nucl. Eng. Des., 24: (2011). 19. Lucas, D.; Krepper, E.; Prasser, H.-M; Manera, A. Investigations on the Stability of the Flow Characteristics in a Bubble Column. Chem. Eng. Technol., 29(9): (2006). 20. Lucas, D.; Krepper, E.; Prassar, H.-M. Use of Models for Lift, Wall and Turbulent Dispersion Forces Acting on Bubbles for Poly-Disperse Flows. Chem. Eng. Sci., 62: (2007). 21. Lucas, D.; Beyer, M.; Kussin, J.; Schütz. P. Benchmark database on the evolution of two-phase flows in a vertical pipe. Nucl. Eng. Des., 240: (2010). 22. Monahan, S.; Vitankar, V.; Fox, R. CFD Predictions for Flow-Regime Transitions in Bubble Columns. AIChE J., 51(7): (2005). 23. Ohnuki, A.; Akimoto, H. Experimental study on transition of flow pattern and phase distribution in upward air-water two-phase flow along a large vertical pipe. Int. J. Mult. Flow, 26: (2000). 24. Prasser, H.-M; Beyer, M.; Carl, H.; Manera, A.; Pietruske, H.; Schütz, P. Experiments on upwards gas/liquid flow in vertical pipes. FZD-482 TOPFLOW Wissenschaftlich- Technische Berichte (2007). 25. Qi, F.S.; Yeoh, G.H.; Cheung, S.C.P.; Tu, J.Y.; Krepper, E.; Lucas, D. Classification of bubbles in vertical gasliquid flow: Part 1 An analysis of experimental data. Int. J. Mult., 39: (2012). 26. Rzehak, R.; Krepper, E. Bubble induced turbulence: Comparison of CFD models. International Symposium on Multiphase flow and Transport Phenomena. Agadir, Morocco. April (2012). 27. Rzehak, R.; Krepper, E.; Lifante, C. Comparative study of wall-force models for the simulation of bubbly flows. Nucl. Eng. Des., 253: (2012). 12 Copyright 2013 by ASME
13 28. Tomiyama, A. Struggle with Computational Bubble Dynamics. Proceedings International Conference on Multiphase Flow, France (1998). 29. Tomiyama, A.; Tamai, H.; Zun, I.; Hosokawa, S. Transverse Migration of Single Bubbles in Simpler Shear Flows. Chem. Eng. Sci., 57: (2002). 30. Tselishcheva, E.; Podowski, M.; Antal, S.; Guillen, D.; Beyer, M; Lucas, D. Analysis of Developing Gas/Liquid Two-Phase Flows. 7th International Conference on Multiphase Flow (ICMF) (2010). 31. Tsuchiya, K.; Furumoto, A.; Fang, L.S.; Zhang, J. Suspension Viscocity and Bubble Rise Velocity in Liquid-Solid Fluidized Beds. Chem. Eng. Sci., 52: 3053 (1997). 32. Wellek, R.M.; Agrawal, A.K.; Skelland, A.H.P. Shapes of Liquid Drops Moving in Liquid Media. AIChE J., 12: (1966). 13 Copyright 2013 by ASME
DEVELOPMENT OF A MULTIPLE VELOCITY MULTIPLE SIZE GROUP MODEL FOR POLY-DISPERSED MULTIPHASE FLOWS
DEVELOPMENT OF A MULTIPLE VELOCITY MULTIPLE SIZE GROUP MODEL FOR POLY-DISPERSED MULTIPHASE FLOWS Jun-Mei Shi, Phil Zwart 1, Thomas Frank 2, Ulrich Rohde, and Horst-Michael Prasser 1. Introduction Poly-dispersed
More informationInvestigation of Three-Dimensional Upward and Downward Directed Gas-Liquid Two-Phase Bubbly Flows in a 180 o -Bent Tube
Investigation of Three-Dimensional Upward and Downward Directed Gas-Liquid Two-Phase Bubbly Flows in a 180 o -Bent Tube Th. Frank, R. Lechner, F. Menter CFX Development, ANSYS Germany GmbH, Staudenfeldweg
More informationEuler-Euler Modeling of Mass-Transfer in Bubbly Flows
Euler-Euler Modeling of Mass-Transfer in Bubbly Flows Roland Rzehak Eckhard Krepper Text optional: Institutsname Prof. Dr. Hans Mustermann www.fzd.de Mitglied der Leibniz-Gemeinschaft Overview Motivation
More informationIHMTC EULER-EULER TWO-FLUID MODEL BASED CODE DEVELOPMENT FOR TWO-PHASE FLOW SYSTEMS
Proceedings of the 24th National and 2nd International ISHMT-ASTFE Heat and Mass Transfer Conference (IHMTC-2017), December 27-30, 2017, BITS-Pilani, Hyderabad, India IHMTC2017-13-0160 EULER-EULER TWO-FLUID
More informationModelling of Gas-Liquid Two-Phase Flows in Vertical Pipes using PHOENICS
Modelling of Gas-Liquid Two-Phase Flows in Vertical Pipes using PHOENICS Vladimir Agranat, Masahiro Kawaji, Albert M.C. Chan* Department of Chemical Engineering and Applied Chemistry University of Toronto,
More informationEulerian model for the prediction of nucleate boiling of refrigerant in heat exchangers
Advanced Computational Methods and Experiments in Heat Transfer XI 51 Eulerian model for the prediction of nucleate boiling of refrigerant in heat exchangers D. Simón, M. C. Paz, A. Eirís&E.Suárez E.T.S.
More informationTurbulence Dispersion Force Physics, Model Derivation and Evaluation
Turbulence Dispersion Force Physics, Model Derivation and Evaluation J.-M. Shi, T. Frank, A. Burns 3 Institute of Safety Research, FZ Rossendorf shi@fz-rossendorf.de ANSYS CFX Germany 3 ANSYS CFX FZR ANSYS
More informationCFD-Modeling of Boiling Processes
CFD-Modeling of Boiling Processes 1 C. Lifante 1, T. Frank 1, A. Burns 2, E. Krepper 3, R. Rzehak 3 conxita.lifante@ansys.com 1 ANSYS Germany, 2 ANSYS UK, 3 HZDR Outline Introduction Motivation Mathematical
More informationDevelopment and Validation of the Wall Boiling Model in ANSYS CFD
Development and Validation of the Wall Boiling Model in ANSYS CFD Th. Frank, C. Lifante, A.D. Burns PBU, Funded CFD Development ANSYS Germany Thomas.Frank@ansys.com E. Krepper, R. Rzehak FZ Dresden-Rossendorf
More informationModelling of Break-up and Coalescence in Bubbly Two-Phase Flows
Modelling of Break-up and Coalescence in Bubbly Two-Phase Flows Simon Lo and Dongsheng Zhang CD-adapco, Trident Park, Didcot OX 7HJ, UK e-mail: simon.lo@uk.cd-adapco.com Abstract Numerical simulations
More informationValidation of Multiphase Flow Modeling in ANSYS CFD
Validation of Multiphase Flow Modeling in ANSYS CFD Th. Frank, C. Lifante, A.D. Burns Head Funded CFD Development ANSYS Germany Thomas.Frank@ansys.com 2009 ANSYS, Inc. All rights reserved. 1 ANSYS, Inc.
More informationModeling Complex Flows! Direct Numerical Simulations! Computational Fluid Dynamics!
http://www.nd.edu/~gtryggva/cfd-course/! Modeling Complex Flows! Grétar Tryggvason! Spring 2011! Direct Numerical Simulations! In direct numerical simulations the full unsteady Navier-Stokes equations
More informationNUMERICAL SIMULATION OF SLUG FLOW REGIME FOR AN AIR-WATER TWO-PHASE FLOW IN HORIZONTAL PIPES
The 11 th International Topical Meeting on Nuclear Reactor Thermal-Hydraulics (NURETH-11) Paper: 038 NUMERICAL SIMULATION OF SLUG FLOW REGIME FOR AN AIR-WATER TWO-PHASE FLOW IN HORIZONTAL PIPES Thomas
More informationExperimental and Numerical Investigation of Two- Phase Flow through Enlarging Singularity
Purdue University Purdue e-pubs International Refrigeration and Air Conditioning Conference School of Mechanical Engineering 212 Experimental and Numerical Investigation of Two- Phase Flow through Enlarging
More informationCFD SIMULATION OF THE DEPARTURE FROM NUCLEATE BOILING
CFD SIMULATION OF THE DEPARTURE FROM NUCLEATE BOILING Ladislav Vyskocil and Jiri Macek UJV Rez a. s., Dept. of Safety Analyses, Hlavni 130, 250 68 Husinec Rez, Czech Republic Ladislav.Vyskocil@ujv.cz;
More informationModeling of Wall-boiling Phenomena from Nucleate Subcooled Boiling up to CHF Conditions
Modeling of Wall-boiling Phenomena from Nucleate Subcooled Boiling up to CHF Conditions Thomas Frank (1), Amine Ben Hadj Ali (1), Conxita Lifante (1), Florian Kaiser (2), Stephan Gabriel (2), Henning Eickenbusch
More informationJ.-M. Laviéville EDF - Fluid Dynamics Power Generation Environment, 6 quai Watier, Chatou, France,
MODELING FREE SURFACE FLOWS RELEVANT TO A PTS SCENARIO: COMPARISON BETWEEN EXPERIMENTAL DATA AND THREE RANS BASED CFD-CODES. COMMENTS ON THE CFD-EXPERIMENT INTEGRATION AND BEST PRACTICE GUIDELINE Yann
More informationOn the validity of the twofluid model for simulations of bubbly flow in nuclear reactors
On the validity of the twofluid model for simulations of bubbly flow in nuclear reactors Henrik Ström 1, Srdjan Sasic 1, Klas Jareteg 2, Christophe Demazière 2 1 Division of Fluid Dynamics, Department
More informationEXPERIMENTS AND NUMERICAL SIMULATIONS OF HORIZONTAL TWO PHASE FLOW REGIMES
Seventh International Conference on CF in the Minerals and Process Industries CSIRO, Melbourne, Australia 9-11 ecember 009 EXPERIMENTS AN NUMERICA SIMUATIONS OF HORIZONTA TWO PHASE FOW REGIMES Thomas HÖHNE
More informationMultiphase Flow and Heat Transfer
Multiphase Flow and Heat Transfer ME546 -Sudheer Siddapureddy sudheer@iitp.ac.in Two Phase Flow Reference: S. Mostafa Ghiaasiaan, Two-Phase Flow, Boiling and Condensation, Cambridge University Press. http://dx.doi.org/10.1017/cbo9780511619410
More informationINVESTIGATION OF THE PWR SUBCHANNEL VOID DISTRIBUTION BENCHMARK (OECD/NRC PSBT BENCHMARK) USING ANSYS CFX
INVESTIGATION OF THE PWR SUBCHANNEL VOID DISTRIBUTION BENCHMARK (OECD/NRC PSBT BENCHMARK) USING ANSYS CFX Th. Frank 1, F. Reiterer 1 and C. Lifante 1 1 ANSYS Germany GmbH, Otterfing, Germany Thomas.Frank@ansys.com,
More informationModelling multiphase flows in the Chemical and Process Industry
Modelling multiphase flows in the Chemical and Process Industry Simon Lo 9/11/09 Contents Breakup and coalescence in bubbly flows Particle flows with the Discrete Element Modelling approach Multiphase
More informationTurbulent Boundary Layers & Turbulence Models. Lecture 09
Turbulent Boundary Layers & Turbulence Models Lecture 09 The turbulent boundary layer In turbulent flow, the boundary layer is defined as the thin region on the surface of a body in which viscous effects
More informationIHTC DRAFT MEASUREMENT OF LIQUID FILM THICKNESS IN MICRO TUBE ANNULAR FLOW
DRAFT Proceedings of the 14 th International Heat Transfer Conference IHTC14 August 8-13, 2010, Washington D.C., USA IHTC14-23176 MEASUREMENT OF LIQUID FILM THICKNESS IN MICRO TUBE ANNULAR FLOW Hiroshi
More informationCFD Simulation of Sodium Boiling in Heated Pipe using RPI Model
Proceedings of the 2 nd World Congress on Momentum, Heat and Mass Transfer (MHMT 17) Barcelona, Spain April 6 8, 2017 Paper No. ICMFHT 114 ISSN: 2371-5316 DOI: 10.11159/icmfht17.114 CFD Simulation of Sodium
More informationThis is a repository copy of Multiphase turbulence in bubbly flows: RANS simulations.
This is a repository copy of Multiphase turbulence in bubbly flows: RANS simulations. White Rose Research Online URL for this paper: http://eprints.whiterose.ac.uk/90676/ Version: Accepted Version Article:
More informationBest Practice Guidelines for Computational Turbulent Dispersed Multiphase Flows. René V.A. Oliemans
Best Practice Guidelines for Computational Turbulent Dispersed Multiphase Flows René V.A. Oliemans ERCOFTAC Seminar, Innventia, Stockholm, June 7-8, 2011 1 Vermelding onderdeel organisatie Department of
More informationCFD analysis of the transient flow in a low-oil concentration hydrocyclone
CFD analysis of the transient flow in a low-oil concentration hydrocyclone Paladino, E. E. (1), Nunes, G. C. () and Schwenk, L. (1) (1) ESSS Engineering Simulation and Scientific Software CELTA - Rod SC-41,
More informationCFD in COMSOL Multiphysics
CFD in COMSOL Multiphysics Mats Nigam Copyright 2016 COMSOL. Any of the images, text, and equations here may be copied and modified for your own internal use. All trademarks are the property of their respective
More informationOMAE FLUID-STRUCTURE INTERACTION MODELING OF SUBSEA JUMPER PIPE
Proceedings of the ASME 2014 33 rd International Conference on Ocean, Offshore and Arctic Engineering OMAE2014 June 8-13, 2014, San Francisco, CA USA OMAE2014-24070 FLUID-STRUCTURE INTERACTION MODELING
More informationCFD SIMULATIONS OF SINGLE AND TWO-PHASE MIXING PROESSES IN STIRRED TANK REACTORS
CFD SIMULATIONS OF SINGLE AND TWO-PHASE MIXING PROESSES IN STIRRED TANK REACTORS Hristo Vesselinov Hristov, Stephan Boden, Günther Hessel, Holger Kryk, Horst-Michael Prasser, and Wilfried Schmitt. Introduction
More informationCFD modelling of multiphase flows
1 Lecture CFD-3 CFD modelling of multiphase flows Simon Lo CD-adapco Trident House, Basil Hill Road Didcot, OX11 7HJ, UK simon.lo@cd-adapco.com 2 VOF Free surface flows LMP Droplet flows Liquid film DEM
More informationA STUDY ON SLUG INDUCED STRESSES USING FILE-BASED COUPLING TECHNIQUE
A STUDY ON SLUG INDUCED STRESSES USING FILE-BASED COUPLING TECHNIQUE Abdalellah O. Mohmmed, Mohammad S. Nasif and Hussain H. Al-Kayiem Department of Mechanical Engineering, Universiti Teknologi Petronas,
More informationInvestigation of slug flow characteristics in inclined pipelines
Computational Methods in Multiphase Flow IV 185 Investigation of slug flow characteristics in inclined pipelines J. N. E. Carneiro & A. O. Nieckele Department of Mechanical Engineering Pontifícia Universidade
More informationHEAT TRANSFER IN A RECIRCULATION ZONE AT STEADY-STATE AND OSCILLATING CONDITIONS - THE BACK FACING STEP TEST CASE
HEAT TRANSFER IN A RECIRCULATION ZONE AT STEADY-STATE AND OSCILLATING CONDITIONS - THE BACK FACING STEP TEST CASE A.K. Pozarlik 1, D. Panara, J.B.W. Kok 1, T.H. van der Meer 1 1 Laboratory of Thermal Engineering,
More informationValidation analyses of advanced turbulence model approaches for stratified two-phase flows
Computational Methods in Multiphase Flow VIII 361 Validation analyses of advanced turbulence model approaches for stratified two-phase flows M. Benz & T. Schulenberg Institute for Nuclear and Energy Technologies,
More informationVALIDATION OF CFD-BWR, A NEW TWO-PHASE COMPUTATIONAL FLUID DYNAMICS MODEL FOR BOILING WATER REACTOR ANALYSIS
VALIDATION OF CFD-BWR, A NEW TWO-PHASE COMPUTATIONAL FLUID DYNAMICS MODEL FOR BOILING WATER REACTOR ANALYSIS V.Ustineno 1, M.Samigulin 1, A.Ioilev 1, S.Lo 2, A.Tentner 3, A.Lychagin 4, A.Razin 4, V.Girin
More informationCFD SIMULATION OF SOLID-LIQUID STIRRED TANKS
CFD SIMULATION OF SOLID-LIQUID STIRRED TANKS Divyamaan Wadnerkar 1, Ranjeet P. Utikar 1, Moses O. Tade 1, Vishnu K. Pareek 1 Department of Chemical Engineering, Curtin University Perth, WA 6102 r.utikar@curtin.edu.au
More informationModel Studies on Slag-Metal Entrainment in Gas Stirred Ladles
Model Studies on Slag-Metal Entrainment in Gas Stirred Ladles Anand Senguttuvan Supervisor Gordon A Irons 1 Approach to Simulate Slag Metal Entrainment using Computational Fluid Dynamics Introduction &
More informationD.R. Rector, M.L. Stewart and A.P. Poloski Pacific Northwest National Laboratory P.O. Box 999, Richland, WA
Modeling of Sediment Bed Behavior for Critical Velocity in Horizontal Piping 9263 D.R. Rector, M.L. Stewart and A.P. Poloski Pacific Northwest National Laboratory P.O. Box 999, Richland, WA ABSTRACT A
More informationMULTIDIMENSIONAL TURBULENCE SPECTRA - STATISTICAL ANALYSIS OF TURBULENT VORTICES
Ninth International Conference on CFD in the Minerals and Process Industries CSIRO, Melbourne, Australia 10-12 December 2012 MULTIDIMENSIONAL TURBULENCE SPECTRA - STATISTICAL ANALYSIS OF TURBULENT VORTICES
More informationA Comparative Analysis of Turbulent Pipe Flow Using k And k Models
A Comparative Analysis of Turbulent Pipe Flow Using k And k Models 1 B. K. Menge, 2 M. Kinyanjui, 3 J. K. Sigey 1 Department of Mathematics and Physics. Technical University of Mombasa, P.O BOX 90420-80100,Mombasa,
More informationNumerical Simulation of the Hagemann Entrainment Experiments
CCC Annual Report UIUC, August 14, 2013 Numerical Simulation of the Hagemann Entrainment Experiments Kenneth Swartz (BSME Student) Lance C. Hibbeler (Ph.D. Student) Department of Mechanical Science & Engineering
More informationThe Reynolds experiment
Chapter 13 The Reynolds experiment 13.1 Laminar and turbulent flows Let us consider a horizontal pipe of circular section of infinite extension subject to a constant pressure gradient (see section [10.4]).
More informationPrediction of Minimum Fluidisation Velocity Using a CFD-PBM Coupled Model in an Industrial Gas Phase Polymerisation Reactor
Journal of Engineering Science, Vol. 10, 95 105, 2014 Prediction of Minimum Fluidisation Velocity Using a CFD-PBM Coupled Model in an Industrial Gas Phase Polymerisation Reactor Vahid Akbari and Mohd.
More informationComparison of the Effects of k-ϵ, k-ω, and Zero Equation Models on Characterization of Turbulent Permeability of Porous Media
Comparison of the Effects of k-ϵ, k-ω, and Zero Equation Models on Characterization of Turbulent Permeability of Porous Media Charlie Matsubara *, Tim Kuo, Helen Wu PIM Consulting, No. 16, Gonyequ 1 st
More informationFlow pattern identification based on a single-wire capacitance probe
Nuclear Science and Techniques 21 (2010) 246 250 Flow pattern identification based on a single-wire capacitance probe HUANG Shanfang 1,2,* ZHANG Bingdong 2 LU Jun 2 WANG Dong 2 YANG Yanhua 3 1 Department
More informationINTERACTION OF AN AIR-BUBBLE DISPERSED PHASE WITH AN INITIALLY ISOTROPIC TURBULENT FLOW FIELD
3rd Workshop on Transport Phenomena in Two-Phase Flow Nessebar, Bulgaria, 2-7 September 1998, p.p. 133-138 INTERACTION OF AN AIR-BUBBLE DISPERSED PHASE WITH AN INITIALLY ISOTROPIC TURBULENT FLOW FIELD
More informationDevelopment of twophaseeulerfoam
ISPRAS OPEN 2016 NUMERICAL STUDY OF SADDLE-SHAPED VOID FRACTION PROFILES EFFECT ON THERMAL HYDRAULIC PARAMETERS OF THE CHANNEL WITH TWO-PHASE FLOW USING OPENFOAM AND COMPARISON WITH EXPERIMENTS Varseev
More informationCENG 501 Examination Problem: Estimation of Viscosity with a Falling - Cylinder Viscometer
CENG 501 Examination Problem: Estimation of Viscosity with a Falling - Cylinder Viscometer You are assigned to design a fallingcylinder viscometer to measure the viscosity of Newtonian liquids. A schematic
More informationNumerical Study of Hydrodynamics in an External-Loop Air-Lift Reactor
Numerical Study of Hydrodynamics in an External-Loop Air-Lift Reactor Joanna Karcz *, Marcelina Bitenc, Marek Domański, Łukasz Kacperski West Pomeranian University of Technology, Szczecin, Dept. of Chem.
More informationNumerical simulations of heat transfer in plane channel flow
Numerical simulations of heat transfer in plane channel flow Najla EL GHARBI 1, 3, a, Rafik ABSI 2, b and Ahmed BENZAOUI 3, c 1 Renewable Energy Development Center, BP 62 Bouzareah 163 Algiers, Algeria
More informationCFD Simulation of Turbulent Flow Structure in Stratified Gas/Liquid Flow and Validation with Experimental Data
SPE-174964-MS CFD Simulation of Turbulent Flow Structure in Stratified Gas/Liquid Flow and Validation with Experimental Data Ramin Dabirian, Amir Mansouri, Ram Mohan, and Ovadia Shoham, The University
More informationNumerical Simulation of Gas-Liquid-Reactors with Bubbly Flows using a Hybrid Multiphase-CFD Approach
Numerical Simulation of Gas-Liquid-Reactors with Bubbly Flows using a Hybrid Multiphase-CFD Approach TFM Hybrid Interface Resolving Two-Fluid Model (HIRES-TFM) by Coupling of the Volume-of-Fluid (VOF)
More informationAnalysis of the Cooling Design in Electrical Transformer
Analysis of the Cooling Design in Electrical Transformer Joel de Almeida Mendes E-mail: joeldealmeidamendes@hotmail.com Abstract This work presents the application of a CFD code Fluent to simulate the
More informationVisualization of flow pattern over or around immersed objects in open channel flow.
EXPERIMENT SEVEN: FLOW VISUALIZATION AND ANALYSIS I OBJECTIVE OF THE EXPERIMENT: Visualization of flow pattern over or around immersed objects in open channel flow. II THEORY AND EQUATION: Open channel:
More informationInlet Diameter and Flow Volume Effects on Separation and Energy Efficiency of Hydrocyclones
IOP Conference Series: Materials Science and Engineering PAPER OPEN ACCESS Inlet Diameter and Flow Volume Effects on Separation and Energy Efficiency of Hydrocyclones To cite this article: Erikli and A
More informationCFD-Modeling of Turbulent Flows in a 3x3 Rod Bundle and Comparison to Experiments
CFD-Modeling of Turbulent Flows in a 3x3 Rod Bundle and Comparison to Experiments C. Lifante 1, B. Krull 1, Th. Frank 1, R. Franz 2, U. Hampel 2 1 PBU, ANSYS Germany, Otterfing 2 Institute of Safety Research,
More informationTHE EFFECT OF TWO PHASE (AIR-WATER) FLOW CHARACTERISTICS ON MOMENTUM FLUX DUE TO FLOW TURNING ELEMENTS AT ATMOSPHERIC CONDITIONS
International Journal of Latest Trends in Engineering and Technology Vol.(8)Issue(1), pp.319-328 DOI: http://dx.doi.org/10.21172/1.81.041 e-issn:2278-621x AN EXPERIMENTAL STUDY OF THE EFFECT OF TWO PHASE
More informationColloquium FLUID DYNAMICS 2012 Institute of Thermomechanics AS CR, v.v.i., Prague, October 24-26, 2012 p.
Colloquium FLUID DYNAMICS 212 Institute of Thermomechanics AS CR, v.v.i., Prague, October 24-26, 212 p. ON A COMPARISON OF NUMERICAL SIMULATIONS OF ATMOSPHERIC FLOW OVER COMPLEX TERRAIN T. Bodnár, L. Beneš
More informationPairwise Interaction Extended Point-Particle (PIEP) Model for droplet-laden flows: Towards application to the mid-field of a spray
Pairwise Interaction Extended Point-Particle (PIEP) Model for droplet-laden flows: Towards application to the mid-field of a spray Georges Akiki, Kai Liu and S. Balachandar * Department of Mechanical &
More informationFluid Structure Interaction Analysis of Two-Phase Flow in an M-shaped Jumper. Star Global Conference University of Houston
Fluid Structure Interaction Analysis of Two-Phase Flow in an M-shaped Jumper Star Global Conference 2012 University of Houston College of Technology Mechanical Engineering Technology Leonardo Chica January
More informationCFD simulation of gas solid bubbling fluidized bed: an extensive assessment of drag models
Computational Methods in Multiphase Flow IV 51 CFD simulation of gas solid bubbling fluidized bed: an extensive assessment of drag models N. Mahinpey 1, F. Vejahati 1 & N. Ellis 2 1 Environmental Systems
More informationMultiphase Flow Modeling & Simulation with Application to Water-Vapor Flows Through Fuel Rod Bundles of Nuclear Reactors
Multiphase Flow Modeling & Simulation with Application to Water-Vapor Flows Through Fuel Rod Bundles of Nuclear Reactors Thomas Frank ANSYS Germany, Otterfing Thomas.Frank@ansys.com 2006 ANSYS, Inc. Th.
More informationTHERMAL HYDRAULIC REACTOR CORE CALCULATIONS BASED ON COUPLING THE CFD CODE ANSYS CFX WITH THE 3D NEUTRON KINETIC CORE MODEL DYN3D
THERMAL HYDRAULIC REACTOR CORE CALCULATIONS BASED ON COUPLING THE CFD CODE ANSYS CFX WITH THE 3D NEUTRON KINETIC CORE MODEL DYN3D A. Grahn, S. Kliem, U. Rohde Forschungszentrum Dresden-Rossendorf, Institute
More informationComputational fluid dynamics study of flow depth in an open Venturi channel for Newtonian fluid
Computational fluid dynamics study of flow depth in an open Venturi channel for Newtonian fluid Prasanna Welahettige 1, Bernt Lie 1, Knut Vaagsaether 1 1 Department of Process, Energy and Environmental
More informationINTERNAL FLOW IN A Y-JET ATOMISER ---NUMERICAL MODELLING---
ILASS-Europe 2002 Zaragoza 9 11 September 2002 INTERNAL FLOW IN A Y-JET ATOMISER ---NUMERICAL MODELLING--- Z. Tapia, A. Chávez e-mail: ztapia@imp.mx Instituto Mexicano del Petróleo Blvd. Adolfo Ruiz Cortines
More informationDetailed 3D modelling of mass transfer processes in two phase flows with dynamic interfaces
Detailed 3D modelling of mass transfer processes in two phase flows with dynamic interfaces D. Darmana, N.G. Deen, J.A.M. Kuipers Fundamentals of Chemical Reaction Engineering, Faculty of Science and Technology,
More informationActive Control of Separated Cascade Flow
Chapter 5 Active Control of Separated Cascade Flow In this chapter, the possibility of active control using a synthetic jet applied to an unconventional axial stator-rotor arrangement is investigated.
More informationINTRODUCTION TO MULTIPHASE FLOW. Mekanika Fluida II -Haryo Tomo-
1 INTRODUCTION TO MULTIPHASE FLOW Mekanika Fluida II -Haryo Tomo- 2 Definitions Multiphase flow is simultaneous flow of Matters with different phases( i.e. gas, liquid or solid). Matters with different
More informationCHAPTER 7 NUMERICAL MODELLING OF A SPIRAL HEAT EXCHANGER USING CFD TECHNIQUE
CHAPTER 7 NUMERICAL MODELLING OF A SPIRAL HEAT EXCHANGER USING CFD TECHNIQUE In this chapter, the governing equations for the proposed numerical model with discretisation methods are presented. Spiral
More informationSimulating Interfacial Tension of a Falling. Drop in a Moving Mesh Framework
Simulating Interfacial Tension of a Falling Drop in a Moving Mesh Framework Anja R. Paschedag a,, Blair Perot b a TU Berlin, Institute of Chemical Engineering, 10623 Berlin, Germany b University of Massachusetts,
More informationComparison of Turbulence Models in the Flow over a Backward-Facing Step Priscila Pires Araujo 1, André Luiz Tenório Rezende 2
Comparison of Turbulence Models in the Flow over a Backward-Facing Step Priscila Pires Araujo 1, André Luiz Tenório Rezende 2 Department of Mechanical and Materials Engineering, Military Engineering Institute,
More informationCFD SIMULATION OF INDUSTRIAL BUBBLE COLUMNS: NUMERICAL AND MODELLING CHALLENGES AND SUCCESSES
Eleventh International Conference on CFD in the Minerals and Process Industries CSIRO, Melbourne, Australia 7-9 December 2015 CFD SIMULATION OF INDUSTRIAL BUBBLE COLUMNS: NUMERICAL AND MODELLING CHALLENGES
More informationA NUMERICAL ANALYSIS OF COMBUSTION PROCESS IN AN AXISYMMETRIC COMBUSTION CHAMBER
SCIENTIFIC RESEARCH AND EDUCATION IN THE AIR FORCE-AFASES 2016 A NUMERICAL ANALYSIS OF COMBUSTION PROCESS IN AN AXISYMMETRIC COMBUSTION CHAMBER Alexandru DUMITRACHE*, Florin FRUNZULICA ** *Institute of
More informationSIEVE TRAY EFFICIENCY USING CFD MODELING AND SIMULATION
SIEVE TRAY EFFICIENCY USING CFD MODELING AND SIMULATION Getye Gesit* School of Chemical and Bio Engineering Addis Ababa Institute of Technology, Addis Ababa University ABSTRACT In this work, computational
More informationEFFECT OF LIQUID PHASE COMPRESSIBILITY ON MODELING OF GAS-LIQUID TWO-PHASE FLOWS USING TWO-FLUID MODEL
EFFECT OF LIQUID PHASE COMPRESSIBILITY ON MODELING OF GAS-LIQUID TWO-PHASE FLOWS USING TWO-FLUID MODEL Vahid SHOKRI 1*,Kazem ESMAEILI 2 1,2 Department of Mechanical Engineering, Sari Branch, Islamic Azad
More informationJournal of NUCLEAR SCIENCE and TECHNOLOGY, Vol. 41, No. 7, p (July 2004)
Journal of NUCLEAR SCIENCE and TECHNOLOGY, Vol. 41, No. 7, p. 765 770 (July 2004) TECHNICAL REPORT Experimental and Operational Verification of the HTR-10 Once-Through Steam Generator (SG) Heat-transfer
More informationNumerical Study of Oil Flow in Tee Junction with Leaks
Advances in Petroleum Exploration and Development Vol. 6, No. 2, 2013, pp. 1-11 DOI:10.3968/j.aped.1925543820130602.1803 ISSN 1925-542X [Print] ISSN 1925-5438 [Online] www.cscanada.net www.cscanada.org
More informationMixing and Evaporation of Liquid Droplets Injected into an Air Stream Flowing at all Speeds
Mixing and Evaporation of Liquid Droplets Injected into an Air Stream Flowing at all Speeds F. Moukalled* and M. Darwish American University of Beirut Faculty of Engineering & Architecture Mechanical Engineering
More informationHEAT TRANSFER CAPABILITY OF A THERMOSYPHON HEAT TRANSPORT DEVICE WITH EXPERIMENTAL AND CFD STUDIES
HEAT TRANSFER CAPABILITY OF A THERMOSYPHON HEAT TRANSPORT DEVICE WITH EXPERIMENTAL AND CFD STUDIES B.M. Lingade a*, Elizabeth Raju b, A Borgohain a, N.K. Maheshwari a, P.K.Vijayan a a Reactor Engineering
More informationThe effect of momentum flux ratio and turbulence model on the numerical prediction of atomization characteristics of air assisted liquid jets
ILASS Americas, 26 th Annual Conference on Liquid Atomization and Spray Systems, Portland, OR, May 204 The effect of momentum flux ratio and turbulence model on the numerical prediction of atomization
More informationThis is a repository copy of CFD Simulation of Boiling Flows with an Eulerian-Eulerian Two-Fluid Model.
This is a repository copy of CFD Simulation of Boiling Flows with an Eulerian-Eulerian Two-Fluid Model. White Rose Research Online URL for this paper: http://eprints.whiterose.ac.uk/128021/ Version: Accepted
More informationApplication of the immersed boundary method to simulate flows inside and outside the nozzles
Application of the immersed boundary method to simulate flows inside and outside the nozzles E. Noël, A. Berlemont, J. Cousin 1, T. Ménard UMR 6614 - CORIA, Université et INSA de Rouen, France emeline.noel@coria.fr,
More informationSimplified Model of WWER-440 Fuel Assembly for ThermoHydraulic Analysis
1 Portál pre odborné publikovanie ISSN 1338-0087 Simplified Model of WWER-440 Fuel Assembly for ThermoHydraulic Analysis Jakubec Jakub Elektrotechnika 13.02.2013 This work deals with thermo-hydraulic processes
More informationLaboratory of Thermal Hydraulics. General Overview
Visit of Nuclear Master Students Laboratory of Thermal Hydraulics General Overview Horst-Michael Prasser December 04, 2009 Paul Scherrer Institut Main Goals Development of analytical and experimental methods
More informationTurbulence - Theory and Modelling GROUP-STUDIES:
Lund Institute of Technology Department of Energy Sciences Division of Fluid Mechanics Robert Szasz, tel 046-0480 Johan Revstedt, tel 046-43 0 Turbulence - Theory and Modelling GROUP-STUDIES: Turbulence
More informationCFD Modelling of the Fluid Flow Characteristics in an External-Loop Air-Lift Reactor
A publication of 145 CHEMICAL ENGINEERING TRANSACTIONS VOL. 2, 21 Chief Editors: Sauro Pierucci, Jiří J. Klemeš Copyright 21, AIDIC Servizi S.r.l., ISBN 978-88-958-2-5; ISSN 1974-9791 The Italian Association
More informationCFD MODEL FOR DETERMINING LOCAL PHASE FRACTION OIL-WATER DISPERSION IN TURBULENT FLOW
CFD MODEL FOR DETERMINING LOCAL PHASE FRACTION OIL-WATER DISPERSION IN TURBULENT FLOW Siti Aslina Hussain 1* and Soo Mee Khuan 1 1 Department of Chemical and Environment Engineering, Faculty of Engineering,
More informationVERTICAL TURBULENT BUOYANT HELIUM JET CFD MODELING AND VALIDATION
VERTICAL TURBULENT BUOYANT HELIUM JET CFD MODELING AND VALIDATION Cheng Z, Agranat V.M. and Tchouvelev A.V. A.V.Tchouvelev & Associates, Inc., 659 Spinnaker Circle, Mississauga, Ontario, Canada L5W R Hydrogenics
More informationLiquid Feed Injection in a High Density Riser
Refereed Proceedings The 12th International Conference on Fluidization - New Horizons in Fluidization Engineering Engineering Conferences International Year 2007 Liquid Feed Injection in a High Density
More informationMultiphase Science and Technology, Vol. 16, Nos. 1-4, pp. 1-20, 2005
Multiphase Science and Technology, Vol. 16, Nos. 1-4, pp. 1-2, 25 EXPERIMENTS ON THE TURBULENT STRUCTURE AND THE VOID FRACTION DISTRIBUTION IN THE TAYLOR BUBBLE WAKE L. Shemer, A. Gulitski and D. Barnea
More informationDEVELOPMENT OF COMPUTATIONAL MULTIFLUID DYNAMICS MODELS FOR NUCLEAR REACTOR APPLICATIONS
DEVELOPMENT OF COMPUTATIONAL MULTIFLUID DYNAMICS MODELS FOR NUCLEAR REACTOR APPLICATIONS Henry Anglart Royal Institute of Technology, Department of Physics Division of Nuclear Reactor Technology Stocholm,
More informationDEVELOPMENT OF INTERFACIAL AREA TRANSPORT EQUATION
DEVELOPMENT OF INTERFACIAL AREA TRANSPORT EQUATION MAMORU ISHII *, SEUNGJIN KIM 1 and JOSEPH KELLY 2 School of Nuclear Engineering Purdue University 400 Central Drive, West Lafaytte, IN 47907-2017, USA
More informationNumerical analysis of fluid flow and heat transfer in 2D sinusoidal wavy channel
Numerical analysis of fluid flow and heat transfer in 2D sinusoidal wavy channel Arunanshu Chakravarty 1* 1 CTU in Prague, Faculty of Mechanical Engineering, Department of Process Engineering,Technická
More informationCalculations on a heated cylinder case
Calculations on a heated cylinder case J. C. Uribe and D. Laurence 1 Introduction In order to evaluate the wall functions in version 1.3 of Code Saturne, a heated cylinder case has been chosen. The case
More informationModeling of dispersed phase by Lagrangian approach in Fluent
Lappeenranta University of Technology From the SelectedWorks of Kari Myöhänen 2008 Modeling of dispersed phase by Lagrangian approach in Fluent Kari Myöhänen Available at: https://works.bepress.com/kari_myohanen/5/
More information150A Review Session 2/13/2014 Fluid Statics. Pressure acts in all directions, normal to the surrounding surfaces
Fluid Statics Pressure acts in all directions, normal to the surrounding surfaces or Whenever a pressure difference is the driving force, use gauge pressure o Bernoulli equation o Momentum balance with
More informationNumerical simulation of fluid flow in a monolithic exchanger related to high temperature and high pressure operating conditions
Advanced Computational Methods in Heat Transfer X 25 Numerical simulation of fluid flow in a monolithic exchanger related to high temperature and high pressure operating conditions F. Selimovic & B. Sundén
More informationNUMERICAL AND EXPERIMENTAL INVESTIGATION OF THE TEMPERATURE DISTRIBUTION INSIDE OIL-COOLED TRANSFORMER WINDINGS
NUMERICAL AND EXPERIMENTAL INVESTIGATION OF THE TEMPERATURE DISTRIBUTION INSIDE OIL-COOLED TRANSFORMER WINDINGS N. Schmidt 1* and S. Tenbohlen 1 and S. Chen 2 and C. Breuer 3 1 University of Stuttgart,
More information