Assessment of RELAP5-3D copyright using data from two-dimensional RPI flow tests

International Nuclear Information System (INIS)

Davis, C.B.

1998-01-01

The capability of the RELAP5-3D copyright computer code to perform multi-dimensional thermal-hydraulic analysis was assessed using data from steady-state flow tests conducted at Rensselaer Polytechnic Institute (RPI). The RPI data were taken in a two-dimensional test section in a low-pressure air/water loop. The test section consisted of a thin vertical channel that simulated a two-dimensional slice through the core of a pressurized water reactor. Single-phase and two-phase flows were supplied to the test section in an asymmetric manner to generate a two-dimensional flow field. A traversing gamma densitometer was used to measure void fraction at many locations in the test section. High speed photographs provided information on the flow patterns and flow regimes. The RPI test section was modeled using the multi-dimensional component in RELAP5-3D Version BF06. Calculations of three RPI experiments were performed. The flow regimes predicted by the base code were in poor agreement with those observed in the tests. The two-phase regions were observed to be in the bubbly and slug flow regimes in the test. However, nearly all of the junctions in the horizontal direction were calculated to be in the stratified flow regime because of the relatively low velocities in that direction. As a result, the void fraction predictions were also in poor agreement with the measured values. Significantly improved results were obtained in sensitivity calculations with a modified version of the code that prevented the horizontal junctions from entering the stratified flow regime. These results indicate that the code's logic in the determination of flow regimes in a multi-dimensional component must be improved. The results of the sensitivity calculations also indicate that RELAP5-3D will provide a significant multi-dimensional hydraulic analysis capability once the flow regime prediction is improved

Multi-dimensional rheology-based two-phase model for sediment transport and applications to sheet flow and pipeline scour

International Nuclear Information System (INIS)

Lee, Cheng-Hsien; Low, Ying Min; Chiew, Yee-Meng

2016-01-01

Sediment transport is fundamentally a two-phase phenomenon involving fluid and sediments; however, many existing numerical models are one-phase approaches, which are unable to capture the complex fluid-particle and inter-particle interactions. In the last decade, two-phase models have gained traction; however, there are still many limitations in these models. For example, several existing two-phase models are confined to one-dimensional problems; in addition, the existing two-dimensional models simulate only the region outside the sand bed. This paper develops a new three-dimensional two-phase model for simulating sediment transport in the sheet flow condition, incorporating recently published rheological characteristics of sediments. The enduring-contact, inertial, and fluid viscosity effects are considered in determining sediment pressure and stresses, enabling the model to be applicable to a wide range of particle Reynolds number. A k − ε turbulence model is adopted to compute the Reynolds stresses. In addition, a novel numerical scheme is proposed, thus avoiding numerical instability caused by high sediment concentration and allowing the sediment dynamics to be computed both within and outside the sand bed. The present model is applied to two classical problems, namely, sheet flow and scour under a pipeline with favorable results. For sheet flow, the computed velocity is consistent with measured data reported in the literature. For pipeline scour, the computed scour rate beneath the pipeline agrees with previous experimental observations. However, the present model is unable to capture vortex shedding; consequently, the sediment deposition behind the pipeline is overestimated. Sensitivity analyses reveal that model parameters associated with turbulence have strong influence on the computed results.

Two-phase flow operational maps for multi-microchannel evaporators

International Nuclear Information System (INIS)

Szczukiewicz, Sylwia; Borhani, Navid; Thome, John Richard

2013-01-01

Highlights: • New operational maps for several different micro-evaporators are presented. • Inlet micro-orifices prevented flow instability, back flow, and flow maldistribution. • Eight different operating regimes were distinguished. • The flashing two-phase flow without back flow operating regime is preferred. -- Abstract: The current paper presents new operational maps for several different multi-microchannel evaporators, with and without any inlet restrictions (micro-orifices), for the two-phase flow of refrigerants R245fa, R236fa, and R1234ze(E). The test fluids flowed in 67 parallel channels, each having a cross-sectional area of 100 × 100 μm 2 . In order to emulate the power dissipated by active components in a 3D CMOS CPU chip, two aluminium microheaters were sputtered onto the back-side of the test section providing a 0.5 cm 2 each. Without any inlet restrictions in the micro-evaporator, significant parallel channel flow instabilities, vapor back flow, and flow maldistribution led to high-amplitude and high-frequency temperature and pressure oscillations. Such undesired phenomena were then prevented by placing restrictions at the inlet of each channel. High-speed flow visualization distinguished eight different operating regimes of the two-phase flow depending on the tested operating conditions. Therefore, the preferred operating regimes can be easily traced. In particular, flashing two-phase flow without back flow appeared to be the best operating regime without any flow and temperature instabilities

Approximate characteristics for one-dimensional two-phase flows

International Nuclear Information System (INIS)

Sarayloo, A.; Peddleson, J.

1985-01-01

An approximate method for determining the characteristics associated with one-dimensional particulate two-phase flow models is presented. The method is based on iteration and is valid for small particulate volume fractions. The method is applied to several special cases involving incompressible particles suspended in a gas. The influences of certain changes in the physical model are investigated

Advanced numerical methods for three dimensional two-phase flow calculations

Energy Technology Data Exchange (ETDEWEB)

Toumi, I. [Laboratoire d`Etudes Thermiques des Reacteurs, Gif sur Yvette (France); Caruge, D. [Institut de Protection et de Surete Nucleaire, Fontenay aux Roses (France)

1997-07-01

This paper is devoted to new numerical methods developed for both one and three dimensional two-phase flow calculations. These methods are finite volume numerical methods and are based on the use of Approximate Riemann Solvers concepts to define convective fluxes versus mean cell quantities. The first part of the paper presents the numerical method for a one dimensional hyperbolic two-fluid model including differential terms as added mass and interface pressure. This numerical solution scheme makes use of the Riemann problem solution to define backward and forward differencing to approximate spatial derivatives. The construction of this approximate Riemann solver uses an extension of Roe`s method that has been successfully used to solve gas dynamic equations. As far as the two-fluid model is hyperbolic, this numerical method seems very efficient for the numerical solution of two-phase flow problems. The scheme was applied both to shock tube problems and to standard tests for two-fluid computer codes. The second part describes the numerical method in the three dimensional case. The authors discuss also some improvements performed to obtain a fully implicit solution method that provides fast running steady state calculations. Such a scheme is not implemented in a thermal-hydraulic computer code devoted to 3-D steady-state and transient computations. Some results obtained for Pressurised Water Reactors concerning upper plenum calculations and a steady state flow in the core with rod bow effect evaluation are presented. In practice these new numerical methods have proved to be stable on non staggered grids and capable of generating accurate non oscillating solutions for two-phase flow calculations.

Advanced numerical methods for three dimensional two-phase flow calculations

International Nuclear Information System (INIS)

Toumi, I.; Caruge, D.

1997-01-01

This paper is devoted to new numerical methods developed for both one and three dimensional two-phase flow calculations. These methods are finite volume numerical methods and are based on the use of Approximate Riemann Solvers concepts to define convective fluxes versus mean cell quantities. The first part of the paper presents the numerical method for a one dimensional hyperbolic two-fluid model including differential terms as added mass and interface pressure. This numerical solution scheme makes use of the Riemann problem solution to define backward and forward differencing to approximate spatial derivatives. The construction of this approximate Riemann solver uses an extension of Roe's method that has been successfully used to solve gas dynamic equations. As far as the two-fluid model is hyperbolic, this numerical method seems very efficient for the numerical solution of two-phase flow problems. The scheme was applied both to shock tube problems and to standard tests for two-fluid computer codes. The second part describes the numerical method in the three dimensional case. The authors discuss also some improvements performed to obtain a fully implicit solution method that provides fast running steady state calculations. Such a scheme is not implemented in a thermal-hydraulic computer code devoted to 3-D steady-state and transient computations. Some results obtained for Pressurised Water Reactors concerning upper plenum calculations and a steady state flow in the core with rod bow effect evaluation are presented. In practice these new numerical methods have proved to be stable on non staggered grids and capable of generating accurate non oscillating solutions for two-phase flow calculations

Multi-scale diffuse interface modeling of multi-component two-phase flow with partial miscibility

Science.gov (United States)

Kou, Jisheng; Sun, Shuyu

2016-08-01

In this paper, we introduce a diffuse interface model to simulate multi-component two-phase flow with partial miscibility based on a realistic equation of state (e.g. Peng-Robinson equation of state). Because of partial miscibility, thermodynamic relations are used to model not only interfacial properties but also bulk properties, including density, composition, pressure, and realistic viscosity. As far as we know, this effort is the first time to use diffuse interface modeling based on equation of state for modeling of multi-component two-phase flow with partial miscibility. In numerical simulation, the key issue is to resolve the high contrast of scales from the microscopic interface composition to macroscale bulk fluid motion since the interface has a nanoscale thickness only. To efficiently solve this challenging problem, we develop a multi-scale simulation method. At the microscopic scale, we deduce a reduced interfacial equation under reasonable assumptions, and then we propose a formulation of capillary pressure, which is consistent with macroscale flow equations. Moreover, we show that Young-Laplace equation is an approximation of this capillarity formulation, and this formulation is also consistent with the concept of Tolman length, which is a correction of Young-Laplace equation. At the macroscopical scale, the interfaces are treated as discontinuous surfaces separating two phases of fluids. Our approach differs from conventional sharp-interface two-phase flow model in that we use the capillary pressure directly instead of a combination of surface tension and Young-Laplace equation because capillarity can be calculated from our proposed capillarity formulation. A compatible condition is also derived for the pressure in flow equations. Furthermore, based on the proposed capillarity formulation, we design an efficient numerical method for directly computing the capillary pressure between two fluids composed of multiple components. Finally, numerical tests

Multi-scale diffuse interface modeling of multi-component two-phase flow with partial miscibility

KAUST Repository

Kou, Jisheng

2016-05-10

In this paper, we introduce a diffuse interface model to simulate multi-component two-phase flow with partial miscibility based on a realistic equation of state (e.g. Peng-Robinson equation of state). Because of partial miscibility, thermodynamic relations are used to model not only interfacial properties but also bulk properties, including density, composition, pressure, and realistic viscosity. As far as we know, this effort is the first time to use diffuse interface modeling based on equation of state for modeling of multi-component two-phase flow with partial miscibility. In numerical simulation, the key issue is to resolve the high contrast of scales from the microscopic interface composition to macroscale bulk fluid motion since the interface has a nanoscale thickness only. To efficiently solve this challenging problem, we develop a multi-scale simulation method. At the microscopic scale, we deduce a reduced interfacial equation under reasonable assumptions, and then we propose a formulation of capillary pressure, which is consistent with macroscale flow equations. Moreover, we show that Young-Laplace equation is an approximation of this capillarity formulation, and this formulation is also consistent with the concept of Tolman length, which is a correction of Young-Laplace equation. At the macroscopical scale, the interfaces are treated as discontinuous surfaces separating two phases of fluids. Our approach differs from conventional sharp-interface two-phase flow model in that we use the capillary pressure directly instead of a combination of surface tension and Young-Laplace equation because capillarity can be calculated from our proposed capillarity formulation. A compatible condition is also derived for the pressure in flow equations. Furthermore, based on the proposed capillarity formulation, we design an efficient numerical method for directly computing the capillary pressure between two fluids composed of multiple components. Finally, numerical tests

Three-dimensional investigation of the two-phase flow structure in a bubbly pipe flow

International Nuclear Information System (INIS)

Hassan, Y.A.; Schmidl, W.D.; Ortiz-Villafuerte, J.

1997-01-01

Particle Image Velocimetry (PIV) is a non-intrusive measurement technique, which can be used to study the structure of various fluid flows. PIV is used to measure the time varying full field velocity data of a particle-seeded flow field within either a two-dimensional plane or three-dimensional volume. PIV is a very efficient measurement technique since it can obtain both qualitative and quantitative spatial information about the flow field being studied. This information can be further processed into information such as vorticity and pathlines. Other flow measurement techniques (Laser Doppler Velocimetry, Hot Wire Anemometry, etc...) only provide quantitative information at a single point. PIV can be used to study turbulence structures if a sufficient amount of data can be acquired and analyzed, and it can also be extended to study two-phase flows if both phases can be distinguished. In this study, the flow structure around a bubble rising in a pipe filled with water was studied in three-dimensions. The velocity of the rising bubble and the velocity field of the surrounding water was measured. Then the turbulence intensities and Reynolds stresses were calculated from the experimental data. (author)

Local measurement and numerical calculation on turbulent two-phase flow in a vertical pipe with sudden expansion

International Nuclear Information System (INIS)

Kondo, Koichi; Yoshida, Kenji; Okawa, Tomio; Kataoka, Isao

2004-01-01

Experiment and numerical calculation were carried out for upward, turbulent bubbly two-phase flow in a vertical pipe with an axisymmetric sudden expansion, which is one of the typical multi-dimensional channel geometries. The void fraction, the liquid velocity and turbulent intensity along the flow direction below and the above the sudden expansion point were measured for various turbulent flow conditions by using a point-electrode resistivity probe and a hot-film anemometry probe. They showed quite complicated behaviors depending upon flow rates of gas and liquid phases and bubble size. In particular, the geometry of sudden expansion affected on the bubble behaviors in multi-dimensional two-phase flow, such as the bubble-stagnation, the bubble-deformation, the enhancement and suppression effects due to the two-phase turbulence etc. Through the measurements, fundamental parameters of the two-phase flow were clarified for the sudden expansion channel. Moreover, a three-dimensional one-way bubble tracking simulation of a single bubble behavior in turbulent flow field along the downstream of the sudden expansion was also demonstrated where equation of motion of bubble was solved by assuming appropriate constitutive models and turbulence model. Based on the trajectories of large number of bubbles, the void fraction distribution was predicted in this calculation. It concretely revealed that the lift force and the two-phase turbulence model were the most important parameters in determining the multi-dimensional void fraction distribution and the calculation should be considered by using the measured experimental data. (author)

Experimental study on two-dimensional film flow with local measurement methods

Energy Technology Data Exchange (ETDEWEB)

Yang, Jin-Hwa, E-mail: evo03@snu.ac.kr [Nuclear Thermal-Hydraulic Engineering Laboratory, Seoul National University, Gwanak 599, Gwanak-ro, Gwanak-gu, Seoul 151-742 (Korea, Republic of); Korea Atomic Energy Research Institute, 989-111, Daedeok-daero, Yuseong-gu, Daejeon 305-600 (Korea, Republic of); Cho, Hyoung-Kyu [Nuclear Thermal-Hydraulic Engineering Laboratory, Seoul National University, Gwanak 599, Gwanak-ro, Gwanak-gu, Seoul 151-742 (Korea, Republic of); Kim, Seok [Korea Atomic Energy Research Institute, 989-111, Daedeok-daero, Yuseong-gu, Daejeon 305-600 (Korea, Republic of); Euh, Dong-Jin, E-mail: djeuh@kaeri.re.kr [Korea Atomic Energy Research Institute, 989-111, Daedeok-daero, Yuseong-gu, Daejeon 305-600 (Korea, Republic of); Park, Goon-Cherl [Nuclear Thermal-Hydraulic Engineering Laboratory, Seoul National University, Gwanak 599, Gwanak-ro, Gwanak-gu, Seoul 151-742 (Korea, Republic of)

2015-12-01

Highlights: • An experimental study on the two-dimensional film flow with lateral air injection was performed. • The ultrasonic thickness gauge was used to measure the local liquid film thickness. • The depth-averaged PIV (Particle Image Velocimetry) method was applied to measure the local liquid film velocity. • The uncertainty of the depth-averaged PIV was quantified with a validation experiment. • Characteristics of two-dimensional film flow were classified following the four different flow patterns. - Abstract: In an accident condition of a nuclear reactor, multidimensional two-phase flows may occur in the reactor vessel downcomer and reactor core. Therefore, those have been regarded as important issues for an advanced thermal-hydraulic safety analysis. In particular, the multi-dimensional two-phase flow in the upper downcomer during the reflood phase of large break loss of coolant accident appears with an interaction between a downward liquid and a transverse gas flow, which determines the bypass flow rate of the emergency core coolant and subsequently, the reflood coolant flow rate. At present, some thermal-hydraulic analysis codes incorporate multidimensional modules for the nuclear reactor safety analysis. However, their prediction capability for the two-phase cross flow in the upper downcomer has not been validated sufficiently against experimental data based on local measurements. For this reason, an experimental study was carried out for the two-phase cross flow to clarify the hydraulic phenomenon and provide local measurement data for the validation of the computational tools. The experiment was performed in a 1/10 scale unfolded downcomer of Advanced Power Reactor 1400 (APR1400). Pitot tubes, a depth-averaged PIV method and ultrasonic thickness gauge were applied for local measurement of the air velocity, the liquid film velocity and the liquid film thickness, respectively. The uncertainty of the depth-averaged PIV method for the averaged

Experimental study on two-dimensional film flow with local measurement methods

International Nuclear Information System (INIS)

Yang, Jin-Hwa; Cho, Hyoung-Kyu; Kim, Seok; Euh, Dong-Jin; Park, Goon-Cherl

2015-01-01

Highlights: • An experimental study on the two-dimensional film flow with lateral air injection was performed. • The ultrasonic thickness gauge was used to measure the local liquid film thickness. • The depth-averaged PIV (Particle Image Velocimetry) method was applied to measure the local liquid film velocity. • The uncertainty of the depth-averaged PIV was quantified with a validation experiment. • Characteristics of two-dimensional film flow were classified following the four different flow patterns. - Abstract: In an accident condition of a nuclear reactor, multidimensional two-phase flows may occur in the reactor vessel downcomer and reactor core. Therefore, those have been regarded as important issues for an advanced thermal-hydraulic safety analysis. In particular, the multi-dimensional two-phase flow in the upper downcomer during the reflood phase of large break loss of coolant accident appears with an interaction between a downward liquid and a transverse gas flow, which determines the bypass flow rate of the emergency core coolant and subsequently, the reflood coolant flow rate. At present, some thermal-hydraulic analysis codes incorporate multidimensional modules for the nuclear reactor safety analysis. However, their prediction capability for the two-phase cross flow in the upper downcomer has not been validated sufficiently against experimental data based on local measurements. For this reason, an experimental study was carried out for the two-phase cross flow to clarify the hydraulic phenomenon and provide local measurement data for the validation of the computational tools. The experiment was performed in a 1/10 scale unfolded downcomer of Advanced Power Reactor 1400 (APR1400). Pitot tubes, a depth-averaged PIV method and ultrasonic thickness gauge were applied for local measurement of the air velocity, the liquid film velocity and the liquid film thickness, respectively. The uncertainty of the depth-averaged PIV method for the averaged

Multi-scale diffuse interface modeling of multi-component two-phase flow with partial miscibility

KAUST Repository

Kou, Jisheng; Sun, Shuyu

2016-01-01

In this paper, we introduce a diffuse interface model to simulate multi-component two-phase flow with partial miscibility based on a realistic equation of state (e.g. Peng-Robinson equation of state). Because of partial miscibility, thermodynamic

Hybrid-dimensional modelling of two-phase flow through fractured porous media with enhanced matrix fracture transmission conditions

Science.gov (United States)

Brenner, Konstantin; Hennicker, Julian; Masson, Roland; Samier, Pierre

2018-03-01

In this work, we extend, to two-phase flow, the single-phase Darcy flow model proposed in [26], [12] in which the (d - 1)-dimensional flow in the fractures is coupled with the d-dimensional flow in the matrix. Three types of so called hybrid-dimensional two-phase Darcy flow models are proposed. They all account for fractures acting either as drains or as barriers, since they allow pressure jumps at the matrix-fracture interfaces. The models also permit to treat gravity dominated flow as well as discontinuous capillary pressure at the material interfaces. The three models differ by their transmission conditions at matrix fracture interfaces: while the first model accounts for the nonlinear two-phase Darcy flux conservations, the second and third ones are based on the linear single phase Darcy flux conservations combined with different approximations of the mobilities. We adapt the Vertex Approximate Gradient (VAG) scheme to this problem, in order to account for anisotropy and heterogeneity aspects as well as for applicability on general meshes. Several test cases are presented to compare our hybrid-dimensional models to the generic equi-dimensional model, in which fractures have the same dimension as the matrix, leading to deep insight about the quality of the proposed reduced models.

Modeling two-phase flow in PEM fuel cell channels

Energy Technology Data Exchange (ETDEWEB)

Wang, Yun; Basu, Suman; Wang, Chao-Yang [Electrochemical Engine Center (ECEC), and Department of Mechanical and Nuclear Engineering, The Pennsylvania State University, University Park, PA 16802 (United States)

2008-05-01

This paper is concerned with the simultaneous flow of liquid water and gaseous reactants in mini-channels of a proton exchange membrane (PEM) fuel cell. Envisaging the mini-channels as structured and ordered porous media, we develop a continuum model of two-phase channel flow based on two-phase Darcy's law and the M{sup 2} formalism, which allow estimate of the parameters key to fuel cell operation such as overall pressure drop and liquid saturation profiles along the axial flow direction. Analytical solutions of liquid water saturation and species concentrations along the channel are derived to explore the dependences of these physical variables vital to cell performance on operating parameters such as flow stoichiometric ratio and relative humility. The two-phase channel model is further implemented for three-dimensional numerical simulations of two-phase, multi-component transport in a single fuel-cell channel. Three issues critical to optimizing channel design and mitigating channel flooding in PEM fuel cells are fully discussed: liquid water buildup towards the fuel cell outlet, saturation spike in the vicinity of flow cross-sectional heterogeneity, and two-phase pressure drop. Both the two-phase model and analytical solutions presented in this paper may be applicable to more general two-phase flow phenomena through mini- and micro-channels. (author)

Analysis of one-dimensional nonequilibrium two-phase flow using control volume method

International Nuclear Information System (INIS)

Minato, Akihiko; Naitoh, Masanori

1987-01-01

A one-dimensional numerical analysis model was developed for prediction of rapid flow transient behavior involving boiling. This model was based on six conservation equations of time averaged parameters of gas and liquid behavior. These equations were solved by using a control volume method with an explicit time integration. This model did not use staggered mesh scheme, which had been commonly used in two-phase flow analysis. Because void fraction and velocity of each phase were defined at the same location in the present model, effects of void fraction on phase velocity calculation were treated directly without interpolation. Though non-staggered mesh scheme was liable to cause numerical instability with zigzag pressure field, stability was achieved by employing the Godunov method. In order to verify the present analytical model, Edwards' pipe blow down and Zaloudek's initially subcooled critical two-phase flow experiments were analyzed. Stable solutions were obtained for rarefaction wave propagation with boiling and transient two-phase flow behavior in a broken pipe by using this model. (author)

Two benchmark cases for the trio two-phase flow module

Energy Technology Data Exchange (ETDEWEB)

Helton, D.; Hassan, Y. [Texas A and M University, Nuclear Engineering Dept., College Station, Texas (United States); Kumbaro, A. [CEA Saclay, 91 - Gif-sur-Yvette (France). Dept. de Mecanique et de Technologie

2001-07-01

This report presents a series of problems that were studied in order to assess the new implementations recently made to the two-phase flow module. Each problem is designed to give insight into a particular area of the code refinement. As such, each problem, and its corresponding results will be discussed individually, with comparisons made to experimental or analytical results whenever possible. TrioU is a thermal hydraulics program created by CEA. It is currently evolving in to a multi-dimensional, multi-fluid, multi-phase program. The purpose of TrioU is to provide a platform for testing of new numerical methods and physical models that are developed by the Nuclear Reactor Division of CEA. TrioU is written in an object-oriented programming language, and maintained by a version-controllable environment, for ease in parallelization and multiple-site development. (author)

Two benchmark cases for the trio two-phase flow module

International Nuclear Information System (INIS)

Helton, D.; Hassan, Y.; Kumbaro, A.

2001-01-01

This report presents a series of problems that were studied in order to assess the new implementations recently made to the two-phase flow module. Each problem is designed to give insight into a particular area of the code refinement. As such, each problem, and its corresponding results will be discussed individually, with comparisons made to experimental or analytical results whenever possible. TrioU is a thermal hydraulics program created by CEA. It is currently evolving in to a multi-dimensional, multi-fluid, multi-phase program. The purpose of TrioU is to provide a platform for testing of new numerical methods and physical models that are developed by the Nuclear Reactor Division of CEA. TrioU is written in an object-oriented programming language, and maintained by a version-controllable environment, for ease in parallelization and multiple-site development. (author)

Application of fast neutron radiography to three-dimensional visualization of steady two-phase flow in a rod bundle

CERN Document Server

Takenaka, N; Fujii, T; Mizubata, M; Yoshii, K

1999-01-01

Three-dimensional void fraction distribution of air-water two-phase flow in a 4x4 rod-bundle near a spacer was visualized by fast neutron radiography using a CT method. One-dimensional cross sectional averaged void fraction distribution was also calculated. The behaviors of low void fraction (thick water) two-phase flow in the rod bundle around the spacer were clearly visualized. It was shown that the void fraction distributions were visualized with a quality similar to those by thermal neutron radiography for low void fraction two-phase flow which is difficult to visualize by thermal neutron radiography. It is concluded that the fast neutron radiography is efficiently applicable to two-phase flow studies.

Surface Tension of Multi-phase Flow with Multiple Junctions Governed by the Variational Principle

International Nuclear Information System (INIS)

Matsutani, Shigeki; Nakano, Kota; Shinjo, Katsuhiko

2011-01-01

We explore a computational model of an incompressible fluid with a multi-phase field in three-dimensional Euclidean space. By investigating an incompressible fluid with a two-phase field geometrically, we reformulate the expression of the surface tension for the two-phase field found by Lafaurie et al. (J Comput Phys 113:134–147, 1994) as a variational problem related to an infinite dimensional Lie group, the volume-preserving diffeomorphism. The variational principle to the action integral with the surface energy reproduces their Euler equation of the two-phase field with the surface tension. Since the surface energy of multiple interfaces even with singularities is not difficult to be evaluated in general and the variational formulation works for every action integral, the new formulation enables us to extend their expression to that of a multi-phase (N-phase, N ≥ 2) flow and to obtain a novel Euler equation with the surface tension of the multi-phase field. The obtained Euler equation governs the equation for motion of the multi-phase field with different surface tension coefficients without any difficulties for the singularities at multiple junctions. In other words, we unify the theory of multi-phase fields which express low dimensional interface geometry and the theory of the incompressible fluid dynamics on the infinite dimensional geometry as a variational problem. We apply the equation to the contact angle problems at triple junctions. We computed the fluid dynamics for a two-phase field with a wall numerically and show the numerical computational results that for given surface tension coefficients, the contact angles are generated by the surface tension as results of balances of the kinematic energy and the surface energy.

Two-phase flow characteristics in BWRs

International Nuclear Information System (INIS)

Katono, Kenichi; Aoyama, Goro; Nagayoshi, Takuji; Yasuda, Kenichi; Nishida, Koji

2014-01-01

Reliable prediction of two-phase flow characteristics is important for safety and economy improvements of BWR plants. We have been developing two-phase flow measurement tools and techniques for BWR thermal hydraulic conditions, such as a 3D time-averaged X-ray CT system, an ultrasonic liquid film sensor and a wire-mesh sensor. We applied the developed items in experiments using the multi-purpose steam-water test facility known as HUSTLE, which can simulate two-phase thermal-hydraulic conditions in a BWR reactor pressure vessel, and we constructed a detailed instrumentation database. We validated a 3D two-phase flow simulator using the database and developed the reactor internal two-phase flow analysis system. (author)

Advanced numerical methods for three dimensional two-phase flow calculations in PWR

International Nuclear Information System (INIS)

Toumi, I.; Gallo, D.; Royer, E.

1997-01-01

This paper is devoted to new numerical methods developed for three dimensional two-phase flow calculations. These methods are finite volume numerical methods. They are based on an extension of Roe's approximate Riemann solver to define convective fluxes versus mean cell quantities. To go forward in time, a linearized conservative implicit integrating step is used, together with a Newton iterative method. We also present here some improvements performed to obtain a fully implicit solution method that provides fast running steady state calculations. This kind of numerical method, which is widely used for fluid dynamic calculations, is proved to be very efficient for the numerical solution to two-phase flow problems. This numerical method has been implemented for the three dimensional thermal-hydraulic code FLICA-4 which is mainly dedicated to core thermal-hydraulic transient and steady-state analysis. Hereafter, we will also find some results obtained for the EPR reactor running in a steady-state at 60% of nominal power with 3 pumps out of 4, and a thermal-hydraulic core analysis for a 1300 MW PWR at low flow steam-line-break conditions. (author)

Multi-phase flow monitoring with electrical impedance tomography using level set based method

International Nuclear Information System (INIS)

Liu, Dong; Khambampati, Anil Kumar; Kim, Sin; Kim, Kyung Youn

2015-01-01

Highlights: • LSM has been used for shape reconstruction to monitor multi-phase flow using EIT. • Multi-phase level set model for conductivity is represented by two level set functions. • LSM handles topological merging and breaking naturally during evolution process. • To reduce the computational time, a narrowband technique was applied. • Use of narrowband and optimization approach results in efficient and fast method. - Abstract: In this paper, a level set-based reconstruction scheme is applied to multi-phase flow monitoring using electrical impedance tomography (EIT). The proposed scheme involves applying a narrowband level set method to solve the inverse problem of finding the interface between the regions having different conductivity values. The multi-phase level set model for the conductivity distribution inside the domain is represented by two level set functions. The key principle of the level set-based method is to implicitly represent the shape of interface as the zero level set of higher dimensional function and then solve a set of partial differential equations. The level set-based scheme handles topological merging and breaking naturally during the evolution process. It also offers several advantages compared to traditional pixel-based approach. Level set-based method for multi-phase flow is tested with numerical and experimental data. It is found that level set-based method has better reconstruction performance when compared to pixel-based method

Two-dimensional phase fraction charts

International Nuclear Information System (INIS)

Morral, J.E.

1984-01-01

A phase fraction chart is a graphical representation of the amount of each phase present in a system as a function of temperature, composition or other variable. Examples are phase fraction versus temperature charts used to characterize specific alloys and as a teaching tool in elementary texts, and Schaeffler diagrams used to predict the amount of ferrite in stainless steel welds. Isothermal-transformation diagrams (TTT diagrams) are examples that give phase (or microconstituent) amount versus temperature and time. The purpose of this communication is to discuss the properties of two-dimensional phase fraction charts in more general terms than have been reported before. It is shown that they can represent multi-component, multiphase equilibria in a way which is easier to read and which contains more information than the isotherms and isopleths of multi-component phase diagrams

POST: a postprocessor computer code for producing three-dimensional movies of two-phase flow in a reactor vessel

International Nuclear Information System (INIS)

Taggart, K.A.; Liles, D.R.

1977-08-01

The development of the TRAC computer code for analysis of LOCAs in light-water reactors involves the use of a three-dimensional (r-theta-z), two-fluid hydrodynamics model to describe the two-phase flow of steam and water through the reactor vessel. One of the major problems involved in interpreting results from this code is the presentation of three-dimensional flow patterns. The purpose of the report is to present a partial solution to this data display problem. A first version of a code which produces three-dimensional movies of flow in the reactor vessel has been written and debugged. This code (POST) is used as a postprocessor in conjunction with a stand alone three-dimensional two-phase hydrodynamics code (CYLTF) which is a test bed for the three-dimensional algorithms to be used in TRAC

One-dimensional transient unequal velocity two-phase flow by the method of characteristics

International Nuclear Information System (INIS)

Rasouli, F.

1981-01-01

An understanding of two-phase flow is important when one is analyzing the accidental loss of coolant or when analyzing industrial processes. If a pipe in the steam generator of a nuclear reactor breaks, the flow will remain critical (or choked) for almost the entire blowdown. For this reason the knowledge of the two-phase maximum (critical) flow rate is important. A six-equation model--consisting of two continuity equations, two energy equations, a mixture momentum equation, and a constitutive relative velocity equation--is solved numerically by the method of characteristics for one-dimensional, transient, two-phase flow systems. The analysis is also extended to the special case of transient critical flow. The six-equation model is used to study the flow of a nonequilibrium sodium-argon system in a horizontal tube in which the nonequilibrium sodium-argon system in a horizontal tube in which the critical flow condition is at the entrance. A four-equation model is used to study the pressure-pulse propagation rate in an isothermal air-water system, and the results that are found are compared with the experimental data. Proper initial and boundary conditions are obtained for the blowdown problem. The energy and mass exchange relations are evaluated by comparing the model predictions with results of void-fraction and heat-transfer experiments. A simplified two-equation model is obtained for the special case of two incompressible phases. This model is used in the preliminary analysis of batch sedimentation. It is also used to predict the shock formation in the gas-solid fluidized bed

Thermo-Fluid Dynamics of Two-Phase Flow

CERN Document Server

Ishii, Mamrou

2011-01-01

"Thermo-fluid Dynamics of Two-Phase Flow, Second Edition" is focused on the fundamental physics of two-phase flow. The authors present the detailed theoretical foundation of multi-phase flow thermo-fluid dynamics as they apply to: Nuclear reactor transient and accident analysis; Energy systems; Power generation systems; Chemical reactors and process systems; Space propulsion; Transport processes. This edition features updates on two-phase flow formulation and constitutive equations and CFD simulation codes such as FLUENT and CFX, new coverage of the lift force model, which is of part

Pressure drop calculation using a one-dimensional mathematical model for two-phase flow through an orifice

DEFF Research Database (Denmark)

Petkov, K.P.; Puton, M; Madsen, Søren Peder

2014-01-01

are accounted for through both friction and acceleration as in a conventional formulation. However, in this analysis the acceleration term is both attributed geometrical effects through the area change and fluid dynamic effects through the expansion of the two-phase flow. The comparison of numerical...... is a one dimensional formulation in space and the equations incorporates the change in tubes and orifice diameter as formulated in (S. Madsen et.al., Dynamic Modeling of Phase Crossings in Two-Phase Flow, Communications in Computational Physics 12 (4), 1129-1147). The pressure changes in the flow...

Multiparticle imaging velocimetry measurements in two-phase flow

International Nuclear Information System (INIS)

Hassan, Y.A.

1998-01-01

The experimental flow visualization tool, Particle Image Velocimetry (PIV), is being extended to determine the velocity fields in two and three-dimensional, two-phase fluid flows. In the past few years, the technique has attracted quite a lot of interest. PIV enables fluid velocities across a region of a flow to be measured at a single instant in time in global domain. This instantaneous velocity profile of a given flow field is determined by digitally recording particle (microspheres or bubbles) images within the flow over multiple successive video frames and then conducting flow pattern identification and analysis of the data. This paper presents instantaneous velocity measurements in various two and three- dimensional, two-phase flow situations. (author)

Numerical analysis for two-dimensional compressible and two-phase flow fields of air-water in Eulerian grid framework

International Nuclear Information System (INIS)

Park, Chan Wook; Lee, Sung Su

2008-01-01

Two-phase compressible flow fields of air-water are investigated numerically in the fixed Eulerian grid framework. The phase interface is captured via volume fractions of ech phase. A way to model two phase compressible flows as a single phase one is found based on an equivalent equation of states of Tait's type for a multiphase cell. The equivalent single phase field is discretized using the Roe's approximate Riemann solver. Two approaches are tried to suppress the pressure oscillation phenomena at the phase interface, a passive advection of volume fraction and a direct pressure relaxation with the compressible form of volume fraction equation. The direct pressure equalizing method suppresses pressure oscillation successfully and generates sharp discontinuities, transmitting and reflecting acoustic waves naturally at the phase interface. In discretizing the compressible form of volume fraction equation, phase interfaces are geometrically reconstructed to minimize the numerical diffusion of volume fraction and relevant variables. The motion of a projectile in a water-filled tube which is fired by the release of highly pressurized air is simulated presuming the flow field as a two dimensional one, and several design factors affecting the projectile movement are investigated

The three-dimensional transient two-phase flow computer programme BACCHUS-3D/TP

International Nuclear Information System (INIS)

Bottoni, M.; Dorr, B.; Homann, C.

1992-04-01

The three-dimensional single-phase flow version of the BACCHUS code, which describes the thermal behaviour of a coolant in hexagonal bundle geometry, developed earlier, provided the basis for the development of the two-phase flow version documented in this report. A detailed description is given of the two-phase Slip Model (SM), and of the Homogeneous Equilibrium Model (HEM) as a subcase, which presents several improvements from both viewpoints of physical modelling and numerical treatment, with respect to usual models found in the literature. The most advanced Separated Phases Model (SPM) is then described in all analytical details necessary to fully understand its implementation in the code. Poblems related to the link between the two above models into an integrated code version are then discussed. The code provides an additional option for modelling of active or passive, permeable or impermeable blockages. This option is documented separately. New numerical methods for solving the algebraic systems of equations derived from the linearization of the fundamental equations have completely superseded previous ones and are explained in detail. Eventually a section is dedicated to an overview of the code verification, made over several years, which goes from steady state single-phase unheated bundle experiments up to fast transient two-phase flow experiments in electrically heated 37-pin bundles. (orig.) [de

Two-phase flow modeling in the rod bundle subchannel analysis

International Nuclear Information System (INIS)

Hisashi, Ninokata

2006-01-01

In order to practice a design-by-analysis of thermohydraulics design of BWR fuel rod bundles, the subchannel analysis would play a major role. There, the immediate concern is improvement in its predictive capability of CHF due in particular to the film dryout (boiling transition phenomena: BT) on the fuel rod surface. Constitutive equations in the subchannel analysis formulation are responsible for the quality of calculated results. The constitutive equations are a result of integration of the local and instantaneous description of two-phase flows over the subchannel control volume. In general, they are expressed in terms of subchannel-control-volume- as well as area-averaged two-phase flow state variables. In principle the information on local and instantaneous physical phenomena taking place inside subchannels must be counted for in the algebraic form of the equations on the basis of a more mechanistic modeling approach. They should include also influences of the multi-dimensional subchannel geometry and fluid material properties. Thermohydraulics phenomena of interests in this deed are: 1) vapor-liquid re-distribution by inter-subchannel exchanges due to the diversion cross flow, turbulent mixing and void drift, 2) liquid film behaviors, 3) transition of two-phase flow regimes, 4) droplet entrainment and deposition and 5) spacer-droplet interactions. These are considered to be five key factors in understanding the BT in BWR fuel rod bundles. In Japan, a university-industry consortium has been formed under the sponsorship of the Ministry of Economics, Trade and Industry. This paper describes an outline of the on-going project and, first, an outline of the current efforts is presented in developing a new two-fluid three field subchannel code NASCA being aimed at predicting onset of BT, and post BT phenomena in advanced BWR fuel rod bundles including those of the tight lattice configuration for a higher conversion. Then the current methodology adopted to improve

Two-phase flow modeling in the rod bundle subchannel analysis

International Nuclear Information System (INIS)

Hisashi, Ninokata

2004-01-01

Full text of publication follows:In order to practice a design-by-analysis of thermohydraulics design of BWR fuel rod bundles, the subchannel analysis would play a major role. There, the immediate concern is improvement in its predictive capability of CHF due in particular to the film dryout (boiling transition phenomena: BT) on the fuel rod surface. Constitutive equations in the subchannel analysis formulation are responsible for the quality of calculated results. The constitutive equations are a result of integration of the local and instantaneous description of two-phase flows over the subchannel control volume. In general, they are expressed in terms of subchannel-control-volume- as well as area-averaged two-phase flow state variables. In principle the information on local and instantaneous physical phenomena taking place inside subchannels must be counted for in the algebraic form of the equations on the basis of a more mechanistic modeling approach. They should include also influences of the multi-dimensional subchannel geometry and fluid material properties. Thermohydraulics phenomena of interests in this deed are: 1) vapor-liquid re-distribution by inter-subchannel exchanges due to the diversion cross flow, turbulent mixing and void drift, 2) liquid film behaviors, 3) transition of two-phase flow regimes, 4) droplet entrainment and deposition and 5) spacer-droplet interactions. These are considered to be five key factors in understanding the BT in BWR fuel rod bundles. In Japan, a university-industry consortium has been formed under the sponsorship of the Ministry of Economics, Trade and Industry. This paper describes an outline of the on-going project and, first, an outline of the current efforts is presented in developing a new two-fluid three field subchannel code NASCA being aimed at predicting onset of BT, and post BT phenomena in advanced BWR fuel rod bundles including those of the tight lattice configuration for a higher conversion. Then the current

Numerical method for three dimensional steady-state two-phase flow calculations

International Nuclear Information System (INIS)

Raymond, P.; Toumi, I.

1992-01-01

This paper presents the numerical scheme which was developed for the FLICA-4 computer code to calculate three dimensional steady state two phase flows. This computer code is devoted to steady state and transient thermal hydraulics analysis of nuclear reactor cores 1,3 . The first section briefly describes the FLICA-4 flow modelling. Then in order to introduce the numerical method for steady state computations, some details are given about the implicit numerical scheme based upon an approximate Riemann solver which was developed for calculation of flow transients. The third section deals with the numerical method for steady state computations, which is derived from this previous general scheme and its optimization. We give some numerical results for steady state calculations and comparisons on required CPU time and memory for various meshing and linear system solvers

A multi-scale network method for two-phase flow in porous media

Energy Technology Data Exchange (ETDEWEB)

Khayrat, Karim, E-mail: khayratk@ifd.mavt.ethz.ch; Jenny, Patrick

2017-08-01

Pore-network models of porous media are useful in the study of pore-scale flow in porous media. In order to extract macroscopic properties from flow simulations in pore-networks, it is crucial the networks are large enough to be considered representative elementary volumes. However, existing two-phase network flow solvers are limited to relatively small domains. For this purpose, a multi-scale pore-network (MSPN) method, which takes into account flow-rate effects and can simulate larger domains compared to existing methods, was developed. In our solution algorithm, a large pore network is partitioned into several smaller sub-networks. The algorithm to advance the fluid interfaces within each subnetwork consists of three steps. First, a global pressure problem on the network is solved approximately using the multiscale finite volume (MSFV) method. Next, the fluxes across the subnetworks are computed. Lastly, using fluxes as boundary conditions, a dynamic two-phase flow solver is used to advance the solution in time. Simulation results of drainage scenarios at different capillary numbers and unfavourable viscosity ratios are presented and used to validate the MSPN method against solutions obtained by an existing dynamic network flow solver.

A multi-scale network method for two-phase flow in porous media

International Nuclear Information System (INIS)

Khayrat, Karim; Jenny, Patrick

2017-01-01

Pore-network models of porous media are useful in the study of pore-scale flow in porous media. In order to extract macroscopic properties from flow simulations in pore-networks, it is crucial the networks are large enough to be considered representative elementary volumes. However, existing two-phase network flow solvers are limited to relatively small domains. For this purpose, a multi-scale pore-network (MSPN) method, which takes into account flow-rate effects and can simulate larger domains compared to existing methods, was developed. In our solution algorithm, a large pore network is partitioned into several smaller sub-networks. The algorithm to advance the fluid interfaces within each subnetwork consists of three steps. First, a global pressure problem on the network is solved approximately using the multiscale finite volume (MSFV) method. Next, the fluxes across the subnetworks are computed. Lastly, using fluxes as boundary conditions, a dynamic two-phase flow solver is used to advance the solution in time. Simulation results of drainage scenarios at different capillary numbers and unfavourable viscosity ratios are presented and used to validate the MSPN method against solutions obtained by an existing dynamic network flow solver.

Numerical method for two-phase flow discontinuity propagation calculation

International Nuclear Information System (INIS)

Toumi, I.; Raymond, P.

1989-01-01

In this paper, we present a class of numerical shock-capturing schemes for hyperbolic systems of conservation laws modelling two-phase flow. First, we solve the Riemann problem for a two-phase flow with unequal velocities. Then, we construct two approximate Riemann solvers: an one intermediate-state Riemann solver and a generalized Roe's approximate Riemann solver. We give some numerical results for one-dimensional shock-tube problems and for a standard two-phase flow heat addition problem involving two-phase flow instabilities

Numerical simulation of transient, adiabatic, two-dimensional two-phase flow using the two-fluid model

International Nuclear Information System (INIS)

Neves Conti, T. das.

1983-01-01

A numerical method is developed to simulate adiabatic, transient, two-dimensional two-phase flow. The two-fluid model is used to obtain the mass and momentum conservation equations. These are solved by an iterative algorithm emphoying a time-marching scheme. Based on the corrective procedure of Hirt and Harlow a poisson equation is derived for the pressure field. This equation is finite-differenced and solved by a suitable matrix inversion technique. In the absence of experiment results several numerical tests were made in order to chec accuracy, convergence and stability of the proposed method. Several tests were also performed to check whether the behavior of void fraction and phasic velocities conforms with previous observations. (Author) [pt

Thermodynamically consistent modeling and simulation of multi-component two-phase flow with partial miscibility

KAUST Repository

Kou, Jisheng; Sun, Shuyu

2017-01-01

A general diffuse interface model with a realistic equation of state (e.g. Peng-Robinson equation of state) is proposed to describe the multi-component two-phase fluid flow based on the principles of the NVT-based framework which is an attractive

Two-phase flow characteristics analysis code: MINCS

International Nuclear Information System (INIS)

Watanabe, Tadashi; Hirano, Masashi; Akimoto, Masayuki; Tanabe, Fumiya; Kohsaka, Atsuo.

1992-03-01

Two-phase flow characteristics analysis code: MINCS (Modularized and INtegrated Code System) has been developed to provide a computational tool for analyzing two-phase flow phenomena in one-dimensional ducts. In MINCS, nine types of two-phase flow models-from a basic two-fluid nonequilibrium (2V2T) model to a simple homogeneous equilibrium (1V1T) model-can be used under the same numerical solution method. The numerical technique is based on the implicit finite difference method to enhance the numerical stability. The code structure is highly modularized, so that new constitutive relations and correlations can be easily implemented into the code and hence evaluated. A flow pattern can be fixed regardless of flow conditions, and state equations or steam tables can be selected. It is, therefore, easy to calculate physical or numerical benchmark problems. (author)

Two-phase flow simulation of scour around a cylindrical pile

Science.gov (United States)

Nagel, T.; Chauchat, J.; Bonamy, C.; Liu, X.; Cheng, Z.; Hsu, T. J.

2017-12-01

Scour around structures is a major engineering issue that requires a detailed description of the flow field but also a consistent description of sediment transport processes that could not only be related to bed shear stress, like Shields parameter based sediment transport formula. In order to address this issue we used a multi-dimensional two-phase flow solver, sedFoam-2.0 (Chauchat et al., GMD 2017) implemented under the open-source CFD toolbox OpenFoam. Three-dimensional simulations have been performed on Roulund et al. (JFM 2005) configurations for clear-water and live bed cases. The k-omega model from Wilcox (AIAA Journal 2006) is used for the turbulent stress and the granular rheology μ(I) is used for the granular stress in the live bed case. The hydrodynamic is validated on the clear water case and the numerical results obtained for the live bed case provide a proof of concept that two-phase flow model is applicable to such problem with quantitative results for the prediction of scour depth upstream and downstream the cylinder at short timescales, up to 300s. Analyzing the simulation results in term of classical dimensionless sediment transport flux versus Shields parameter allows to get more insight into the fine scale sediment transport mechanisms involved in the scour process.

Thermodynamically consistent modeling and simulation of multi-component two-phase flow model with partial miscibility

KAUST Repository

Kou, Jisheng; Sun, Shuyu

2016-01-01

A general diffuse interface model with a realistic equation of state (e.g. Peng-Robinson equation of state) is proposed to describe the multi-component two-phase fluid flow based on the principles of the NVT-based framework which is a latest

Mechanistic multidimensional analysis of horizontal two-phase flows

International Nuclear Information System (INIS)

Tselishcheva, Elena A.; Antal, Steven P.; Podowski, Michael Z.

2010-01-01

The purpose of this paper is to discuss the results of analysis of two-phase flow in horizontal tubes. Two flow situations have been considered: gas/liquid flow in a long straight pipe, and similar flow conditions in a pipe with 90 deg. elbow. The theoretical approach utilizes a multifield modeling concept. A complete three-dimensional two-phase flow model has been implemented in a state-of-the-art computational multiphase fluid dynamics (CMFD) computer code, NPHASE. The overall model has been tested parametrically. Also, the results of NPHASE simulations have been compared against experimental data for a pipe with 90 deg. elbow.

Multi-dimensional modeling of gas-liquid two-phase flows. Application to the simulation of ascending bubble flows in vertical ducts

International Nuclear Information System (INIS)

Morel, Ch.

1997-01-01

The aim of this thesis is the 3-D modeling and numerical simulation of liquid/gas (water/vapor or water/air) two-phase flows in cooling circuits of nuclear power plants during normal and accidental situations. The development of a multidimensional dual-fluid model encounters two problems: the statistical effects of turbulence and the interface mass, momentum and energy transfers. The models developed in this study were introduced in the 3-D module of the CATHARE code developed by the CEA and the results were compared to experimental results available in the literature. The first chapter describes the equations of the local dual-fluid model for the 3-D description of two-phase flows. Closing relations adapted to dispersed flows with isothermal bubbles and without phase transformation are proposed and focus on the momentum transfer at the interfaces. The theoretical study of turbulence in the liquid phase of a bubble flow is modelled in chapter 2. Chapter 3 deals with the voluminal interface area used in the interface mass, momentum and energy transfers, and chapters 4 and 5 concern the application of the developed models to concrete situations. Chapter 4 describes in details the 3-D module of the CATHARE code while chapter 5 gives a comparison of numerical results obtained using the CATHARE code with other experimental results obtained at EdF. (J.S.)

Multidimensional analysis of developing two-phase flows in an ESBWR chimney with and without riser channels

International Nuclear Information System (INIS)

Murakawa, H.; Antal, S.P.; Lahey, R T.

2008-01-01

The object of this work was to simulate developing multidimensional velocity and void fraction distributions in bubbly and churn turbulent two-phase flows. An advanced Computational Multiphase Fluid Dynamics (CMFD) code, NPHASE, was used to perform three-dimensional, multi-field simulations of the developing phasic velocity and phase distributions in vertical adiabatic conduits. The NPHASE code employed a multi-field two-fluid model, in which, for churn turbulent flow, the vapor phase was divided into small and large, cap bubble fields. In addition, state-of-the-art interfacial area density and field-to-field mass transfer models were used for both the small and large, cap bubbles. In particular, the bubble breakup and coalescence processes were quantified using a two-group interfacial area density transport equation. This allowed the CMFD simulation of developing churn turbulent flows in an ESBWR with and without vertical riser channels in the chimney region above the core. Based on these simulations it was concluded that riser channels have little adverse effect on the induced natural circulation flow through the core and the stability characteristics of an ESBWR. (authors)

A Two-Temperature Open-Source CFD Model for Hypersonic Reacting Flows, Part Two: Multi-Dimensional Analysis †

Directory of Open Access Journals (Sweden)

Vincent Casseau

2016-12-01

Full Text Available hy2Foam is a newly-coded open-source two-temperature computational fluid dynamics (CFD solver that has previously been validated for zero-dimensional test cases. It aims at (1 giving open-source access to a state-of-the-art hypersonic CFD solver to students and researchers; and (2 providing a foundation for a future hybrid CFD-DSMC (direct simulation Monte Carlo code within the OpenFOAM framework. This paper focuses on the multi-dimensional verification of hy2Foam and firstly describes the different models implemented. In conjunction with employing the coupled vibration-dissociation-vibration (CVDV chemistry–vibration model, novel use is made of the quantum-kinetic (QK rates in a CFD solver. hy2Foam has been shown to produce results in good agreement with previously published data for a Mach 11 nitrogen flow over a blunted cone and with the dsmcFoam code for a Mach 20 cylinder flow for a binary reacting mixture. This latter case scenario provides a useful basis for other codes to compare against.

Application of Light Reflection Visualization for Measuring Organic-Liquid Saturation for Two-Phase Systems in Two-Dimensional Flow Cells.

Science.gov (United States)

DiFilippo, Erica L; Brusseau, Mark L

2011-11-01

A simple, noninvasive imaging technique was used to obtain in situ measurements of organic-liquid saturation in a two-phase system under dynamic conditions. Efficacy of the light reflection visualization (LRV) imaging method was tested through comparison of measured and known volumes of organic liquid for experiments conducted with a two-dimensional flow cell. Two sets of experiments were conducted, with source-zone configurations representing two archetypical residual-and-pool architectures. LRV measurements were collected during the injection of organic liquid and during a dissolution phase induced by water flushing. There was a strong correlation between measured and known organic-liquid volumes, with the LRV-measured values generally somewhat lower than the known volumes. Errors were greater for the system wherein organic liquid was present in multiple zones comprised of porous media of different permeabilities, and for conditions of multiphase flow. This method proved effective at determining organic-liquid distribution in a two-phase system using minimal specialized equipment.

Transient two-phase flow

International Nuclear Information System (INIS)

Hsu, Y.Y.

1974-01-01

The following papers related to two-phase flow are summarized: current assumptions made in two-phase flow modeling; two-phase unsteady blowdown from pipes, flow pattern in Laval nozzle and two-phase flow dynamics; dependence of radial heat and momentum diffusion; transient behavior of the liquid film around the expanding gas slug in a vertical tube; flooding phenomena in BWR fuel bundles; and transient effects in bubble two-phase flow. (U.S.)

Multi-Use Non-Intrusive Flow Characterization System (FCS), Phase I

Data.gov (United States)

National Aeronautics and Space Administration — The innovation is a Multi-Use Non-Intrusive Flow Characterization System (FCS) for densified, normal boiling point, and two-phase cryogenic flows, capable of...

Stability of interfacial waves in two-phase flows

Energy Technology Data Exchange (ETDEWEB)

Liu, W S [Ontario Hydro, Toronto, ON (Canada)

1996-12-31

The influence of the interfacial pressure and the flow distribution in the one-dimensional two-fluid model on the stability problems of interfacial waves is discussed. With a proper formulation of the interfacial pressure, the following two-phase phenomena can be predicted from the stability and stationary criteria of the interfacial waves: onset of slug flow, stationary hydraulic jump in a stratified flow, flooding in a vertical pipe, and the critical void fraction of a bubbly flow. It can be concluded that the interfacial pressure plays an important role in the interfacial wave propagation of the two-fluid model. The flow distribution parameter may enhance the flow stability range, but only plays a minor role in the two-phase characteristics. (author). 20 refs., 3 tabs., 4 figs.

Measurement of two phase flow properties using the nuclear reactor instruments

International Nuclear Information System (INIS)

Albrecht, R.W.; Washington Univ., Seattle; Crowe, R.D.; Dailey, D.J.; Kosaly, G.; Damborg, M.J.

1982-01-01

A procedure is introduced for characterizing one dimensional, two phase flow in terms of three properties; propagation, structure, and dynamics. It is shown that all of these properties can be measured by analyzing the response of the reactor neutron field to a two phase flow perturbation. Therefore, a nuclear reactor can be regarded as a two phase flow instrument. (author)

Multi-GPU hybrid programming accelerated three-dimensional phase-field model in binary alloy

Directory of Open Access Journals (Sweden)

Changsheng Zhu

2018-03-01

Full Text Available In the process of dendritic growth simulation, the computational efficiency and the problem scales have extremely important influence on simulation efficiency of three-dimensional phase-field model. Thus, seeking for high performance calculation method to improve the computational efficiency and to expand the problem scales has a great significance to the research of microstructure of the material. A high performance calculation method based on MPI+CUDA hybrid programming model is introduced. Multi-GPU is used to implement quantitative numerical simulations of three-dimensional phase-field model in binary alloy under the condition of multi-physical processes coupling. The acceleration effect of different GPU nodes on different calculation scales is explored. On the foundation of multi-GPU calculation model that has been introduced, two optimization schemes, Non-blocking communication optimization and overlap of MPI and GPU computing optimization, are proposed. The results of two optimization schemes and basic multi-GPU model are compared. The calculation results show that the use of multi-GPU calculation model can improve the computational efficiency of three-dimensional phase-field obviously, which is 13 times to single GPU, and the problem scales have been expanded to 8193. The feasibility of two optimization schemes is shown, and the overlap of MPI and GPU computing optimization has better performance, which is 1.7 times to basic multi-GPU model, when 21 GPUs are used.

Mathematical modeling of disperse two-phase flows

CERN Document Server

Morel, Christophe

2015-01-01

This book develops the theoretical foundations of disperse two-phase flows, which are characterized by the existence of bubbles, droplets or solid particles finely dispersed in a carrier fluid, which can be a liquid or a gas. Chapters clarify many difficult subjects, including modeling of the interfacial area concentration. Basic knowledge of the subjects treated in this book is essential to practitioners of Computational Fluid Dynamics for two-phase flows in a variety of industrial and environmental settings. The author provides a complete derivation of the basic equations, followed by more advanced subjects like turbulence equations for the two phases (continuous and disperse) and multi-size particulate flow modeling. As well as theoretical material, readers will discover chapters concerned with closure relations and numerical issues. Many physical models are presented, covering key subjects including heat and mass transfers between phases, interfacial forces and fluid particles coalescence and breakup, a...

The measurement of gas–liquid two-phase flows in a small diameter pipe using a dual-sensor multi-electrode conductance probe

International Nuclear Information System (INIS)

Zhai, Lu-Sheng; Bian, Peng; Han, Yun-Feng; Gao, Zhong-Ke; Jin, Ning-De

2016-01-01

We design a dual-sensor multi-electrode conductance probe to measure the flow parameters of gas–liquid two-phase flows in a vertical pipe with an inner diameter of 20 mm. The designed conductance probe consists of a phase volume fraction sensor (PVFS) and a cross-correlation velocity sensor (CCVS). Through inserting an insulated flow deflector in the central part of the pipe, the gas–liquid two-phase flows are forced to pass through an annual space. The multiple electrodes of the PVFS and the CCVS are flush-mounted on the inside of the pipe wall and the outside of the flow deflector, respectively. The geometry dimension of the PVFS is optimized based on the distribution characteristics of the sensor sensitivity field. In the flow loop test of vertical upward gas–liquid two-phase flows, the output signals from the dual-sensor multi-electrode conductance probe are collected by a data acquisition device from the National Instruments (NI) Corporation. The information transferring characteristics of local flow structures in the annular space are investigated using the transfer entropy theory. Additionally, the kinematic wave velocity is measured based on the drift velocity model to investigate the propagation behavior of the stable kinematic wave in the annular space. Finally, according to the motion characteristics of the gas–liquid two-phase flows, the drift velocity model based on the flow patterns is constructed to measure the individual phase flow rate with higher accuracy. (paper)

Numerical simulation of complex multi-dimensional two-phase flows in nuclear power plant coolant circuits by means of the best-estimate thermal-hydraulic code BAGIRA

International Nuclear Information System (INIS)

Kalinichenko, S.D.; Kroshilin, A.E.; Kroshilin, V.E.; Smirnov, A.V.

2009-01-01

Recent results are exposed, obtained by applying the best-estimate thermal hydraulic code BAGIRA for three-dimensional modeling complex two-phase flow dynamics inside the vessel of the horizontal steam generator PGV-1000 used in reactor units with VVER-1000. Spatial volumetric void fraction and velocity distributions are calculated and compared with available experimental data. (author)

A Three-Dimensional Model of Two-Phase Flows in a Porous Medium Accounting for Motion of the Liquid–Liquid Interface

DEFF Research Database (Denmark)

Shapiro, Alexander A.

2018-01-01

A new three-dimensional hydrodynamic model for unsteady two-phase flows in a porous medium, accounting for the motion of the interface between the flowing liquids, is developed. In a minimum number of interpretable geometrical assumptions, a complete system of macroscale flow equations is derived......, their expansion or contraction is also described, while rotation has been proven negligible. A detailed comparison with the previous studies for the two-phase flows accounting for propagation of the interface on micro- and macroscale has been carried out. A numerical algorithm has been developed allowing...

Study of two-phase flow redistribution between two passes of a heat exchanger

International Nuclear Information System (INIS)

Mendes de Moura, L.F.

1989-04-01

The object of the present thesis deals with the study of two-phase flow redistribution between two passes of a heat exchanger. Mass flow rate measurements of each component performed at each channel outlet of the second pass allowed us to determine the influence of mass flow, gas quality, flow direction (upward or downward) and common header geometry upon flow redistribution. Local void fraction inside common header was measured with an optical probe. A two-dimensional two-phase flow computational code was developed from a two-fluid model. Modelling of interfacial momentum transfer was used in order to take into account twp-phase flow patterns in common headers. Numerical simulation results show qualitative agreement with experimental results. Present theoretical model limitations are analysed and future improvements are proposed [fr

Geometric analysis of the solutions of two-phase flows: two-fluid model

International Nuclear Information System (INIS)

Kestin, J.; Zeng, D.L.

1984-01-01

This report contains a lightly edited draft of a study of the two-fluid model in two-phase flow. The motivation for the study stems from the authors' conviction that the construction of a computer code for any model should be preceded by a geometrical analysis of the pattern of trajectories in the phase space appropriate for the model. Such a study greatly facilitates the understanding of the phenomenon of choking and anticipates the computational difficulties which arise from the existence of singularities. The report contains a derivation of the six conservation equations of the model which includes a consideration of the simplifications imposed on a one-dimensional treatment by the presence of boundary layers at the wall and between the phases. The model is restricted to one-dimensional adiabatic flows of a single substance present in two phases, but thermodynamic equilibrium between the phases is not assumed. The role of closure conditions is defined but no specific closure conditions, or explicit equations of state, are introduced

A Simple Free Surface Tracking Model for Multi-dimensional Two-Fluid Approaches

International Nuclear Information System (INIS)

Lee, Seungjun; Yoon, Han Young

2014-01-01

The development in two-phase experiments devoted to find unknown phenomenological relationships modified conventional flow pattern maps into a sophisticated one and even extended to the multi-dimensional usage. However, for a system including a large void fraction gradient, such as a pool with the free surface, the flow patterns varies spatially throughout small number of cells and sometimes results in an unstable and unrealistic prediction of flows at the large gradient void fraction cells. Then, the numerical stability problem arising from the free surface is the major interest in the analyses of a passive cooling pool convecting the decay heat naturally, which has become a design issue to increase the safety level of nuclear reactors recently. In this research, a new and simple free surface tracking method combined with a simplified topology map is presented. The method modified the interfacial drag coefficient only for the cells defined as the free surface. The performance is shown by comparing the natural convection analysis of a small scale pool with respect to single- and two-phase condition. A simple free surface tracking model with a simplified topology map is developed

Two-phase flow models in unbounded two-phase critical flows

International Nuclear Information System (INIS)

Celata, G.P.; Cumo, M.; Farello, G.E.

1985-01-01

With reference to a Loss-of-Coolant Accident in Light Water Reactors, an analysis of the unbounded two-phase critical flow (i.e. the issuing two-phase jet) has been accomplished. Considering jets external shape, obtained by means of photographic pictures; pressure profiles inside the jet, obtained by means of a movable ''Pitot;'' and jet phases distribution information, obtained by means of X-rays pictures; a characterization of the flow pattern in the unbounded region of a two-phase critical flow is given. Jets X-ray pictures show the existence of a central high density ''core'' gradually evaporating all around, which gives place to a characteristic ''dartflow'' the length of which depends on stagnation thermodynamic conditions

One-Dimensional, Two-Phase Flow Modeling Toward Interpreting Motor Slag Expulsion Phenomena

Science.gov (United States)

Kibbey, Timothy P.

2012-01-01

Aluminum oxide slag accumulation and expulsion was previously shown to be a player in various solid rocket motor phenomena, including the Space Shuttle's Reusable Solid Rocket Motor (RSRM) pressure perturbation, or "blip," and phantom moment. In the latter case, such un ]commanded side accelerations near the end of burn have also been identified in several other motor systems. However, efforts to estimate the mass expelled during a given event have come up short. Either bulk calculations are performed without enough physics present, or multiphase, multidimensional Computational Fluid Dynamic analyses are performed that give a snapshot in time and space but do not always aid in grasping the general principle. One ]dimensional, two ]phase compressible flow calculations yield an analytical result for nozzle flow under certain assumptions. This can be carried further to relate the bulk motor parameters of pressure, thrust, and mass flow rate under the different exhaust conditions driven by the addition of condensed phase mass flow. An unknown parameter is correlated to airflow testing with water injection where mass flow rates and pressure are known. Comparison is also made to full ]scale static test motor data where thrust and pressure changes are known and similar behavior is shown. The end goal is to be able to include the accumulation and flow of slag in internal ballistics predictions. This will allow better prediction of the tailoff when much slag is ejected and of mass retained versus time, believed to be a contributor to the widely-observed "flight knockdown" parameter.

Macroscopic balance equations for two-phase flow models

International Nuclear Information System (INIS)

Hughes, E.D.

1979-01-01

The macroscopic, or overall, balance equations of mass, momentum, and energy are derived for a two-fluid model of two-phase flows in complex geometries. These equations provide a base for investigating methods of incorporating improved analysis methods into computer programs, such as RETRAN, which are used for transient and steady-state thermal-hydraulic analyses of nuclear steam supply systems. The equations are derived in a very general manner so that three-dimensional, compressible flows can be analysed. The equations obtained supplement the various partial differential equation two-fluid models of two-phase flow which have recently appeared in the literature. The primary objective of the investigation is the macroscopic balance equations. (Auth.)

Multi-level adaptive simulation of transient two-phase flow in heterogeneous porous media

KAUST Repository

Chueh, C.C.

2010-10-01

An implicit pressure and explicit saturation (IMPES) finite element method (FEM) incorporating a multi-level shock-type adaptive refinement technique is presented and applied to investigate transient two-phase flow in porous media. Local adaptive mesh refinement is implemented seamlessly with state-of-the-art artificial diffusion stabilization allowing simulations that achieve both high resolution and high accuracy. Two benchmark problems, modelling a single crack and a random porous medium, are used to demonstrate the robustness of the method and illustrate the capabilities of the adaptive refinement technique in resolving the saturation field and the complex interaction (transport phenomena) between two fluids in heterogeneous media. © 2010 Elsevier Ltd.

Modelling of two-phase flow based on separation of the flow according to velocity

International Nuclear Information System (INIS)

Narumo, T.

1997-01-01

The thesis concentrates on the development work of a physical one-dimensional two-fluid model that is based on Separation of the Flow According to Velocity (SFAV). The conventional way to model one-dimensional two-phase flow is to derive conservation equations for mass, momentum and energy over the regions occupied by the phases. In the SFAV approach, the two-phase mixture is divided into two subflows, with as distinct average velocities as possible, and momentum conservation equations are derived over their domains. Mass and energy conservation are treated equally with the conventional model because they are distributed very accurately according to the phases, but momentum fluctuations follow better the flow velocity. Submodels for non-uniform transverse profile of velocity and density, slip between the phases within each subflow and turbulence between the subflows have been derived. The model system is hyperbolic in any sensible flow conditions over the whole range of void fraction. Thus, it can be solved with accurate numerical methods utilizing the characteristics. The characteristics agree well with the used experimental data on two-phase flow wave phenomena Furthermore, the characteristics of the SFAV model are as well in accordance with their physical counterparts as of the best virtual-mass models that are typically optimized for special flow regimes like bubbly flow. The SFAV model has proved to be applicable in describing two-phase flow physically correctly because both the dynamics and steady-state behaviour of the model has been considered and found to agree well with experimental data This makes the SFAV model especially suitable for the calculation of fast transients, taking place in versatile form e.g. in nuclear reactors

One-dimensional two-phase thermal hydraulics (ENSTA course)

International Nuclear Information System (INIS)

Olive, J.

1995-11-01

This course is part of the ENSTA 3rd year thermal hydraulics program (nuclear power option). Its purpose is to provide the theoretical basis and main physical notions pertaining to two-phase flow, mainly focussed on water-steam flows. The introduction describes the physical specificities of these flows, emphasizing their complexity. The mathematical bases are then presented (partial derivative equations), leading to a one-dimensional type, simplified description. Balances drawn up for a pipe length volume are used to introduce the mass conservation. motion and energy equations for each phase. Various postulates used to simplify two-phase models are presented, culminating in homogeneous model definitions and equations, several common examples of which are given. The model is then applied to the calculation of pressure drops in two-phase flows. This involves presenting the models most frequently used to represent pressure drops by friction or due to pipe irregularities, without giving details (numerical values of parameters). This chapter terminates with a brief description of static and dynamic instabilities in two-phase flows. Finally, heat transfer conditions frequently encountered in liquid-steam flows are described, still in the context of a 1D model. This chapter notably includes reference to under-saturated boiling conditions and the various forms of DNB. The empirical heat transfer laws are not discussed in detail. Additional material is appended, some of which is in the form of corrected exercises. (author). 6 appends

Linearly decoupled energy-stable numerical methods for multi-component two-phase compressible flow

KAUST Repository

Kou, Jisheng

2017-12-06

In this paper, for the first time we propose two linear, decoupled, energy-stable numerical schemes for multi-component two-phase compressible flow with a realistic equation of state (e.g. Peng-Robinson equation of state). The methods are constructed based on the scalar auxiliary variable (SAV) approaches for Helmholtz free energy and the intermediate velocities that are designed to decouple the tight relationship between velocity and molar densities. The intermediate velocities are also involved in the discrete momentum equation to ensure a consistency relationship with the mass balance equations. Moreover, we propose a component-wise SAV approach for a multi-component fluid, which requires solving a sequence of linear, separate mass balance equations. We prove that the methods have the unconditional energy-dissipation feature. Numerical results are presented to verify the effectiveness of the proposed methods.

Analyses of liquid-gas two-phase flow in fermentation tanks

International Nuclear Information System (INIS)

Toi, Takashi; Serizawa, Akimi; Takahashi, Osamu; Kawara, Zensaku; Gofuku, Akio; Kataoka, Isao.

1993-01-01

The understanding of two-phase flow is one of the important problems for both design and safety analyses of various engineering systems. For example, the flow conditions in beer fermentation tanks have an influence on the quality of production and productivity of tank. In this study, a two-dimensional numerical calculation code based on the one-pressure two-fluid model is developed to understand the circulation structure of low quality liquid-gas two-phase flows induced by bubble plume in a tank. (author)

Analysis of water hammer in two-component two-phase flows

International Nuclear Information System (INIS)

Warde, H.; Marzouk, E.; Ibrahim, S.

1989-01-01

The water hammer phenomena caused by a sudden valve closure in air-water two-phase flows must be clarified for the safety analysis of LOCA in reactors and further for the safety of boilers, chemical plants, pipe transport of fluids such as petroleum and natural gas. In the present work water hammer phenomena caused by sudden valve closure in two-component two-phase flows are investigated theoretically and experimentally. The phenomena are more complicated than in single phase-flows due to the fact of the presence of compressible component. Basic partial differential equations based on a one-dimensional homogeneous flow model are solved by the method of characteristic. The analysis is extended to include friction in a two-phase mixture depending on the local flow pattern. The profiles of the pressure transients, the propagation velocity of pressure waves and the effect of valve closure on the transient pressure are found. Different two-phase flow pattern and frictional pressure drop correlations were used including Baker, Chesholm and Beggs and Bril correlations. The effect of the flow pattern on the characteristic of wave propagation is discussed primarily to indicate the effect of void fraction on the velocity of wave propagation and on the attenuation of pressure waves. Transient pressure in the mixture were recorded at different air void fractions, rates of uniform valve closure and liquid flow velocities with the aid of pressure transducers, transient wave form recorders interfaced with an on-line pc computer. The results are compared with computation, and good agreement was obtained within experimental accuracy

K-FIX: a computer program for transient, two-dimensional, two-fluid flow

International Nuclear Information System (INIS)

Rivard, W.C.; Torrey, M.D.

1976-11-01

The transient dynamics of two-dimensional, two-phase flow with interfacial exchange are calculated at all flow speeds using the K-FIX program. Each phase is described in terms of its own density, velocity, and temperature. The six field equations for the two phases couple through mass, momentum, and energy exchange. The equations are solved using an Eulerian finite difference technique that implicitly couples the rates of phase transitions, momentum, and energy exchange to determination of the pressure, density, and velocity fields. The implicit solution is accomplished iteratively without linearizing the equations, thus eliminating the need for numerous derivative terms. K-FIX is written in a highly modular form to be easily adaptable to a variety of problems. It is applied to growth of an isolated steam bubble in a superheated water pool

Waves in separated two-phase flow

International Nuclear Information System (INIS)

Pols, R.M.

1998-06-01

This dissertation presents an integral approach to the modelling of co-current flow of liquid and gas for a class of non-linear wave problems. Typically the liquid phase and the gas phase are decoupled and the liquid is depth averaged. The resulting non-linear shallow water equations are solved to predict the behaviour of the finite amplitude waves. The integral approach is applied to the modelling of two-dimensional waves in a horizontal and slightly inclined rectangular channel, two-dimensional waves in a vertical pipe and three-dimensional waves in a horizontal tube. For flow in a horizontal or slightly inclined channel the liquid is driven by the interfacial shear from the gas phase and the surface is subject to extensive wave action. For thin liquid films the pressure in the liquid may be taken as hydrostatic and gravity acts as a restoring force on the liquid. Roll wave solutions to the non-linear shallow water equations are sought corresponding to an interfacial shear stress dependent on the liquid film height. Wave solutions are shown to exist but only for parameters within a defined range dependent on the channel inclination, interfacial roughness and linear dependence on the liquid film height of the shear stresses. Such solutions are discontinuous and are pieced together by a jump where mass and momentum are conserved. The model calculations on wave height and wave velocity are compared with experimental data. The essentially two-dimensional analysis developed for stratified horizontal flow can be extended to quasi three-dimensional flow in the case of shallow liquid depth for a circular pipe. In this case the liquid depth changes with circumferential position and consequently modifies the interfacial shear exerted on the liquid surface creating a wave spreading mechanism alongside changing the wave profile across the pipe. The wave spreading mechanism supposes a wave moving in axial direction at a velocity faster than the liquid thereby sweeping liquid

Measurement of multi-dimensional flow structure for flow boiling in a tube

International Nuclear Information System (INIS)

Adachi, Yu; Ito, Daisuke; Saito, Yasushi

2014-01-01

With an aim of the measurement of multi-dimensional flow structure of in-tube boiling two-phase flow, the authors built their own wire mesh measurement system based on electrical conductivity measurement, and examined the relationship between the electrical conductivity obtained by the wire mesh sensor and the void fraction. In addition, the authors measured the void fraction using neutron radiography, and compared the result with the measured value using the wire mesh sensor. From the comparison with neutron radiography, it was found that the new method underestimated the void fraction in the flow in the vicinity of the void fraction of 0.2-0.5, similarly to the conventional result. In addition, since the wire mesh sensor cannot measure dispersed droplets, it tends to overestimate the void fraction in the high void fraction region, such as churn flow accompanied by droplet generation. In the electrical conductivity wire-mesh sensor method, it is necessary to correctly take into account the effect of liquid film or droplets. The authors also built a measurement system based on the capacitance wire mesh sensor method using the difference in dielectric constant, performed the confirmation of transmission and reception signals using deionized water as a medium, and showed the validity of the system. As for the dispersed droplets, the capacitance method has a potential to be able to measure them. (A.O.)

A two-dimensional, two-phase mass transport model for liquid-feed DMFCs

International Nuclear Information System (INIS)

Yang, W.W.; Zhao, T.S.

2007-01-01

A two-dimensional, isothermal two-phase mass transport model for a liquid-feed direct methanol fuel cell (DMFC) is presented in this paper. The two-phase mass transport in the anode and cathode porous regions is formulated based on the classical multiphase flow in porous media without invoking the assumption of constant gas pressure in the unsaturated porous medium flow theory. The two-phase flow behavior in the anode flow channel is modeled by utilizing the drift-flux model, while in the cathode flow channel the homogeneous mist-flow model is used. In addition, a micro-agglomerate model is developed for the cathode catalyst layer. The model also accounts for the effects of both methanol and water crossover through the membrane. The comprehensive model formed by integrating those in the different regions is solved numerically using a home-written computer code and validated against the experimental data in the literature. The model is then used to investigate the effects of various operating and structural parameters, such as methanol concentration, anode flow rate, porosities of both anode and cathode electrodes, the rate of methanol crossover, and the agglomerate size, on cell performance

Present status of numerical analysis on transient two-phase flow

International Nuclear Information System (INIS)

Akimoto, Masayuki; Hirano, Masashi; Nariai, Hideki.

1987-01-01

The Special Committee for Numerical Analysis of Thermal Flow has recently been established under the Japan Atomic Energy Association. Here, some methods currently used for numerical analysis of transient two-phase flow are described citing some information given in the first report of the above-mentioned committee. Many analytical models for transient two-phase flow have been proposed, each of which is designed to describe a flow by using differential equations associated with conservation of mass, momentum and energy in a continuous two-phase flow system together with constructive equations that represent transportation of mass, momentum and energy though a gas-liquid interface or between a liquid flow and the channel wall. The author has developed an analysis code, called MINCS, that serves for systematic examination of conservation equation and constructive equations for two-phase flow models. A one-dimensional, non-equilibrium two-liquid flow model that is used as the basic model for the code is described. Actual procedures for numerical analysis is shown and some problems concerning transient two-phase analysis are described. (Nogami, K.)

Modelling of two-phase flow based on separation of the flow according to velocity

Energy Technology Data Exchange (ETDEWEB)

Narumo, T. [VTT Energy, Espoo (Finland). Nuclear Energy

1997-12-31

The thesis concentrates on the development work of a physical one-dimensional two-fluid model that is based on Separation of the Flow According to Velocity (SFAV). The conventional way to model one-dimensional two-phase flow is to derive conservation equations for mass, momentum and energy over the regions occupied by the phases. In the SFAV approach, the two-phase mixture is divided into two subflows, with as distinct average velocities as possible, and momentum conservation equations are derived over their domains. Mass and energy conservation are treated equally with the conventional model because they are distributed very accurately according to the phases, but momentum fluctuations follow better the flow velocity. Submodels for non-uniform transverse profile of velocity and density, slip between the phases within each subflow and turbulence between the subflows have been derived. The model system is hyperbolic in any sensible flow conditions over the whole range of void fraction. Thus, it can be solved with accurate numerical methods utilizing the characteristics. The characteristics agree well with the used experimental data on two-phase flow wave phenomena Furthermore, the characteristics of the SFAV model are as well in accordance with their physical counterparts as of the best virtual-mass models that are typically optimized for special flow regimes like bubbly flow. The SFAV model has proved to be applicable in describing two-phase flow physically correctly because both the dynamics and steady-state behaviour of the model has been considered and found to agree well with experimental data This makes the SFAV model especially suitable for the calculation of fast transients, taking place in versatile form e.g. in nuclear reactors. 45 refs. The thesis includes also five previous publications by author.

Characterizing the correlations between local phase fractions of gas–liquid two-phase flow with wire-mesh sensor

Science.gov (United States)

Liu, W. L.; Dong, F.

2016-01-01

Understanding of flow patterns and their transitions is significant to uncover the flow mechanics of two-phase flow. The local phase distribution and its fluctuations contain rich information regarding the flow structures. A wire-mesh sensor (WMS) was used to study the local phase fluctuations of horizontal gas–liquid two-phase flow, which was verified through comparing the reconstructed three-dimensional flow structure with photographs taken during the experiments. Each crossing point of the WMS is treated as a node, so the measurement on each node is the phase fraction in this local area. An undirected and unweighted flow pattern network was established based on connections that are formed by cross-correlating the time series of each node under different flow patterns. The structure of the flow pattern network reveals the relationship of the phase fluctuations at each node during flow pattern transition, which is then quantified by introducing the topological index of the complex network. The proposed analysis method using the WMS not only provides three-dimensional visualizations of the gas–liquid two-phase flow, but is also a thorough analysis for the structure of flow patterns and the characteristics of flow pattern transition. This article is part of the themed issue ‘Supersensing through industrial process tomography’. PMID:27185959

Characterizing the correlations between local phase fractions of gas-liquid two-phase flow with wire-mesh sensor.

Science.gov (United States)

Tan, C; Liu, W L; Dong, F

2016-06-28

Understanding of flow patterns and their transitions is significant to uncover the flow mechanics of two-phase flow. The local phase distribution and its fluctuations contain rich information regarding the flow structures. A wire-mesh sensor (WMS) was used to study the local phase fluctuations of horizontal gas-liquid two-phase flow, which was verified through comparing the reconstructed three-dimensional flow structure with photographs taken during the experiments. Each crossing point of the WMS is treated as a node, so the measurement on each node is the phase fraction in this local area. An undirected and unweighted flow pattern network was established based on connections that are formed by cross-correlating the time series of each node under different flow patterns. The structure of the flow pattern network reveals the relationship of the phase fluctuations at each node during flow pattern transition, which is then quantified by introducing the topological index of the complex network. The proposed analysis method using the WMS not only provides three-dimensional visualizations of the gas-liquid two-phase flow, but is also a thorough analysis for the structure of flow patterns and the characteristics of flow pattern transition. This article is part of the themed issue 'Supersensing through industrial process tomography'. © 2016 The Author(s).

Robust second-order scheme for multi-phase flow computations

Science.gov (United States)

Shahbazi, Khosro

2017-06-01

A robust high-order scheme for the multi-phase flow computations featuring jumps and discontinuities due to shock waves and phase interfaces is presented. The scheme is based on high-order weighted-essentially non-oscillatory (WENO) finite volume schemes and high-order limiters to ensure the maximum principle or positivity of the various field variables including the density, pressure, and order parameters identifying each phase. The two-phase flow model considered besides the Euler equations of gas dynamics consists of advection of two parameters of the stiffened-gas equation of states, characterizing each phase. The design of the high-order limiter is guided by the findings of Zhang and Shu (2011) [36], and is based on limiting the quadrature values of the density, pressure and order parameters reconstructed using a high-order WENO scheme. The proof of positivity-preserving and accuracy is given, and the convergence and the robustness of the scheme are illustrated using the smooth isentropic vortex problem with very small density and pressure. The effectiveness and robustness of the scheme in computing the challenging problem of shock wave interaction with a cluster of tightly packed air or helium bubbles placed in a body of liquid water is also demonstrated. The superior performance of the high-order schemes over the first-order Lax-Friedrichs scheme for computations of shock-bubble interaction is also shown. The scheme is implemented in two-dimensional space on parallel computers using message passing interface (MPI). The proposed scheme with limiter features approximately 50% higher number of inter-processor message communications compared to the corresponding scheme without limiter, but with only 10% higher total CPU time. The scheme is provably second-order accurate in regions requiring positivity enforcement and higher order in the rest of domain.

A method for identifying gas-liquid two-phase flow patterns on the basis of wavelet packet multi-scale information entropy and HMM

International Nuclear Information System (INIS)

Zhou Yunlong; Zhang Xueqing; Gao Yunpeng; Cheng Yue

2009-01-01

For studying flow regimes of gas/liquid two-phase in a vertical upward pipe, the conductance fluctuation information of four typical flow regimes was collected by a measuring the system with self-made multiple conductivity probes. Owing to the non-stationarity of conductance fluctuation signals of gas-liquid two-phase flow, a kind of' flow regime identification method based on wavelet packet Multi-scale Information Entropy and Hidden Markov Model (HMM) was put forward. First of all, the collected conductance fluctuation signals were decomposed into eight different frequency bands signals. Secondly, the wavelet packet multi-scale information entropy of different frequency bands signals were regarded as the input characteristic vectors of all states HMM which had been trained. In the end the regime identification of' the gas-liquid two-phase flow could be performed. The study showed that the method that HMM was applied to identify the flow regime was superior to the one that BP neural network was used, and the results proved that the method was efficient and feasible. (authors)

Two-dimensional numerical experiments with DRIX-2D on two-phase-water-flows referring to the HDR-blowdown-experiments

International Nuclear Information System (INIS)

Moesinger, H.

1979-08-01

The computer program DRIX-2D has been developed from SOLA-DF. The essential elements of the program structure are described. In order to verify DRIX-2D an Edwards-Blowdown-Experiment is calculated and other numerical results are compared with steady state experiments and models. Numerical experiments on transient two-phase flow, occurring in the broken pipe of a PWR in the case of a hypothetic LOCA, are performed. The essential results of the two-dimensional calculations are: 1. The appearance of a radial profile of void-fraction, velocity, sound speed and mass flow-rate inside the blowdown nozzle. The reason for this is the flow contraction at the nozzle inlet leading to more vapour production in the vicinity of the pipe wall. 2. A comparison between modelling in axisymmetric and Cartesian coordinates and calculations with and without the core barrel show the following: a) The three-dimensional flow pattern at the nozzle inlet is poorly described using Cartesian coordinates. In consequence a considerable difference in pressure history results. b) The core barrel alters the reflection behaviour of the pressure waves oscillating in the blowdown-nozzle. Therefore, the core barrel should be modelled as a wall normal to the nozzle axis. (orig./HP) [de

Two-phase flow in porous media: power-law scaling of effective permeability

Energy Technology Data Exchange (ETDEWEB)

Groeva, Morten; Hansen, Alex, E-mail: Morten.Grova@ntnu.no, E-mail: Alex.Hansen@ntnu.no [Department of Physics, NTNU, NO-7491 Trondheim (Norway)

2011-09-15

A recent experiment has reported power-law scaling of effective permeability of two-phase flow with respect to capillary number for a two-dimensional model porous medium. In this paper, we consider the simultaneous flow of two phases through a porous medium under steady-state conditions, fixed total flow-rate and saturation, using a two-dimensional network simulator. We obtain power-law exponents for the scaling of effective permeability with respect to capillary number. The simulations are performed both for viscosity matched fluids and for a high viscosity ratio resembling that of air and water. Good power-law behaviour is found for both cases. Different exponents are found, depending on saturation.

Hydrodynamics of single- and two-phase flow in inclined rod arrays

International Nuclear Information System (INIS)

Ebeling-Koning, D.B.; Todreas, N.E.

1983-09-01

Required inputs for thermal-hydraulic codes are constitutive relations for fluid-solid flow resistance, in single-phase flow, and interfacial momentum exchange (relative phase motion), in two-phase flow. An inclined rod array air-water experiment was constructed to study the hydrodynamics of multidimensional porous medium flow in rod arrays. Velocities, pressures, and bubble distributions were measured in square rod arrays of P/d = 1.5, at 0, 30, 45, and 90 degree inclinations to the vertical flow direction. Constitutive models for single-phase flow resistance are reviewed, new comprehensive models developed, and an assessment with previously published and new data made. The principle of superimposing one-dimensional correlations proves successful for turbulent single-phase inclined flow. For bubbly two-phase incline flow a new flow separation phenomena was observed and modeled. A two-region liquid velocity model is developed to explain the experimentally observed phenomena. Fundamental data for bubbles rising in rod arrays were also taken

Characteristics and stability analyses of transient one-dimensional two-phase flow equations and their finite difference approximations

International Nuclear Information System (INIS)

Lyczkowski, R.W.; Gidaspow, D.; Solbrig, C.W.; Hughes, E.D.

1975-01-01

Equation systems describing one-dimensional, transient, two-phase flow with separate continuity, momentum, and energy equations for each phase are classified by use of the method of characteristics. Little attempt is made to justify the physics of these equations. Many of the equation systems possess complex-valued characteristics and hence, according to well-known mathematical theorems, are not well-posed as initial-value problems (IVPs). Real-valued characteristics are necessary but not sufficient to insure well-posedness. In the absence of lower order source or sink terms (potential type flows), which can affect the well-posedness of IVPs, the complex characteristics associated with these two-phase flow equations imply unbounded exponential growth for disturbances of all wavelengths. Analytical and numerical examples show that the ill-posedness of IVPs for the two-phase flow partial differential equations which possess complex characteristics produce unstable numerical schemes. These unstable numerical schemes can produce apparently stable and even accurate results if the growth rate resulting from the complex characteristics remains small throughout the time span of the numerical experiment or if sufficient numerical damping is present for the increment size used. Other examples show that clearly nonphysical numerical instabilities resulting from the complex characteristics can be produced. These latter types of numerical instabilities are shown to be removed by the addition of physically motivated differential terms which eliminate the complex characteristics. (auth)

Numerical modeling and investigation of two-phase reactive flow in a high-low pressure chambers system

International Nuclear Information System (INIS)

Cheng, Cheng; Zhang, Xiaobing

2016-01-01

Highlights: • A novel two-dimensional two-phase flow model is established for the high-low pressure chambers system. • A strong packing of particles is observed at the projectile base and will cause the pressure to rise faster. • Different length–diameter ratios can affect the flow behavior through the vent-holes obviously. • The muzzle velocity decreases with the length–diameter ratio of the high-pressure chamber. - Abstract: A high-low pressure chambers system is proposed to meet the demands of low launch acceleration for informative equipment in many special fields such as Aeronautics, Astronautics and Weaponry. A two-dimensional two-phase flow numerical model is established to describe the complex physical process based on a modified two-fluid theory, which takes into account gas production, interphase drag, intergranular stress, and heat transfer between two phases. In order to reduce the computational cost, the parameters in the high-pressure chamber at the instant the vent-holes open are calculated by the zero-dimensional model as the initial conditions for the two-phase flow simulation in the high-low pressure chambers system. The simulation results reveal good agreement with the experiments and the launch acceleration of a projectile can be improved by this system. The propellant particles can be tracked clearly in both chambers and a strong packing of particles at the base of projectile will cause the pressure to rise faster than at other areas both in the axis and radial directions. The length–diameter ratio of the high-pressure chamber (a typical multi-dimensional parameter) is investigated. Different length–diameter ratios can affect the maximum pressure drop and the loss of total pressure impulse through the vent-hole, then the muzzle velocity and the launch acceleration of projectiles can be influenced directly. This article puts forward a new prediction tool for the understanding and design of transient processes in high-low pressure

An investigation of two-dimensional, two-phase flow of steam in a cascade of turbine blading by the time-marching method

International Nuclear Information System (INIS)

Teymourtash, A. R.; Mahpeykar, M. R.

2003-01-01

During the course of expansion in turbines, the steam at first super cools and then nucleated to become a two-phase mixture. This is an area where greater understanding can lead to improved design. This paper describes a numerical method for the solution of two-dimensional two-phase flow of steam in a cascade of turbine blading; the unsteady euler equations governing the overall behaviour of the fluid are combined with equations describing droplet behaviour and treated by Jasmine fourth order runge Kutta time marching scheme which modified to allow for two-phase effects. The theoretical surface pressure distributions, droplet radii and contours of constant wetness fraction are presented and results are discussed in the light of knowledge of actual surface pressure distributions

Study on flow instabilities in two-phase mixtures

International Nuclear Information System (INIS)

Ishii, M.

1976-03-01

Various mechanisms that can induce flow instabilities in two-phase flow systems are reviewed and their relative importance discussed. In view of their practical importance, the density-wave instabilities have been analyzed in detail based on the one-dimensional two-phase flow formulation. The dynamic response of the system to the inlet flow perturbations has been derived from the model; thus the characteristic equation that predicts the onset of instabilities has been obtained. The effects of various system parameters, such as the heat flux, subcooling, pressure, inlet velocity, inlet orificing, and exit orificing on the stability boundary have been analyzed. In addition to numerical solutions, some simple stability criteria under particular conditions have been obtained. Both results have been compared with various experimental data, and a satisfactory agreement has been demonstrated

Local wettability reversal during steady-state two-phase flow in porous media.

Science.gov (United States)

Sinha, Santanu; Grøva, Morten; Ødegården, Torgeir Bryge; Skjetne, Erik; Hansen, Alex

2011-09-01

We study the effect of local wettability reversal on remobilizing immobile fluid clusters in steady-state two-phase flow in porous media. We consider a two-dimensional network model for a porous medium and introduce a wettability alteration mechanism. A qualitative change in the steady-state flow patterns, destabilizing the percolating and trapped clusters, is observed as the system wettability is varied. When capillary forces are strong, a finite wettability alteration is necessary to move the system from a single-phase to a two-phase flow regime. When both phases are mobile, we find a linear relationship between fractional flow and wettability alteration.

HRL Aespoe - two-phase flow experiment - gas and water flow in fractured crystalline rock

International Nuclear Information System (INIS)

Kull, H.; Liedtke, L.

1998-01-01

(The full text of the contribution follows:) Gas generated from radioactive waste may influence the hydraulic and mechanical properties of the man-made barriers and the immediate surroundings of the repository. Prediction of alteration in fractured crystalline rock is difficult. There is a lack of experimental data, and calibrated models are not yet available. Because of the general importance of this matter the German Federal Ministry for Education, Science, Research and Technology decided to conduct a two-phase flow study at HRL Aespoe within the scope of the co-operation agreement with SKB. Within the presentation an overview of field experiments and modelling studies scheduled until end of '99 are given. Conceptual models for one- and two-phase flow, methodologies and with respect to numerical calculations necessary parameter set-ups are discussed. Common objective of in-situ experiments is to calibrate flow models to improve the reliability of predictions for gas migration through fractured rock mass. Hence, in a defined dipole flow field in niche 2/715 at HRL Aespoe effective hydraulic parameters are evaluated. Numerical modelling of non-isothermal, two-phase, two-component processes is feasible only for two-dimensional representation of a porous medium. To overcome this restriction a computer program will be developed to model three-dimensional, fractured, porous media. Rational aspects of two-phase flow studies are for the designing of geotechnical barriers and for the long-term safety analysis of potential radionuclide transport in a future repository required for the licensing process

Effective Rheology of Two-Phase Flow in Three-Dimensional Porous Media: Experiment and Simulation.

Science.gov (United States)

Sinha, Santanu; Bender, Andrew T; Danczyk, Matthew; Keepseagle, Kayla; Prather, Cody A; Bray, Joshua M; Thrane, Linn W; Seymour, Joseph D; Codd, Sarah L; Hansen, Alex

2017-01-01

We present an experimental and numerical study of immiscible two-phase flow of Newtonian fluids in three-dimensional (3D) porous media to find the relationship between the volumetric flow rate ( Q ) and the total pressure difference ([Formula: see text]) in the steady state. We show that in the regime where capillary forces compete with the viscous forces, the distribution of capillary barriers at the interfaces effectively creates a yield threshold ([Formula: see text]), making the fluids reminiscent of a Bingham viscoplastic fluid in the porous medium. In this regime, Q depends quadratically on an excess pressure drop ([Formula: see text]). While increasing the flow rate, there is a transition, beyond which the overall flow is Newtonian and the relationship is linear. In our experiments, we build a model porous medium using a column of glass beads transporting two fluids, deionized water and air. For the numerical study, reconstructed 3D pore networks from real core samples are considered and the transport of wetting and non-wetting fluids through the network is modeled by tracking the fluid interfaces with time. We find agreement between our numerical and experimental results. Our results match with the mean-field results reported earlier.

Stochastic Rotation Dynamics simulations of wetting multi-phase flows

Science.gov (United States)

Hiller, Thomas; Sanchez de La Lama, Marta; Brinkmann, Martin

2016-06-01

Multi-color Stochastic Rotation Dynamics (SRDmc) has been introduced by Inoue et al. [1,2] as a particle based simulation method to study the flow of emulsion droplets in non-wetting microchannels. In this work, we extend the multi-color method to also account for different wetting conditions. This is achieved by assigning the color information not only to fluid particles but also to virtual wall particles that are required to enforce proper no-slip boundary conditions. To extend the scope of the original SRDmc algorithm to e.g. immiscible two-phase flow with viscosity contrast we implement an angular momentum conserving scheme (SRD+mc). We perform extensive benchmark simulations to show that a mono-phase SRDmc fluid exhibits bulk properties identical to a standard SRD fluid and that SRDmc fluids are applicable to a wide range of immiscible two-phase flows. To quantify the adhesion of a SRD+mc fluid in contact to the walls we measure the apparent contact angle from sessile droplets in mechanical equilibrium. For a further verification of our wettability implementation we compare the dewetting of a liquid film from a wetting stripe to experimental and numerical studies of interfacial morphologies on chemically structured surfaces.

Two multi-dimensional uncertainty relations

International Nuclear Information System (INIS)

Skala, L; Kapsa, V

2008-01-01

Two multi-dimensional uncertainty relations, one related to the probability density and the other one related to the probability density current, are derived and discussed. Both relations are stronger than the usual uncertainty relations for the coordinates and momentum

Generalized modeling of multi-component vaporization/condensation phenomena for multi-phase-flow analysis

International Nuclear Information System (INIS)

Morita, K.; Fukuda, K.; Tobita, Y.; Kondo, Sa.; Suzuki, T.; Maschek, W.

2003-01-01

A new multi-component vaporization/condensation (V/C) model was developed to provide a generalized model for safety analysis codes of liquid metal cooled reactors (LMRs). These codes simulate thermal-hydraulic phenomena of multi-phase, multi-component flows, which is essential to investigate core disruptive accidents of LMRs such as fast breeder reactors and accelerator driven systems. The developed model characterizes the V/C processes associated with phase transition by employing heat transfer and mass-diffusion limited models for analyses of relatively short-time-scale multi-phase, multi-component hydraulic problems, among which vaporization and condensation, or simultaneous heat and mass transfer, play an important role. The heat transfer limited model describes the non-equilibrium phase transition processes occurring at interfaces, while the mass-diffusion limited model is employed to represent effects of non-condensable gases and multi-component mixture on V/C processes. Verification of the model and method employed in the multi-component V/C model of a multi-phase flow code was performed successfully by analyzing a series of multi-bubble condensation experiments. The applicability of the model to the accident analysis of LMRs is also discussed by comparison between steam and metallic vapor systems. (orig.)

Characterizing the correlations between local phase fractions of gas���liquid two-phase flow with wire-mesh sensor

OpenAIRE

Tan, C.; Liu, W. L.; Dong, F.

2016-01-01

Understanding of flow patterns and their transitions is significant to uncover the flow mechanics of two-phase flow. The local phase distribution and its fluctuations contain rich information regarding the flow structures. A wire-mesh sensor (WMS) was used to study the local phase fluctuations of horizontal gas���liquid two-phase flow, which was verified through comparing the reconstructed three-dimensional flow structure with photographs taken during the experiments. Each crossing point of t...

Two dimensional numerical model for steam--water flow in a sudden contraction

International Nuclear Information System (INIS)

Crowe, C.T.; Choi, H.N.

1976-01-01

A computational model developed for two-dimensional dispersed two-phase flows is applied to steam--water flow in a sudden contraction. The calculational scheme utilizes the cellular approach in which each cell is regarded as a control volume and the droplets are regarded as sources of mass, momentum and energy to the conveying (steam) phase. The predictions show how droplets channel in the entry region and affect the velocity and pressure distributions along the duct

Internal structure and interfacial velocity development for bubbly two-phase flow

International Nuclear Information System (INIS)

Kocamustafaogullari, G.; Huang, W.D.

1994-01-01

This paper describes an experimental study of the internal structure of air-water flowing horizontally. The double-sensor resistivity probe technique was applied for measurements of local interfacial parameters, including void fraction, interfacial area concentration, bubble size distributions, bubble passing frequency and bubble interface velocity. Bubbly flow patterns at several flow conditions were examined at three axial locations, L/D=25, 148 and 253, in which the first measurement represents the entrance region where the flow develops, and the second and third may represent near fully developed bubbly flow patterns. The experimental results are presented in three-dimensional perspective plots of the interfacial parameters over the cross-section. These multi-dimensional presentations showed that the local values of the void fraction, interfacial area concentration and bubble passing frequency were nearly constant over the cross-section at L/D=25, with slight local peaking close to the channel wall. Although similar local peakings were observed at the second and third locations, the internal flow structure segregation due to buoyancy appeared to be very strong in the axial direction. A simple comparison of profiles of the interfacial parameters at the three locations indicated that the flow pattern development was a continuous process. Finally, it was shown that the so-called ''fully developed'' bubbly two-phase flow pattern cannot be established in a horizontal pipe and that there was no strong correspondence between void fraction and interface velocity profiles. ((orig.))

Development of a Two-Phase Flow Analysis Code based on a Unstructured-Mesh SIMPLE Algorithm

Energy Technology Data Exchange (ETDEWEB)

Kim, Jong Tae; Park, Ik Kyu; Cho, Heong Kyu; Yoon, Han Young; Kim, Kyung Doo; Jeong, Jae Jun

2008-09-15

For analyses of multi-phase flows in a water-cooled nuclear power plant, a three-dimensional SIMPLE-algorithm based hydrodynamic solver CUPID-S has been developed. As governing equations, it adopts a two-fluid three-field model for the two-phase flows. The three fields represent a continuous liquid, a dispersed droplets, and a vapour field. The governing equations are discretized by a finite volume method on an unstructured grid to handle the geometrical complexity of the nuclear reactors. The phasic momentum equations are coupled and solved with a sparse block Gauss-Seidel matrix solver to increase a numerical stability. The pressure correction equation derived by summing the phasic volume fraction equations is applied on the unstructured mesh in the context of a cell-centered co-located scheme. This paper presents the numerical method and the preliminary results of the calculations.

Comparison of differential pressure model based on flow regime for gas/liquid two-phase flow

International Nuclear Information System (INIS)

Dong, F; Zhang, F S; Li, W; Tan, C

2009-01-01

Gas/liquid two-phase flow in horizontal pipe is very common in many industry processes, because of the complexity and variability, the real-time parameter measurement of two-phase flow, such as the measurement of flow regime and flow rate, becomes a difficult issue in the field of engineering and science. The flow regime recognition plays a fundamental role in gas/liquid two-phase flow measurement, other parameters of two-phase flow can be measured more easily and correctly based on the correct flow regime recognition result. A multi-sensor system is introduced to make the flow regime recognition and the mass flow rate measurement. The fusion system is consisted of temperature sensor, pressure sensor, cross-section information system and v-cone flow meter. After the flow regime recognition by cross-section information system, comparison of four typical differential pressure (DP) models is discussed based on the DP signal of v-cone flow meter. Eventually, an optimum DP model has been chosen for each flow regime. The experiment result of mass flow rate measurement shows it is efficient to classify the DP models by flow regime.

Heating limits of boiling downward two-phase flow in parallel channels

International Nuclear Information System (INIS)

Fukuda, Kenji; Kondoh, Tetsuya; Hasegawa, Shu; Sakai, Takaaki.

1989-01-01

Flow characteristics and heating limits of downward two-phase flow in single or parallel multi-channels are investigated experimentally and analytically. The heating section used is made of glass tube, in which the heater tube is inserted, and the flow regime inside it is observed. In single channel experiments with low flow rate conditions, it is found that, initially, gas phase which flows upward against the downward liquid phase flow condenses and diminishes as it flows up being cooled by inflowing liquid. However, as the heating power is increased, some portion of the gas phase reaches the top and accumulates to form an liquid level, which eventually causes the dryout. On the other hand, for high flow rate condition, the flooding at the bottom of the heated section is the cause of the dryout. In parallel multi-channels experiments, reversed (upward) flow which leads to the dryout is observed in some of these channels for low flow rate conditions, while the situation is the same to the single channel case for high flow rate conditions. Analyses are carried out to predict the onset of dryout in single channel using the drift flux model as well as the Wallis' flooding correlation. Above-mentioned two types of the dryout and their boundary are predicted which agree well with the experimental results. (author)

Coupling two-phase fluid flow with two-phase darcy flow in anisotropic porous media

KAUST Repository

Chen, J.

2014-06-03

This paper reports a numerical study of coupling two-phase fluid flow in a free fluid region with two-phase Darcy flow in a homogeneous and anisotropic porous medium region. The model consists of coupled Cahn-Hilliard and Navier-Stokes equations in the free fluid region and the two-phase Darcy law in the anisotropic porous medium region. A Robin-Robin domain decomposition method is used for the coupled Navier-Stokes and Darcy system with the generalized Beavers-Joseph-Saffman condition on the interface between the free flow and the porous media regions. Obtained results have shown the anisotropic properties effect on the velocity and pressure of the two-phase flow. 2014 Jie Chen et al.

Coupling Two-Phase Fluid Flow with Two-Phase Darcy Flow in Anisotropic Porous Media

Directory of Open Access Journals (Sweden)

Jie Chen

2014-06-01

Full Text Available This paper reports a numerical study of coupling two-phase fluid flow in a free fluid region with two-phase Darcy flow in a homogeneous and anisotropic porous medium region. The model consists of coupled Cahn-Hilliard and Navier-Stokes equations in the free fluid region and the two-phase Darcy law in the anisotropic porous medium region. A Robin-Robin domain decomposition method is used for the coupled Navier-Stokes and Darcy system with the generalized Beavers-Joseph-Saffman condition on the interface between the free flow and the porous media regions. Obtained results have shown the anisotropic properties effect on the velocity and pressure of the two-phase flow.

Measurement of Liquid-Metal Two-Phase Flow with a Dynamic Neutron Radiography

International Nuclear Information System (INIS)

Cha, J. E.; Lim, I. C.; Kim, H. R.; Kim, C. M.; Nam, H. Y.; Saito, Y.

2005-01-01

The dynamic neutron radiography(DNR) has complementary characteristics to X-ray radiography and is suitable to visualization and measurement of a multi-phase flow research in a metallic duct and liquid metal flow. The flow-field information of liquid metal system is very important for the safety analysis of fast breeder reactor and the design of the spallation target of accelerator driven system. A DNR technique was applied to visualize the flow field in the gas-liquid metal two-phase flow with the HANARO-beam facility. The lead bismuth eutectic and the nitrogen gas were used to construct the two-phase flow field in the natural circulation U-channel. The two-phase flow images in the riser were taken at various combinations of the liquid flow and gas flow with high frame-rate neutron radiography at 1000 fps

Modulating patterns of two-phase flow with electric fields.

Science.gov (United States)

Liu, Dingsheng; Hakimi, Bejan; Volny, Michael; Rolfs, Joelle; Anand, Robbyn K; Turecek, Frantisek; Chiu, Daniel T

2014-07-01

This paper describes the use of electro-hydrodynamic actuation to control the transition between three major flow patterns of an aqueous-oil Newtonian flow in a microchannel: droplets, beads-on-a-string (BOAS), and multi-stream laminar flow. We observed interesting transitional flow patterns between droplets and BOAS as the electric field was modulated. The ability to control flow patterns of a two-phase fluid in a microchannel adds to the microfluidic tool box and improves our understanding of this interesting fluid behavior.

Numerical simulation of two-phase flow with front-capturing

International Nuclear Information System (INIS)

Tzanos, C.P.; Weber, D.P.

2000-01-01

Because of the complexity of two-phase flow phenomena, two-phase flow codes rely heavily on empirical correlations. This approach has a number of serious shortcomings. Advances in parallel computing and continuing improvements in computer speed and memory have stimulated the development of numerical simulation tools that rely less on empirical correlations and more on fundamental physics. The objective of this work is to take advantage of developments in massively parallel computing, single-phase computational fluid dynamics of complex systems, and numerical methods for front capturing in two-phase flows to develop a computer code for direct numerical simulation of two-phase flow. This includes bubble/droplet transport, interface deformation and topology change, bubble-droplet interactions, interface mass, momentum, and energy transfer. In this work, the Navier-Stokes and energy equations are solved by treating both phases as a single fluid with interfaces between the two phases, and a discontinuity in material properties across the moving interfaces. The evolution of the interfaces is simulated by using the front capturing technique of the level-set methods. In these methods, the boundary of a two-fluid interface is modeled as the zero level set of a smooth function φ. The level-set function φ is defined as the signed distance from the interface (φ is negative inside a droplet/bubble and positive outside). Compared to other front-capturing or front-tracking methods, the level-set approach is relatively easy to implement even in three-dimensional flows, and it has been shown to simulate well the coalescence and breakup of droplets/bubbles

Dynamic modelling for two-phase flow systems

International Nuclear Information System (INIS)

Guerra, M.A.

1991-06-01

Several models for two-phase flow have been studied, developing a thermal-hydraulic analysis code with one of these models. The program calculates, for one-dimensional cases with variable flow area, the transient behaviour of system process variables, when the boundary conditions (heat flux, flow rate, enthalpy and pressure) are functions of time. The modular structure of the code, eases the program growth. In fact, the present work is the basis for a general purpose accident and transient analysis code in nuclear reactors. Code verification has been made against RETRAN-02 results. Satisfactory results have been achieved with the present version of the code. (Author) [es

Is it possible to design universal multi-phase flow analyzer?

International Nuclear Information System (INIS)

Ivanov Kolev, N.

2005-01-01

Transient 3D-multiphase flows consisting of many chemical constituents in nature and technology (Figs. 1 and 2) are the common case of flows. In many technical applications we have to do with particular realization of the multi-phase flows like steady state flows, or single component flows or single phase flows etc. Engineers and scientists created hundreds of computer codes for description of more or less specific realizations of multi-phase flows. If one compares the structure of these codes one is astonished by the waste of the human resources for programming repeating model elements like equations of state, friction lows in variety of geometry, heat transfer coefficients, mathematical equation solvers, data handling procedures, graphical environment etc. It is hardly to expect, that the best solution for the specific sub-phenomenon is available in all codes. Looking in other branches of the technology like computer chips production we realize that the revolutionary idea of having common ''chips'' within complex applications is very far from its practical realization in the computational multi-phase flow dynamics. Following this line of arguments I expressed several times in my publications explicitly or implicitly the idea, that it is possible to create a universal multi-phase flow analyzer in the sense of computer architecture, that is capable to absorb the adequate multi-phase knowledge data base specified in Appendix 1. The subject of this paper is to summarize some of the main ideas, some of them already realized by this author, on the way of creating such computer code architecture, to illustrate haw they work, and to make an outlook regarding what are the challenges in the future developments. We confine deliberately our attention to the solution of the so called local volume and time averaged system of PDE's for a simple reason: Direct numerical resolution of interacting fluids is possible as demonstrated for small scales by many researchers, but for

A real two-phase submarine debris flow and tsunami

International Nuclear Information System (INIS)

Pudasaini, Shiva P.; Miller, Stephen A.

2012-01-01

The general two-phase debris flow model proposed by Pudasaini is employed to study subaerial and submarine debris flows, and the tsunami generated by the debris impact at lakes and oceans. The model, which includes three fundamentally new and dominant physical aspects such as enhanced viscous stress, virtual mass, and generalized drag (in addition to buoyancy), constitutes the most generalized two-phase flow model to date. The advantage of this two-phase debris flow model over classical single-phase, or quasi-two-phase models, is that the initial mass can be divided into several parts by appropriately considering the solid volume fraction. These parts include a dry (landslide or rock slide), a fluid (water or muddy water; e.g., dams, rivers), and a general debris mixture material as needed in real flow simulations. This innovative formulation provides an opportunity, within a single framework, to simultaneously simulate the sliding debris (or landslide), the water lake or ocean, the debris impact at the lake or ocean, the tsunami generation and propagation, the mixing and separation between the solid and fluid phases, and the sediment transport and deposition process in the bathymetric surface. Applications of this model include (a) sediment transport on hill slopes, river streams, hydraulic channels (e.g., hydropower dams and plants); lakes, fjords, coastal lines, and aquatic ecology; and (b) submarine debris impact and the rupture of fiber optic, submarine cables and pipelines along the ocean floor, and damage to offshore drilling platforms. Numerical simulations reveal that the dynamics of debris impact induced tsunamis in mountain lakes or oceans are fundamentally different than the tsunami generated by pure rock avalanches and landslides. The analysis includes the generation, amplification and propagation of super tsunami waves and run-ups along coastlines, debris slide and deposition at the bottom floor, and debris shock waves. It is observed that the

A real two-phase submarine debris flow and tsunami

Energy Technology Data Exchange (ETDEWEB)

Pudasaini, Shiva P.; Miller, Stephen A. [Department of Geodynamics and Geophysics, Steinmann Institute, University of Bonn Nussallee 8, D-53115, Bonn (Germany)

2012-09-26

The general two-phase debris flow model proposed by Pudasaini is employed to study subaerial and submarine debris flows, and the tsunami generated by the debris impact at lakes and oceans. The model, which includes three fundamentally new and dominant physical aspects such as enhanced viscous stress, virtual mass, and generalized drag (in addition to buoyancy), constitutes the most generalized two-phase flow model to date. The advantage of this two-phase debris flow model over classical single-phase, or quasi-two-phase models, is that the initial mass can be divided into several parts by appropriately considering the solid volume fraction. These parts include a dry (landslide or rock slide), a fluid (water or muddy water; e.g., dams, rivers), and a general debris mixture material as needed in real flow simulations. This innovative formulation provides an opportunity, within a single framework, to simultaneously simulate the sliding debris (or landslide), the water lake or ocean, the debris impact at the lake or ocean, the tsunami generation and propagation, the mixing and separation between the solid and fluid phases, and the sediment transport and deposition process in the bathymetric surface. Applications of this model include (a) sediment transport on hill slopes, river streams, hydraulic channels (e.g., hydropower dams and plants); lakes, fjords, coastal lines, and aquatic ecology; and (b) submarine debris impact and the rupture of fiber optic, submarine cables and pipelines along the ocean floor, and damage to offshore drilling platforms. Numerical simulations reveal that the dynamics of debris impact induced tsunamis in mountain lakes or oceans are fundamentally different than the tsunami generated by pure rock avalanches and landslides. The analysis includes the generation, amplification and propagation of super tsunami waves and run-ups along coastlines, debris slide and deposition at the bottom floor, and debris shock waves. It is observed that the

Introduction to investigations of the negative corona and EHD flow in gaseous two-phase fluids

Science.gov (United States)

Jerzy, MIZERACZYK; Artur, BERENDT

2018-05-01

Research interests have recently been directed towards electrical discharges in multi-phase environments. Natural electrical discharges, such as lightning and coronas, occur in the Earth’s atmosphere, which is actually a mixture of gaseous phase (air) and suspended solid and liquid particulate matters (PMs). An example of an anthropogenic gaseous multi-phase environment is the flow of flue gas through electrostatic precipitators (ESPs), which are generally regarded as a mixture of a post-combustion gas with solid PM and microdroplets suspended in it. Electrical discharges in multi-phase environments, the knowledge of which is scarce, are becoming an attractive research subject, offering a wide variety of possible discharges and multi-phase environments to be studied. This paper is an introduction to electrical discharges in multi-phase environments. It is focused on DC negative coronas and accompanying electrohydrodynamic (EHD) flows in a gaseous two-phase fluid formed by air (a gaseous phase) and solid PM (a solid phase), run under laboratory conditions. The introduction is based on a review of the relevant literature. Two cases will be considered: the first case is of a gaseous two-phase fluid, initially motionless in a closed chamber before being subjected to a negative corona (with the needle-to-plate electrode arrangement), which afterwards induces an EHD flow in the chamber, and the second, of a gaseous two-phase fluid flowing transversely with respect to the needle-to-plate electrode axis along a chamber with a corona discharge running between the electrodes. This review-based introductory paper should be of interest to theoretical researchers and modellers in the field of negative corona discharges in single- or two-phase fluids, and for engineers who work on designing EHD devices (such as ESPs, EHD pumps, and smoke detectors).

Visualization of velocity field and phase distribution in gas-liquid two-phase flow by NMR imaging

International Nuclear Information System (INIS)

Matsui, G.; Monji, H.; Obata, J.

2004-01-01

NMR imaging has been applied in the field of fluid mechanics, mainly single phase flow, to visualize the instantaneous flow velocity field. In the present study, NMR imaging was used to visualize simultaneously both the instantaneous phase structure and velocity field of gas-liquid two-phase flow. Two methods of NMR imaging were applied. One is useful to visualize both the one component of liquid velocity and the phase distribution. This method was applied to horizontal two-phase flow and a bubble rising in stagnant oil. It was successful in obtaining some pictures of velocity field and phase distribution on the cross section of the pipe. The other is used to visualize a two-dimensional velocity field. This method was applied to a bubble rising in a stagnant water. The velocity field was visualized after and before the passage of a bubble at the measuring cross section. Furthermore, the distribution of liquid velocity was obtained. (author)

System identification on two-phase flow stability

International Nuclear Information System (INIS)

Wu Shaorong; Zhang Youjie; Wang Dazhong; Bo Jinghai; Wang Fei

1996-01-01

The theoretical principle, experimental method and results of interrelation analysis identification for the instability of two-phase flow are described. A completely new concept of test technology and method on two-phase flow stability was developed by using he theory of information science on system stability and system identification for two-phase flow stability in thermo-physics field. Application of this method would make it possible to identify instability boundary of two-phase flow under stable operation conditions of two-phase flow system. The experiment was carried out on the thermohydraulic test system HRTL-5. Using reverse repeated pseudo-random sequences of heating power as input signal sources and flow rate as response function in the test, the two-phase flow stability and stability margin of the natural circulation system are investigated. The effectiveness and feasibility of identifying two-phase flow stability by using this system identification method were experimentally demonstrated. Basic data required for mathematics modeling of two-phase flow and analysis of two-phase flow stability were obtained, which are useful for analyzing, monitoring of the system operation condition, and forecasting of two-phase flow stability in engineering system

Two-dimensional turbulent flows on a bounded domain

NARCIS (Netherlands)

Kramer, W.

2006-01-01

Large-scale flows in the oceans and the atmosphere reveal strong similarities with purely two-dimensional flows. One of the most typical features is the cascade of energy from smaller flow scales towards larger scales. This is opposed to three-dimensional turbulence where larger flow structures

Tube Radial Distribution Flow Separation in a Microchannel Using an Ionic Liquid Aqueous Two-Phase System Based on Phase Separation Multi-Phase Flow.

Science.gov (United States)

Nagatani, Kosuke; Shihata, Yoshinori; Matsushita, Takahiro; Tsukagoshi, Kazuhiko

2016-01-01

Ionic liquid aqueous two-phase systems were delivered into a capillary tube to achieve tube radial distribution flow (TRDF) or annular flow in a microspace. The phase diagram, viscosity of the phases, and TRDF image of the 1-butyl-3-methylimidazolium chloride and NaOH system were examined. The TRDF was formed with inner ionic liquid-rich and outer ionic liquid-poor phases in the capillary tube. The phase configuration was explained using the viscous dissipation principle. We also examined the distribution of rhodamine B in a three-branched microchannel on a microchip with ionic liquid aqueous two-phase systems for the first time.

Stochastic modelling of two-phase flows including phase change

International Nuclear Information System (INIS)

Hurisse, O.; Minier, J.P.

2011-01-01

Stochastic modelling has already been developed and applied for single-phase flows and incompressible two-phase flows. In this article, we propose an extension of this modelling approach to two-phase flows including phase change (e.g. for steam-water flows). Two aspects are emphasised: a stochastic model accounting for phase transition and a modelling constraint which arises from volume conservation. To illustrate the whole approach, some remarks are eventually proposed for two-fluid models. (authors)

Turbine flow meter response in two-phase flows

International Nuclear Information System (INIS)

Shim, W.J.; Dougherty, T.J.; Cheh, H.Y.

1996-01-01

The purpose of this paper is to suggest a simple method of calibrating turbine flow meters to measure the flow rates of each phase in a two-phase flow. The response of two 50.8 mm (2 inch) turbine flow meters to air-water, two-phase mixtures flowing vertically in a 57 mm I.D. (2.25 inch) polycarbonate tube has been investigated for both upflow and downflow. The flow meters were connected in series with an intervening valve to provide an adjustable pressure difference between them. Void fractions were measured by two gamma densitometers, one upstream of the flow meters and the other downstream. The output signal of the turbine flow meters was found to depend only on the actual volumetric flow rate of the gas, F G , and liquid, F L , at the location of the flow meter

Single phase and two-phase flow pressure losses through restrictions, expansions and inserts

International Nuclear Information System (INIS)

Glenat, P.; Solignac, P.

1984-11-01

We give a selection of methods to predict pressure losses through retrictions, expansions and inserts. In single phase flow, we give the classical method based on the one-dimensional momentum and mass balances. In two-phase flow, we propose the method given by Harshe et al. and an empirical approach suggested by Chisholm. We notice the distinction between long and short inserts depends upon wether or not the vena contracta lies within insert. Finally, we propose three correlations to calculate void fraction through the singularities which have been considered [fr

Hydrodynamics of single- and two-phase flow in inclined rod arrays

International Nuclear Information System (INIS)

Todreas, N.E.

1984-01-01

Required inputs for thermal-hydraulic codes are constitutive relations for fluid-solid flow resistance, in single-phase flow, and interfacial momentum exchange (relative phase motion), in two-phase flow. An inclined rod array air-water experiment was constructed to study the hydrodynamics of multidimensional porous medium flow in rod arrays. Velocities, pressures, bubble distributions, and void fractions were measured in inline and rotational square rod arrays of P/d = 1.5, at 0, 30, 45, and 90 degree inclinations to the vertical flow direction. Constitutive models for single-phase flow resistance are reviewed, new comprehensive models developed, and an assessment with previously published and new data made. The principle of superimposing one-dimensional correlations proves successful for turbulent single-phase inclined flow. For bubbly two-phase yawed flow through incline rod arrays a new flow separation phenomena was observed and modeled. Bubbles of diameters significantly smaller than the rod diameter travel along the rod axis, while larger diameter bubbles move through the rod array gaps. The outcome is a flow separation not predictable with current interfacial momentum exchange models. This phenomenon was not observed in rotated square rod arrays. Current interfacial momentum exchange models were confirmed for this rod arrangement. Models for the two phase flow resistance multiplier for cross flow were reviewed and compared with data from cross and yawed flow rod arrays. Both drag and lift components of the multiplier were well predicted by the homogenous model. Other models reviewed overpredicted the data by a factor of two

Domain decomposition parallel computing for transient two-phase flow of nuclear reactors

Energy Technology Data Exchange (ETDEWEB)

Lee, Jae Ryong; Yoon, Han Young [KAERI, Daejeon (Korea, Republic of); Choi, Hyoung Gwon [Seoul National University, Seoul (Korea, Republic of)

2016-05-15

KAERI (Korea Atomic Energy Research Institute) has been developing a multi-dimensional two-phase flow code named CUPID for multi-physics and multi-scale thermal hydraulics analysis of Light water reactors (LWRs). The CUPID code has been validated against a set of conceptual problems and experimental data. In this work, the CUPID code has been parallelized based on the domain decomposition method with Message passing interface (MPI) library. For domain decomposition, the CUPID code provides both manual and automatic methods with METIS library. For the effective memory management, the Compressed sparse row (CSR) format is adopted, which is one of the methods to represent the sparse asymmetric matrix. CSR format saves only non-zero value and its position (row and column). By performing the verification for the fundamental problem set, the parallelization of the CUPID has been successfully confirmed. Since the scalability of a parallel simulation is generally known to be better for fine mesh system, three different scales of mesh system are considered: 40000 meshes for coarse mesh system, 320000 meshes for mid-size mesh system, and 2560000 meshes for fine mesh system. In the given geometry, both single- and two-phase calculations were conducted. In addition, two types of preconditioners for a matrix solver were compared: Diagonal and incomplete LU preconditioner. In terms of enhancement of the parallel performance, the OpenMP and MPI hybrid parallel computing for a pressure solver was examined. It is revealed that the scalability of hybrid calculation was enhanced for the multi-core parallel computation.

Generalized network modeling of capillary-dominated two-phase flow.

Science.gov (United States)

Raeini, Ali Q; Bijeljic, Branko; Blunt, Martin J

2018-02-01

We present a generalized network model for simulating capillary-dominated two-phase flow through porous media at the pore scale. Three-dimensional images of the pore space are discretized using a generalized network-described in a companion paper [A. Q. Raeini, B. Bijeljic, and M. J. Blunt, Phys. Rev. E 96, 013312 (2017)2470-004510.1103/PhysRevE.96.013312]-which comprises pores that are divided into smaller elements called half-throats and subsequently into corners. Half-throats define the connectivity of the network at the coarsest level, connecting each pore to half-throats of its neighboring pores from their narrower ends, while corners define the connectivity of pore crevices. The corners are discretized at different levels for accurate calculation of entry pressures, fluid volumes, and flow conductivities that are obtained using direct simulation of flow on the underlying image. This paper discusses the two-phase flow model that is used to compute the averaged flow properties of the generalized network, including relative permeability and capillary pressure. We validate the model using direct finite-volume two-phase flow simulations on synthetic geometries, and then present a comparison of the model predictions with a conventional pore-network model and experimental measurements of relative permeability in the literature.

Generalized network modeling of capillary-dominated two-phase flow

Science.gov (United States)

Raeini, Ali Q.; Bijeljic, Branko; Blunt, Martin J.

2018-02-01

We present a generalized network model for simulating capillary-dominated two-phase flow through porous media at the pore scale. Three-dimensional images of the pore space are discretized using a generalized network—described in a companion paper [A. Q. Raeini, B. Bijeljic, and M. J. Blunt, Phys. Rev. E 96, 013312 (2017), 10.1103/PhysRevE.96.013312]—which comprises pores that are divided into smaller elements called half-throats and subsequently into corners. Half-throats define the connectivity of the network at the coarsest level, connecting each pore to half-throats of its neighboring pores from their narrower ends, while corners define the connectivity of pore crevices. The corners are discretized at different levels for accurate calculation of entry pressures, fluid volumes, and flow conductivities that are obtained using direct simulation of flow on the underlying image. This paper discusses the two-phase flow model that is used to compute the averaged flow properties of the generalized network, including relative permeability and capillary pressure. We validate the model using direct finite-volume two-phase flow simulations on synthetic geometries, and then present a comparison of the model predictions with a conventional pore-network model and experimental measurements of relative permeability in the literature.

Energy law preserving C0 finite element schemes for phase field models in two-phase flow computations

International Nuclear Information System (INIS)

Hua Jinsong; Lin Ping; Liu Chun; Wang Qi

2011-01-01

Highlights: → We study phase-field models for multi-phase flow computation. → We develop an energy-law preserving C0 FEM. → We show that the energy-law preserving method work better. → We overcome unphysical oscillation associated with the Cahn-Hilliard model. - Abstract: We use the idea in to develop the energy law preserving method and compute the diffusive interface (phase-field) models of Allen-Cahn and Cahn-Hilliard type, respectively, governing the motion of two-phase incompressible flows. We discretize these two models using a C 0 finite element in space and a modified midpoint scheme in time. To increase the stability in the pressure variable we treat the divergence free condition by a penalty formulation, under which the discrete energy law can still be derived for these diffusive interface models. Through an example we demonstrate that the energy law preserving method is beneficial for computing these multi-phase flow models. We also demonstrate that when applying the energy law preserving method to the model of Cahn-Hilliard type, un-physical interfacial oscillations may occur. We examine the source of such oscillations and a remedy is presented to eliminate the oscillations. A few two-phase incompressible flow examples are computed to show the good performance of our method.

Measurements of liquid-phase turbulence in gas–liquid two-phase flows using particle image velocimetry

International Nuclear Information System (INIS)

Zhou, Xinquan; Doup, Benjamin; Sun, Xiaodong

2013-01-01

Liquid-phase turbulence measurements were performed in an air–water two-phase flow loop with a circular test section of 50 mm inner diameter using a particle image velocimetry (PIV) system. An optical phase separation method-–planar laser-induced fluorescence (PLIF) technique—which uses fluorescent particles and an optical filtration technique, was employed to separate the signals of the fluorescent seeding particles from those due to bubbles and other noises. An image pre-processing scheme was applied to the raw PIV images to remove the noise residuals that are not removed by the PLIF technique. In addition, four-sensor conductivity probes were adopted to measure the radial distribution of the void fraction. Two benchmark tests were performed: the first was a comparison of the PIV measurement results with those of similar flow conditions using thermal anemometry from previous studies; the second quantitatively compared the superficial liquid velocities calculated from the local liquid velocity and void fraction measurements with the global liquid flow rate measurements. The differences of the superficial liquid velocity obtained from the two measurements were bounded within ±7% for single-phase flows and two-phase bubbly flows with the area-average void fraction up to 18%. Furthermore, a preliminary uncertainty analysis was conducted to investigate the accuracy of the two-phase PIV measurements. The systematic uncertainties due to the circular pipe curvature effects, bubble surface reflection effects and other potential uncertainty sources of the PIV measurements were discussed. The purpose of this work is to facilitate the development of a measurement technique (PIV-PLIF) combined with image pre-processing for the liquid-phase turbulence in gas–liquid two-phase flows of relatively high void fractions. The high-resolution data set can be used to more thoroughly understand two-phase flow behavior, develop liquid-phase turbulence models, and assess high

A Two-Temperature Open-Source CFD Model for Hypersonic Reacting Flows, Part Two: Multi-Dimensional Analysis †

OpenAIRE

Vincent Casseau; Daniel E. R. Espinoza; Thomas J. Scanlon; Richard E. Brown

2016-01-01

hy2Foam is a newly-coded open-source two-temperature computational fluid dynamics (CFD) solver that has previously been validated for zero-dimensional test cases. It aims at (1) giving open-source access to a state-of-the-art hypersonic CFD solver to students and researchers; and (2) providing a foundation for a future hybrid CFD-DSMC (direct simulation Monte Carlo) code within the OpenFOAM framework. This paper focuses on the multi-dimensional verification of hy2Foam and firstly describes th...

Three-dimensional two-phase mass transport model for direct methanol fuel cells

International Nuclear Information System (INIS)

Yang, W.W.; Zhao, T.S.; Xu, C.

2007-01-01

A three-dimensional (3D) steady-state model for liquid feed direct methanol fuel cells (DMFC) is presented in this paper. This 3D mass transport model is formed by integrating five sub-models, including a modified drift-flux model for the anode flow field, a two-phase mass transport model for the porous anode, a single-phase model for the polymer electrolyte membrane, a two-phase mass transport model for the porous cathode, and a homogeneous mist-flow model for the cathode flow field. The two-phase mass transport models take account the effect of non-equilibrium evaporation/ condensation at the gas-liquid interface. A 3D computer code is then developed based on the integrated model. After being validated against the experimental data reported in the literature, the code was used to investigate numerically transport behaviors at the DMFC anode and their effects on cell performance

A multi-phase ferrofluid flow model with equation of state for thermomagnetic pumping and heat transfer

Energy Technology Data Exchange (ETDEWEB)

Aursand, Eskil, E-mail: eskil.aursand@sintef.no; Gjennestad, Magnus Aa.; Yngve Lervåg, Karl; Lund, Halvor

2016-03-15

A one-dimensional multi-phase flow model for thermomagnetically pumped ferrofluid with heat transfer is proposed. The thermodynamic model is a combination of a simplified particle model and thermodynamic equations of state for the base fluid. The magnetization model is based on statistical mechanics, taking into account non-uniform particle size distributions. An implementation of the proposed model is validated against experiments from the literature, and found to give good predictions for the thermomagnetic pumping performance. However, the results reveal a very large sensitivity to uncertainties in heat transfer coefficient predictions. - Highlights: • A multi-phase flow model for thermomagnetically pumped ferrofluid is proposed. • An implementation is validated against experiments from the literature. • Predicted thermomagnetic pumping effect agrees with experiments. • However, a very large sensitivity to heat transfer coefficient is revealed.

Development of an Efficient Meso- scale Multi-phase Flow Solver in Nuclear Applications

Energy Technology Data Exchange (ETDEWEB)

Lee, Taehun [City Univ. (CUNY), NY (United States)

2015-10-20

The proposed research aims at formulating a predictive high-order Lattice Boltzmann Equation for multi-phase flows relevant to nuclear energy related application - namely, saturated and sub-cooled boiling in reactors, and liquid- liquid mixing and extraction for fuel cycle separation. An efficient flow solver will be developed based on the Finite Element based Lattice Boltzmann Method (FE- LBM), accounting for phase-change heat transfer and capable of treating multiple phases over length scales from the submicron to the meter. A thermal LBM will be developed in order to handle adjustable Prandtl number, arbitrary specific heat ratio, a wide range of temperature variations, better numerical stability during liquid-vapor phase change, and full thermo-hydrodynamic consistency. Two-phase FE-LBM will be extended to liquid–liquid–gas multi-phase flows for application to high-fidelity simulations building up from the meso-scale up to the equipment sub-component scale. While several relevant applications exist, the initial applications for demonstration of the efficient methods to be developed as part of this project include numerical investigations of Critical Heat Flux (CHF) phenomena in nuclear reactor fuel bundles, and liquid-liquid mixing and interfacial area generation for liquid-liquid separations. In addition, targeted experiments will be conducted for validation of this advanced multi-phase model.

Large eddy simulation of a two-phase reacting swirl flow inside a cement cyclone

International Nuclear Information System (INIS)

Mikulčić, Hrvoje; Vujanović, Milan; Ashhab, Moh'd Sami; Duić, Neven

2014-01-01

This work presents a numerical study of the highly swirled gas–solid flow inside a cement cyclone. The computational fluid dynamics – CFD simulation for continuum fluid flow and heat exchange was used for the investigation. The Eulearian–Lagrangian approach was used to describe the two-phase flow, and the large eddy simulation – LES method was used for correctly obtaining the turbulent fluctuations of the gas phase. A model describing the reaction of the solid phase, e.g. the calcination process, has been developed and implemented within the commercial finite volume CFD code FIRE. Due to the fact that the calcination process has a direct influence on the overall energy efficiency of the cement production, it is of great importance to have a certain degree of limestone degradation at the cyclone's outlet. The heat exchange between the gas and solid phase is of particular importance when studying cement cyclones, as it has a direct effect on the calcination process. In order to study the heat exchange phenomena and the flow characteristics, a three dimensional geometry of a real industrial scroll type cyclone was used for the CFD simulation. The gained numerical results, characteristic for cyclones, such as the pressure drop, and concentration of particles can thus be used for better understanding of the complex swirled two-phase flow inside the cement cyclone and also for improving the heat exchange phenomena. - Highlights: • CFD (computational fluid dynamics) is being increasingly used to enhance efficiency of reacting multi-phase flows. • Numerical model of calcination process was presented. • A detailed industrial geometry was used for the CFD simulation. • Presented model and measurement data are in good agreement

Stratified steady and unsteady two-phase flows between two parallel plates

International Nuclear Information System (INIS)

Sim, Woo Gun

2006-01-01

To understand fluid dynamic forces acting on a structure subjected to two-phase flow, it is essential to get detailed information about the characteristics of two-phase flow. Stratified steady and unsteady two-phase flows between two parallel plates have been studied to investigate the general characteristics of the flow related to flow-induced vibration. Based on the spectral collocation method, a numerical approach has been developed for the unsteady two-phase flow. The method is validated by comparing numerical result to analytical one given for a simple harmonic two-phase flow. The flow parameters for the steady two-phase flow, such as void fraction and two-phase frictional multiplier, are evaluated. The dynamic characteristics of the unsteady two-phase flow, including the void fraction effect on the complex unsteady pressure, are illustrated

Multi-phase flow modeling of soil contamination and soil remediation

NARCIS (Netherlands)

Dijke, van M.I.J.

1997-01-01

In this thesis multi-phase flow models are used to study the flow behavior of liquid contaminants in aquifers and of gases that are injected below the groundwater table for remediation purposes. Considered problems are redistribution of a lens of light nonaqueous phase

Nonlinear analysis of gas-water/oil-water two-phase flow in complex networks

CERN Document Server

Gao, Zhong-Ke; Wang, Wen-Xu

2014-01-01

Understanding the dynamics of multi-phase flows has been a challenge in the fields of nonlinear dynamics and fluid mechanics. This chapter reviews our work on two-phase flow dynamics in combination with complex network theory. We systematically carried out gas-water/oil-water two-phase flow experiments for measuring the time series of flow signals which is studied in terms of the mapping from time series to complex networks. Three network mapping methods were proposed for the analysis and identification of flow patterns, i.e. Flow Pattern Complex Network (FPCN), Fluid Dynamic Complex Network (FDCN) and Fluid Structure Complex Network (FSCN). Through detecting the community structure of FPCN based on K-means clustering, distinct flow patterns can be successfully distinguished and identified. A number of FDCN’s under different flow conditions were constructed in order to reveal the dynamical characteristics of two-phase flows. The FDCNs exhibit universal power-law degree distributions. The power-law exponent ...

Multi dimensional analysis of Design Basis Events using MARS-LMR

International Nuclear Information System (INIS)

Woo, Seung Min; Chang, Soon Heung

2012-01-01

the multi-dimensional analysis. In particular, we can see the reverse and mixing flow phenomena in LOHS better than TOP and LOF, which affect the peak temperature in the one dimensional analysis at 330 s but is not observed in the multi-dimensional analysis. In addition, the different temperatures at the central line of fuel and cladding are shown due to different mass flow in the core. The two PDRCs show the almost same results in one and multi-dimensional analysis, respectively. However, in the comparison between the one and multi-dimension, the heat removal capability of PDRC in the multi-dimensional analysis is better than the one dimensional analysis at the beginning of calculation. The aspects change at the ending of calculation due to the different mass flow in PDRC and over flow sodium temperature. In wrap up, the sodium temperatures of the core in the one dimensional analysis are generally higher than the multidimensional analysis. The temperature peak point in the one dimensional analysis under LOHS does not show in the multi-dimensional analysis, because the reverse and mixing flow pattern between the hot pool and the core outlet region. Even thought there are some different results between the two different dimensional analyses, we can say that the KALIMER-600 system has the safety margins.

Experimental study on liquid velocity in upward and downward two-phase flows

International Nuclear Information System (INIS)

Sun, X.; Paranjape, S.; Kim, S.; Ozar, B.; Ishii, M.

2003-01-01

Local characteristics of the liquid phase in upward and downward air-water two-phase flows were experimentally investigated in a 50.8-mm inner-diameter round pipe. An integral Laser Doppler Anemometry (LDA) system was used to measure the axial liquid velocity and its fluctuations. No effect of the flow direction on the liquid velocity radial profile was observed in single-phase liquid benchmark experiments. Local multi-sensor conductivity probes were used to measure the radial profiles of the bubble velocity and the void fraction. The measurement results in the upward and downward two-phase flows are compared and discussed. The results in the downward flow demonstrated that the presence of the bubbles tended to flatten the liquid velocity radial profile, and the maximum liquid velocity could occur off the pipe centerline, in particular at relatively low flow rates. However, the maximum liquid velocity always occurred at the pipe center in the upward flow. Also, noticeable turbulence enhancement due to the bubbles in the two-phase flows was observed in the current experimental flow conditions. Furthermore, the distribution parameter and the void weighted area-averaged drift velocity were obtained based on the definitions

An incompressible two-dimensional multiphase particle-in-cell model for dense particle flows

Energy Technology Data Exchange (ETDEWEB)

Snider, D.M. [SAIC, Albuquerque, NM (United States); O`Rourke, P.J. [Los Alamos National Lab., NM (United States); Andrews, M.J. [Texas A and M Univ., College Station, TX (United States). Dept. of Mechanical Engineering

1997-06-01

A two-dimensional, incompressible, multiphase particle-in-cell (MP-PIC) method is presented for dense particle flows. The numerical technique solves the governing equations of the fluid phase using a continuum model and those of the particle phase using a Lagrangian model. Difficulties associated with calculating interparticle interactions for dense particle flows with volume fractions above 5% have been eliminated by mapping particle properties to a Eulerian grid and then mapping back computed stress tensors to particle positions. This approach utilizes the best of Eulerian/Eulerian continuum models and Eulerian/Lagrangian discrete models. The solution scheme allows for distributions of types, sizes, and density of particles, with no numerical diffusion from the Lagrangian particle calculations. The computational method is implicit with respect to pressure, velocity, and volume fraction in the continuum solution thus avoiding courant limits on computational time advancement. MP-PIC simulations are compared with one-dimensional problems that have analytical solutions and with two-dimensional problems for which there are experimental data.

One-dimensional two-fluid model for wavy flow beyond the Kelvin–Helmholtz instability: Limit cycles and chaos

Energy Technology Data Exchange (ETDEWEB)

Lopez de Bertodano, Martín, E-mail: bertodan@purdue.edu [School of Nuclear Engineering, Purdue University, West Lafayette, IN 47907 (United States); Fullmer, William D. [Department of Chemical and Biological Engineering, U. of Colorado, Boulder, CO 80309 (United States); Clausse, Alejandro [CNEA-CONICET and Universidad Nacional del Centro, 7000 Tandil (Argentina)

2016-12-15

A 1D TFM numerical simulation of near horizontal stratified two-phase flow is performed where the TFM, including surface tension and viscous stresses, is simplified to a two-equation model using the fixed-flux approximation. As the angle of inclination of the channel increases so does the driving body force, so the flow becomes KH unstable, and waves grow and develop nonlinearities. It is shown that these waves grow until they reach a limit cycle due to viscous dissipation at wave fronts. Upon further inclination of the channel, chaos is observed. The appearance of chaos in a 1D TFM implies a nonlinear process that transfers energy intermittently from long wavelengths where energy is produced to short wavelengths where energy is dissipated by viscosity, so that an averaged energy equilibrium in frequency space is attained. This is comparable to the well-known turbulent stability mechanism of the multi-dimensional Navier–Stokes equations, i.e., chaos implies Lyapunov stability, but in this case it is strictly a two-phase phenomenon.

Optical measuring techniques and their application to two-phase and three-phase flows

International Nuclear Information System (INIS)

Liu Xiaozhi.

1992-01-01

First of all it is shown that by an optical system based on the Laser-Doppler technology, which uses a pair of cylindrical waves and two optical detectors, the particle size, speed and refractive index can be measured by means of the signal frequencies. The second optical method to characterize spherical particles in a multi-phase flow is an extended phase-Doppler system. By means of an additional pair of photodetectors it has been possible for the first time to measure the refractive index in addition to speed and particle size. The last part of the paper shows that by a special phase-Doppler anemometry system with only two detectors it is also possible to distinguish between reflecting and refractive particles. By means of such PDA system measurements were made in a gas-fluid-solid three-phase flow directed vertically upwards. (orig./DG) [de

An analytical model for prediction of two-phase (noncondensable) flow pump performance

International Nuclear Information System (INIS)

Furuya, O.

1985-01-01

During operational transients or a hypothetical LOCA (loss of coolant accident) condition, the recirculating coolant of PWR (pressurized water reactor) may flash into steam due to a loss of line pressure. Under such two-phase flow conditions, it is well known that the recirculation pump becomes unable to generate the same head as that of the single-phase flow case. Similar situations also exist in oil well submersible pumps where a fair amount of gas is contained in oil. Based on the one dimensional control volume method, an analytical method has been developed to determine the performance of pumps operating under two-phase flow conditions. The analytical method has incorporated pump geometry, void fraction, flow slippage and flow regime into the basic formula, but neglected the compressibility and condensation effects. During the course of model development, it has been found that the head degradation is mainly caused by higher acceleration on liquid phase and deceleration on gas phase than in the case of single-phase flows. The numerical results for head degradations and torques obtained with the model favorably compared with the air/water two-phase flow test data of Babcock and Wilcox (1/3 scale) and Creare (1/20 scale) pumps

Numerical and dimensional analysis of nanoparticles transport with two-phase flow in porous media

KAUST Repository

El-Amin, Mohamed

2015-04-01

In this paper, a mathematical model and numerical simulation are developed to describe the imbibition of nanoparticles-water suspension into two-phase flow in a porous medium. The flow system may be changed from oil-wet to water-wet due to nanoparticles (which are also water-wet) deposition on surface of the pores. So, the model is extended to include the negative capillary pressure and mixed-wet relative permeability correlations to fit with the mixed-wet system. Moreover, buoyancy and capillary forces as well as Brownian diffusion and mechanical dispersion are considered in the mathematical model. An example of countercurrent imbibition in a core of small scale is considered. A dimensional analysis of the governing equations is introduced to examine contributions of each term of the model. Several important dimensionless numbers appear in the dimensionless equations, such as Darcy number Da, capillary number Ca, and Bond number Bo. Throughout this investigation, we monitor the changing of the fluids and solid properties due to addition of the nanoparticles using numerical experiments.

Electric capacitance tomography and two-phase flow for the nuclear reactor safety analysis

International Nuclear Information System (INIS)

Lee, Jae Young

2008-01-01

Recently electric capacitance tomography has been developed to be used in the analysis of two-phase flow. Although its electric field is not focused as the hard ray tomography such as the X-ray or gamma ray, its convenience of easy access to the system and easy maintenance due to no requirement of radiation shielding benefits us in its application in the two-phase flow study, one of important area in the nuclear safety analysis. In the present paper, the practical technologies in the electric capacitance tomography are represented in both parts of hardware and software. In the software part, both forward problem and inverse problem are discussed and the method of regularization. In the hardware part, the brief discussion of the electronics circuits is made which provides femto farad resolution with a reasonable speed (150 frame/sec for 16 electrodes). Some representative ideal cases are studied to demonstrate its potential capability for the two-phase flow analysis. Also, some variations of the tomography such as axial tomography, and three dimensional tomography are discussed. It was found that the present ECT is expected to become a useful tool to understand the complicated three dimensional two-phase flow which may be an important feature to be equipped by the safety analysis codes. (author)

Implicit approximate Riemann solver for two fluid two phase flow models

International Nuclear Information System (INIS)

Raymond, P.; Toumi, I.; Kumbaro, A.

1993-01-01

This paper is devoted to the description of new numerical methods developed for the numerical treatment of two phase flow models with two velocity fields which are now widely used in nuclear engineering for design or safety calculations. These methods are finite volumes numerical methods and are based on the use of Approximate Riemann Solver's concepts in order to define convective flux versus mean cell quantities. The first part of the communication will describe the numerical method for a three dimensional drift flux model and the extensions which were performed to make the numerical scheme implicit and to have fast running calculations of steady states. Such a scheme is now implemented in the FLICA-4 computer code devoted to 3-D steady state and transient core computations. We will present results obtained for a steady state flow with rod bow effect evaluation and for a Steam Line Break calculation were the 3-D core thermal computation was coupled with a 3-D kinetic calculation and a thermal-hydraulic transient calculation for the four loops of a Pressurized Water Reactor. The second part of the paper will detail the development of an equivalent numerical method based on an approximate Riemann Solver for a two fluid model with two momentum balance equations for the liquid and the gas phases. The main difficulty for these models is due to the existence of differential modelling terms such as added mass effects or interfacial pressure terms which make hyperbolic the model. These terms does not permit to write the balance equations system in a conservative form, and the classical theory for discontinuity propagation for non-linear systems cannot be applied. Meanwhile, the use of non-conservative products theory allows the study of discontinuity propagation for a non conservative model and this will permit the construction of a numerical scheme for two fluid two phase flow model. These different points will be detailed in that section which will be illustrated by

CFD prediction of flow and phase distribution in fuel assemblies with spacers

Energy Technology Data Exchange (ETDEWEB)

Anglart, H.; Nylund, O. [ABB Atom AB, Vasteras (Switzerland); Kurul, N. [Rensselaer Polytechnic Institute, Troy, NY (United States)] [and others

1995-09-01

This paper is concerned with the modeling and computation of multi-dimensional two-phase flows in BWR fuel assemblies. The modeling principles are presented based on using a two-fluid model in which lateral interfacial effects are accounted for. This model has been used to evaluate the velocity fields of both vapor and liquid phases, as well as phase distribution, between fuel elements in geometries similar to BWR fuel bundles. Furthermore, this model has been used to predict, in a detailed mechanistic manner, the effects of spacers on flow and phase distribution between, and pressure drop along, fuel elements. The related numerical simulations have been performed using a CFD computer code, CFDS-FLOW3D.

Sensitivity studies on the multi-sensor conductivity probe measurement technique for two-phase flows

Energy Technology Data Exchange (ETDEWEB)

Worosz, Ted [Department of Mechanical and Nuclear Engineering, The Pennsylvania State University, 230 Reber Building, University Park, PA 16802 (United States); Bernard, Matt [The United States Nuclear Regulatory Commission, 11545 Rockville Pike, Rockville, MD 20852 (United States); Kong, Ran; Toptan, Aysenur [Department of Mechanical and Nuclear Engineering, The Pennsylvania State University, 230 Reber Building, University Park, PA 16802 (United States); Kim, Seungjin, E-mail: skim@psu.edu [Department of Mechanical and Nuclear Engineering, The Pennsylvania State University, 230 Reber Building, University Park, PA 16802 (United States); Hoxie, Chris [The United States Nuclear Regulatory Commission, 11545 Rockville Pike, Rockville, MD 20852 (United States)

2016-12-15

Highlights: • Revised conductivity probe circuit to eliminate signal “ghosting” among sensors. • Higher sampling frequencies suggested for bubble number frequency and a{sub i} measurements. • Two-phase parameter sensitivity to measurement duration and bubble number investigated. • Sensors parallel to pipe wall recommended for symmetric bubble velocity measurements. • Sensor separation distance ratio (s/d) greater than four minimizes bubble velocity error. - Abstract: The objective of this study is to advance the local multi-sensor conductivity probe measurement technique through systematic investigation into several practical aspects of a conductivity probe measurement system. Firstly, signal “ghosting” among probe sensors is found to cause artificially high bubble velocity measurements and low interfacial area concentration (a{sub i}) measurements that depend on sampling frequency and sensor impedance. A revised electrical circuit is suggested to eliminate this artificial variability. Secondly, the sensitivity of the probe measurements to sampling frequency is investigated in 13 two-phase flow conditions with superficial liquid and gas velocities ranging from 1.00–5.00 m/s and 0.17–2.0 m/s, respectively. With increasing gas flow rate, higher sampling frequencies, greater than 100 kHz in some cases, are required to adequately capture the bubble number frequency and a{sub i} measurements. This trend is due to the increase in gas velocity and the transition to the slug flow regime. Thirdly, the sensitivity of the probe measurements to the measurement duration as well as the sample number is investigated for the same flow conditions. Measurements of both group-I (spherical/distorted) and group-II (cap/slug/churn-turbulent) bubbles are found to be relatively insensitive to both the measurement duration and the number of bubbles, as long as the measurements are made for a duration long enough to capture a collection of samples characteristic to a

A turbulent two-phase flow model for nebula flows

International Nuclear Information System (INIS)

Champney, J.M.; Cuzzi, J.N.

1990-01-01

A new and very efficient turbulent two-phase flow numericaly model is described to analyze the environment of a protoplanetary nebula at a stage prior to the formation of planets. Focus is on settling processes of dust particles in flattened gaseous nebulae. The model employs a perturbation technique to improve the accuracy of the numerical simulations of such flows where small variations of physical quantities occur over large distance ranges. The particles are allowed to be diffused by gas turbulence in addition to settling under gravity. Their diffusion coefficients is related to the gas turbulent viscosity by the non-dimensional Schmidt number. The gas turbulent viscosity is determined by the means of the eddy viscosity hypothesis that assumes the Reynolds stress tensor proportional to the mean strain rate tensor. Zero- and two-equation turbulence models are employed. Modeling assumptions are detailed and discussed. The numerical model is shown to reproduce an existing analytical solution for the settling process of particles in an inviscid nebula. Results of nebula flows are presented taking into account turbulence effects of nebula flows. Diffusion processes are found to control the settling of particles. 24 refs

Three-dimensional fluid mechanics of particulate two-phase flows in U-bend and helical conduits

Science.gov (United States)

Tiwari, Prashant; Antal, Steven P.; Podowski, Michael Z.

2006-04-01

The results of numerous studies performed to date have shown that the performance of various hydraulic systems can be significantly improved by using curved conduit geometries instead of straight tubes. In particular, the formation of Dean vortices, which enhance the development of centrifugal instabilities, has been identified as a factor behind reducing the near-wall concentration buildup in particulate flow devices (e.g., in membrane filtration modules). Still, several issues regarding the effect of conduit curvature on local multidimensional phenomena governing fluid flow still remain open. A related issue is concerned with the impact that conduit geometry makes on the concentration distribution of a dispersed phase in two-phase flows in general, and in particulate flows (solid/liquid or solid/gas suspensions) in particular. It turns out that only very limited efforts have been made in the past to understand the fluid mechanics of such flows via advanced computer simulations. The purpose of this paper is to present the results of full three-dimensional (3D) theoretical and numerical analyses of single- and two-phase dilute particle/liquid flows in U-bend and helical curved conduits. The numerical analysis is based on computational fluid dynamics (CFD) simulations performed using a state-of-the-art multiphase flow computer code, NPHASE. The major issues discussed in the first part of the paper are concerned with the effect of curved/coiled geometry on the evolution of flow field and the associated wall shear. It has been demonstrated that the primary curvature (a common factor for both the U-bend and helix geometries) may cause a substantial asymmetry in the radial distribution of the main flow velocity. This, in turn, leads to a significant, albeit highly nonuniform, increase in the wall shear stress. Specifically, the wall shear around the outer half of tube circumference may become twice the corresponding value for a straight tube, and gradually decrease to

Quasi-3D Modelling of Two-Phase Slug Flow in Pipes

Directory of Open Access Journals (Sweden)

S. Mo

2014-03-01

Full Text Available In this paper, we present progress obtained by the Quasi 3-Dimensional (Q3D model for pipe flows. This model is based on a multi-fluid multi-field formulation with construction and tracking of the large-scale interfaces (LSIs. The computational time is significantly reduced compared to full 3D by using a specially adopted slice-averaging technique. However, the slice-averaging generates new terms in the model equations. These terms are related to important mechanisms such as wall shear stress and turbulence production at side walls. We present some basic performance tests of the Q3D model, including single phase wall friction and the velocities of single Taylor bubbles at inclinations ranging from horizontal to vertical. Finally we report the performance of the model for slug flow in horizontal and 10° inclined pipes. The model reproduces the experimental data satisfactorily for both cases in a very short simulation time compared to full 3D.

History independence of steady state in simultaneous two-phase flow through two-dimensional porous media.

Science.gov (United States)

Erpelding, Marion; Sinha, Santanu; Tallakstad, Ken Tore; Hansen, Alex; Flekkøy, Eirik Grude; Måløy, Knut Jørgen

2013-11-01

It is well known that the transient behavior during drainage or imbibition in multiphase flow in porous media strongly depends on the history and initial condition of the system. However, when the steady-state regime is reached and both drainage and imbibition take place at the pore level, the influence of the evolution history and initial preparation is an open question. Here, we present an extensive experimental and numerical work investigating the history dependence of simultaneous steady-state two-phase flow through porous media. Our experimental system consists of a Hele-Shaw cell filled with glass beads which we model numerically by a network of disordered pores transporting two immiscible fluids. From measurements of global pressure evolution, histograms of saturation, and cluster-size distributions, we find that when both phases are flowing through the porous medium, the steady state does not depend on the initial preparation of the system or on the way it has been reached.

Interaction between local parameters of two-phase flow and random forces on a cylinder

International Nuclear Information System (INIS)

Sylviane Pascal-Ribot; Yves Blanchet; Franck Baj; Phillippe Piteau

2005-01-01

Full text of publication follows: In the frame of assessments of steam generator tube bundle vibrations, a study was conducted in order to investigate the effects of an air/water flow on turbulent buffeting forces induced on a cylinder. The main purpose is to relate the physical parameters characterizing an air/water two-phase crossflow with the structural loading of a fixed cylindrical tube. In this first approach, the experiments are carried out in a rectangular acrylic test section supplied with a vertical upward bubbly flow. This flow is transversally impeded by a fixed rigid 12,15 mm diameter cylinder. Different turbulence grids are used in order to modify two-phase characteristics such as bubble diameter, void fraction profile, fluctuation parameters. Preliminarily, a dimensional analysis of fluid-structure interaction under two-phase turbulent solicitations has enabled to identify a list of physically relevant variables which must be measured to evaluate the random forces. The meaning of these relevant parameters as well as the effect of flow patterns are discussed. Direct measurements of two-phase flow parameters are performed simultaneously with measurements of forces exerted on the cylinder. The main descriptive parameters of a two-phase flow are measured using a bi-optical probe, in particular void fraction profiles, interfacial velocities, bubble diameters, void fraction fluctuations. In the same time, the magnitude of random forces caused by two-phase flow is measured with a force transducer. A thorough analysis of the experimental data is then undertaken in order to correlate physical two-phase mechanisms with the random forces exerted on the cylinder. The hypotheses made while applying the dimensional analysis are verified and their pertinence is discussed. Finally, physical parameters involved in random buffeting forces applied on a transverse tube are proposed to scale the spectral magnitude of these forces and comparisons with other authors

Surface representations of two- and three-dimensional fluid flow topology

Science.gov (United States)

Helman, James L.; Hesselink, Lambertus

1990-01-01

We discuss our work using critical point analysis to generate representations of the vector field topology of numerical flow data sets. Critical points are located and characterized in a two-dimensional domain, which may be either a two-dimensional flow field or the tangential velocity field near a three-dimensional body. Tangent curves are then integrated out along the principal directions of certain classes of critical points. The points and curves are linked to form a skeleton representing the two-dimensional vector field topology. When generated from the tangential velocity field near a body in a three-dimensional flow, the skeleton includes the critical points and curves which provide a basis for analyzing the three-dimensional structure of the flow separation. The points along the separation curves in the skeleton are used to start tangent curve integrations to generate surfaces representing the topology of the associated flow separations.

Numerical simulation of compressible two-phase flow using a diffuse interface method

International Nuclear Information System (INIS)

Ansari, M.R.; Daramizadeh, A.

2013-01-01

Highlights: ► Compressible two-phase gas–gas and gas–liquid flows simulation are conducted. ► Interface conditions contain shock wave and cavitations. ► A high-resolution diffuse interface method is investigated. ► The numerical results exhibit very good agreement with experimental results. -- Abstract: In this article, a high-resolution diffuse interface method is investigated for simulation of compressible two-phase gas–gas and gas–liquid flows, both in the presence of shock wave and in flows with strong rarefaction waves similar to cavitations. A Godunov method and HLLC Riemann solver is used for discretization of the Kapila five-equation model and a modified Schmidt equation of state (EOS) is used to simulate the cavitation regions. This method is applied successfully to some one- and two-dimensional compressible two-phase flows with interface conditions that contain shock wave and cavitations. The numerical results obtained in this attempt exhibit very good agreement with experimental results, as well as previous numerical results presented by other researchers based on other numerical methods. In particular, the algorithm can capture the complex flow features of transient shocks, such as the material discontinuities and interfacial instabilities, without any oscillation and additional diffusion. Numerical examples show that the results of the method presented here compare well with other sophisticated modeling methods like adaptive mesh refinement (AMR) and local mesh refinement (LMR) for one- and two-dimensional problems

Three dimensional turbulence structure measurements in air/water two phase flow

International Nuclear Information System (INIS)

Wang, S.K.L.

1986-01-01

The phenomena of turbulent air/water two phase upward and downward flows in a circular test section were investigated. Important flow quantities such as void fraction, liquid velocity, and Reynolds stresses were measured by using both single sensor and three sensor hot film probes. A digital data processing technique based on combined derivative and level thresholding was developed to determine the local void fraction from hot-film anemometer signals. The measured local void fraction was integrated and the result was compared with the chordal averaged void fraction measured by a gamma ray densitometer. It was found that the local measurement underestimated local void fraction due to surface tension effects and bubble deflection by the probe. A correlation based on local parameters characterizing probe/bubble interaction was developed, and it corrected the measured void fraction successfully. The measured void fraction profiles in upward flow and downward flow showed two distinct patterns. In upward flow, bubbles tend to migrate toward the wall and the void fraction profile shows a sharp peak near the wall. In downward flow, as the liquid velocity increases, the wall peaking phenomenon fades out and bubbles tend to migrate toward the center of the pipe

Simulation of horizontal pipe two-phase slug flows using the two-fluid model

Energy Technology Data Exchange (ETDEWEB)

Ortega Malca, Arturo J. [Pontificia Univ. Catolica do Rio de Janeiro, RJ (Brazil). Dept. de Engenharia Mecanica. Nucleo de Simulacao Termohidraulica de Dutos (SIMDUT); Nieckele, Angela O. [Pontificia Univ. Catolica do Rio de Janeiro, RJ (Brazil). Dept. de Engenharia Mecanica

2005-07-01

Slug flow occurs in many engineering applications, mainly in the transport of hydrocarbon fluids in pipelines. The intermittency of slug flow causes severe unsteady loading on the pipelines carrying the fluids, which gives rise to design problems. Therefore, it is important to be able to predict the onset and development of slug flow as well as slug characteristics. The present work consists in the simulation of two-phase flow in slug pattern through horizontal pipes using the two-fluid model in its transient and one-dimensional form. The advantage of this model is that the flow field is allowed to develop naturally from a given initial conditions as part of the transient calculation; the slug evolves automatically as a product of the computed flow development. Simulations are then carried out for a large number of flow conditions that lead a slug flow. (author)

Rheology of dense granular flows in two dimensions: Comparison of fully two-dimensional flows to unidirectional shear flow

Science.gov (United States)

Bhateja, Ashish; Khakhar, Devang V.

2018-06-01

We consider the rheology of steady two-dimensional granular flows, in different geometries, using discrete element method-based simulations of soft spheres. The flow classification parameter (ψ ), which defines the local flow type (ranging from pure rotation to simple shear to pure extension), varies spatially, to a significant extent, in the flows. We find that the material behaves as a generalized Newtonian fluid. The μ -I scaling proposed by Jop et al. [Nature (London) 441, 727 (2006), 10.1038/nature04801] is found to be valid in both two-dimensional and unidirectional flows, as observed in previous studies; however, the data for each flow geometry fall on a different curve. The results for the two-dimensional silo flow indicate that the viscosity does not depend directly on the flow type parameter, ψ . We find that the scaling based on "granular fluidity" [Zhang and Kamrin, Phys. Rev. Lett. 118, 058001 (2017), 10.1103/PhysRevLett.118.058001] gives good collapse of the data to a single curve for all the geometries. The data for the variation of the solid faction with inertial number show a reasonable collapse for the different geometries.

Multi-dimensional Mixing Behavior of Steam-Water Flow in a Downcomer Annulus during LBLOCA Reflood Phase with a DVI Injection Mode

International Nuclear Information System (INIS)

Kwon, T.S.; Yun, B.J.; Euh, D.J.; Chu, I.C.; Song, C.H.

2002-01-01

Multi-dimensional thermal-hydraulic behavior in the downcomer annulus of a pressurized water reactor vessel with a Direct Vessel Injection (DVI) mode is presented based on the experimental observation in the MIDAS (Multi-dimensional Investigation in Downcomer Annulus Simulation) steam-water test facility. From the steady-state test results to simulate the late reflood phase of a Large Break Loss-of-Coolant Accidents(LBLOCA), isothermal lines show the multidimensional phenomena of a phasic interaction between steam and water in the downcomer annulus very well. MIDAS is a steam-water separate effect test facility, which is 1/4.93 linearly scaled-down of 1400 MWe PWR type of a nuclear reactor, focused on understanding multi-dimensional thermalhydraulic phenomena in downcomer annulus with various types of safety injection during the refill or reflood phase of a LBLOCA. The initial and the boundary conditions are scaled from the pre-test analysis based on the preliminary calculation using the TRAC code. The superheated steam with a superheating degree of 80 K at a given downcomer pressure of 180 kPa is injected equally through three intact cold legs into the downcomer. (authors)

Modeling of annular two-phase flow using a unified CFD approach

Energy Technology Data Exchange (ETDEWEB)

Li, Haipeng, E-mail: haipengl@kth.se; Anglart, Henryk, E-mail: henryk@kth.se

2016-07-15

Highlights: • Annular two-phase flow has been modeled using a unified CFD approach. • Liquid film was modeled based on a two-dimensional thin film assumption. • Both Eulerian and Lagrangian methods were employed for the gas core flow modeling. - Abstract: A mechanistic model of annular flow with evaporating liquid film has been developed using computational fluid dynamics (CFD). The model is employing a separate solver with two-dimensional conservation equations to predict propagation of a thin boiling liquid film on solid walls. The liquid film model is coupled to a solver of three-dimensional conservation equations describing the gas core, which is assumed to contain a saturated mixture of vapor and liquid droplets. Both the Eulerian–Eulerian and the Eulerian–Lagrangian approach are used to describe the droplet and vapor motion in the gas core. All the major interaction phenomena between the liquid film and the gas core flow have been accounted for, including the liquid film evaporation as well as the droplet deposition and entrainment. The resultant unified framework for annular flow has been applied to the steam-water flow with conditions typical for a Boiling Water Reactor (BWR). The simulation results for the liquid film flow rate show good agreement with the experimental data, with the potential to predict the dryout occurrence based on criteria of critical film thickness or critical film flow rate.

Modeling of annular two-phase flow using a unified CFD approach

International Nuclear Information System (INIS)

Li, Haipeng; Anglart, Henryk

2016-01-01

Highlights: • Annular two-phase flow has been modeled using a unified CFD approach. • Liquid film was modeled based on a two-dimensional thin film assumption. • Both Eulerian and Lagrangian methods were employed for the gas core flow modeling. - Abstract: A mechanistic model of annular flow with evaporating liquid film has been developed using computational fluid dynamics (CFD). The model is employing a separate solver with two-dimensional conservation equations to predict propagation of a thin boiling liquid film on solid walls. The liquid film model is coupled to a solver of three-dimensional conservation equations describing the gas core, which is assumed to contain a saturated mixture of vapor and liquid droplets. Both the Eulerian–Eulerian and the Eulerian–Lagrangian approach are used to describe the droplet and vapor motion in the gas core. All the major interaction phenomena between the liquid film and the gas core flow have been accounted for, including the liquid film evaporation as well as the droplet deposition and entrainment. The resultant unified framework for annular flow has been applied to the steam-water flow with conditions typical for a Boiling Water Reactor (BWR). The simulation results for the liquid film flow rate show good agreement with the experimental data, with the potential to predict the dryout occurrence based on criteria of critical film thickness or critical film flow rate.

Thermo-fluid dynamics of two-phase flow

CERN Document Server

Ishii, Mamoru; Ishii, Mamoru; Ishii, M

2006-01-01

Provides a very systematic treatment of two phase flow problems from a theoretical perspectiveProvides an easy to follow treatment of modeling and code devlopemnt of two phase flow related phenomenaCovers new results of two phase flow research such as coverage of fuel cells technology.

Development of a large-scale general purpose two-phase flow analysis code

International Nuclear Information System (INIS)

Terasaka, Haruo; Shimizu, Sensuke

2001-01-01

A general purpose three-dimensional two-phase flow analysis code has been developed for solving large-scale problems in industrial fields. The code uses a two-fluid model to describe the conservation equations for two-phase flow in order to be applicable to various phenomena. Complicated geometrical conditions are modeled by FAVOR method in structured grid systems, and the discretization equations are solved by a modified SIMPLEST scheme. To reduce computing time a matrix solver for the pressure correction equation is parallelized with OpenMP. Results of numerical examples show that the accurate solutions can be obtained efficiently and stably. (author)

Shroud leakage flow models and a multi-dimensional coupling CFD (computational fluid dynamics) method for shrouded turbines

International Nuclear Information System (INIS)

Zou, Zhengping; Liu, Jingyuan; Zhang, Weihao; Wang, Peng

2016-01-01

Multi-dimensional coupling simulation is an effective approach for evaluating the flow and aero-thermal performance of shrouded turbines, which can balance the simulation accuracy and computing cost effectively. In this paper, 1D leakage models are proposed based on classical jet theories and dynamics equations, which can be used to evaluate most of the main features of shroud leakage flow, including the mass flow rate, radial and circumferential momentum, temperature and the jet width. Then, the 1D models are expanded to 2D distributions on the interface by using a multi-dimensional scaling method. Based on the models and multi-dimensional scaling, a multi-dimensional coupling simulation method for shrouded turbines is developed, in which, some boundary source and sink are set on the interface between the shroud and the main flow passage. To verify the precision, some simulations on the design point and off design points of a 1.5 stage turbine are conducted. It is indicated that the models and methods can give predictions with sufficient accuracy for most of the flow field features and will contribute to pursue deeper understanding and better design methods of shrouded axial turbines, which are the important devices in energy engineering. - Highlights: • Free and wall attached jet theories are used to model the leakage flow in shrouds. • Leakage flow rate is modeled by virtual labyrinth number and residual-energy factor. • A scaling method is applied to 1D model to obtain 2D distributions on interfaces. • A multi-dimensional coupling CFD method for shrouded turbines is proposed. • The proposed coupling method can give accurate predictions with low computing cost.

Two-phase flow in refrigeration systems

CERN Document Server

Gu, Junjie; Gan, Zhongxue

2013-01-01

Two-Phase Flow in Refrigeration Systems presents recent developments from the authors' extensive research programs on two-phase flow in refrigeration systems. This book covers advanced mass and heat transfer and vapor compression refrigeration systems and shows how the performance of an automotive air-conditioning system is affected through results obtained experimentally and theoretically, specifically with consideration of two-phase flow and oil concentration. The book is ideal for university postgraduate students as a textbook, researchers and professors as an academic reference book, and b

Two-group interfacial area concentration correlations of two-phase flows in large diameter pipes

International Nuclear Information System (INIS)

Shen, Xiuzhong; Hibiki, Takashi

2015-01-01

The reliable empirical correlations and models are one of the important ways to predict the interfacial area concentration (IAC) in two-phase flows. However, up to now, no correlation or model is available for the prediction of the IAC in the two-phase flows in large diameter pipes. This study collected an IAC experimental database of two-phase flows taken under various flow conditions in large diameter pipes and presented a systematic way to predict the IAC for two-phase flows from bubbly, cap-bubbly to churn flow in large diameter pipes by categorizing bubbles into two groups (group-1: spherical and distorted bubble, group-2: cap bubble). Correlations were developed to predict the group-1 void fraction from the void fraction of all bubble. The IAC contribution from group-1 bubbles was modeled by using the dominant parameters of group-1 bubble void fraction and Reynolds number based on the parameter-dependent analysis of Hibiki and Ishii (2001, 2002) using one-dimensional bubble number density and interfacial area transport equations. A new drift velocity correlation for two-phase flow with large cap bubbles in large diameter pipes was derived in this study. By comparing the newly-derived drift velocity correlation with the existing drift velocity correlation of Kataoka and Ishii (1987) for large diameter pipes and using the characteristics of the representative bubbles among the group 2 bubbles, we developed the model of IAC and bubble size for group 2 cap bubbles. The developed models for estimating the IAC are compared with the entire collected database. A reasonable agreement was obtained with average relative errors of ±28.1%, ±54.4% and ±29.6% for group 1, group 2 and all bubbles respectively. (author)

Refrigeration. Two-Phase Flow. Flow Regimes and Pressure Drop

DEFF Research Database (Denmark)

Knudsen, Hans-Jørgen Høgaard

2002-01-01

The note gives the basic definitions used in two-phase flow. Flow regimes and flow regimes map are introduced. The different contributions to the pressure drop are stated together with an imperical correlation from the litterature.......The note gives the basic definitions used in two-phase flow. Flow regimes and flow regimes map are introduced. The different contributions to the pressure drop are stated together with an imperical correlation from the litterature....

Theoretical investigations on two-phase flow instability in parallel channels under axial non-uniform heating

International Nuclear Information System (INIS)

Lu, Xiaodong; Wu, Yingwei; Zhou, Linglan; Tian, Wenxi; Su, Guanghui; Qiu, Suizheng; Zhang, Hong

2014-01-01

Highlights: • We developed a model based on homogeneous flow model to analyze two-phase flow instability in parallel channels. • The influence of axial non-uniform heating on the system stability has been investigated. • Influences of various factors on system instability under cosine heat flux have been studied. • The system under top-peaked heat flux is the most stable system. - Abstract: Two-phase flow instability in parallel channels heated by axial non-uniform heat flux has been theoretically studied in this paper. The system control equations of parallel channels were established based on the homogeneous flow model in two-phase region. Semi-implicit finite-difference scheme and staggered mesh method were used to discretize the equations, and the difference equations were solved by chasing method. Cosine, bottom-peaked and top-peaked heat fluxes were used to study the influence of non-uniform heating on two-phase flow instability of the parallel channels system. The marginal stability boundaries (MSB) of parallel channels and three-dimensional instability spaces (or instability reefs) under different heat flux conditions have been obtained. Compared with axial uniform heating, axial non-uniform heating will affect the system stability. Cosine and bottom-peaked heat fluxes can destabilize the system stability in high inlet subcooling region, while the opposite effect can be found in low inlet subcooling region. However, top-peaked heat flux can enhance the system stability in the whole region. In addition, for cosine heat flux, increasing the system pressure or inlet resistance coefficient can strengthen the system stability, and increasing the heating power will destabilize the system stability. The influence of inlet subcooling number on the system stability is multi-valued under cosine heat flux

Three-dimensional lattice Boltzmann model for immiscible two-phase flow simulations.

Science.gov (United States)

Liu, Haihu; Valocchi, Albert J; Kang, Qinjun

2012-04-01

We present an improved three-dimensional 19-velocity lattice Boltzmann model for immisicible binary fluids with variable viscosity and density ratios. This model uses a perturbation step to generate the interfacial tension and a recoloring step to promote phase segregation and maintain surfaces. A generalized perturbation operator is derived using the concept of a continuum surface force together with the constraints of mass and momentum conservation. A theoretical expression for the interfacial tension is determined directly without any additional analysis and assumptions. The recoloring algorithm proposed by Latva-Kokko and Rothman is applied for phase segregation, which minimizes the spurious velocities and removes lattice pinning. This model is first validated against the Laplace law for a stationary bubble. It is found that the interfacial tension is predicted well for density ratios up to 1000. The model is then used to simulate droplet deformation and breakup in simple shear flow. We compute droplet deformation at small capillary numbers in the Stokes regime and find excellent agreement with the theoretical Taylor relation for the segregation parameter β=0.7. In the limit of creeping flow, droplet breakup occurs at a critical capillary number 0.35

Industrial aspects of gas-liquid two-phase flow

International Nuclear Information System (INIS)

Hewitt, G.F.

1977-01-01

The lecture begins by reviewing the various types of plant in which two phase flow occurs. Specifically, boiling plant, condensing plant and pipelines are reviewed, and the various two phase flow problems occurring in them are described. Of course, many other kinds of chemical engineering plant involve two phase flow, but are somewhat outside the scope of this lecture. This would include distillation columns, vapor-liquid separators, absorption towers etc. Other areas of industrial two phase flow which have been omitted for space reasons from this lecture are those concerned with gas/solids, liquid/solid and liquid/liquid flows. There then follows a description of some of the two phase flow processes which are relevant in industrial equipment and where special problems occur. The topics chosen are as follows: (1) pressure drop; (2) horizontal tubes - separation effects non-uniformites in heat transfer coefficient, effect of bends on dryout; (3) multicomponent mixtures - effects in pool boiling, mass transfer effects in condensation and Marangoni effects; (4) flow distribution - manifold problems in single phase flow, separation effects at a single T-junction in two phase flow and distribution in manifolds in two phase flow; (5) instability - oscillatory instability, special forms of instability in cryogenic systems; (6) nucleate boiling - effect of variability of surface, unresolved problems in forced convective nucleate boiling; and (7) shell side flows - flow patterns, cross flow boiling, condensation in cross flow

Theoretical modelling and experimental investigation of single-phase and two-phase flow division at a tee-junction

International Nuclear Information System (INIS)

Lemonnier, H.; Hervieu, E.

1991-01-01

Phase separation in a tee-junction is modelled in the particular case of bubbly-flow. The model is based on a two-dimensional approach and hence, uses local equations. The first step consists in modelling the single-phase flow in the tee-junction. The free streamline theory is used to predict the flow of the continuous phase. The two recirculation zones which are presented in this case are predicted by the model. The second step consists in predicting the gas bubble paths as a result of the actions of the single-phase flow. Finally, the trajectories of gas bubbles are used to predict the separation characteristics of the tee-junction. Each step of the modelling procedure has been carefully tested by an in-depth experimental investigation. Excellent quantitative agreement is obtained between experimental results and model predictions. Moreover, the phase separation phenomenon is found to be clearly described by the model. (orig.)

The use of magnetic resonance imaging to quantify multi-phase flow patterns and transitions

International Nuclear Information System (INIS)

Reyes, J.N. Jr.; Lafi, A.Y.; Saloner, D.

1998-01-01

Conventional measurement techniques have given limited insights into the complex structure of multi-phase flows. This has led to highly subjective flow pattern classifications which have been cast in terms of flow regime maps. Rather than using static flow regime maps, some of the next generation of multi-phase flow analysis codes will implement interfacial area transport equations that would calculate the flow patterns that evolve spatially and temporally. To assess these new codes, a large data base needs to be established to quantify the essential characteristics of multi-phase flow structure. One such characteristic is the interfacial area concentration. In this paper, we discuss the current benefits and limitations of using Magnetic Resonance Imaging (MRI) to examine multi- phase flow patterns and transitions. Of particular interest, are the MRI measurements of interfacial area concentration for slug flow in an air-water system. These tests were performed at the University of California, San Francisco (UCSF) School of Medicine MRI Center as a collaborative research effort with Oregon State University (OSU). The special scanning sequences designed by UCSF were capable of imaging at repetition intervals as fast as 7 milliseconds. (author)

The use of magnetic resonance imaging to quantify multi-phase flow patterns and transitions

Energy Technology Data Exchange (ETDEWEB)

Reyes, Jr, J N; Lafi, A Y [Department of Nuclear Engineering, Oregon State University, Corvallis, OR (United States); Saloner, D [University of California, San Francisco School of Medicine, Veterans Administration Medical Center, San Francisco, CA (United States)

1998-09-01

Conventional measurement techniques have given limited insights into the complex structure of multi-phase flows. This has led to highly subjective flow pattern classifications which have been cast in terms of flow regime maps. Rather than using static flow regime maps, some of the next generation of multi-phase flow analysis codes will implement interfacial area transport equations that would calculate the flow patterns that evolve spatially and temporally. To assess these new codes, a large data base needs to be established to quantify the essential characteristics of multi-phase flow structure. One such characteristic is the interfacial area concentration. In this paper, we discuss the current benefits and limitations of using Magnetic Resonance Imaging (MRI) to examine multi- phase flow patterns and transitions. Of particular interest, are the MRI measurements of interfacial area concentration for slug flow in an air-water system. These tests were performed at the University of California, San Francisco (UCSF) School of Medicine MRI Center as a collaborative research effort with Oregon State University (OSU). The special scanning sequences designed by UCSF were capable of imaging at repetition intervals as fast as 7 milliseconds. (author)

Transient analysis of air-water two-phase flow in channels and bends

International Nuclear Information System (INIS)

Khan, H.J.; Ye, W.; Pertmer, G.A.

1992-01-01

The algorithm used in this paper is the Newton Block Gauss Seidel method, which has been applied to both simple and complex flow conditions in two-phase flow. This paper contains a description of difference techniques and an iterative solution algorithm that is used to solve the field and constitutive equations of the two-fluid model. In practice, this solution procedure has been proven to be stable and capable of generating solutions in problems where other schemes have failed. The method converges rapidly for reasonable error tolerances and is easily extended to three-dimensional geometries. Using air-water as the two-phase medium, transient flow behavior in several geometries of interest are shown. Flow through a vertical channel with flow obstruction, large U bends, and 90-deg bends are being demonstrated with variation of inlet void fraction and slip ratio. Significant changes in the velocity and void distribution profiles have been observed. Various regions of flow recirculation are obtained in the flow domain for each phase. The phasic velocity and void distributions are dominated by gravity-induced phase separation causing air to accumulate in the upper region. The influence of inlet slip ratio and interfacial momentum transfer on the transient flow profile has been demonstrated in detail

Application of the level set method for multi-phase flow computation in fusion engineering

International Nuclear Information System (INIS)

Luo, X-Y.; Ni, M-J.; Ying, A.; Abdou, M.

2006-01-01

Numerical simulation of multi-phase flow is essential to evaluate the feasibility of a liquid protection scheme for the power plant chamber. The level set method is one of the best methods for computing and analyzing the motion of interface among the multi-phase flow. This paper presents a general formula for the second-order projection method combined with the level set method to simulate unsteady incompressible multi-phase flow with/out phase change flow encountered in fusion science and engineering. The third-order ENO scheme and second-order semi-implicit Crank-Nicholson scheme is used to update the convective and diffusion term. The numerical results show this method can handle the complex deformation of the interface and the effect of liquid-vapor phase change will be included in the future work

Development of multidimensional two-phase flow measurement sensor in rod bundle

International Nuclear Information System (INIS)

Arai, Takahiro; Furuya, Masahiro; Shirakawa, Kenetsu; Kanai, Taizo

2011-01-01

In order to acquire multidimensional two-phase flow in 10x10 bundle, SubChannel Void Sensor (SCVC) consisting of 11-wire by 11-wire and 10-rod by 10-rod electrodes is developed. A conductance value in a proximity region of one wire and another gives void fraction in the center of subchannel region. A phasic velocity can be estimated by using two layers of wire meshes, like as so-called wire mesh sensor. 121 points (=11x11) of void fraction as well as those of phasic velocity are acquired. It is peculiarity of the devised sensor that void fraction near rod surface can be estimated by a conductance value in a proximity region of one wire and one rod. 400 additional points of void fraction in 10x10 bundle can be, therefore, acquired. The time resolution of measurement is up to 1250 frames (cross sections) per second. We capability in a 10x10 bundle with o.d. 10 mm and 3110 mm long is demonstrated. The devised sensor is installed in 8 height levels to acquire the two-phase flow dynamics along axial direction. A pair of sensor layers is mounted in each level and is placed by 30 mm apart with each other to estimate a phasic velocity distribution on the basis of cross-correlation function of the two layers. Air bubbles are injected through sintered metal nozzles from the bottom end of 10x10 rods. Air flow rate distribution can vary with a controlled valves connected to each nozzle. The devised sensor exhibited the quasi three-dimensional flow structures, i.e. void fraction, phasic velocity and bubble chord length distributions. These quasi three-dimensional structures explorer complexity of two-phase flow dynamics such as coalescence and breakup of bubbles in the transient phasic velocity distributions. (author)

Development and assessment of Multi-dimensional flow models in the thermal-hydraulic system analysis code MARS

Energy Technology Data Exchange (ETDEWEB)

Chung, B. D.; Bae, S. W.; Jeong, J. J.; Lee, S. M

2005-04-15

A new multi-dimensional component has been developed to allow for more flexible 3D capabilities in the system code, MARS. This component can be applied in the Cartesian and cylindrical coordinates. For the development of this model, the 3D convection and diffusion terms are implemented in the momentum and energy equation. And a simple Prandtl's mixing length model is applied for the turbulent viscosity. The developed multi-dimensional component was assessed against five conceptual problems with analytic solution. And some SETs are calculated and compared with experimental data. With this newly developed multi-dimensional flow module, the MARS code can realistic calculate the flow fields in pools such as those occurring in the core, steam generators and IRWST.

Development and assessment of Multi-dimensional flow models in the thermal-hydraulic system analysis code MARS

International Nuclear Information System (INIS)

Chung, B. D.; Bae, S. W.; Jeong, J. J.; Lee, S. M.

2005-04-01

A new multi-dimensional component has been developed to allow for more flexible 3D capabilities in the system code, MARS. This component can be applied in the Cartesian and cylindrical coordinates. For the development of this model, the 3D convection and diffusion terms are implemented in the momentum and energy equation. And a simple Prandtl's mixing length model is applied for the turbulent viscosity. The developed multi-dimensional component was assessed against five conceptual problems with analytic solution. And some SETs are calculated and compared with experimental data. With this newly developed multi-dimensional flow module, the MARS code can realistic calculate the flow fields in pools such as those occurring in the core, steam generators and IRWST

Some issues in the simulation of two-phase flows: The relative velocity

International Nuclear Information System (INIS)

Gräbel, J.; Hensel, S.; Ueberholz, P.; Farber, P.; Zeidan, D.

2016-01-01

In this paper we compare numerical approximations for solving the Riemann problem for a hyperbolic two-phase flow model in two-dimensional space. The model is based on mixture parameters of state where the relative velocity between the two-phase systems is taken into account. This relative velocity appears as a main discontinuous flow variable through the complete wave structure and cannot be recovered correctly by some numerical techniques when simulating the associated Riemann problem. Simulations are validated by comparing the results of the numerical calculation qualitatively with OpenFOAM software. Simulations also indicate that OpenFOAM is unable to resolve the relative velocity associated with the Riemann problem.

Some issues in the simulation of two-phase flows: The relative velocity

Energy Technology Data Exchange (ETDEWEB)

Gräbel, J.; Hensel, S.; Ueberholz, P.; Farber, P. [Niederrhein University of Applied Sciences, Institute for Modelling and High Performance Computing, Reinarzstraße 49, 47805 Krefeld (Germany); Zeidan, D. [School of Basic Sciences and Humanities, German Jordanian University, Amman (Jordan)

2016-06-08

In this paper we compare numerical approximations for solving the Riemann problem for a hyperbolic two-phase flow model in two-dimensional space. The model is based on mixture parameters of state where the relative velocity between the two-phase systems is taken into account. This relative velocity appears as a main discontinuous flow variable through the complete wave structure and cannot be recovered correctly by some numerical techniques when simulating the associated Riemann problem. Simulations are validated by comparing the results of the numerical calculation qualitatively with OpenFOAM software. Simulations also indicate that OpenFOAM is unable to resolve the relative velocity associated with the Riemann problem.

A review of damping of two-phase flows

International Nuclear Information System (INIS)

Hara, Fumio

1993-01-01

Damping of two-phase flows has been recognized as one of the most unknown parameters in analyzing vibrational characteristics of structures subjected to two-phase flows since it seems to be influenced by many physical parameters involved in the physics of dynamic energy dissipation of a vibrating structure, for example, liquid viscosity, surface tension, flow velocity, mass ratio, frequency, void fraction, flow regime and so forth. This paper deals with a review of scientific works done to date on the damping of two phase flows and discussions about what has been clarified and what has not been known to us, or what kinds of research are needed about two-phase flow damping. The emphasis is put on the definition of two-phase fluid damping, damping measurement techniques, damping characteristics in relation to two phase flow configurations, and damping generation mechanisms

Extension of CFD Codes Application to Two-Phase Flow Safety Problems - Phase 3

International Nuclear Information System (INIS)

Bestion, D.; Anglart, H.; Mahaffy, J.; Lucas, D.; Song, C.H.; Scheuerer, M.; Zigh, G.; Andreani, M.; Kasahara, F.; Heitsch, M.; Komen, E.; Moretti, F.; Morii, T.; Muehlbauer, P.; Smith, B.L.; Watanabe, T.

2014-11-01

. Identification of needs for additional experimental validation; - identification of a matrix of numerical benchmarks of special interest for the selected NRS problems; - establish the foundation of Best Practice Guidelines for two-phase CFD application to the selected NRS problems. This document is an extension of the first report produced by the Writing Group 3. A few NRS problems where two-phase CFD may bring real benefit are first selected to be further analysed in more detail according to some criteria. They must be high priority issues from the point of view of nuclear safety with some investigations going on and CFD investigations must have a reasonable chance to be successful in a reasonable period of time. They must address both the present generation of PWR and BWR and the Generation 3 water reactors and should possibly address all flow regimes so that they may, to some extent, envelop many other issues. A general multi-step methodology for application of two-phase CFD to nuclear safety issues is proposed. Many options are possible when using two phase CFD, for the basic model (one-fluid, two-fluid, multi-field,...), for the averaging or filtering of turbulent and two-phase scales (using RANS, LES,...), for the treatment of the interface either by an Interface Tracking Method or statistically by a volume fraction, an interfacial area equation,... The choices have to be justified after an in-depth analysis of the issue and an identification of all basic flow processes. Then closure relations have to be selected or developed for interfacial transfers, turbulent transfers and wall transfers, and a test matrix has to be established to validate all the models in a separate effect way. The method has been applied to the selected issues resulting in an updated state of the art and gaps were identified in the modelling. Available data for validation were reviewed and needs of additional data were identified. Verification tests were also identified. A few benchmarks are proposed

Study of the Riemann problem and construction of multidimensional Godunov-type schemes for two-phase flow models

International Nuclear Information System (INIS)

Toumi, I.

1990-04-01

This thesis is devoted to the study of the Riemann problem and the construction of Godunov type numerical schemes for one or two dimensional two-phase flow models. In the first part, we study the Riemann problem for the well-known Drift-Flux, model which has been widely used for the analysis of thermal hydraulics transients. Then we use this study to construct approximate Riemann solvers and we describe the corresponding Godunov type schemes for simplified equation of state. For computation of complex two-phase flows, a weak formulation of Roe's approximate Riemann solver, which gives a method to construct a Roe-averaged jacobian matrix with a general equation of state, is proposed. For two-dimensional flows, the developed methods are based upon an approximate solver for a two-dimensional Riemann problem, according to Harten-Lax-Van Leer principles. The numerical results for standard test problems show the good behaviour of these numerical schemes for a wide range of flow conditions [fr

Two-phase flow in fractured rock

International Nuclear Information System (INIS)

Davies, P.; Long, J.; Zuidema, P.

1993-11-01

This report gives the results of a three-day workshop on two-phase flow in fractured rock. The workshop focused on two-phase flow processes that are important in geologic disposal of nuclear waste as experienced in a variety of repository settings. The goals and objectives of the workshop were threefold: exchange information; describe the current state of understanding; and identify research needs. The participants were divided into four subgroups. Each group was asked to address a series of two-phase flow processes. The following groups were defined to address these processes: basic flow processes; fracture/matrix interactions; complex flow processes; and coupled processes. For each process, the groups were asked to address these four issues: (1) describe the two-phase flow processes that are important with respect to repository performance; (2) describe how this process relates to the specific driving programmatic issues given above for nuclear waste storage; (3) evaluate the state of understanding for these processes; and (4) suggest additional research to address poorly understood processes relevant to repository performance. The reports from each of the four working groups are given here

Characterization of two-phase flow regimes in horizontal tubes using 81mKr tracer experiments.

Science.gov (United States)

Oriol, Jean; Leclerc, Jean Pierre; Berne, Philippe; Gousseau, Georges; Jallut, Christian; Tochon, Patrice; Clement, Patrice

2008-10-01

The diagnosis of heat exchangers on duty with respect to flow mal-distributions needs the development of non-intrusive inlet-outlet experimental techniques in order to perform an online fault diagnosis. Tracer experiments are an example of such techniques. They can be applied to mono-phase heat exchangers but also to multi-phase ones. In this case, the tracer experiments are more difficult to perform. In order to check for the capabilities of tracer experiments to be used for the flow mal-distribution diagnosis in the case of multi-phase heat exchangers, we present here a preliminary study on the simplest possible system: two-phase flows in a horizontal tube. (81m)Kr is used as gas tracer and properly collimated NaI (TI) crystal scintillators as detectors. The specific shape of the tracer response allows two-phase flow regimes to be characterized. Signal analysis allows the estimation of the gas phase real average velocity and consequently of the liquid phase real average velocity as well as of the volumetric void fraction. These results are compared successfully to those obtained with liquid phase tracer experiments previously presented by Oriol et al. 2007. Characterization of the two-phase flow regimes and liquid dispersion in horizontal and vertical tubes using coloured tracer and no intrusive optical detector. Chem. Eng. Sci. 63(1), 24-34, as well as to those given by correlations from literature.

Characterization of two-phase flow regimes in horizontal tubes using 81mKr tracer experiments

International Nuclear Information System (INIS)

Oriol, Jean; Leclerc, Jean Pierre; Berne, Philippe; Gousseau, Georges; Jallut, Christian; Tochon, Patrice; Clement, Patrice

2008-01-01

The diagnosis of heat exchangers on duty with respect to flow mal-distributions needs the development of non-intrusive inlet-outlet experimental techniques in order to perform an online fault diagnosis. Tracer experiments are an example of such techniques. They can be applied to mono-phase heat exchangers but also to multi-phase ones. In this case, the tracer experiments are more difficult to perform. In order to check for the capabilities of tracer experiments to be used for the flow mal-distribution diagnosis in the case of multi-phase heat exchangers, we present here a preliminary study on the simplest possible system: two-phase flows in a horizontal tube. 81m Kr is used as gas tracer and properly collimated NaI (TI) crystal scintillators as detectors. The specific shape of the tracer response allows two-phase flow regimes to be characterized. Signal analysis allows the estimation of the gas phase real average velocity and consequently of the liquid phase real average velocity as well as of the volumetric void fraction. These results are compared successfully to those obtained with liquid phase tracer experiments previously presented by Oriol et al. 2007. Characterization of the two-phase flow regimes and liquid dispersion in horizontal and vertical tubes using coloured tracer and no intrusive optical detector. Chem. Eng. Sci. 63(1), 24-34, as well as to those given by correlations from literature

Two-fluid model of two-phase flow in a pin bundle of a nuclear reactor

International Nuclear Information System (INIS)

Chawla, T.C.; Ishii, M.

1980-01-01

By considering two-phase flow as a field which is subdivided into two turbulent single-phase regions with moving boundaries separating the two constituent phases, such that the differential balances for three-dimensional turbulent flow hold for each subregion and for the interface, we perform the Eulerian area averaging over the cross-sectional area of each phase in a given channel and segment averaging of transverse momentum equation along the phase intercepts at the interchannel boundaries. To simplify the governing equations obtained as a result of these operations, we invoke the assumption that the motion of the fluid in each phase is dominantly in axial direction, that is the transverse components of velocity are small compared to axial components. We further assume that the variation of axial component of velocity within a channel is much stronger than the variation along the axial direction. We also assume that similar arguments can also be applied to the variation of enthalpy in a channel. As a result of these considerations, we obtain two sets of continuity, momentum, and energy equations describing motion of each phase in the axial direction. The phasic interaction terms which appear in these equations are governed by interfacial transfer conditions obtained from interface balances. The segment-averaged transverse-momentum equation for each phase provides the governing equation for cross flow. (author)

Experimental observation of a multi-dimensional mixing behavior of steam-water flow in the MIDAS test facility

International Nuclear Information System (INIS)

Kweon, T. S.; Yun, B. J.; Ah, D. J.; Ju, I. C.; Song, C. H.; Park, J. K.

2001-01-01

Multi-dimensional thermal-hydraulic hehavior, such as ECC (Emergency Core Cooling) bypass, ECC penetration, steam-water condensation and accumulated water level, in an annular downcomer of a PWR (Pressurized Water Reactor) reactor vessel with a DVI(Direct Vessel Injection) injection mode is presented based on the experimental observations in the MIDAS (Multi-dimensional Investigation in Downcomer Annulus Simulation) steam-water facility. From the steady-state tests to similate a late reflood phase of LBLOCA (Large Break Loss-of-Coolant Accidents), major thermal-hydraulic phenomena in the downcomer are quantified under a wide range of test conditions. Especially, isothermal lines show well multi-dimensional phenomena of phase interaction between steam and water in the annulus downcomer. Overall test results show that multi-dimensional thermal-hydraulic behaviors occur in the downcomer annulus region as expected. The MIDAS test facility is a steam-water separate effect test facility, which is 1/4.93 linearly scaled-down of a 1400 MWe PWR type of nuclear reactor, with focusing on understanding multi-dimensional thermal-hydraulic phenomena in annulus downcomer with various types of safety injection location during refill or reflood phase of a LBLOCA in PWR

Unsteady interfacial coupling of two-phase flow models

International Nuclear Information System (INIS)

Hurisse, O.

2006-01-01

The primary coolant circuit in a nuclear power plant contains several distinct components (vessel, core, pipes,...). For all components, specific codes based on the discretization of partial differential equations have already been developed. In order to obtain simulations for the whole circuit, the interfacial coupling of these codes is required. The approach examined within this work consists in coupling codes by providing unsteady information through the coupling interface. The numerical technique relies on the use of an interface model, which is combined with the basic strategy that was introduced by Greenberg and Leroux in order to compute approximations of steady solutions of non-homogeneous hyperbolic systems. Three different coupling cases have been examined: (i) the coupling of a one-dimensional Euler system with a two-dimensional Euler system; (ii) the coupling of two distinct homogeneous two-phase flow models; (iii) the coupling of a four-equation homogeneous model with the standard two-fluid model. (author)

Microgravity Two-Phase Flow Transition

Science.gov (United States)

Parang, M.; Chao, D.

1999-01-01

Two-phase flows under microgravity condition find a large number of important applications in fluid handling and storage, and spacecraft thermal management. Specifically, under microgravity condition heat transfer between heat exchanger surfaces and fluids depend critically on the distribution and interaction between different fluid phases which are often qualitatively different from the gravity-based systems. Heat transfer and flow analysis in two-phase flows under these conditions require a clear understanding of the flow pattern transition and development of appropriate dimensionless scales for its modeling and prediction. The physics of this flow is however very complex and remains poorly understood. This has led to various inadequacies in flow and heat transfer modeling and has made prediction of flow transition difficult in engineering design of efficient thermal and flow systems. In the present study the available published data for flow transition under microgravity condition are considered for mapping. The transition from slug to annular flow and from bubbly to slug flow are mapped using dimensionless variable combination developed in a previous study by the authors. The result indicate that the new maps describe the flow transitions reasonably well over the range of the data available. The transition maps are examined and the results are discussed in relation to the presumed balance of forces and flow dynamics. It is suggested that further evaluation of the proposed flow and transition mapping will require a wider range of microgravity data expected to be made available in future studies.

Modeling of Two-Phase Flow in Rough-Walled Fracture Using Level Set Method

Directory of Open Access Journals (Sweden)

Yunfeng Dai

2017-01-01

Full Text Available To describe accurately the flow characteristic of fracture scale displacements of immiscible fluids, an incompressible two-phase (crude oil and water flow model incorporating interfacial forces and nonzero contact angles is developed. The roughness of the two-dimensional synthetic rough-walled fractures is controlled with different fractal dimension parameters. Described by the Navier–Stokes equations, the moving interface between crude oil and water is tracked using level set method. The method accounts for differences in densities and viscosities of crude oil and water and includes the effect of interfacial force. The wettability of the rough fracture wall is taken into account by defining the contact angle and slip length. The curve of the invasion pressure-water volume fraction is generated by modeling two-phase flow during a sudden drainage. The volume fraction of water restricted in the rough-walled fracture is calculated by integrating the water volume and dividing by the total cavity volume of the fracture while the two-phase flow is quasistatic. The effect of invasion pressure of crude oil, roughness of fracture wall, and wettability of the wall on two-phase flow in rough-walled fracture is evaluated.

Numerical study for two phase flow in the near nozzle region of turbine combustors

International Nuclear Information System (INIS)

Pervez, K.; Mushtaq, S.

1999-01-01

In the present study flow conditions in the near nozzle region of the combustion chamber have been investigated. There exists two-phase flow in this region. The overall performance and pollutant formation in the combustion chamber have been investigated. There exists two-phase flow in this region. The overall performance and pollutant formation in the combustion zone largely depends on the spray field in the near nozzle region the studies are conducted to determined the effects of multi jets on the flow pattern in the near nozzle region The phase doppler particle analyzer (PDPA) has been used to measure the velocities and sizes of the droplets. The flow field of two-phase liquid drop-air jets is formed from three injectors arranged in t line. Furthermore the two-phase flow field has been analyzed numerically also. The numerical analysis consists of two computational models, namely (i) 3 non-evaporating two-phase jets, (II) 3 evaporating two phase jets. The Eulerian-Eulerian approach in incorporated in both the numerical models. Since the flow is turbulent, a two-equation model (k-Epsilon) is implemented in the numerical analysis. Numerical solution of the conservation equation is obtained using PHOENICS computer code. Boundary conditions are provided from the experimental measurements. Numerical domain for the two models of the analysis starts at some distance (about 10 diameters of the injector orifice) where the atomization process is complete and droplet size and velocity could be measured experimentally. (author)

Comparison of Two-Phase Pipe Flow in OpenFOAM with a Mechanistic Model

Science.gov (United States)

Shuard, Adrian M.; Mahmud, Hisham B.; King, Andrew J.

2016-03-01

Two-phase pipe flow is a common occurrence in many industrial applications such as power generation and oil and gas transportation. Accurate prediction of liquid holdup and pressure drop is of vast importance to ensure effective design and operation of fluid transport systems. In this paper, a Computational Fluid Dynamics (CFD) study of a two-phase flow of air and water is performed using OpenFOAM. The two-phase solver, interFoam is used to identify flow patterns and generate values of liquid holdup and pressure drop, which are compared to results obtained from a two-phase mechanistic model developed by Petalas and Aziz (2002). A total of 60 simulations have been performed at three separate pipe inclinations of 0°, +10° and -10° respectively. A three dimensional, 0.052m diameter pipe of 4m length is used with the Shear Stress Transport (SST) k - ɷ turbulence model to solve the turbulent mixtures of air and water. Results show that the flow pattern behaviour and numerical values of liquid holdup and pressure drop compare reasonably well to the mechanistic model.

Transient two-dimensional flow in porous media

International Nuclear Information System (INIS)

Sharpe, L. Jr.

1979-01-01

The transient flow of an isothermal ideal gas from the cavity formed by an underground nuclear explosion is investigated. A two-dimensional finite element method is used in analyzing the gas flow. Numerical results of the pressure distribution are obtained for both the stemming column and the surrounding porous media

SedFoam-2.0: a 3-D two-phase flow numerical model for sediment transport

Directory of Open Access Journals (Sweden)

J. Chauchat

2017-11-01

Full Text Available In this paper, a three-dimensional two-phase flow solver, SedFoam-2.0, is presented for sediment transport applications. The solver is extended from twoPhaseEulerFoam available in the 2.1.0 release of the open-source CFD (computational fluid dynamics toolbox OpenFOAM. In this approach the sediment phase is modeled as a continuum, and constitutive laws have to be prescribed for the sediment stresses. In the proposed solver, two different intergranular stress models are implemented: the kinetic theory of granular flows and the dense granular flow rheology μ(I. For the fluid stress, laminar or turbulent flow regimes can be simulated and three different turbulence models are available for sediment transport: a simple mixing length model (one-dimensional configuration only, a k − ε, and a k − ω model. The numerical implementation is demonstrated on four test cases: sedimentation of suspended particles, laminar bed load, sheet flow, and scour at an apron. These test cases illustrate the capabilities of SedFoam-2.0 to deal with complex turbulent sediment transport problems with different combinations of intergranular stress and turbulence models.

SedFoam-2.0: a 3-D two-phase flow numerical model for sediment transport

Science.gov (United States)

Chauchat, Julien; Cheng, Zhen; Nagel, Tim; Bonamy, Cyrille; Hsu, Tian-Jian

2017-11-01

In this paper, a three-dimensional two-phase flow solver, SedFoam-2.0, is presented for sediment transport applications. The solver is extended from twoPhaseEulerFoam available in the 2.1.0 release of the open-source CFD (computational fluid dynamics) toolbox OpenFOAM. In this approach the sediment phase is modeled as a continuum, and constitutive laws have to be prescribed for the sediment stresses. In the proposed solver, two different intergranular stress models are implemented: the kinetic theory of granular flows and the dense granular flow rheology μ(I). For the fluid stress, laminar or turbulent flow regimes can be simulated and three different turbulence models are available for sediment transport: a simple mixing length model (one-dimensional configuration only), a k - ɛ, and a k - ω model. The numerical implementation is demonstrated on four test cases: sedimentation of suspended particles, laminar bed load, sheet flow, and scour at an apron. These test cases illustrate the capabilities of SedFoam-2.0 to deal with complex turbulent sediment transport problems with different combinations of intergranular stress and turbulence models.

Lattice Boltzmann model for simulating immiscible two-phase flows

International Nuclear Information System (INIS)

Reis, T; Phillips, T N

2007-01-01

The lattice Boltzmann equation is often promoted as a numerical simulation tool that is particularly suitable for predicting the flow of complex fluids. This paper develops a two-dimensional 9-velocity (D2Q9) lattice Boltzmann model for immiscible binary fluids with variable viscosities and density ratio using a single relaxation time for each fluid. In the macroscopic limit, this model is shown to recover the Navier-Stokes equations for two-phase flows. This is achieved by constructing a two-phase component of the collision operator that induces the appropriate surface tension term in the macroscopic equations. A theoretical expression for surface tension is determined. The validity of this analysis is confirmed by comparing numerical and theoretical predictions of surface tension as a function of density. The model is also shown to predict Laplace's law for surface tension and Poiseuille flow of layered immiscible binary fluids. The spinodal decomposition of two fluids of equal density but different viscosity is then studied. At equilibrium, the system comprises one large low viscosity bubble enclosed by the more viscous fluid in agreement with theoretical arguments of Renardy and Joseph (1993 Fundamentals of Two-Fluid Dynamics (New York: Springer)). Two other simulations, namely the non-equilibrium rod rest and the coalescence of two bubbles, are performed to show that this model can be used to simulate two fluids with a large density ratio

A new two dimensional spectral/spatial multi-diagonal code for noncoherent optical code division multiple access (OCDMA) systems

Science.gov (United States)

Kadhim, Rasim Azeez; Fadhil, Hilal Adnan; Aljunid, S. A.; Razalli, Mohamad Shahrazel

2014-10-01

A new two dimensional codes family, namely two dimensional multi-diagonal (2D-MD) codes, is proposed for spectral/spatial non-coherent OCDMA systems based on the one dimensional MD code. Since the MD code has the property of zero cross correlation, the proposed 2D-MD code also has this property. So that, the multi-access interference (MAI) is fully eliminated and the phase induced intensity noise (PIIN) is suppressed with the proposed code. Code performance is analyzed in terms of bit error rate (BER) while considering the effect of shot noise, PIIN, and thermal noise. The performance of the proposed code is compared with the related MD, modified quadratic congruence (MQC), two dimensional perfect difference (2D-PD) and two dimensional diluted perfect difference (2D-DPD) codes. The analytical and the simulation results reveal that the proposed 2D-MD code outperforms the other codes. Moreover, a large number of simultaneous users can be accommodated at low BER and high data rate.

Non-isothermal effects on multi-phase flow in porous medium

DEFF Research Database (Denmark)

Singh, Ashok; Wang, W; Park, C. H.

2010-01-01

In this paper a ppT -formulation for non-isothermal multi-phase flow is given including diffusion and latent heat effects. Temperature and pressure dependencies of governing parameters are considered, in particular surface tension variation on phase interfaces along with temperature changes. A we...

Zero-G two phase flow regime modeling in adiabatic flow

International Nuclear Information System (INIS)

Reinarts, T.R.; Best, F.R.; Wheeler, M.; Miller, K.M.

1993-01-01

Two-phase flow, thermal management systems are currently being considered as an alternative to conventional, single phase systems for future space missions because of their potential to reduce overall system mass, size, and pumping power requirements. Knowledge of flow regime transitions, heat transfer characteristics, and pressure drop correlations is necessary to design and develop two-phase systems. This work is concerned with microgravity, two-phase flow regime analysis. The data come from a recent sets of experiments. The experiments were funded by NASA Johnson Space Center (JSC) and conducted by NASA JSC with Texas A ampersand M University. The experiment was on loan to NASA JSC from Foster-Miller, Inc., who constructed it with funding from the Air Force Phillips Laboratory. The experiment used R12 as the working fluid. A Foster-Miller two phase pump was used to circulate the two phase mixture and allow separate measurements of the vapor and liquid flow streams. The experimental package was flown 19 times for 577 parabolas aboard the NASA KC-135 aircraft which simulates zero-G conditions by its parabolic flight trajectory. Test conditions included bubbly, slug and annular flow regimes in 0-G. The superficial velocities of liquid and vapor have been obtained from the measured flow rates and are presented along with the observed flow regimes and several flow regime transition predictions. None of the predictions completely describe the transitions as indicated by the data

Hamiltonian formalism of two-dimensional Vlasov kinetic equation.

Science.gov (United States)

Pavlov, Maxim V

2014-12-08

In this paper, the two-dimensional Benney system describing long wave propagation of a finite depth fluid motion and the multi-dimensional Russo-Smereka kinetic equation describing a bubbly flow are considered. The Hamiltonian approach established by J. Gibbons for the one-dimensional Vlasov kinetic equation is extended to a multi-dimensional case. A local Hamiltonian structure associated with the hydrodynamic lattice of moments derived by D. J. Benney is constructed. A relationship between this hydrodynamic lattice of moments and the two-dimensional Vlasov kinetic equation is found. In the two-dimensional case, a Hamiltonian hydrodynamic lattice for the Russo-Smereka kinetic model is constructed. Simple hydrodynamic reductions are presented.

Flow patterns in vertical two-phase flow

International Nuclear Information System (INIS)

McQuillan, K.W.; Whalley, P.B.

1985-01-01

This paper is concerned with the flow patterns which occur in upwards gas-liquid two-phase flow in vertical tubes. The basic flow patterns are described and the use of flow patter maps is discussed. The transition between plug flow and churn flow is modelled under the assumption that flooding of the falling liquid film limits the stability of plug flow. The resulting equation is combined with other flow pattern transition equations to produce theoretical flow pattern maps, which are then tested against experimental flow pattern data. Encouraging agreement is obtained

Two-dimensional macroscopic quantum tunneling in multi-gap superconductor Josephson junctions

International Nuclear Information System (INIS)

Asai, Hidehiro; Kawabata, Shiro; Ota, Yukihiro; Machida, Masahiko

2014-01-01

Low-temperature characters of superconducting devices yield definite probes for different superconducting phenomena. We study the macroscopic quantum tunneling (MQT) in a Josephson junction, composed of a single-gap superconductor and a two-gap superconductor. Since this junction has two kinds to the superconducting phase differences, calculating the MQT escape rate requires the analysis of quantum tunneling in a multi-dimensional configuration space. Our approach is the semi-classical approximation along a 1D curve in a 2D potential- energy landscape, connecting two adjacent potential (local) minimums through a saddle point. We find that this system has two plausible tunneling paths; an in-phase path and an out-of-phase path. The former is characterized by the Josephson-plasma frequency, whereas the latter is by the frequency of the characteristic collective mode in a two-band superconductor, Josephson- Leggett mode. Depending on external bias current and inter-band Josephson-coupling energy, one of them mainly contributes to the MQT. Our numerical calculations show that the difference between the in-phase path and the out-of-phase path is manifest, with respect to the bias- current-dependence of the MQT escape rate. This result suggests that our MQT setting be an indicator of the Josephson-Leggett mode

Complex network analysis of phase dynamics underlying oil-water two-phase flows

Science.gov (United States)

Gao, Zhong-Ke; Zhang, Shan-Shan; Cai, Qing; Yang, Yu-Xuan; Jin, Ning-De

2016-01-01

Characterizing the complicated flow behaviors arising from high water cut and low velocity oil-water flows is an important problem of significant challenge. We design a high-speed cycle motivation conductance sensor and carry out experiments for measuring the local flow information from different oil-in-water flow patterns. We first use multivariate time-frequency analysis to probe the typical features of three flow patterns from the perspective of energy and frequency. Then we infer complex networks from multi-channel measurements in terms of phase lag index, aiming to uncovering the phase dynamics governing the transition and evolution of different oil-in-water flow patterns. In particular, we employ spectral radius and weighted clustering coefficient entropy to characterize the derived unweighted and weighted networks and the results indicate that our approach yields quantitative insights into the phase dynamics underlying the high water cut and low velocity oil-water flows. PMID:27306101

Encyclopedia of two-phase heat transfer and flow I fundamentals and methods

CERN Document Server

2015-01-01

The aim of the two–set series is to present a very detailed and up–to–date reference for researchers and practicing engineers in the fields of mechanical, refrigeration, chemical, nuclear and electronics engineering on the important topic of two-phase heat transfer and two-phase flow. The scope of the first set of 4 volumes presents the fundamentals of the two-phase flows and heat transfer mechanisms, and describes in detail the most important prediction methods, while the scope of the second set of 4 volumes presents numerous special topics and numerous applications, also including numerical simulation methods. Practicing engineers will find extensive coverage to applications involving: multi-microchannel evaporator cold plates for electronics cooling, boiling on enhanced tubes and tube bundles, flow pattern based methods for predicting boiling and condensation inside horizontal tubes, pressure drop methods for singularies (U-bends and contractions), boiling in multiport tubes, and boiling and condens...

Two-phase-flow models and their limitations

International Nuclear Information System (INIS)

Ishii, M.; Kocamustafaogullari, G.

1982-01-01

An accurate prediction of transient two-phase flow is essential to safety analyses of nuclear reactors under accident conditions. The fluid flow and heat transfer encountered are often extremely complex due to the reactor geometry and occurrence of transient two-phase flow. Recently considerable progresses in understanding and predicting these phenomena have been made by a combination of rigorous model development, advanced computational techniques, and a number of small and large scale supporting experiments. In view of their essential importance, the foundation of various two-phase-flow models and their limitations are discussed in this paper

Geometric effects of 90-degree vertical elbows on local two-phase flow parameters

International Nuclear Information System (INIS)

Yadav, M.; Worosz, T.; Kim, S.

2011-01-01

This study presents the geometric effects of 90-degree vertical elbows on the development of the local two-phase flow parameters. A multi-sensor conductivity probe is used to measure local two-phase flow parameters. It is found that immediately downstream of the vertical-upward elbow, the bubbles have a bimodal distribution along the horizontal radius of the pipe cross-section causing a dual-peak in the profiles of local void fraction and local interfacial area concentration. Immediately downstream of the vertical-downward elbow it is observed that the bubbles tend to migrate towards the inside of the elbow's curvature. The axial transport of void fraction and interfacial area concentration indicates that the elbows promote bubble disintegration. Preliminary predictions are obtained from group-one interfacial area transport equation (IATE) model for vertical-upward and vertical-downward two-phase flow. (author)

Influence of the Fin on Two-Dimensional Characteristics of Dispersed Flow With Wall Liquid Film in the Vicinity of Obstacle

International Nuclear Information System (INIS)

Stosic, Zoran V.; Stevanovic, Vladimir D.; Serizawa, Akimi

2002-01-01

Spacers have positive effects on the heat transfer enhancement and critical heat flux (CHF) increase downstream of their location in the boiling channel. These effects are further increased by the inclusion of the fin on the spacer rear edge. Numerical simulation of a separation in a high void gas phase and dispersed droplets flow around a spacer, with a liquid film flowing on the wall, is performed. Mechanisms leading to the CHF increase due to the two-phase flow separation and liquid film thickening downstream the spacer are demonstrated. Numerical simulations of gas phase, entrained droplets and wall liquid film flows were performed with the three-fluid model and with the application of the high order numerical scheme for the liquid film surface interface tracking. Predicted is a separation of gas and entrained droplets streams around the spacer without and with a fin inclined 30 and 60 degrees to the wall, as well as a change of wall liquid film thickness in the vicinity of spacer. Results of liquid film dynamic behaviour are compared with the recently obtained experimental results. Multi-dimensional characteristics of surface waves on the liquid film were measured with newly developed ultrasonic transmission technique in a 3 3 rod bundle test section with air-water flow under atmospheric conditions. Obtained numerical results are in good agreement with experimental observations. The presented investigation gives insight into the complex mechanisms of separated two-phase flow with wall liquid film around the spacer and support thermal-hydraulic design and optimisation of flow obstacles in various thermal equipment. (authors)

Universality Results for Multi-phase Hele-Shaw Flows

Science.gov (United States)

Daripa, Prabir

2013-03-01

Saffman-Taylor instability is a well known viscosity driven instability of an interface separating two immiscible fluids. We study linear stability of displacement processes in a Hele-Shaw cell involving an arbitrary number of immiscible fluid phases. This is a problem involving many interfaces. Universal stability results have been obtained for this multi-phase immiscible flow in the sense that the results hold for arbitrary number of interfaces. These stability results have been applied to design displacement processes that are considerably less unstable than the pure Saffman-Taylor case. In particular, we derive universal formula which gives specific values of the viscosities of the fluid layers corresponding to smallest unstable band. Other similar universal results will also be presented. The talk is based on the following paper. This work was supported by the Qatar National Research Fund (a member of The Qatar Foundation).

Characterization of two-phase flow regimes in horizontal tubes using {sup 81m}Kr tracer experiments

Energy Technology Data Exchange (ETDEWEB)

Oriol, Jean [LPAC, CEA Grenoble, 17, rue des Martyrs, 38054 Grenoble Cedex 9 (France); Leclerc, Jean Pierre [Laboratoire des Sciences du Genie Chimique (LSGC), Nancy-Universite, CNRS, BP 20451, F-54001 Nancy (France)], E-mail: leclerc@ensic.inpl-nancy.fr; Berne, Philippe; Gousseau, Georges [L2T, CEA Grenoble, 17, rue des Martyrs, 38054 Grenoble Cedex 9 (France); Jallut, Christian [Universite de Lyon, Universite Lyon 1, LAGEP, UMR CNRS 5007, ESCPE, 43 Bd du 11 novembre 1918, 69622 Villeurbanne Cedex (France); Tochon, Patrice; Clement, Patrice [GRETh, CEA Grenoble, 17, rue des Martyrs, 38054 Grenoble Cedex 9 (France)

2008-10-15

The diagnosis of heat exchangers on duty with respect to flow mal-distributions needs the development of non-intrusive inlet-outlet experimental techniques in order to perform an online fault diagnosis. Tracer experiments are an example of such techniques. They can be applied to mono-phase heat exchangers but also to multi-phase ones. In this case, the tracer experiments are more difficult to perform. In order to check for the capabilities of tracer experiments to be used for the flow mal-distribution diagnosis in the case of multi-phase heat exchangers, we present here a preliminary study on the simplest possible system: two-phase flows in a horizontal tube. {sup 81m}Kr is used as gas tracer and properly collimated NaI (TI) crystal scintillators as detectors. The specific shape of the tracer response allows two-phase flow regimes to be characterized. Signal analysis allows the estimation of the gas phase real average velocity and consequently of the liquid phase real average velocity as well as of the volumetric void fraction. These results are compared successfully to those obtained with liquid phase tracer experiments previously presented by Oriol et al. 2007. Characterization of the two-phase flow regimes and liquid dispersion in horizontal and vertical tubes using coloured tracer and no intrusive optical detector. Chem. Eng. Sci. 63(1), 24-34, as well as to those given by correlations from literature.

Lennard-Jones fluids in two-dimensional nano-pores. Multi-phase coexistence and fluid structure

Science.gov (United States)

Yatsyshin, Petr; Savva, Nikos; Kalliadasis, Serafim

2014-03-01

We present a number of fundamental findings on the wetting behaviour of nano-pores. A popular model for fluid confinement is a one-dimensional (1D) slit pore formed by two parallel planar walls and it exhibits capillary condensation (CC): a first-order phase transition from vapour to capillary-liquid (Kelvin shift). Capping such a pore at one end by a third orthogonal wall forms a prototypical two-dimensional (2D) pore. We show that 2D pores possess a wetting temperature such that below this temperature CC remains of first order, above it becomes a continuous phase transition manifested by a slab of capillary-liquid filling the pore from the capping wall. Continuous CC exhibits hysteresis and can be preceded by a first-order capillary prewetting transition. Additionally, liquid drops can form in the corners of the 2D pore (remnant of 2D wedge prewetting). The three fluid phases, vapour, capillary-liquid slab and corner drops, can coexist at the pore triple point. Our model is based on the statistical mechanics of fluids in the density functional formulation. The fluid-fluid and fluid-substrate interactions are dispersive. We analyze in detail the microscopic fluid structure, isotherms and full phase diagrams. Our findings also suggest novel ways to control wetting of nano-pores. We are grateful to the European Research Council via Advanced Grant No. 247031 for support.

Fluid-elastic vibration in two-phase cross flow

International Nuclear Information System (INIS)

Sasakawa, T.; Serizawa, A.; Kawara, Z.

2003-01-01

The present work aims at clarifying the mechanisms of fluid elastic vibration of tube bundles in two-phase cross flow. The experiment is conducted using air-water two-phase flow under atmospheric pressure. The test section is a 1.03m long transparent acrylic square duct with 128 x 128 mm 2 cross section, which consists of 3 rod-rows with 5 rods in each row. The rods are 125mm long aluminum rods with 22 mm in diameter (p/D=1.45). The natural frequency of rod vibration is about 30Hz. The result indicated a diversion of observed trend in vibration behavior depending on two-phase flow patterns either bubbly flow or churn flow. Specifically, in churn flow, the fluid elastic vibration has been observed to occur when the frequency in void fraction fluctuation approached to the natural frequency of the rods, but this was not the case in fluid elastic vibration in bubbly flow. This fact suggests the existence of mechanisms closely coupled with two-phase flow structures depending on the flow patterns, that is, static two-phase character-controlled mechanism in bubbly flow and dynamic character- controlled in churn flow

An interfacial shear term evaluation study for adiabatic dispersed air–water two-phase flow with the two-fluid model using CFD

Energy Technology Data Exchange (ETDEWEB)

Sharma, S.L., E-mail: sharma55@purdue.edu [School of Nuclear Engineering, Purdue University, West Lafayette, IN (United States); Hibiki, T.; Ishii, M. [School of Nuclear Engineering, Purdue University, West Lafayette, IN (United States); Schlegel, J.P. [Department of Mining and Nuclear Engineering, Missouri University of Science and Technology, Rolla, MO (United States); Buchanan, J.R.; Hogan, K.J. [Bettis Laboratory, Naval Nuclear Laboratory, West Mifflin, PA (United States); Guilbert, P.W. [ANSYS UK Ltd, Oxfordshire (United Kingdom)

2017-02-15

Highlights: • Closure form of the interfacial shear term in three-dimensional form is investigated. • Assessment against adiabatic upward bubbly air–water flow data using CFD. • Effect of addition of the interfacial shear term on the phase distribution. - Abstract: In commercially available Computational Fluid Dynamics (CFD) codes such as ANSYS CFX and Fluent, the interfacial shear term is missing in the field momentum equations. The derivation of the two-fluid model (Ishii and Hibiki, 2011) indicates the presence of this term as a momentum source in the right hand side of the field momentum equation. The inclusion of this term is considered important for proper modeling of the interfacial momentum coupling between phases. For separated flows, such as annular flow, the importance of the shear term is understood in the one-dimensional (1-D) form as the major mechanism by which the wall shear is transferred to the gas phase (Ishii and Mishima, 1984). For gas dispersed two-phase flow CFD simulations, it is important to assess the significance of this term in the prediction of phase distributions. In the first part of this work, the closure of this term in three-dimensional (3-D) form in a CFD code is investigated. For dispersed gas–liquid flow, such as bubbly or churn-turbulent flow, bubbles are dispersed in the shear layer of the continuous phase. The continuous phase shear stress is mainly due to the presence of the wall and the modeling of turbulence through the Boussinesq hypothesis. In a 3-D simulation, the continuous phase shear stress can be calculated from the continuous fluid velocity gradient, so that the interfacial shear term can be closed using the local values of the volume fraction and the total stress of liquid phase. This form also assures that the term acts as an action-reaction force for multiple phases. In the second part of this work, the effect of this term on the volume fraction distribution is investigated. For testing the model two-phase

Equivalence of two models in single-phase multicomponent flow simulations

KAUST Repository

Wu, Yuanqing

2016-02-28

In this work, two models to simulate the single-phase multicomponent flow in reservoirs are introduced: single-phase multicomponent flow model and two-phase compositional flow model. Because the single-phase multicomponent flow is a special case of the two-phase compositional flow, the two-phase compositional flow model can also simulate the case. We compare and analyze the two models when simulating the single-phase multicomponent flow, and then demonstrate the equivalence of the two models mathematically. An experiment is also carried out to verify the equivalence of the two models.

Equivalence of two models in single-phase multicomponent flow simulations

KAUST Repository

Wu, Yuanqing; Sun, Shuyu

2016-01-01

In this work, two models to simulate the single-phase multicomponent flow in reservoirs are introduced: single-phase multicomponent flow model and two-phase compositional flow model. Because the single-phase multicomponent flow is a special case of the two-phase compositional flow, the two-phase compositional flow model can also simulate the case. We compare and analyze the two models when simulating the single-phase multicomponent flow, and then demonstrate the equivalence of the two models mathematically. An experiment is also carried out to verify the equivalence of the two models.

Statistical mechanics and correlation properties of a rotating two-dimensional flow of like-sign vortices

International Nuclear Information System (INIS)

Viecelli, J.A.

1993-01-01

The Hamiltonian flow of a set of point vortices of like sign and strength has a low-temperature phase consisting of a rotating triangular lattice of vortices, and a normal temperature turbulent phase consisting of random clusters of vorticity that orbit about a common center along random tracks. The mean-field flow in the normal temperature phase has similarities with turbulent quasi-two-dimensional rotating laboratory and geophysical flows, whereas the low-temperature phase displays effects associated with quantum fluids. In the normal temperature phase the vortices follow power-law clustering distributions, while in the time domain random interval modulation of the vortex orbit radii fluctuations produces singular fractional exponent power-law low-frequency spectra corresponding to time autocorrelation functions with fractional exponent power-law tails. Enhanced diffusion is present in the turbulent state, whereas in the solid-body rotation state vortices thermally diffuse across the lattice. Over the entire temperature range the interaction energy of a single vortex in the field of the rest of the vortices follows positive temperature Fermi--Dirac statistics, with the zero temperature limit corresponding to the rotating crystal phase, and the infinite temperature limit corresponding to a Maxwellian distribution. Analyses of weather records dependent on the large-scale quasi-two-dimensional atmospheric circulation suggest the presence of singular fractional exponent power-law spectra and fractional exponent power-law autocorrelation tails, consistent with the theory

Spectral Cascade-Transport Turbulence Model Development for Two-Phase Flows

Science.gov (United States)

Brown, Cameron Scott

Turbulence modeling remains a challenging problem in nuclear reactor applications, particularly for the turbulent multiphase flow conditions in nuclear reactor subchannels. Understanding the fundamental physics of turbulent multiphase flows is crucial for the improvement and further development of multiphase flow models used in reactor operation and safety calculations. Reactor calculations with Reynolds-averaged Navier-Stokes (RANS) approach continue to become viable tools for reactor analysis. The on-going increase in available computational resources allows for turbulence models that are more complex than the traditional two-equation models to become practical choices for nuclear reactor computational fluid dynamic (CFD) and multiphase computational fluid dynamic (M-CFD) simulations. Similarly, increased computational capabilities continue to allow for higher Reynolds numbers and more complex geometries to be evaluated using direct numerical simulation (DNS), thus providing more validation and verification data for turbulence model development. Spectral turbulence models are a promising approach to M-CFD simulations. These models resolve mean flow parameters as well as the turbulent kinetic energy spectrum, reproducing more physical details of the turbulence than traditional two-equation type models. Previously, work performed by other researchers on a spectral cascade-transport model has shown that the model behaves well for single and bubbly twophase decay of isotropic turbulence, single and two-phase uniform shear flow, and single-phase flow in a channel without resolving the near-wall boundary layer for relatively low Reynolds number. Spectral models are great candidates for multiphase RANS modeling since bubble source terms can be modeled as contributions to specific turbulence scales. This work focuses on the improvement and further development of the spectral cascadetransport model (SCTM) to become a three-dimensional (3D) turbulence model for use in M

Nonlinear dynamics of two-phase flow

International Nuclear Information System (INIS)

Rizwan-uddin

1986-01-01

Unstable flow conditions can occur in a wide variety of laboratory and industry equipment that involve two-phase flow. Instabilities in industrial equipment, which include boiling water reactor (BWR) cores, steam generators, heated channels, cryogenic fluid heaters, heat exchangers, etc., are related to their nonlinear dynamics. These instabilities can be of static (Ledinegg instability) or dynamic (density wave oscillations) type. Determination of regions in parameters space where these instabilities can occur and knowledge of system dynamics in or near these regions is essential for the safe operation of such equipment. Many two-phase flow engineering components can be modeled as heated channels. The set of partial differential equations that describes the dynamics of single- and two-phase flow, for the special case of uniform heat flux along the length of the channel, can be reduced to a set of two coupled ordinary differential equations [in inlet velocity v/sub i/(t) and two-phase residence time tau(t)] involving history integrals: a nonlinear ordinary functional differential equation and an integral equation. Hence, to solve these equations, the dependent variables must be specified for -(nu + tau) ≤ t ≤ 0, where nu is the single-phase residence time. This system of nonlinear equations has been solved analytically using asymptotic expansion series for finite but small perturbations and numerically using finite difference techniques

Smart multi-channel two-dimensional micro-gas chromatography for rapid workplace hazardous volatile organic compounds measurement.

Science.gov (United States)

Liu, Jing; Seo, Jung Hwan; Li, Yubo; Chen, Di; Kurabayashi, Katsuo; Fan, Xudong

2013-03-07

We developed a novel smart multi-channel two-dimensional (2-D) micro-gas chromatography (μGC) architecture that shows promise to significantly improve 2-D μGC performance. In the smart μGC design, a non-destructive on-column gas detector and a flow routing system are installed between the first dimensional separation column and multiple second dimensional separation columns. The effluent from the first dimensional column is monitored in real-time and decision is then made to route the effluent to one of the second dimensional columns for further separation. As compared to the conventional 2-D μGC, the greatest benefit of the smart multi-channel 2-D μGC architecture is the enhanced separation capability of the second dimensional column and hence the overall 2-D GC performance. All the second dimensional columns are independent of each other, and their coating, length, flow rate and temperature can be customized for best separation results. In particular, there is no more constraint on the upper limit of the second dimensional column length and separation time in our architecture. Such flexibility is critical when long second dimensional separation is needed for optimal gas analysis. In addition, the smart μGC is advantageous in terms of elimination of the power intensive thermal modulator, higher peak amplitude enhancement, simplified 2-D chromatogram re-construction and potential scalability to higher dimensional separation. In this paper, we first constructed a complete smart 1 × 2 channel 2-D μGC system, along with an algorithm for automated control/operation of the system. We then characterized and optimized this μGC system, and finally employed it in two important applications that highlight its uniqueness and advantages, i.e., analysis of 31 workplace hazardous volatile organic compounds, and rapid detection and identification of target gas analytes from interference background.

Two phase pressure drop across abrupt area changes in oscillatory flow

International Nuclear Information System (INIS)

Weisman, J.; Ake, T.; Knott, R.

1975-05-01

The ability of one-dimensional momentum balances to predict behavior during two-phase oscillatory flow has been examined. Flow oscillations in a Freon-Freon vapor system were induced by cycling a three-way valve so as to divert a varying portion of the total flow from the test section. The data taken included the phase shift between pressure and flow curves and the head fluctuation to flow fluctuation ratio. These data were compared to theoretical predictions. For void fractions below about 0.55, reasonable agreement between predictions and observations was obtained. At the highest void fractions examined, agreement between prediction and measurement was poor. However, at these high void fractions measurement errors may have been very substantial. Further, the simplifying assumptions required by the numerical calculation procedure chosen could no longer be considered valid. No real conclusion could therefore be drawn from the data at high voids

One-dimensional two-phase thermal hydraulics (ENSTA course); Thermo-hydraulique diphasique monodimensionnelle. Cours ENSTA

Energy Technology Data Exchange (ETDEWEB)

Olive, J

1995-11-01

This course is part of the ENSTA 3rd year thermal hydraulics program (nuclear power option). Its purpose is to provide the theoretical basis and main physical notions pertaining to two-phase flow, mainly focussed on water-steam flows. The introduction describes the physical specificities of these flows, emphasizing their complexity. The mathematical bases are then presented (partial derivative equations), leading to a one-dimensional type, simplified description. Balances drawn up for a pipe length volume are used to introduce the mass conservation. motion and energy equations for each phase. Various postulates used to simplify two-phase models are presented, culminating in homogeneous model definitions and equations, several common examples of which are given. The model is then applied to the calculation of pressure drops in two-phase flows. This involves presenting the models most frequently used to represent pressure drops by friction or due to pipe irregularities, without giving details (numerical values of parameters). This chapter terminates with a brief description of static and dynamic instabilities in two-phase flows. Finally, heat transfer conditions frequently encountered in liquid-steam flows are described, still in the context of a 1D model. This chapter notably includes reference to under-saturated boiling conditions and the various forms of DNB. The empirical heat transfer laws are not discussed in detail. Additional material is appended, some of which is in the form of corrected exercises. (author). 6 appends.

Identification of dominant structures and their flow dynamics in the turbulent two-phase flow using POD technique

Energy Technology Data Exchange (ETDEWEB)

Munir, Shahzad; Siddiqui, Muhammad Israr; Heikal, Morgan; Aziz, Abdul Rashid Abdul [Universiti Teknologi PETRONAS, Bander Seri Iskandar (Malaysia); Sercey, Guillaume de [University of Brighton, Brighton (United Kingdom)

2015-11-15

The Proper orthogonal decomposition (POD) method has seen increasingly used in the last two decades and has a lot of applications for the comparison of experimental and numerically simulated data. The POD technique is often used to extract information about coherent structures dominating the flow. The two-dimensional and two-component instantaneous velocity fields of both liquid and gas phases of a slug flow were obtained by Particle image velocimetry (PIV) combined with Laser induced fluorescence (LIF). POD was applied to the velocity fields of both phases separately to identify the coherent flow structures. We focused on POD eigenmodes and their corresponding energy contents of both liquid and gas phases. The sum of first few eigenmodes that contain maximum turbulent kinetic energy of the flow represents the coherent structures. In the case of liquid phase the first eigenmode contained 42% of the total energy, while in the gas phase the decaying energy distribution was flat. The POD results showed that the coefficient of mode 1 for the liquid phase oscillated between positive and negative values and had the highest amplitude. For the visualization of coherent motion different linear combinations of eigenmodes for liquid and gas phases were used. The phenomena of turbulent bursting events associated with Q2 events (low momentum fluid moving away from the wall) and Q4 events (high momentum flow moving towards the wall) were also discussed to assess its contribution in turbulence production.

Comparison of Two-Phase Pipe Flow in OpenFOAM with a Mechanistic Model

International Nuclear Information System (INIS)

Shuard, Adrian M; Mahmud, Hisham B; King, Andrew J

2016-01-01

Two-phase pipe flow is a common occurrence in many industrial applications such as power generation and oil and gas transportation. Accurate prediction of liquid holdup and pressure drop is of vast importance to ensure effective design and operation of fluid transport systems. In this paper, a Computational Fluid Dynamics (CFD) study of a two-phase flow of air and water is performed using OpenFOAM. The two-phase solver, interFoam is used to identify flow patterns and generate values of liquid holdup and pressure drop, which are compared to results obtained from a two-phase mechanistic model developed by Petalas and Aziz (2002). A total of 60 simulations have been performed at three separate pipe inclinations of 0°, +10° and -10° respectively. A three dimensional, 0.052m diameter pipe of 4m length is used with the Shear Stress Transport (SST) k - ω turbulence model to solve the turbulent mixtures of air and water. Results show that the flow pattern behaviour and numerical values of liquid holdup and pressure drop compare reasonably well to the mechanistic model. (paper)

Analysis of phase dynamics in two-phase flow using latticegas automata

International Nuclear Information System (INIS)

Ohashi, H.; Hashimoto, Y.; Tsumaya, A.; Chen, Y.; Akiyama, M.

1998-01-01

In this paper, we describe lattice gas automaton models appropriate for two-phase flow simulation and their applications to study various phase dynamics of two-fluid mixtures. Several algorithms are added to the original immiscible Lattice Gas model to adjust surface tension and to introduce density difference between two fluids. Surface tension is controlled by the collision rules an difference in density is due to nonlocal forces between automaton particles. We simulate the relative motion of the dispersed phase in another continuous fluid. Deformation and disintegration of rising drops are reproduced. The interaction between multiple drops is also observed in calculations. Furutre, we obtain the transition of the two-phase flow pattern from bubbly, slug to annular flow. Density difference of two phase is one of the key ingredients to generate the annular flow pattern

Multi-Band Light Curves from Two-Dimensional Simulations of Gamma-Ray Burst Afterglows

Science.gov (United States)

MacFadyen, Andrew

2010-01-01

The dynamics of gamma-ray burst outflows is inherently multi-dimensional. 1.) We present high resolution two-dimensional relativistic hydrodynamics simulations of GRBs in the afterglow phase using adaptive mesh refinement (AMR). Using standard synchrotron radiation models, we compute multi-band light curves, from the radio to X-ray, directly from the 2D hydrodynamics simulation data. We will present on-axis light curves for both constant density and wind media. We will also present off-axis light curves relevant for searches for orphan afterglows. We find that jet breaks are smoothed due to both off-axis viewing and wind media effects. 2.) Non-thermal radiation mechanisms in GRB afterglows require substantial magnetic field strengths. In turbulence driven by shear instabilities in relativistic magnetized gas, we demonstrate that magnetic field is naturally amplified to half a percent of the total energy (epsilon B = 0.005). We will show high resolution three dimensional relativistic MHD simulations of this process as well as particle in cell (PIC) simulations of mildly relativistic collisionless shocks.

Development of subchannel void measurement sensor and multidimensional two-phase flow dynamics in rod bundle

International Nuclear Information System (INIS)

Arai, T.; Furuya, M.; Kanai, T.; Shirakawa, K.

2011-01-01

An accurate subchannel database is crucial for modeling the multidimensional two-phase flow in a rod bundle and for validating subchannel analysis codes. Based on available reference, it can be said that a point-measurement sensor for acquiring void fractions and bubble velocity distributions do not infer interactions of the subchannel flow dynamics, such as a cross flow and flow distribution, etc. In order to acquire multidimensional two-phase flow in a 10×10 rod bundle with an o.d. of 10 mm and 3110 mm length, a new sensor consisting of 11-wire by 11-wire and 10-rod by 10-rod electrodes was developed. Electric potential in the proximity region between two wires creates a void fraction in the center subchannel region, like a so-called wire mesh sensor. A unique aspect of the devised sensor is that the void fraction near the rod surface can be estimated from the electric potential in the proximity region between one wire and one rod. The additional 400 points of void fraction and phasic velocity in 10×10 bundle can therefore be acquired. The devised sensor exhibits the quasi three-dimensional flow structures, i.e. void fraction, phasic velocity and bubble chord length distributions. These quasi three-dimensional structures exhibit the complexity of two-phase flow dynamics, such as coalescence and the breakup of bubbles in transient phasic velocity distributions. (author)

A positivity preserving and conservative variational scheme for phase-field modeling of two-phase flows

Science.gov (United States)

Joshi, Vaibhav; Jaiman, Rajeev K.

2018-05-01

We present a positivity preserving variational scheme for the phase-field modeling of incompressible two-phase flows with high density ratio. The variational finite element technique relies on the Allen-Cahn phase-field equation for capturing the phase interface on a fixed Eulerian mesh with mass conservative and energy-stable discretization. The mass conservation is achieved by enforcing a Lagrange multiplier which has both temporal and spatial dependence on the underlying solution of the phase-field equation. To make the scheme energy-stable in a variational sense, we discretize the spatial part of the Lagrange multiplier in the phase-field equation by the mid-point approximation. The proposed variational technique is designed to reduce the spurious and unphysical oscillations in the solution while maintaining the second-order accuracy of both spatial and temporal discretizations. We integrate the Allen-Cahn phase-field equation with the incompressible Navier-Stokes equations for modeling a broad range of two-phase flow and fluid-fluid interface problems. The coupling of the implicit discretizations corresponding to the phase-field and the incompressible flow equations is achieved via nonlinear partitioned iterative procedure. Comparison of results between the standard linear stabilized finite element method and the present variational formulation shows a remarkable reduction of oscillations in the solution while retaining the boundedness of the phase-indicator field. We perform a standalone test to verify the accuracy and stability of the Allen-Cahn two-phase solver. We examine the convergence and accuracy properties of the coupled phase-field solver through the standard benchmarks of the Laplace-Young law and a sloshing tank problem. Two- and three-dimensional dam break problems are simulated to assess the capability of the phase-field solver for complex air-water interfaces involving topological changes on unstructured meshes. Finally, we demonstrate the phase

Two-phase flow

International Nuclear Information System (INIS)

Olive, J.

1990-01-01

The design, operation and safety of nuclear components requires increasingly accurate knowledge of two-phase flows. This knowledge is also necessary for some studies related to electricity applications. The author presents some concrete examples showing the range of problems and the complexity of the phenomena involved in these types of flows. Then, the basic principles of their numerical modelling are explained, as well as the new tendency to use increasingly local and refined models. The newest computer codes developed at EDF are briefly presented. Experimental studies dealing with twophase flow are also referred to, and their connections to numerical modelling are explained. Emphasis is placed on the major efforts devoted to the development of new test rigs and instrumentation [fr

Bubbly flows around a two-dimensional circular cylinder

Science.gov (United States)

Lee, Jubeom; Park, Hyungmin

2016-11-01

Two-phase cross flows around a bluff body occur in many thermal-fluid systems like steam generators, heat exchangers and nuclear reactors. However, our current knowledge on the interactions among bubbles, bubble-induced flows and the bluff body are limited. In the present study, the gas-liquid bubbly flows around a solid circular cylinder are experimentally investigated while varying the mean void fraction from 5 to 27%. The surrounding liquid (water) is initially static and the liquid flow is only induced by the air bubbles. For the measurements, we use the high-speed two-phase particle image velocimetry techniques. First, depending on the mean void fraction, two regimes are classified with different preferential concentration of bubbles in the cylinder wake, which are explained in terms of hydrodynamic force balances acting on rising bubbles. Second, the differences between the two-phase and single-phase flows (while matching their Reynolds numbers) around a circular cylinder will be discussed in relation to effects of bubble dynamics and the bubble-induced turbulence on the cylinder wake. Supported by a Grant (MPSS-CG-2016-02) through the Disaster and Safety Management Institute funded by Ministry of Public Safety and Security of Korean government.

Parallel Computing Characteristics of Two-Phase Thermal-Hydraulics code, CUPID

International Nuclear Information System (INIS)

Lee, Jae Ryong; Yoon, Han Young

2013-01-01

Parallelized CUPID code has proved to be able to reproduce multi-dimensional thermal hydraulic analysis by validating with various conceptual problems and experimental data. In this paper, the characteristics of the parallelized CUPID code were investigated. Both single- and two phase simulation are taken into account. Since the scalability of a parallel simulation is known to be better for fine mesh system, two types of mesh system are considered. In addition, the dependency of the preconditioner for matrix solver was also compared. The scalability for the single-phase flow is better than that for two-phase flow due to the less numbers of iterations for solving pressure matrix. The CUPID code was investigated the parallel performance in terms of scalability. The CUPID code was parallelized with domain decomposition method. The MPI library was adopted to communicate the information at the interface cells. As increasing the number of mesh, the scalability is improved. For a given mesh, single-phase flow simulation with diagonal preconditioner shows the best speedup. However, for the two-phase flow simulation, the ILU preconditioner is recommended since it reduces the overall simulation time

CFD model of diabatic annular two-phase flow using the Eulerian–Lagrangian approach

International Nuclear Information System (INIS)

Li, Haipeng; Anglart, Henryk

2015-01-01

Highlights: • A CFD model of annular two-phase flow with evaporating liquid film has been developed. • A two-dimensional liquid film model is developed assuming that the liquid film is sufficiently thin. • The liquid film model is coupled to the gas core flow, which is represented using the Eulerian–Lagrangian approach. - Abstract: A computational fluid dynamics (CFD) model of annular two-phase flow with evaporating liquid film has been developed based on the Eulerian–Lagrangian approach, with the objective to predict the dryout occurrence. Due to the fact that the liquid film is sufficiently thin in the diabatic annular flow and at the pre-dryout conditions, it is assumed that the flow in the wall normal direction can be neglected, and the spatial gradients of the dependent variables tangential to the wall are negligible compared to those in the wall normal direction. Subsequently the transport equations of mass, momentum and energy for liquid film are integrated in the wall normal direction to obtain two-dimensional equations, with all the liquid film properties depth-averaged. The liquid film model is coupled to the gas core flow, which currently is represented using the Eulerian–Lagrangian technique. The mass, momentum and energy transfers between the liquid film, gas, and entrained droplets have been taken into account. The resultant unified model for annular flow has been applied to the steam–water flow with conditions typical for a Boiling Water Reactor (BWR). The simulation results for the liquid film flow rate show favorable agreement with the experimental data, with the potential to predict the dryout occurrence based on criteria of critical film thickness or critical film flow rate

Drift flux model as approximation of two fluid model for two phase dispersed and slug flow in tube

International Nuclear Information System (INIS)

Nigmatulin, R.I.

1995-01-01

The analysis of one-dimensional schematizing for non-steady two-phase dispersed and slug flow in tube is presented. Quasi-static approximation, when inertia forces because of the accelerations of the phases may be neglected, is considered. Gas-liquid bubbly and slug vertical upward flows are analyzed. Non-trivial theoretical equations for slip velocity for these flows are derived. Juxtaposition of the derived equations for slip velocity with the famous Zuber-Findlay correlation as cross correlation coefficients is criticized. The generalization of non-steady drift flux Wallis theory taking into account influence of wall friction on the bubbly or slug flows for kinematical waves is considered

Drift flux model as approximation of two fluid model for two phase dispersed and slug flow in tube

Energy Technology Data Exchange (ETDEWEB)

Nigmatulin, R.I.

1995-09-01

The analysis of one-dimensional schematizing for non-steady two-phase dispersed and slug flow in tube is presented. Quasi-static approximation, when inertia forces because of the accelerations of the phases may be neglected, is considered. Gas-liquid bubbly and slug vertical upward flows are analyzed. Non-trivial theoretical equations for slip velocity for these flows are derived. Juxtaposition of the derived equations for slip velocity with the famous Zuber-Findlay correlation as cross correlation coefficients is criticized. The generalization of non-steady drift flux Wallis theory taking into account influence of wall friction on the bubbly or slug flows for kinematical waves is considered.

Fluid dynamics of cryogenic two-phase flows

International Nuclear Information System (INIS)

Verfondern, K.; Jahn, W.

2004-01-01

The objective of this study was to examine the flow behavior of a methane hydrate/methane-liquid hydrogen dispersed two-phase fluid through a given design of a moderator chamber for the ESS target system. The calculations under simplified conditions, e.g., taking no account of heat input from outside, have shown that the computer code used, CFX, was able to simulate the behavior of the two-phase flow through the moderator chamber, producing reasonable results up to a certain level of the solid phase fraction, that allowed a continuous flow process through the chamber. Inlet flows with larger solid phase fractions than 40 vol% were found to be a ''problem'' for the computer code. From the computer runs based on fractions between 20 and 40 vol%, it was observed that with increasing solid phase fraction at the inlet, the resulting flow pattern revealed a strong tendency for blockage within the chamber, supported by the ''heavy weight'' of the pellets compared to the carrying liquid. Locations which are prone to the development of such uneven flow behavior are the areas around the turning points in the semispheres and near the exit of the moderator. The considered moderator chamber with horizontal inlet and outlet flow for a solid-liquid two-phase fluid does not seem to be an appropriate design. (orig.)

Development of gas-liquid two-phase flow measurement technique in narrow channel. Application of micro wire-mesh sensor to the flow between parallel plates

International Nuclear Information System (INIS)

Ito, Daisuke; Kikura, Hiroshige; Aritomi, Masanori

2009-01-01

A novel two-phase flow measuring technique based on local electrical conductivity measurement was developed for clarifications of three-dimensional flow structure in gas-liquid two-phase flow in a narrow channel. The measuring method applies the principle of conventional wire-mesh tomography, which can measure the instantaneous void fraction distributions in a cross-section of a flow channel. In this technique, the electrodes are fixed on the inside of the walls facing each other, and the local void fractions were obtained by the electrical conductivity measurement between electrodes arranged on each wall. Therefore, the flow structure and the bubble behavior can be investigated by three-dimensional void fraction distributions in the channel with narrow gap. In this paper, a micro Wire-Mesh Sensor (μWMS) which has the gap of 3 mm was developed, and the instantaneous void fraction distributions were measured. From the measured distributions, three-dimensional bubble distributions were reconstructed, and bubble volumes and bubble velocities were estimated. (author)

A Well-Posed Two Phase Flow Model and its Numerical Solutions for Reactor Thermal-Fluids Analysis

Energy Technology Data Exchange (ETDEWEB)

Kadioglu, Samet Y. [Idaho National Lab. (INL), Idaho Falls, ID (United States); Berry, Ray [Idaho National Lab. (INL), Idaho Falls, ID (United States); Martineau, Richard [Idaho National Lab. (INL), Idaho Falls, ID (United States)

2016-08-01

A 7-equation two-phase flow model and its numerical implementation is presented for reactor thermal-fluids applications. The equation system is well-posed and treats both phases as compressible flows. The numerical discretization of the equation system is based on the finite element formalism. The numerical algorithm is implemented in the next generation RELAP-7 code (Idaho National Laboratory (INL)’s thermal-fluids code) built on top of an other INL’s product, the massively parallel multi-implicit multi-physics object oriented code environment (MOOSE). Some preliminary thermal-fluids computations are presented.

A Well-Posed Two Phase Flow Model and its Numerical Solutions for Reactor Thermal-Fluids Analysis

International Nuclear Information System (INIS)

Kadioglu, Samet Y.; Berry, Ray; Martineau, Richard

2016-01-01

A 7-equation two-phase flow model and its numerical implementation is presented for reactor thermal-fluids applications. The equation system is well-posed and treats both phases as compressible flows. The numerical discretization of the equation system is based on the finite element formalism. The numerical algorithm is implemented in the next generation RELAP-7 code (Idaho National Laboratory (INL)'s thermal-fluids code) built on top of an other INL's product, the massively parallel multi-implicit multi-physics object oriented code environment (MOOSE). Some preliminary thermal-fluids computations are presented.

Vortex-Induced Waves in Two-Phase Liquid-Liquid Flows past Bluff Body

Science.gov (United States)

Zainal Abidin, M. I. I.; Park, Kyeong H.; Angeli, Panagiota; Xie, Zhihua; Kahouadji, Lyes; Matar, Omar K.

2017-11-01

Transverse cylinders of various sizes are used to generate vortex-induced interfacial waves in two-phase oil-water flows and to influence flow pattern transitions. The vortex shedding properties at different cylinder sizes and the resulting induced waves are studied experimentally with Particle Image Velocimetry (PIV) and high-speed imaging. The system consists of a 7 m long horizontal 37 mm ID acrylic pipe and different cylinders with diameters of 2, 5 and 8 mm, located in the water phase, 460 mm after the two phases come into contact. The cylinder generates waves with frequencies similar to the von Karman vortices and changes in vortex shedding properties at different cylinder size are reflected on the resulting interfacial wave characteristics. The presence of the transverse cylinder actuates the transition from stratified to dispersed flows; the boundary between the two patterns is shifted to lower mixture velocity with increasing cylinder size. Three-dimensional numerical simulation of the system is developed to assist in designing new system. Project funded by EPSRC UK and Memphis Grant.

Apparatus for monitoring two-phase flow

Science.gov (United States)

Sheppard, John D.; Tong, Long S.

1977-03-01

A method and apparatus for monitoring two-phase flow is provided that is particularly related to the monitoring of transient two-phase (liquid-vapor) flow rates such as may occur during a pressurized water reactor core blow-down. The present invention essentially comprises the use of flanged wire screens or similar devices, such as perforated plates, to produce certain desirable effects in the flow regime for monitoring purposes. One desirable effect is a measurable and reproducible pressure drop across the screen. The pressure drop can be characterized for various known flow rates and then used to monitor nonhomogeneous flow regimes. Another useful effect of the use of screens or plates in nonhomogeneous flow is that such apparatus tends to create a uniformly dispersed flow regime in the immediate downstream vicinity. This is a desirable effect because it usually increases the accuracy of flow rate measurements determined by conventional methods.

A generalized volumetric dispersion model for a class of two-phase separation/reaction: finite difference solutions

Science.gov (United States)

Siripatana, Chairat; Thongpan, Hathaikarn; Promraksa, Arwut

2017-03-01

This article explores a volumetric approach in formulating differential equations for a class of engineering flow problems involving component transfer within or between two phases. In contrast to conventional formulation which is based on linear velocities, this work proposed a slightly different approach based on volumetric flow-rate which is essentially constant in many industrial processes. In effect, many multi-dimensional flow problems found industrially can be simplified into multi-component or multi-phase but one-dimensional flow problems. The formulation is largely generic, covering counter-current, concurrent or batch, fixed and fluidized bed arrangement. It was also intended to use for start-up, shut-down, control and steady state simulation. Since many realistic and industrial operation are dynamic with variable velocity and porosity in relation to position, analytical solutions are rare and limited to only very simple cases. Thus we also provide a numerical solution using Crank-Nicolson finite difference scheme. This solution is inherently stable as tested against a few cases published in the literature. However, it is anticipated that, for unconfined flow or non-constant flow-rate, traditional formulation should be applied.

Two-phase flow characterisation by nuclear magnetic resonance

International Nuclear Information System (INIS)

Leblond, J.; Javelot, S.; Lebrun, D.; Lebon, L.

1998-01-01

The results presented in this paper demonstrate the performance of the PFGSE-NMR to obtain a complete characterisation of two-phase flows. Different methods are proposed to characterise air-water flows in different regimes: stationary two-phase flows and flows in transient condition. Finally a modified PFGSE is proposed to analyse the turbulence of air-water bubbly flow. (author)

A two-dimensional, transient, compressible isothermal and two-phase model for the air-side electrode of PEM fuel cells

International Nuclear Information System (INIS)

Khakbaz Baboli, M.; Kermani, M.J.

2008-01-01

A two-dimensional, transient, compressible, isothermal and two-phase flow of reactant-product mixture in the air-side electrode of proton exchange membrane fuel cells (PEMFC) are numerically studied in the present paper. The mixture is composed of four species: oxygen, nitrogen, liquid water and water vapor. The governing PDE's are conservation of the water vapor and oxygen species, momentum equation of the mixture (gas+liquid), mass conservation of the liquid phase, and mass conservation of the mixture. In this study, a separate PDE for the mass conservation of the liquid water is solved to calculate the saturation levels. The capillary pressure was used to determine the slip velocity between the phases. A full compressible form of the momentum equation was used, with the ∇.V preserved in the equation. The Maxwell-Stefan equation was used to model the diffusive fluxes of the multi-component gas mixture. The strongly coupled equations are solved based on a recently developed finite volume SIMPLER scheme of S.V. Patankar, Numerical Heat Transfer and Fluid Flow, Hemisphere Publishing Corp., McGraw-Hill Book Company, 1984. The computational domain consists of two regions; an open area (gas delivery channel) linked to a porous gas diffusion layer (GDL). A single (unified) set of the PDE's are used for the whole domain with the corresponding properties of each sub-domain. A polarization curve for the whole spectrum of the dry and wet regions were obtained. The results were compared with the experiments of E.A. Ticianelli, C.R. Derouin, A. Redondo, S. Srinivasan, J. Electrochem. Soc. 135 (1988) 2209, and good agreements were achieved

Measurement system of bubbly flow using ultrasonic velocity profile monitor and video data processing unit. 2. Flow characteristics of bubbly countercurrent flow

International Nuclear Information System (INIS)

Aritomi, Masanori; Zhou, Shirong; Nakajima, Makoto; Takeda, Yasushi; Mori, Michitsugu.

1997-01-01

The authors have developed a measurement system which is composed of an ultrasonic velocity profile monitor and a video data processing unit in order to clarify its multi-dimensional flow characteristics in bubbly flows and to offer a data base to validate numerical codes for multi-dimensional two-phase flow. In this paper, the measurement system was applied for bubbly countercurrent flows in a vertical rectangular channel. At first, both bubble and water velocity profiles and void fraction profiles in the channel were investigated statistically. Next, turbulence intensity in a continuous liquid phase was defined as a standard deviation of velocity fluctuation, and the two-phase multiplier profile of turbulence intensity in the channel was clarified as a ratio of the standard deviation of flow fluctuation in a bubbly countercurrent flow to that in a water single phase flow. Finally, the distribution parameter and drift velocity used in the drift flux model for bubbly countercurrent flows were calculated from the obtained velocity profiles of both phases and void fraction profile, and were compared with the correlation proposed for bubbly countercurrent flows. (author)

Studies on turbulence structure and liquid film behavior in annular two-phase flow flowing in a throat section

International Nuclear Information System (INIS)

Yoshida, Kenji; Miyabe, Masaya; Matsumoto, Tadayoshi; Kataoka, Isao; Ohmori, Shuichi; Mori, Michitsugu

2004-01-01

Experimental studies on turbulence structure and liquid film behavior in annular two-phase flow were carried out concerned with the steam injector systems for a next-generation nuclear reactor. In the steam injector, steam/water annular two-phase flow is formed at the mixing nozzle. To make an appropriate design for high-performance steam injector system, it is very important to accumulate the fundamental data of thermo-hydro dynamic characteristics of annular flow in the steam injector. Especially, the turbulence modification in multi-phase flow due to the phase interaction is one of the most important phenomena and has attracted research attention. In this study, the liquid film behavior and the resultant turbulence modification due to the phase interaction were investigated. The behavior of the interfacial waves on liquid film flow such as the ripple or disturbance waves were observed to make clear the interfacial velocity and the special structure of the interfacial waves by using the high-speed video camera and the digital camera. The measurements for gas-phase velocity profiles and turbulent intensity in annular flow passing through the throat section were precisely performed to investigate quantitatively the turbulent modification in annular flow by using the constant temperature hot-wire anemometer. The measurements for liquid film thickness by the electrode needle method were also carried out. (author)

Encyclopedia of two-phase heat transfer and flow II special topics and applications

CERN Document Server

Kim, Jungho

2015-01-01

The aim of the two–set series is to present a very detailed and up–to–date reference for researchers and practicing engineers in the fields of mechanical, refrigeration, chemical, nuclear and electronics engineering on the important topic of two-phase heat transfer and two-phase flow. The scope of the first set of 4 volumes presents the fundamentals of the two-phase flows and heat transfer mechanisms, and describes in detail the most important prediction methods, while the scope of the second set of 4 volumes presents numerous special topics and numerous applications, also including numerical simulation methods. Practicing engineers will find extensive coverage to applications involving: multi-microchannel evaporator cold plates for electronics cooling, boiling on enhanced tubes and tube bundles, flow pattern based methods for predicting boiling and condensation inside horizontal tubes, pressure drop methods for singularies (U-bends and contractions), boiling in multiport tubes, and boiling and condens...

Application of a method for comparing one-dimensional and two-dimensional models of a ground-water flow system

International Nuclear Information System (INIS)

Naymik, T.G.

1978-01-01

To evaluate the inability of a one-dimensional ground-water model to interact continuously with surrounding hydraulic head gradients, simulations using one-dimensional and two-dimensional ground-water flow models were compared. This approach used two types of models: flow-conserving one-and-two dimensional models, and one-dimensional and two-dimensional models designed to yield two-dimensional solutions. The hydraulic conductivities of controlling features were varied and model comparison was based on the travel times of marker particles. The solutions within each of the two model types compare reasonably well, but a three-dimensional solution is required to quantify the comparison

Dynamic simulation of dispersed gas-liquid two-phase flow using a discrete bubble model.

NARCIS (Netherlands)

Delnoij, E.; Lammers, F.A.; Kuipers, J.A.M.; van Swaaij, Willibrordus Petrus Maria

1997-01-01

In this paper a detailed hydrodynamic model for gas-liquid two-phase flow will be presented. The model is based on a mixed Eulerian-Lagrangian approach and describes the time-dependent two-dimensional motion of small, spherical gas bubbles in a bubble column operating in the homogeneous regime. The

Apparatus for monitoring two-phase flow

International Nuclear Information System (INIS)

Sheppard, J.D.; Tong, L.S.

1977-01-01

A method and apparatus for monitoring two-phase flow is provided that is particularly related to the monitoring of transient two-phase (liquid-vapor) flow rates such as may occur during a pressurized water reactor core blow-down. The present invention essentially comprises the use of flanged wire screens or similar devices, such as perforated plates, to produce certain desirable effects in the flow regime for monitoring purposes. One desirable effect is a measurable and reproducible pressure drop across the screen. The pressure drop can be characterized for various known flow rates and then used to monitor nonhomogeneous flow regimes. Another useful effect of the use of screens or plates in nonhomogeneous flow is that such apparatus tends to create a uniformly dispersed flow regime in the immediate downstream vicinity. This is a desirable effect because it usually increases the accuracy of flow rate measurements determined by conventional methods. 3 claims, 9 figures

Two-phase flow patterns and their relationship to two-phase heat transfer

International Nuclear Information System (INIS)

Hewitt, G.F.

1977-01-01

The objective of this lecture was to discuss the general nature of two phase flows, to define the various regimes of flow and to discuss the influence of these regimes on the heat transfer processes taking place. The methods of regime delineation are briefly described and regime descriptions introduced for both vertical and horizontal flows in tubes. ''Flow regime maps'' have been widely used as an aid to determination of the regime which occurs in a given situation. Some of the more widely used maps are described and the limitations of this approach discussed. There have been many attempts to obtain a better phenomenological description of two phase flow patterns. In this lecture, these attempts will be reviewed in the context of the bubble/plug, plug/churn and churn/annular flow transitions in vertical flow. The latter two transitions are related to the flooding/flow reversal phenomena. For horizontal flows, recent work on the onset of slugging will be reviewed. In flows with evaporation or condensation, the situation is influenced by departures from thermodynamic equilibrium and the types of departure observed are discuss briefly. Flow patterns and their relationships with heat transfer regimes are then reviewed for the case of condensation in horizontal tubes and evaporation in vertical tubes

Two-dimensional heat flow apparatus

Science.gov (United States)

McDougall, Patrick; Ayars, Eric

2014-06-01

We have created an apparatus to quantitatively measure two-dimensional heat flow in a metal plate using a grid of temperature sensors read by a microcontroller. Real-time temperature data are collected from the microcontroller by a computer for comparison with a computational model of the heat equation. The microcontroller-based sensor array allows previously unavailable levels of precision at very low cost, and the combination of measurement and modeling makes for an excellent apparatus for the advanced undergraduate laboratory course.

Development of a Laser Dopper Anemometer technique for the measurement of two phase dispersed flow

International Nuclear Information System (INIS)

Srinivasan, J.

1978-05-01

A new optical technique using Laser-Doppler Anemometry is presented for the measurement of the local number densities and two-dimensional velocity probability densities of a turbulent dilute two-phase dispersion which has a distribution of particle size and a predominant direction of flow. This technique establishes that by a suitable scheme of discrimination on the signal amplitude, residence time and frequency of the Doppler signals caused by the scattered light from individual particles in the probing volume, the size distribution of moderately large particles in a dilute dispersed flow can be determined. The newly developed Laser-Doppler Anemometer (LDA) technique was applied to a solid particle-water two-phase flow and a water droplet-air two-phase flow. Particular emphasis was placed on turbulent two-phase water droplet-air flow inside a vertical rectangular channel. At each of nine different measuring locations along the transverse axis (starting at 250μ from the channel wall), over 20,000 Doppler signals were individually examined. The particle size and number density distributions, and the axial and lateral velocity distributions of both phases are reported. The analysis reveals some interesting features of two-phase dispersed flow. A film of water on the channel wall was formed due to the deposition of droplets from the flow. The water droplet entrainment from the wall film and the subsequent breakup of some of these into the flow are discussed. A discussion of the relationship between the particle distributions and turbulent flow characteristics is presented

Approximate Riemann solvers and flux vector splitting schemes for two-phase flow

International Nuclear Information System (INIS)

Toumi, I.; Kumbaro, A.; Paillere, H.

1999-01-01

These course notes, presented at the 30. Von Karman Institute Lecture Series in Computational Fluid Dynamics, give a detailed and through review of upwind differencing methods for two-phase flow models. After recalling some fundamental aspects of two-phase flow modelling, from mixture model to two-fluid models, the mathematical properties of the general 6-equation model are analysed by examining the Eigen-structure of the system, and deriving conditions under which the model can be made hyperbolic. The following chapters are devoted to extensions of state-of-the-art upwind differencing schemes such as Roe's Approximate Riemann Solver or the Characteristic Flux Splitting method to two-phase flow. Non-trivial steps in the construction of such solvers include the linearization, the treatment of non-conservative terms and the construction of a Roe-type matrix on which the numerical dissipation of the schemes is based. Extension of the 1-D models to multi-dimensions in an unstructured finite volume formulation is also described; Finally, numerical results for a variety of test-cases are shown to illustrate the accuracy and robustness of the methods. (authors)

Entrainment in vertical annular two-phase flow

International Nuclear Information System (INIS)

Sawant, Pravin; Ishii, Mamoru; Mori, Michitsugu

2009-01-01

Prediction of amount of entrained droplets or entrainment fraction in annular two-phase flow is essential for the estimation of dryout condition and analysis of post dryout heat transfer in light water nuclear reactors and steam boilers. In this study, air-water and organic fluid (Freon-113) annular flow entrainment experiments have been carried out in 9.4 and 10.2 mm diameter test sections, respectively. Both the experiments covered three distinct pressure conditions and wide range of liquid and gas flow conditions. The organic fluid experiments simulated high pressure steam-water annular flow conditions. In each of the experiments, measurements of entrainment fraction, droplet entrainment rate and droplet deposition rate have been performed by using a liquid film extraction method. A simple, explicit and non-dimensional correlation developed by Sawant et al. (2008a) for the prediction of entrainment fraction is further improved in this study in order to account for the existence of critical gas and liquid flow rates below which no entrainment is possible. Additionally, a new correlation is proposed for the estimation of minimum liquid film flow rate at the maximum entrainment fraction condition. The improved correlation successfully predicted the newly collected air-water and Freon-113 entrainment fraction data. Furthermore, the correlations satisfactorily compared with the air-water, helium-water and air-genklene experimental data measured by Willetts (1987). (author)

Ferrofluid-in-oil two-phase flow patterns in a flow-focusing microchannel

Science.gov (United States)

Sheu, T. S.; Chen, Y. T.; Lih, F. L.; Miao, J. M.

This study investigates the two-phase flow formation process of water-based Fe3O4 ferrofluid (dispersed phase) in a silicon oil (continuous phase) flow in the microfluidic flow-focusing microchannel under various operational conditions. With transparent PDMS chip and optical microscope, four main two-phase flow patterns as droplet flow, slug flow, ring flow and churn flow are observed. The droplet shape, size, and formation mechanism were also investigated under different Ca numbers and intended to find out the empirical relations. The paper marks an original flow pattern map of the ferrofluid-in-oil flows in the microfluidic flow-focusing microchannels. The flow pattern transiting from droplet flow to slug flow appears for an operational conditions of QR < 1 and Lf / W < 1. The power law index that related Lf / W to QR was 0.36 in present device.

Sugarcane molasses-based bio-ethanol wastewater treatment by two-phase multi-staged up-flow anaerobic sludge blanket (UASB) combination with up-flow UASB and down-flow hanging sponge.

Science.gov (United States)

Choeisai, P; Jitkam, N; Silapanoraset, K; Yubolsai, C; Yoochatchaval, W; Yamaguchi, T; Onodera, T; Syutsubo, K

2014-01-01

This study was designed to evaluate a treatment system for high strength wastewater (vinasse) from a sugarcane molasses-based bio-ethanol plant in Thailand. A laboratory-scale two-phase treatment system composed of a sulfate reducing (SR) tank and multi-staged up-flow anaerobic sludge blanket (MS-UASB) reactor was used as the pre-treatment unit. Conventional UASB and down-flow hanging sponge (DHS) reactors were used as the post-treatment unit. The treatment system was operated for 300 days under ambient temperature conditions (24.6-29.6 °C). The hydraulic retention time (HRT) in each unit was kept at 25 h for the two-phase system and 23 h for the UASB&DHS. The influent concentration was allowed to reach up to 15,000 mg chemical oxygen demand (COD)/L. COD removal efficiency (based on influent COD) of the two-phase MS-UASB and the UASB&DHS was 54.9 and 18.7%, respectively. Due to the effective removal of sulfide in the SR tank, the MS-UASB achieved a high methane conversion ratio of up to 97%. In DHS, nitrification occurred at the outside portion of the sponge media while denitrification occurred at the inside. Consequently, 27% of the total nitrogen (TN) was removed. An amount of 32% of residual nitrogen (28 mgN/L) was in the form of nitrate, a better nitrogen state for fertilizer.

DISTRIBUTION OF TWO-PHASE FLOW IN A DISTRIBUTOR

Directory of Open Access Journals (Sweden)

AZRIDJAL AZIZ

2012-02-01

Full Text Available The flow configuration and distribution behavior of two-phase flow in a distributor made of acrylic resin have been investigated experimentally. In this study, air and water were used as two-phase flow working fluids. The distributor consists of one inlet and two outlets, which are set as upper and lower, respectively. The flow visualization at the distributor was made by using a high–speed camera. The flow rates of air and water flowing out from the upper and lower outlet branches were measured. Effects of inclination angle of the distributor were investigated. By changing the inclination angle from vertical to horizontal, uneven distributions were also observed. The distribution of two-phase flow through distributor tends even flow distribution on the vertical position and tends uneven distribution on inclined and horizontal positions. It is shown that even distribution could be achieved at high superficial velocities of both air and water.

Coherent structures and flow topology of transitional separated-reattached flow over two and three dimensional geometrical shapes

Science.gov (United States)

Diabil, Hayder Azeez; Li, Xin Kai; Abdalla, Ibrahim Elrayah

2017-09-01

Large-scale organized motions (commonly referred to coherent structures) and flow topology of a transitional separated-reattached flow have been visualised and investigated using flow visualisation techniques. Two geometrical shapes including two-dimensional flat plate with rectangular leading edge and three-dimensional square cylinder are chosen to shed a light on the flow topology and present coherent structures of the flow over these shapes. For both geometries and in the early stage of the transition, two-dimensional Kelvin-Helmholtz rolls are formed downstream of the leading edge. They are observed to be twisting around the square cylinder while they stay flat in the case of the two-dimensional flat plate. For both geometrical shapes, the two-dimensional Kelvin-Helmholtz rolls move downstream of the leading edge and they are subjected to distortion to form three-dimensional hairpin structures. The flow topology in the flat plate is different from that in the square cylinder. For the flat plate, there is a merging process by a pairing of the Kelvin-Helmholtz rolls to form a large structure that breaks down directly into many hairpin structures. For the squire cylinder case, the Kelvin-Helmholtz roll evolves topologically to form a hairpin structure. In the squire cylinder case, the reattachment length is much shorter and a forming of the three-dimensional structures is closer to the leading edge than that in the flat plate case.

Pore-scale Simulation and Imaging of Multi-phase Flow and Transport in Porous Media (Invited)

Science.gov (United States)

Crawshaw, J.; Welch, N.; Daher, I.; Yang, J.; Shah, S.; Grey, F.; Boek, E.

2013-12-01

-NMR experiments. We then use our preferred multi-phase model to directly calculate flow in pore space images of two different sandstones and observe excellent agreement with experimental relative permeabilities. Also we calculate cluster size distributions in good agreement with experimental studies. Our analysis shows that the simulations are able to predict both multi-phase flow and transport properties directly on large 3D pore space images of real rocks. Pore space images, left and velocity distributions, right (Yang and Boek, 2013)

Review on two-phase flow instabilities in narrow spaces

International Nuclear Information System (INIS)

Tadrist, L.

2007-01-01

Instabilities in two-phase flow have been studied since the 1950s. These phenomena may appear in power generation and heat transfer systems where two-phase flow is involved. Because of thermal management in small size systems, micro-fluidics plays an important role. Typical processes must be considered when the channel hydraulic diameter becomes very small. In this paper, a brief review of two-phase flow instabilities encountered in channels having hydraulic diameters greater than 10 mm are presented. The main instability types are discussed according to the existing experimental results and models. The second part of the paper examines two-phase flow instabilities in narrow spaces. Pool and flow boiling cases are considered. Experiments as well as theoretical models existing in the literature are examined. It was found that several experimental works evidenced these instabilities meanwhile only limited theoretical developments exist in the literature. In the last part of the paper an interpretation of the two-phase flow instabilities linked to narrow spaces are presented. This approach is based on characteristic time scales of the two-phase flow and bubble growth in the capillaries

Parallel two-phase-flow-induced vibrations in fuel pin model

International Nuclear Information System (INIS)

Hara, Fumio; Yamashita, Tadashi

1978-01-01

This paper reports the experimental results of vibrations of a fuel pin model -herein meaning the essential form of a fuel pin from the standpoint of vibration- in a parallel air-and-water two-phase flow. The essential part of the experimental apparatus consisted of a flat elastic strip made of stainless steel, both ends of which were firmly supported in a circular channel conveying the two-phase fluid. Vibrational strain of the fuel pin model, pressure fluctuation of the two-phase flow and two-phase-flow void signals were measured. Statistical measures such as power spectral density, variance and correlation function were calculated. The authors obtained (1) the relation between variance of vibrational strain and two-phase-flow velocity, (2) the relation between variance of vibrational strain and two-phase-flow pressure fluctuation, (3) frequency characteristics of variance of vibrational strain against the dominant frequency of the two-phase-flow pressure fluctuation, and (4) frequency characteristics of variance of vibrational strain against the dominant frequency of two-phase-flow void signals. The authors conclude that there exist two kinds of excitation mechanisms in vibrations of a fuel pin model inserted in a parallel air-and-water two-phase flow; namely, (1) parametric excitation, which occurs when the fundamental natural frequency of the fuel pin model is related to the dominant travelling frequency of water slugs in the two-phase flow by the ratio 1/2, 1/1, 3/2 and so on; and (2) vibrational resonance, which occurs when the fundamental frequency coincides with the dominant frequency of the two-phase-flow pressure fluctuation. (auth.)

Thermoplastic Micromodel Investigation of Two-Phase Flows in a Fractured Porous Medium

Directory of Open Access Journals (Sweden)

Shao-Yiu Hsu

2017-01-01

Full Text Available In the past few years, micromodels have become a useful tool for visualizing flow phenomena in porous media with pore structures, e.g., the multifluid dynamics in soils or rocks with fractures in natural geomaterials. Micromodels fabricated using glass or silicon substrates incur high material cost; in particular, the microfabrication-facility cost for making a glass or silicon-based micromold is usually high. This may be an obstacle for researchers investigating the two-phase-flow behavior of porous media. A rigid thermoplastic material is a preferable polymer material for microfluidic models because of its high resistance to infiltration and deformation. In this study, cyclic olefin copolymer (COC was selected as the substrate for the micromodel because of its excellent chemical, optical, and mechanical properties. A delicate micromodel with a complex pore geometry that represents a two-dimensional (2D cross-section profile of a fractured rock in a natural oil or groundwater reservoir was developed for two-phase-flow experiments. Using an optical visualization system, we visualized the flow behavior in the micromodel during the processes of imbibition and drainage. The results show that the flow resistance in the main channel (fracture with a large radius was higher than that in the surrounding area with small pore channels when the injection or extraction rates were low. When we increased the flow rates, the extraction efficiency of the water and oil in the mainstream channel (fracture did not increase monotonically because of the complex two-phase-flow dynamics. These findings provide a new mechanism of residual trapping in porous media.

Two-phase flow model with nonequilibrium and critical flow

International Nuclear Information System (INIS)

Sureau, H.; Houdayer, G.

1976-01-01

The model proposed includes the three conservation equations (mass, momentum, energy) applied to the two phase flows and a fourth partial derivative equation which takes into account the nonequilibriums and describes the mass transfer process. With this model, the two phase critical flow tests performed on the Moby-Dick loop (CENG) with several geometries, are interpreted by a unique law. Extrapolations to industrial dimension problems show that geometry and size effects are different from those obtained with earlier models (Zaloudek, Moody, Fauske) [fr

Two-phase flow regimes and mechanisms of critical heat flux under subcooled flow boiling conditions

International Nuclear Information System (INIS)

Le Corre, Jean-Marie; Yao, Shi-Chune; Amon, Cristina H.

2010-01-01

A literature review of critical heat flux (CHF) experimental visualizations under subcooled flow boiling conditions was performed and systematically analyzed. Three major types of CHF flow regimes were identified (bubbly, vapor clot and slug flow regime) and a CHF flow regime map was developed, based on a dimensional analysis of the phenomena and available experimental information. It was found that for similar geometric characteristics and pressure, a Weber number (We)/thermodynamic quality (x) map can be used to predict the CHF flow regime. Based on the experimental observations and the review of the available CHF mechanistic models under subcooled flow boiling conditions, hypothetical CHF mechanisms were selected for each CHF flow regime, all based on a concept of wall dry spot overheating, rewetting prevention and subsequent dry spot spreading. Even though the selected concept has not received much attention (in term or theoretical developments and applications) as compared to other more popular DNB models, its basis have often been cited by experimental investigators and is considered by the authors as the 'most-likely' mechanism based on the literature review and analysis performed in this work. The selected modeling concept has the potential to span the CHF conditions from highly subcooled bubbly flow to early stage of annular flow and has been numerically implemented and validated in bubbly flow and coupled with one- and three-dimensional (CFD) two-phase flow codes, in a companion paper. [Le Corre, J.M., Yao, S.C., Amon, C.H., in this issue. A mechanistic model of critical heat flux under subcooled flow boiling conditions for application to one and three-dimensional computer codes. Nucl. Eng. Des.].

Modeling analyses of two-phase flow instabilities for straight and helical tubes in nuclear power plants

International Nuclear Information System (INIS)

Dong, Ruiting; Niu, Fenglei; Zhou, Yuan; Yu, Yu; Guo, Zhangpeng

2016-01-01

Highlights: • Two-phase flow instabilities in straight and helical tubes were studied. • The effects of system pressure, mass flux, inlet subcooling on DWO were studied. • The simulation results are consistent with the experimental results. • The RELAP5 results are consistent with frequency domain method results. - Abstract: The effects of system pressure, mass flux and inlet subcooling on two-phase flow instability for the test section consisted of two heated straight channels or two helical channels are studied by means of RELAP5/MOD3.3 and multi-variable frequency domain control theory. The experimental data in two straight channels are used to verify the RELAP5 and multi-variable frequency domain control theory results. The thermal hydraulic behaviors and parametric effects are simulated and compared with the experimental data. The RELAP5 results show that the flow stability increases with the system pressure, mass velocity, and inlet subcooling at high subcoolings. The frequency domain theory presents the same results as those given by the time domain theory (RELAP5). The effects of system pressure, mass velocity and inlet subcooling are simulated to find the difference between the straight and the helical tube flows. The RELAP5 and the multi-variable frequency domain control theory are used in modeling and simulating density wave oscillation to study their advantages and disadvantages in straight and helical tubes.

Modeling analyses of two-phase flow instabilities for straight and helical tubes in nuclear power plants

Energy Technology Data Exchange (ETDEWEB)

Dong, Ruiting [Beijing Key Laboratory of Passive Nuclear Power Safety and Technology, North China Electric Power University, Beijing 102206 (China); Niu, Fenglei, E-mail: niufenglei@ncepu.edu.cn [Beijing Key Laboratory of Passive Nuclear Power Safety and Technology, North China Electric Power University, Beijing 102206 (China); Zhou, Yuan [School of Nuclear Science and Engineering, Shanghai Jiaotong University, Shanghai 200240 (China); Yu, Yu; Guo, Zhangpeng [Beijing Key Laboratory of Passive Nuclear Power Safety and Technology, North China Electric Power University, Beijing 102206 (China)

2016-10-15

Highlights: • Two-phase flow instabilities in straight and helical tubes were studied. • The effects of system pressure, mass flux, inlet subcooling on DWO were studied. • The simulation results are consistent with the experimental results. • The RELAP5 results are consistent with frequency domain method results. - Abstract: The effects of system pressure, mass flux and inlet subcooling on two-phase flow instability for the test section consisted of two heated straight channels or two helical channels are studied by means of RELAP5/MOD3.3 and multi-variable frequency domain control theory. The experimental data in two straight channels are used to verify the RELAP5 and multi-variable frequency domain control theory results. The thermal hydraulic behaviors and parametric effects are simulated and compared with the experimental data. The RELAP5 results show that the flow stability increases with the system pressure, mass velocity, and inlet subcooling at high subcoolings. The frequency domain theory presents the same results as those given by the time domain theory (RELAP5). The effects of system pressure, mass velocity and inlet subcooling are simulated to find the difference between the straight and the helical tube flows. The RELAP5 and the multi-variable frequency domain control theory are used in modeling and simulating density wave oscillation to study their advantages and disadvantages in straight and helical tubes.

New modeling and experimental approaches for characterization of two-phase flow interfacial structure

International Nuclear Information System (INIS)

Ishii, Mamoru; Sun, Xiaodong

2004-01-01

This paper presents new experimental and modeling approaches in characterizing interfacial structures in gas-liquid two-phase flow. For the experiments, two objective approaches are developed to identify flow regimes and to obtain local interfacial structure data. First, a global measurement technique using a non-intrusive ring-type impedance void-meter and a self-organizing neural network is presented to identify the one-dimensional'' flow regimes. In the application of this measurement technique, two methods are discussed, namely, one based on the probability density function of the impedance probe measurement (PDF input method) and the other based on the sorted impedance signals, which is essentially the cumulative probability distribution function of the impedance signals (instantaneous direct signal input method). In the latter method, the identification can be made close to instantaneously since the required signals can be acquired over a very short time period. In addition, a double-sensor conductivity probe can also be used to obtain ''local'' flow regimes by using the instantaneous direct signal input method with the bubble chord length information. Furthermore, a newly designed conductivity probe with multiple double-sensor heads is proposed to obtain ''two-dimensional'' flow regimes across the flow channel. Secondly, a state-of-the-art four-sensor conductivity probe technique has been developed to obtain detailed local interfacial structure information. The four-sensor conductivity probe accommodates the double-sensor probe capability and can be applied in a wide range of flow regimes spanning from bubbly to churn-turbulent flows. The signal processing scheme is developed such that it categorizes the acquired parameters into two groups based on bubble cord length information. Furthermore, for the modeling of the interfacial structure characterization, the interfacial area transport equation proposed earlier has been studied to provide a dynamic and

On intermittent flow characteristics of gas–liquid two-phase flow

Energy Technology Data Exchange (ETDEWEB)

Thaker, Jignesh; Banerjee, Jyotirmay, E-mail: jbaner@gmail.com

2016-12-15

Highlights: • Unified correlations for intermittent flow characteristics are developed. • Influence of inflow conditions on intermittent flow characteristics is analysed. • Developed correlations can be used for effective design of piping components. - Abstract: Flow visualisation experiments are reported for intermittent regime of gas–liquid two-phase flow. Intermittent flow characteristics, which include plug/slug frequency, liquid plug/slug velocity, liquid plug/slug length, and plug/slug bubble length are determined by image processing of flow patterns captured at a rate of 1600 frames per second (FPS). Flow characteristics are established as a function of inlet superficial velocity of both the phases (in terms of Re{sub SL} and Re{sub SG}). The experimental results are first validated with the existing correlations for slug flow available in literature. It is observed that the correlations proposed in literature for slug flow do not accurately predict the flow characteristics in the plug flow regime. The differences are clearly highlighted in this paper. Based on the measured database for both plug and slug flow regime, modified correlations for the intermittent flow regime are proposed. The correlations reported in the present paper, which also include plug flow characteristics will aid immensely to the effective design and optimization of operating conditions for safer operation of two-phase flow piping systems.

A Variational Model for Two-Phase Immiscible Electroosmotic Flow at Solid Surfaces

KAUST Repository

Shao, Sihong

2012-01-01

We develop a continuum hydrodynamic model for two-phase immiscible flows that involve electroosmotic effect in an electrolyte and moving contact line at solid surfaces. The model is derived through a variational approach based on the Onsager principle of minimum energy dissipation. This approach was first presented in the derivation of a continuum hydrodynamic model for moving contact line in neutral two-phase immiscible flows (Qian, Wang, and Sheng, J. Fluid Mech. 564, 333-360 (2006)). Physically, the electroosmotic effect can be formulated by the Onsager principle as well in the linear response regime. Therefore, the same variational approach is applied here to the derivation of the continuum hydrodynamic model for charged two-phase immiscible flows where one fluid component is an electrolyte exhibiting electroosmotic effect on a charged surface. A phase field is employed to model the diffuse interface between two immiscible fluid components, one being the electrolyte and the other a nonconductive fluid, both allowed to slip at solid surfaces. Our model consists of the incompressible Navier-Stokes equation for momentum transport, the Nernst-Planck equation for ion transport, the Cahn-Hilliard phase-field equation for interface motion, and the Poisson equation for electric potential, along with all the necessary boundary conditions. In particular, all the dynamic boundary conditions at solid surfaces, including the generalized Navier boundary condition for slip, are derived together with the equations of motion in the bulk region. Numerical examples in two-dimensional space, which involve overlapped electric double layer fields, have been presented to demonstrate the validity and applicability of the model, and a few salient features of the two-phase immiscible electroosmotic flows at solid surface. The wall slip in the vicinity of moving contact line and the Smoluchowski slip in the electric double layer are both investigated. © 2012 Global-Science Press.

Numerical flow analyses of a two-phase hydraulic coupling

Energy Technology Data Exchange (ETDEWEB)

Hur, N.; Kwak, M.; Moshfeghi, M. [Sogang University, Seoul (Korea, Republic of); Chang, C.-S.; Kang, N.-W. [VS Engineering, Seoul (Korea, Republic of)

2017-05-15

We investigated flow characteristics in a hydraulic coupling at different charged water conditions and speed ratios. Hence, simulations were performed for three-dimensional two-phase flow by using the VOF method. The realizable k-ε turbulence model was adopted. To resolve the interaction of passing blades of the primary and secondary wheels, simulations were conducted in the unsteady framework using a sliding grid technique. The results show that the water-air distribution inside the wheel is strongly dependent upon both amount of charged water and speed ratio. Generally, air is accumulated in the center of the wheel, forming a toroidal shape wrapped by the circulating water. The results also show that at high speed ratios, the solid-body-like rotation causes dry areas on the periphery of the wheels and, hence, considerably decreases the circulating flow rate and the transmitted torque. Furthermore, the momentum transfer was investigated through the concept of a mass flux triangle based on the local velocity multiplied by the local mixture density instead of the velocity triangle commonly used in a single-phase turbomachine analysis. Also, the mass fluxes along the radius of the coupling in the partially charged and fully charged cases were found to be completely different. It is shown that the flow rate at the interfacial plane and also the transmitted torque are closely related and are strongly dependent upon both the amount of charged water and speed ratio. Finally, a conceptual categorization together with two comprehensive maps was provided for the torque transmission and also circulating flow rates. These two maps in turn exhibit valuable engineering information and can serve as bases for an optimal design of a hydraulic coupling.

Mechanistic multidimensional analysis of two-phase flow in horizontal tube with 90 deg elbow

International Nuclear Information System (INIS)

Tselishcheva, E.A.; Antal, St.P.; Podowski, M.Z.; Marshall, S.

2007-01-01

The development of modeling and simulation capabilities of two-phase flow and heat transfer is very important for the design, operation and safety of nuclear reactors. Whereas a significant progress in this field has been made over the recent years, further advancements are clearly needed for new concepts of advanced (Generation-IV in particular) reactors. Difficulties in analyzing gas/liquid flows are due to the fact that such two-phase mixtures can assume several different flow patterns, each characterized by flow-regime specific interfacial phenomena of mass, momentum and energy transfer. The level of difficulty increases even further in the case of a complex tube geometries and spatial orientations. The purpose of this paper is to discuss the results of the analysis of a two-phase flow in a horizontal pipe with a 90-degree elbow. The overall objective of the present work is the development of a 3-dimensional computational model of a two-phase high-Reynolds number turbulent flow. The overall new model has been encoded in the next-generation Computational Multiphase Fluid Dynamics (CMFD) computer code, NPHASE. The model has been tested parametrically and the results of NPHASE calculations have been compared against experimental data. It has been demonstrated that the proposed model is consistent both physically and numerically, the predictions are in a reasonable agreement with the measurements

Numerical simulation of gas-liquid two-phase flow and convective heat transfer in a micro tube

International Nuclear Information System (INIS)

Fukagata, Koji; Kasagi, Nobuhide; Ua-arayaporn, Poychat; Himeno, Takehiro

2007-01-01

Numerical simulation of an air and water two-phase flow in a 20 μm ID tube is carried out. A focus is laid upon the flow and heat transfer characteristics in bubble-train flows. An axisymmetric two-dimensional flow is assumed. The finite difference method is used to solve the governing equations, while the level set method is adopted for capturing the interface of gas and liquid. In each simulation, the mean pressure gradient and the wall heat flux are kept constant. The simulation is repeated under different conditions of pressure gradient and void fraction. The superficial Reynolds numbers of gas and liquid phases studied are 0.34-13 and 16-490, respectively, and the capillary number is 0.0087-0.27. Regardless of the flow conditions, the gas-phase velocity is found approximately 1.2 times higher than the liquid-phase velocity. This is in accordance with the Armand correlation valid for two-phase flows in macro-sized tubes. The two-phase friction coefficient is found to be scaled with the Reynolds number based on the effective viscosity of the Einstein type. The computed wall temperature distribution is qualitatively similar to that observed experimentally in a mini channel. The local Nusselt number beneath the bubble is found notably higher than that of single-phase flow

Inlet effects on vertical-downward air–water two-phase flow

Energy Technology Data Exchange (ETDEWEB)

Qiao, Shouxu; Mena, Daniel; Kim, Seungjin, E-mail: skim@psu.edu

2017-02-15

Highlights: • Inlet effects on two-phase flow parameters in vertical-downward flow are studied. • Flow regimes in the vertical-downward two-phase flow are defined. • Vertical-downward flow regime maps for three inlet configurations are developed. • Frictional pressure loss analysis for three different inlets is performed. • Database of local two-phase flow parameters for each inlet configuration. - Abstract: This paper focuses on investigating the geometric effects of inlets on global and local two-phase flow parameters in vertical-downward air–water two-phase flow. Flow visualization, frictional pressure loss analysis, and local experiments are performed in a test facility constructed from 50.8 mm inner diameter acrylic pipes. Three types of inlets of interest are studied: (1) two-phase flow injector without a flow straightener (Type A), (2) two-phase flow injector with a flow straightener (Type B), and (3) injection through a horizontal-to-vertical-downward 90° vertical elbow (Type C). A detailed flow visualization study is performed to characterize flow regimes including bubbly, slug, churn-turbulent, and annular flow. Flow regime maps for each inlet are developed and compared to identify the effects of each inlet. Frictional pressure loss analysis shows that the Lockhart–Martinelli method is capable of correlating the frictional loss data acquired for Type B and Type C inlets with a coefficient value of C = 25, but additional data may be needed to model the Type A inlet. Local two-phase flow parameters measured by a four-sensor conductivity probe in four bubbly and near bubbly flow conditions are analyzed. It is observed that vertical-downward two-phase flow has a characteristic center-peaked void profile as opposed to a wall-peaked profile as seen in vertical-upward flow. Furthermore, it is shown that the Type A inlet results in the most pronounced center-peaked void fraction profile, due to the coring phenomenon. Type B and Type C inlets

Low-diffusion rotated upwind schemes, multigrid and defect correction for steady, multi-dimensional Euler flows

NARCIS (Netherlands)

Koren, B.; Hackbusch, W.; Trottenberg, U.

1991-01-01

Two simple, multi-dimensional upwind discretizations for the steady Euler equations are derived, with the emphasis Iying on bath a good accuracy and a good solvability. The multi-dimensional upwinding consists of applying a one-dimensional Riemann solver with a locally rotated left and right state,

Numerical simulation for two-phase jet problem

International Nuclear Information System (INIS)

Lee, W.H.; Shah, V.L.

1981-01-01

A computer program TWOP was developed for obtaining the numerical solutions of three-dimensional, transient, two-phase flow system with nonequilibrium and nonhomogeneous conditions. TWOP employs two-fluid model and a set of the conservation equations formulated by Harlow and Amsden along with their Implicit Multi-Field (IMF) numerical technique that allows all degrees of couplings between the two fields. We have further extended the procedure of Harlow and Amsden by incorporating the implicit couplings of phase transition and interfacial heat transfer terms in the energy equations. Numerical results of two tested problems are presented to demonstrate the capabilities of the TWOP code. The first problem is the separation of vapor and liquid, showing that the code can handle the computational difficulties such as liquid packing and sharp interface phenomena. The second problem is the high pressure two-phase jet impinged on vertical plate, demonstrating the important role of the interfacial mass and momentum exchange

Modeling two-phase flow in three-dimensional complex flow-fields of proton exchange membrane fuel cells

Science.gov (United States)

Kim, Jinyong; Luo, Gang; Wang, Chao-Yang

2017-10-01

3D fine-mesh flow-fields recently developed by Toyota Mirai improved water management and mass transport in proton exchange membrane (PEM) fuel cell stacks, suggesting their potential value for robust and high-power PEM fuel cell stack performance. In such complex flow-fields, Forchheimer's inertial effect is dominant at high current density. In this work, a two-phase flow model of 3D complex flow-fields of PEMFCs is developed by accounting for Forchheimer's inertial effect, for the first time, to elucidate the underlying mechanism of liquid water behavior and mass transport inside 3D complex flow-fields and their adjacent gas diffusion layers (GDL). It is found that Forchheimer's inertial effect enhances liquid water removal from flow-fields and adds additional flow resistance around baffles, which improves interfacial liquid water and mass transport. As a result, substantial improvements in high current density cell performance and operational stability are expected in PEMFCs with 3D complex flow-fields, compared to PEMFCs with conventional flow-fields. Higher current density operation required to further reduce PEMFC stack cost per kW in the future will necessitate optimizing complex flow-field designs using the present model, in order to efficiently remove a large amount of product water and hence minimize the mass transport voltage loss.

Novel image reconstruction algorithm for multi-phase flow tomography system using γ ray method

International Nuclear Information System (INIS)

Hao Kuihong; Wang Huaxiang; Gao Mei

2007-01-01

After analyzing the reason of image reconstructed algorithm by using the conventional back projection (IBP) is prone to produce spurious line, and considering the characteristic of multi-phase flow tomography, a novel image reconstruction algorithm is proposed, which carries out the intersection calculation using back projection data. This algorithm can obtain a perfect system point spread function, and can eliminate spurious line better. Simulating results show that the algorithm is effective for identifying multi-phase flow pattern. (authors)

Numerical simulation of 3-D incompressible, multi-phase flows over cavitating projectiles

Energy Technology Data Exchange (ETDEWEB)

Owis, F.M.; Nayfeh, A.H. [Blacksburg State University, Dept. of Engineering Science and Mechanics, MC 0219, Virginia Polytechnic Institute, VA (United States)

2004-04-01

The hydrodynamic cavitation over axisymmetric projectiles is computed using the unsteady incompressible Navier-Stokes equations for multi-fluid elements. The governing equations are discretized on a structured grid using an upwind difference scheme with flux limits. A preconditioning dual-time stepping method is used for the unsteady computations. The Eigen-system is derived for the Jacobian matrices. This Eigen-system is suitable for high-density ratio multi-fluid flows and it provides high numerical stability and fast convergence. This method can be used to compute single- as well as multi-phase flows. Cavitating flows over projectiles with different geometries are computed and the results are in good agreement with available experimental data and other published computations. (authors)

Thermal hydraulics-I. 1. Phasic Discrimination in Two-Phase-Flow Measurements Using Particle Image Velocimetry

International Nuclear Information System (INIS)

Todd, D.R.; Ortiz-Villafuerte, J.; Schmidl, W.D.; Hassan, Y.A.; Sanchez-Silva, F.

2001-01-01

Information about the dispersed phase parameters -such as location, displacement, and interfacial area -are very important in the analysis of two-phase flows. Local flow disturbances in the continuous phase can be quite significant when the dispersed phase (i.e., a particle, drop, or bubble) passes through the medium. Application of point-wise measurement methods such as hot wire anemometry and laser anemometry suffer significant limitations in two-phase-flow measurements when these local disturbances are strong. Also, these two methods typically lack the ability to quantify the dispersed phase. Previous work has shown that meaningful analysis of the instantaneous continuous phase velocity field requires knowledge of the dispersed phase parameters, especially location and trajectory. Continuous phase parameters such as the local instantaneous vorticity and local turbulence fluctuations are influenced by the passage of the dispersed phase. Thus, development of two-phase-flow models (such as a bubble wake model) requires knowledge of the relative location of a local continuous phase parameter to the dispersed flow object (i.e., directly behind or off the side of the object). Also, conditional sampling must be performed using a meaningful parameter as the sampling point, i.e., the passage of a specific size of bubble. A system has been developed at Texas A and M University to quantify the dispersed phase parameters for two-phase bubbly flow in a vertical pipe with co-current upward flow. This system uses an orthogonal shadow particle image velocimetry (SPIV) technique, which instantaneously measures three-dimensional bubble locations, volumes, and interfacial areas -while measuring the three-dimensional bubble velocities and accelerations over a sequence of discrete measurements. The SPIV system is capable of analyzing flows with a large number of bubbles in close proximity. A set of sample images has been collected as part of the preliminary testing and development

A model for non-equilibrium, non-homogeneous two-phase critical flow

International Nuclear Information System (INIS)

Bassel, Wageeh Sidrak; Ting, Daniel Kao Sun

1999-01-01

Critical two phase flow is a very important phenomena in nuclear reactor technology for the analysis of loss of coolant accident. Several recent papers, Lee and Shrock (1990), Dagan (1993) and Downar (1996) , among others, treat the phenomena using complex models which require heuristic parameters such as relaxation constants or interfacial transfer models. In this paper a mathematical model for one dimensional non equilibrium and non homogeneous two phase flow in constant area duct is developed. The model is constituted of three conservation equations type mass ,momentum and energy. Two important variables are defined in the model: equilibrium constant in the energy equation and the impulse function in the momentum equation. In the energy equation, the enthalpy of the liquid phase is determined by a linear interpolation function between the liquid phase enthalpy at inlet condition and the saturated liquid enthalpy at local pressure. The interpolation coefficient is the equilibrium constant. The momentum equation is expressed in terms of the impulse function. It is considered that there is slip between the liquid and vapor phases, the liquid phase is in metastable state and the vapor phase is in saturated stable state. The model is not heuristic in nature and does not require complex interface transfer models. It is proved numerically that for the critical condition the partial derivative of two phase pressure drop with respect to the local pressure or to phase velocity must be zero.This criteria is demonstrated by numerical examples. The experimental work of Fauske (1962) and Jeandey (1982) were analyzed resulting in estimated numerical values for important parameters like slip ratio, equilibrium constant and two phase frictional drop. (author)

Slug flow transitions in horizontal gas/liquid two-phase flows. Dependence on channel height and system pressure for air/water and steam/water two-phase flows

International Nuclear Information System (INIS)

Nakamura, Hideo

1996-05-01

The slug flow transitions and related phenomena for horizontal two-phase flows were studied for a better prediction of two-phase flows that typically appear during the reactor loss-of-coolant accidents (LOCAs). For better representation of the flow conditions experimentally, two large-scaled facility: TPTF for high-pressure steam/water two-phase flows and large duct test facility for air/water two-phase flows, were used. The visual observation of the flow using a video-probe was performed in the TPTF experiments for good understanding of the phenomena. The currently-used models and correlations based mostly on the small-scale low-pressure experiments were reviewed and improved based on these experimental results. The modified Taitel-Dukler model for prediction of transition into slug flow from wavy flow and the modified Steen-Wallis correlation for prediction of onset of liquid entrainment from the interfacial waves were obtained. An empirical correlation for the gas-liquid interfacial friction factor was obtained further for prediction of liquid levels at wavy flow. The region of slug flow regime that is generally under influences of the channel height and system pressure was predicted well when these models and correlations were applied together. (author). 90 refs

Magnetic liquid metal two-phase flow research. Phase 1. Final report

International Nuclear Information System (INIS)

Graves, R.D.

1983-04-01

The Phase I research demonstrates the feasibility of the magnetic liquid metal (MLM) two-phase flow concept. A dispersion analysis is presented based on a complete set of two-phase-flow equations augmented to include stresses due to magnetic polarization of the fluid. The analysis shows that the stability of the MLM two-phase flow is determined by the magnetic Mach number, the slip ratio, geometry of the flow relative to the applied magnetic field, and by the voidage dependence of the interfacial forces. Results of a set of experiments concerned with magnetic effects on the dynamics of single bubble motion in an aqueous-based, viscous, conducting magnetic fluid are presented. Predictions in the theoretical literature are qualitatively verified using a bench-top experimental apparatus. In particular, applied magnetic fields are seen to lead to reduced bubble size at fixed generating orifice pressure

One- and two-dimensional Stirling machine simulation using experimentally generated reversing flow turbuulence models

International Nuclear Information System (INIS)

Goldberg, L.F.

1990-08-01

The activities described in this report do not constitute a continuum but rather a series of linked smaller investigations in the general area of one- and two-dimensional Stirling machine simulation. The initial impetus for these investigations was the development and construction of the Mechanical Engineering Test Rig (METR) under a grant awarded by NASA to Dr. Terry Simon at the Department of Mechanical Engineering, University of Minnesota. The purpose of the METR is to provide experimental data on oscillating turbulent flows in Stirling machine working fluid flow path components (heater, cooler, regenerator, etc.) with particular emphasis on laminar/turbulent flow transitions. Hence, the initial goals for the grant awarded by NASA were, broadly, to provide computer simulation backup for the design of the METR and to analyze the results produced. This was envisaged in two phases: First, to apply an existing one-dimensional Stirling machine simulation code to the METR and second, to adapt a two-dimensional fluid mechanics code which had been developed for simulating high Rayleigh number buoyant cavity flows to the METR. The key aspect of this latter component was the development of an appropriate turbulence model suitable for generalized application to Stirling simulation. A final-step was then to apply the two-dimensional code to an existing Stirling machine for which adequate experimental data exist. The work described herein was carried out over a period of three years on a part-time basis. Forty percent of the first year's funding was provided as a match to the NASA funds by the Underground Space Center, University of Minnesota, which also made its computing facilities available to the project at no charge

Dynamic Modeling Strategy for Flow Regime Transition in Gas-Liquid Two-Phase Flows

Directory of Open Access Journals (Sweden)

Xia Wang

2012-12-01

Full Text Available In modeling gas-liquid two-phase flows, the concept of flow regimes has been widely used to characterize the global interfacial structure of the flows. Nearly all constitutive relations that provide closures to the interfacial transfers in two-phase flow models, such as the two-fluid model, are flow regime dependent. Current nuclear reactor safety analysis codes, such as RELAP5, classify flow regimes using flow regime maps or transition criteria that were developed for steady-state, fully-developed flows. As two-phase flows are dynamic in nature, it is important to model the flow regime transitions dynamically to more accurately predict the two-phase flows. The present work aims to develop a dynamic modeling strategy to determine flow regimes in gas-liquid two-phase flows through introduction of interfacial area transport equations (IATEs within the framework of a two-fluid model. The IATE is a transport equation that models the interfacial area concentration by considering the creation of the interfacial area, fluid particle (bubble or liquid droplet disintegration, boiling and evaporation, and the destruction of the interfacial area, fluid particle coalescence and condensation. For flow regimes beyond bubbly flows, a two-group IATE has been proposed, in which bubbles are divided into two groups based on their size and shapes, namely group-1 and group-2 bubbles. A preliminary approach to dynamically identify the flow regimes is discussed, in which discriminators are based on the predicted information, such as the void fraction and interfacial area concentration. The flow regime predicted with this method shows good agreement with the experimental observations.

Experimental on two sensors combination used in horizontal pipe gas-water two-phase flow

International Nuclear Information System (INIS)

Wu, Hao; Dong, Feng

2014-01-01

Gas-water two phase flow phenomenon widely exists in production and living and the measurement of it is meaningful. A new type of long-waist cone flow sensor has been designed to measure two-phase mass flow rate. Six rings structure of conductance probe is used to measure volume fraction and axial velocity. The calibration of them have been made. Two sensors have been combined in horizontal pipeline experiment to measure two-phase flow mass flow rate. Several model of gas-water two-phase flow has been discussed. The calculation errors of total mass flow rate measurement is less than 5% based on the revised homogeneous flow model

Prediction of two-phase choked-flow through safety valves

International Nuclear Information System (INIS)

Arnulfo, G; Bertani, C; De Salve, M

2014-01-01

Different models of two-phase choked flow through safety valves are applied in order to evaluate their capabilities of prediction in different thermal-hydraulic conditions. Experimental data available in the literature for two-phase fluid and subcooled liquid upstream the safety valve have been compared with the models predictions. Both flashing flows and non-flashing flows of liquid and incondensable gases have been considered. The present paper shows that for flashing flows good predictions are obtained by using the two-phase valve discharge coefficient defined by Lenzing and multiplying it by the critical flow rate in an ideal nozzle evaluated by either Omega Method or the Homogeneous Non-equilibrium Direct Integration. In case of non-flashing flows of water and air, Leung/Darby formulation of the two-phase valve discharge coefficient together with the Omega Method is more suitable to the prediction of flow rate.

TWO-DIMENSIONAL CORE-COLLAPSE SUPERNOVA MODELS WITH MULTI-DIMENSIONAL TRANSPORT

International Nuclear Information System (INIS)

Dolence, Joshua C.; Burrows, Adam; Zhang, Weiqun

2015-01-01

We present new two-dimensional (2D) axisymmetric neutrino radiation/hydrodynamic models of core-collapse supernova (CCSN) cores. We use the CASTRO code, which incorporates truly multi-dimensional, multi-group, flux-limited diffusion (MGFLD) neutrino transport, including all relevant O(v/c) terms. Our main motivation for carrying out this study is to compare with recent 2D models produced by other groups who have obtained explosions for some progenitor stars and with recent 2D VULCAN results that did not incorporate O(v/c) terms. We follow the evolution of 12, 15, 20, and 25 solar-mass progenitors to approximately 600 ms after bounce and do not obtain an explosion in any of these models. Though the reason for the qualitative disagreement among the groups engaged in CCSN modeling remains unclear, we speculate that the simplifying ''ray-by-ray'' approach employed by all other groups may be compromising their results. We show that ''ray-by-ray'' calculations greatly exaggerate the angular and temporal variations of the neutrino fluxes, which we argue are better captured by our multi-dimensional MGFLD approach. On the other hand, our 2D models also make approximations, making it difficult to draw definitive conclusions concerning the root of the differences between groups. We discuss some of the diagnostics often employed in the analyses of CCSN simulations and highlight the intimate relationship between the various explosion conditions that have been proposed. Finally, we explore the ingredients that may be missing in current calculations that may be important in reproducing the properties of the average CCSNe, should the delayed neutrino-heating mechanism be the correct mechanism of explosion

Experimental and numerical investigation on two-phase flow instabilities

Energy Technology Data Exchange (ETDEWEB)

Ruspini, Leonardo Carlos

2013-03-01

Two-phase flow instabilities are experimentally and numerically studied within this thesis. In particular, the phenomena called Ledinegg instability, density wave oscillations and pressure drop oscillations are investigated. The most important investigations regarding the occurrence of two-phase flow instabilities are reviewed. An extensive description of the main contributions in the experimental and analytical research is presented. In addition, a critical discussion and recommendations for future investigations are presented. A numerical framework using a hp-adaptive method is developed in order to solve the conservation equations modelling general thermo-hydraulic systems. A natural convection problem is analysed numerically in order to test the numerical solver. Moreover, the description of an adaptive strategy to solve thermo-hydraulic problems is presented. In the second part of this dissertation, a homogeneous model is used to study Ledinegg, density wave and pressure drop oscillations phenomena numerically. The dynamic characteristics of the Ledinegg (flow excursion) phenomenon are analysed through the simulation of several transient examples. In addition, density wave instabilities in boiling and condensing systems are investigated. The effects of several parameters, such as the fluid inertia and compressibility volumes, on the stability limits of Ledinegg and density wave instabilities are studied, showing a strong influence of these parameters. Moreover, the phenomenon called pressure drop oscillations is numerically investigated. A discussion of the physical representation of several models is presented with reference to the obtained numerical results. Finally, the influence of different parameters on these phenomena is analysed. In the last part, an experimental investigation of these phenomena is presented. The designing methodology used for the construction of the experimental facility is described. Several simulations and a non-dimensional

Visual Analysis of Inclusion Dynamics in Two-Phase Flow.

Science.gov (United States)

Karch, Grzegorz Karol; Beck, Fabian; Ertl, Moritz; Meister, Christian; Schulte, Kathrin; Weigand, Bernhard; Ertl, Thomas; Sadlo, Filip

2018-05-01

In single-phase flow visualization, research focuses on the analysis of vector field properties. In two-phase flow, in contrast, analysis of the phase components is typically of major interest. So far, visualization research of two-phase flow concentrated on proper interface reconstruction and the analysis thereof. In this paper, we present a novel visualization technique that enables the investigation of complex two-phase flow phenomena with respect to the physics of breakup and coalescence of inclusions. On the one hand, we adapt dimensionless quantities for a localized analysis of phase instability and breakup, and provide detailed inspection of breakup dynamics with emphasis on oscillation and its interplay with rotational motion. On the other hand, we present a parametric tightly linked space-time visualization approach for an effective interactive representation of the overall dynamics. We demonstrate the utility of our approach using several two-phase CFD datasets.

Predicting transition in two- and three-dimensional separated flows

International Nuclear Information System (INIS)

Cutrone, L.; De Palma, P.; Pascazio, G.; Napolitano, M.

2008-01-01

This paper is concerned with the numerical prediction of two- and three-dimensional transitional separated flows of turbomachinery interest. The recently proposed single-point transition model based on the use of a laminar kinetic energy transport equation is considered, insofar as it does not require to evaluate any integral parameter, such as boundary-layer thickness, and is thus directly applicable to three-dimensional flows. A well established model, combining a transition-onset correlation with an intermittency transport equation, is also used for comparison. Both models are implemented within a Reynolds-averaged Navier-Stokes solver employing a low-Reynolds-number k-ω turbulence model. The performance of the transition models have been evaluated and tested versus well-documented incompressible flows past a flat plate with semi-circular leading edge, namely: tests T3L2, T3L3, T3L5, and T3LA1 of ERCOFTAC, with different Reynolds numbers and free-stream conditions, the last one being characterized by a non-zero pressure gradient. In all computations, the first model has proven as adequate as or superior to the second one and has been then applied with success to two more complex test cases, for which detailed experimental data are available in the literature, namely: the two- and three-dimensional flows through the T106 linear turbine cascade

A Semi-implicit Numerical Scheme for a Two-dimensional, Three-field Thermo-Hydraulic Modeling

International Nuclear Information System (INIS)

Hwang, Moonkyu; Jeong, Jaejoon

2007-07-01

The behavior of two-phase flow is modeled, depending on the purpose, by either homogeneous model, drift flux model, or separated flow model, Among these model, in the separated flow model, the behavior of each flow phase is modeled by its own governing equation, together with the interphase models which describe the thermal and mechanical interactions between the phases involved. In this study, a semi-implicit numerical scheme for two-dimensional, transient, two-fluid, three-field is derived. The work is an extension to the previous study for the staggered, semi-implicit numerical scheme in one-dimensional geometry (KAERI/TR-3239/2006). The two-dimensional extension is performed by specifying a relevant governing equation set and applying the related finite differencing method. The procedure for employing the semi-implicit scheme is also described in detail. Verifications are performed for a 2-dimensional vertical plate for a single-phase and two-phase flows. The calculations verify the mass and energy conservations. The symmetric flow behavior, for the verification problem, also confirms the momentum conservation of the numerical scheme

Modelling aspects of two phase flow

International Nuclear Information System (INIS)

Mayinger, F.

1977-01-01

In two phase flow scaling is much more limited to very narrowly defined physical phenomena than in single phase fluids. For complex and combined phenomena it can be achieved not by using dimensionless numbers alone but in addition a detailed mathematical description of the physical problem - usually in the form of a computer program - must be available. An important role plays the scaling of the thermodynamic data of the modelling fluid. From a literature survey and from own scaling experiments the conclusion can be drawn that Freon is a quite suitable modelling fluid for scaling steam-water mixtures. However, whithout a theoretical description of the phenomena nondimensional numbers for scaling two phase flow must be handled very carefully. (orig.) [de

Development of a two-phase SPH model for sediment laden flows

Science.gov (United States)

Shi, Huabin; Yu, Xiping; Dalrymple, Robert A.

2017-12-01

A SPH model based on a general formulation for solid-fluid two-phase flows is proposed for suspended sediment motion in free surface flows. The water and the sediment are treated as two miscible fluids, and the multi-fluid system is discretized by a single set of SPH particles, which move with the water velocity and carry properties of the two phases. Large eddy simulation (LES) is introduced to deal with the turbulence effect, and the widely used Smagorinsky model is modified to take into account the influence of sediment particles on the turbulence. The drag force is accurately formulated by including the hindered settling effect. In the model, the water is assumed to be weakly compressible while the sediment is incompressible, and a new equation of state is proposed for the pressure in the sediment-water mixture. Dynamic boundary condition is employed to treat wall boundaries, and a new strategy of Shepard filtering is adopted to damp the pressure oscillation. The developed two-phase SPH model is validated by comparing the numerical results with analytical solutions for idealized cases of still water containing both neutrally buoyant and naturally settling sand and for plane Poiseuille flows carrying neutrally buoyant particles, and is then applied to sand dumping from a line source into a water tank, where the sand cloud settles with a response of the free water surface. It is shown that the numerical results are in good agreement with the experimental data as well as the empirical formulas. The characteristics of the settling sand cloud, the pressure field, and the flow vortices are studied. The motion of the free water surface is also discussed. The proposed two-phase SPH model is proven to be effective for numerical simulation of sand dumping into waters.

Interfacial structures in downward two-phase bubbly flow

International Nuclear Information System (INIS)

Paranjape, S.S.; Kim, S.; Ishii, M.; Kelly, J.

2003-01-01

Downward two-phase flow was studied considering its significance in view of Light Water Reactor Accidents (LWR) such as Loss of Heat Sink (LOHS) by feed water loss or secondary pipe break. The flow studied, was an adiabatic, air-water, co-current, vertically downward two-phase flow. The experimental test sections had internal hydraulic diameters of 25.4 mm and 50.8 mm. Flow regime map was obtained using the characteristic signals obtained from an impedance void meter, employing neural network based identification methodology to minimize the subjective judgment in determining the flow regimes. A four sensor conductivity probe was used to measure the local two phase flow parameters, which characterize the interfacial structures. The local time averaged two-phase flow parameters measured were: void fraction (α), interfacial area concentration (a i ), bubble velocity (v g ), and Sauter mean diameter (D Sm ). The flow conditions were from the bubbly flow regime. The local profiles of these parameters as well as their axial development revealed the nature of the interfacial structures and the bubble interaction mechanisms occurring in the flow. Furthermore, this study provided a good database for the development of the interfacial area transport equation, which dynamically models the changes in the interfacial area along the flow field. An interfacial area transport equation was developed for downward flow based on that developed for the upward flow, with certain modifications in the bubble interaction terms. The area averaged values of the interfacial area concentration were compared with those predicted by the interfacial area transport model. (author)

Mechanics of occurrence of critical flow in compressible two-phase flow

International Nuclear Information System (INIS)

Katto, Yoshiro; Sudo, Yukio

1976-01-01

Fundamental framework of mechanics for the occurrence of critical flow is investigated, following the principle that the critical flow appears as a limit in a continuous change of state of flow along a nozzle (or a pipe) and should be derived only from simultaneous mechanical equations concerned with the flow. Mathematical procedures with which the critical flow: (i) the single phase flow of an arbitrary fluid, unrestricted by the equation of state of ideal gas, where the number of simultaneous equations is equal to the number of independent variables, and (ii) the one-component, separated two-phase flow under saturated condition, where the number of equations exceeds that of variables. In each case, interesting mechanism of leading to the occurrence of a limiting state of flow at a definite cross-section in a nozzle (incl. a pipe) is clarified, and a definite state of flow at the critical cross-section is also determined. Then, the analysis is extended to the critical flow which should appear in the completely isolated and the homogeneously dispersed, two-component, two-phase flow (composed of a compressible and an incompressible substance). It is found that the analyses of these special flow patterns provide several supplementary information to the mechanics of critical flow. (auth.)

Violent flows in aqueous foams III: physical multi-phase model comparison with aqueous foam shock tube experiments

Science.gov (United States)

Redford, J. A.; Ghidaglia, J.-M.; Faure, S.

2018-06-01

Mitigation of blast waves in aqueous foams is a problem that has a strong dependence on multi-phase effects. Here, a simplified model is developed from the previous articles treating violent flows (D'Alesio et al. in Eur J Mech B Fluids 54:105-124, 2015; Faure and Ghidaglia in Eur J Mech B Fluids 30:341-359, 2011) to capture the essential phenomena. The key is to have two fluids with separate velocities to represent the liquid and gas phases. This allows for the interaction between the two phases, which may include terms for drag, heat transfer, mass transfer due to phase change, added mass effects, to be included explicitly in the model. A good test for the proposed model is provided by two experimental data sets that use a specially designed shock tube. The first experiment has a test section filled with spray droplets, and the second has a range of aqueous foams in the test section. A substantial attenuation of the shock wave is seen in both cases, but a large difference is observed in the sound speeds. The droplets cause no observable change from the air sound speed, while the foams have a reduced sound speed of approximately 50-75 m/s . In the model given here, an added mass term is introduced in the governing equations to capture the low sound speed. The match between simulation and experiment is found to be satisfactory for both droplets and the foam. This is especially good when considering the complexity of the physics and the effects that are unaccounted for, such as three-dimensionality and droplet atomisation. The resulting statistics illuminate the processes occurring in such flows.

Multiscale time-splitting strategy for multiscale multiphysics processes of two-phase flow in fractured media

KAUST Repository

Sun, S.; Kou, J.; Yu, B.

2011-01-01

The temporal discretization scheme is one important ingredient of efficient simulator for two-phase flow in the fractured porous media. The application of single-scale temporal scheme is restricted by the rapid changes of the pressure and saturation in the fractured system with capillarity. In this paper, we propose a multi-scale time splitting strategy to simulate multi-scale multi-physics processes of two-phase flow in fractured porous media. We use the multi-scale time schemes for both the pressure and saturation equations; that is, a large time-step size is employed for the matrix domain, along with a small time-step size being applied in the fractures. The total time interval is partitioned into four temporal levels: the first level is used for the pressure in the entire domain, the second level matching rapid changes of the pressure in the fractures, the third level treating the response gap between the pressure and the saturation, and the fourth level applied for the saturation in the fractures. This method can reduce the computational cost arisen from the implicit solution of the pressure equation. Numerical examples are provided to demonstrate the efficiency of the proposed method.

A numerical model of two-phase flow at the micro-scale using the volume-of-fluid method

Science.gov (United States)

Shams, Mosayeb; Raeini, Ali Q.; Blunt, Martin J.; Bijeljic, Branko

2018-03-01

This study presents a simple and robust numerical scheme to model two-phase flow in porous media where capillary forces dominate over viscous effects. The volume-of-fluid method is employed to capture the fluid-fluid interface whose dynamics is explicitly described based on a finite volume discretization of the Navier-Stokes equations. Interfacial forces are calculated directly on reconstructed interface elements such that the total curvature is preserved. The computed interfacial forces are explicitly added to the Navier-Stokes equations using a sharp formulation which effectively eliminates spurious currents. The stability and accuracy of the implemented scheme is validated on several two- and three-dimensional test cases, which indicate the capability of the method to model two-phase flow processes at the micro-scale. In particular we show how the co-current flow of two viscous fluids leads to greatly enhanced flow conductance for the wetting phase in corners of the pore space, compared to a case where the non-wetting phase is an inviscid gas.

The effect of the virtual mass force term on the stability of transient two-phase flow analysis

International Nuclear Information System (INIS)

Watanabe, Tadashi; Hirano, Masashi; Tanabe, Fumiya

1989-08-01

The effect of the virtual mass force term on the stability of transient two-phase flow analysis is studied. The objective form of the virtual mass acceleration is used. The virtual mass coefficient is determined from the stability condition of basic equations against infinitesimal high wave-number perturbations. The parameter is chosen so that a reasonable agreement between the analytical and experimental sound speed in two-phase flows can be obtained. A one-dimensional sedimentation problem is simulated by the MINCS code which is a tool for transient two-phase flow analysis. The stability analysis is performed for the numerical procedure. It is shown that calculated results are stabilized so long as the virtual mass coefficient satisfies the stability condition of differential equations. (author)

Dynamics of lava flow - Thickness growth characteristics of steady two-dimensional flow

Science.gov (United States)

Park, S.; Iversen, J. D.

1984-01-01

The thickness growth characteristics of flowing lava are investigated using a heat balance model and a two-dimensional model for flow of a Bingham plastic fluid down an inclined plane. It is found that yield strength plays a crucial role in the thickening of a lava flow of given flow rate. To illustrate this point, downstream thickness profiles and yield strength distributions were calculated for flows with mass flow rates of 10,000 and 100,000 kg/m-sec. Higher flow rates led to slow cooling rates which resulted in slow rate of increase of yield strength and thus greater flow lengths.

Regimes of Two-Phase Flow in Short Rectangular Channel

Science.gov (United States)

Chinnov, Evgeny A.; Guzanov, Vladimir V.; Cheverda, Vyacheslav; Markovich, Dmitry M.; Kabov, Oleg A.

2009-08-01

Experimental study of two-phase flow in the short rectangular horizontal channel with height 440 μm has been performed. Characteristics of liquid motion inside the channel have been registered and measured by the Laser Induced Fluorescence technique. New information has allowed determining more precisely the characteristics of churn regime and boundaries between different regimes of two-phase flow. It was shown that formation of some two-phase flow regimes and transitions between them are determined by instability of the flow in the lateral parts of the channel.

Forced convection flow boiling and two-phase flow phenomena in a microchannel

Science.gov (United States)

Na, Yun Whan

2008-07-01

) technique. The effects of different constant heat fluxes and different channel heights on the boiling mechanisms were investigated. The effects of liquid velocity on the bubble departure diameter were analyzed. The obtained results showed that the wall superheats at the position of nucleate boiling are relatively independent of the mass flow rates at the same channel height. The obtained results, however, showed that the heat flux at the onset of nucleate boiling strongly depends on the channel height. With a decrease of the channel height and an increase of the liquid velocity at the channel inlet, the departure diameter of a bubble was smaller. The periodic flow patterns, such as the bubbly flow, elongated slug flow, and churn flow were observed in the microchannel. Flow instabilities of two-phase flow boiling in a trapezoidal microchannel using a three-dimensional model were investigated. Fluctuation behaviors of flow boiling parameters such as wall temperature and inlet pressure caused by periodic flow patterns were studied at different heat fluxes and mass fluxes. The numerical results showed large amplitude and short period oscillations for wall temperature and inlet pressure fluctuations. Stable and unstable flow boiling regime with short period oscillations were investigated. Those flow boiling regimes were not listed in stable and unstable boiling regime map proposed by Wang et al. (2007).

Unsteady State Two Phase Flow Pressure Drop Calculations

OpenAIRE

Ayatollahi, Shahaboddin

1992-01-01

A method is presented to calculate unsteady state two phase flow in a gas-liquid line based on a quasi-steady state approach. A computer program for numerical solution of this method was prepared. Results of calculations using the computer program are presented for several unsteady state two phase flow systems

Two-Phase Annular Flow in Helical Coil Flow Channels in a Reduced Gravity Environment

Science.gov (United States)

Keshock, Edward G.; Lin, Chin S.

1996-01-01

A brief review of both single- and two-phase flow studies in curved and coiled flow geometries is first presented. Some of the complexities of two-phase liquid-vapor flow in curved and coiled geometries are discussed, and serve as an introduction to the advantages of observing such flows under a low-gravity environment. The studies proposed -- annular two-phase air-water flow in helical coil flow channels are described. Objectives of the studies are summarized.

Experimental study of single- and two-phase flow fields around PWR steam generator tube support plates

International Nuclear Information System (INIS)

Bates, J.M.; Stewart, C.W.

1979-08-01

Laser-Doppler anemometry (LDA) was used to measure local mean axial velocities and turbulence intnsities at selected locations within a study model dimensionally protypic of an existing PWR steam generator design. The model tube bundle with support plate was installed in a special flow housing that formed part of an isothermal recirculating water flow loop. Flow conditions for this experiment were intended to simulate only typical single-phase flow velocities and were not an attempt to completely model actual steam generator, boiling, two-phase flow conditions. The measurements were performed in water at approximately 85 0 F with test section average velocities of approximately 0.55 and 1.1 fps. These conditions corresponded to Reynolds numbers of approximately 7,000 and approximately 14,000, respectively. Normalized velocity and turbulence intensity ratios are graphically reported. Additional qualitative, photographic investigations of air-water two-phase flows in a PWR steam generator study model were also performed

Qualitative behaviour of incompressible two-phase flows with phase ...

Indian Academy of Sciences (India)

Jan Prüss

2017-11-07

Nov 7, 2017 ... Qualitative behaviour of incompressible two-phase flows with phase ... Germany. 2Graduate School of Human and Environmental Studies, Kyoto University, ... Note that j is a dummy variable as it can be eliminated from the ...

Characterization of horizontal air–water two-phase flow

Energy Technology Data Exchange (ETDEWEB)

Kong, Ran; Kim, Seungjin, E-mail: skim@psu.edu

2017-02-15

Highlights: • A visualization study is performed to develop flow regime map in horizontal flow. • Database in horizontal bubbly flow is extended using a local conductivity probe. • Frictional pressure drop analysis is performed in horizontal bubbly flow. • Drift flux analysis is performed in horizontal bubbly flow. - Abstract: This paper presents experimental studies performed to characterize horizontal air–water two-phase flow in a round pipe with an inner diameter of 3.81 cm. A detailed flow visualization study is performed using a high-speed video camera in a wide range of two-phase flow conditions to verify previous flow regime maps. Two-phase flows are classified into bubbly, plug, slug, stratified, stratified-wavy, and annular flow regimes. While the transition boundaries identified in the present study compare well with the existing ones (Mandhane et al., 1974) in general, some discrepancies are observed for bubbly-to-plug/slug, and plug-to-slug transition boundaries. Based on the new transition boundaries, three additional test conditions are determined in horizontal bubbly flow to extend the database by Talley et al. (2015a). Various local two-phase flow parameters including void fraction, interfacial area concentration, bubble velocity, and bubble Sauter mean diameter are obtained. The effects of increasing gas flow rate on void fraction, bubble Sauter mean diameter, and bubble velocity are discussed. Bubbles begin to coalesce near the gas–liquid layer instead of in the highly packed region when gas flow rate increases. Using all the current experimental data, two-phase frictional pressure loss analysis is performed using the Lockhart–Martinelli method. It is found that the coefficient C = 24 yields the best agreement with the data with the minimum average difference. Moreover, drift flux analysis is performed to predict void-weighted area-averaged bubble velocity and area-averaged void fraction. Based on the current database, functional

Multi-perspective views of students’ difficulties with one-dimensional vector and two-dimensional vector

Science.gov (United States)

Fauzi, Ahmad; Ratna Kawuri, Kunthi; Pratiwi, Retno

2017-01-01

Researchers of students’ conceptual change usually collects data from written tests and interviews. Moreover, reports of conceptual change often simply refer to changes in concepts, such as on a test, without any identification of the learning processes that have taken place. Research has shown that students have difficulties with vectors in university introductory physics courses and high school physics courses. In this study, we intended to explore students’ understanding of one-dimensional and two-dimensional vector in multi perspective views. In this research, we explore students’ understanding through test perspective and interviews perspective. Our research study adopted the mixed-methodology design. The participants of this research were sixty students of third semester of physics education department. The data of this research were collected by testand interviews. In this study, we divided the students’ understanding of one-dimensional vector and two-dimensional vector in two categories, namely vector skills of the addition of one-dimensionaland two-dimensional vector and the relation between vector skills and conceptual understanding. From the investigation, only 44% of students provided correct answer for vector skills of the addition of one-dimensional and two-dimensional vector and only 27% students provided correct answer for the relation between vector skills and conceptual understanding.

Film boiling from spheres in single- and two-phase flow

International Nuclear Information System (INIS)

Liu, C.; Theofanous, T.G.; Yuen, W.W.

1992-01-01

Experimental data on film boiling heat transfer from single, inductively heated, spheres in single- and two-phase flow (saturated water and steam, respectively) are presented. In the single-phase-flow experiments water velocities ranged from 0.1 to 2.0 m/s; in the two-phase-flow experiments superficial water and steam velocities covered 0.1 to 0.6 m/s and 4 to 10 m/s, respectively. All experiments were run at atmospheric pressure and with sphere temperatures from 900C down to quenching. Limited interpretations of the single-phase- flow data are possible, but the two-phase-flow data are new and unique

Two-phase flow patterns in horizontal rectangular minichannel

Directory of Open Access Journals (Sweden)

Ron’shin Fedor

2016-01-01

Full Text Available The two-phase flow in a short horizontal channel of rectangular cross-section of 1 × 19 mm2 has been studied experimentally. Five conventional two-phase flow patterns have been detected (bubble, churn, stratified, annular and jet and transitions between them have been determined. It is shown that a change in the width of the horizontal channels has a substantial effect on the boundaries between the flow regimes.

Turbulence in two-phase flows

International Nuclear Information System (INIS)

Sullivan, J.P.; Houze, R.N.; Buenger, D.E.; Theofanous, T.G.

1981-01-01

Hot film Anemometry and Laser Doppler Velocimetry have been employed in this work to study the turbulence characteristics of Bubbly and Stratified two-phase flows, respectively. Extensive consistency checks were made to establish the reliability and hence the utility of these experimental techniques for the measurement of turbulence in two-phase flows. Buoyancy-driven turbulence in vertical bubbly flows has been identified experimentally and correlated in terms of a shear velocity superposition approach. This approach provides a criterion for the demarcation of the buoyancy-driven turbulence region from the wall shear-generated turbulence region. Our data confirm the roughly isotropic behavior expected for buoyancy-driven turbulence. Upgrading of our experimental system will permit investigations of the wall-shear dominated regime (i.e., isotropy, superposition approach, etc.). The stratified flow data demonstrate clearly that the maximum in the mean velocity profile does not coincide with the zero shear plane, indicating the existence of a negative eddy viscosity region. Previous studies do not take into account this difference and thus they yield incorrect friction factor data in addition to certain puzzling behavior in the upper wall region. The conditioned turbulence data in the wavy region indicate interesting trends and that an appropriate normalization of intensities must take into account the shear velocity at the interfacial (wavy) region

Experimental Study about Two-phase Damping Ratio on a Tube Bundle Subjected to Homogeneous Two-phase Flow

Energy Technology Data Exchange (ETDEWEB)

Sim, Woo Gun; Dagdan, Banzragch [Hannam Univ., Daejeon (Korea, Republic of)

2017-03-15

Two-phase cross flow exists in many shell-and-tube heat exchangers such as condensers, evaporators, and nuclear steam generators. The drag force acting on a tube bundle subjected to air/water flow is evaluated experimentally. The cylinders subjected to two-phase flow are arranged in a normal square array. The ratio of pitch to diameter is 1.35, and the diameter of the cylinder is 18 mm. The drag force along the flow direction on the tube bundles is measured to calculate the drag coefficient and the two-phase damping ratio. The two-phase damping ratios, given by the analytical model for a homogeneous two-phase flow, are compared with experimental results. The correlation factor between the frictional pressure drop and the hydraulic drag coefficient is determined from the experimental results. The factor is used to calculate the drag force analytically. It is found that with an increase in the mass flux, the drag force, and the drag coefficients are close to the results given by the homogeneous model. The result shows that the damping ratio can be calculated using the homogeneous model for bubbly flow of sufficiently large mass flux.

DSMC simulation of two-phase plume flow with UV radiation

Energy Technology Data Exchange (ETDEWEB)

Li, Jie; Liu, Ying; Wang, Ning; Jin, Ling [College of Aerospace Science and Engineering, National University of Defense Technology, Changsha, Hunan, 410073 (China)

2014-12-09

Rarefied gas-particle two-phase plume in which the phase of particles is liquid or solid flows from a solid propellant rocket of hypersonic vehicle flying at high altitudes, the aluminum oxide particulates not only impact the rarefied gas flow properties, but also make a great difference to plume radiation signature, so the radiation prediction of the rarefied gas-particle two-phase plume flow is very important for space target detection of hypersonic vehicles. Accordingly, this project aims to study the rarefied gas-particle two-phase flow and ultraviolet radiation (UV) characteristics. Considering a two-way interphase coupling of momentum and energy, the direct simulation Monte Carlo (DSMC) method is developed for particle phase change and the particle flow, including particulate collision, coalescence as well as separation, and a Monte Carlo ray trace model is implemented for the particulate UV radiation. A program for the numerical simulation of the gas-particle two-phase flow and radiation in which the gas flow nonequilibrium is strong is implemented as well. Ultraviolet radiation characteristics of the particle phase is studied based on the calculation of the flow field coupled with the radiation calculation, the radiation model for different size particles is analyzed, focusing on the effects of particle emission, absorption, scattering as well as the searchlight emission of the nozzle. A new approach may be proposed to describe the rarefied gas-particle two-phase plume flow and radiation transfer characteristics in this project.

DSMC simulation of two-phase plume flow with UV radiation

Science.gov (United States)

Li, Jie; Liu, Ying; Wang, Ning; Jin, Ling

2014-12-01

Rarefied gas-particle two-phase plume in which the phase of particles is liquid or solid flows from a solid propellant rocket of hypersonic vehicle flying at high altitudes, the aluminum oxide particulates not only impact the rarefied gas flow properties, but also make a great difference to plume radiation signature, so the radiation prediction of the rarefied gas-particle two-phase plume flow is very important for space target detection of hypersonic vehicles. Accordingly, this project aims to study the rarefied gas-particle two-phase flow and ultraviolet radiation (UV) characteristics. Considering a two-way interphase coupling of momentum and energy, the direct simulation Monte Carlo (DSMC) method is developed for particle phase change and the particle flow, including particulate collision, coalescence as well as separation, and a Monte Carlo ray trace model is implemented for the particulate UV radiation. A program for the numerical simulation of the gas-particle two-phase flow and radiation in which the gas flow nonequilibrium is strong is implemented as well. Ultraviolet radiation characteristics of the particle phase is studied based on the calculation of the flow field coupled with the radiation calculation, the radiation model for different size particles is analyzed, focusing on the effects of particle emission, absorption, scattering as well as the searchlight emission of the nozzle. A new approach may be proposed to describe the rarefied gas-particle two-phase plume flow and radiation transfer characteristics in this project.

Two-Phase Quality/Flow Meter

Science.gov (United States)

Moerk, J. Steven (Inventor); Youngquist, Robert C. (Inventor); Werlink, Rudy J. (Inventor)

1999-01-01

A quality and/or flow meter employs a capacitance probe assembly for measuring the dielectric constant of flow stream, particularly a two-phase flow stream including liquid and gas components.ne dielectric constant of the flow stream varies depending upon the volume ratios of its liquid and gas components, and capacitance measurements can therefore be employed to calculate the quality of the flow, which is defined as the volume ratio of liquid in the flow to the total volume ratio of gas and liquid in the flow. By using two spaced capacitance sensors, and cross-correlating the time varying capacitance values of each, the velocity of the flow stream can also be determined. A microcontroller-based processing circuit is employed to measure the capacitance of the probe sensors.The circuit employs high speed timer and counter circuits to provide a high resolution measurement of the time interval required to charge each capacitor in the probe assembly. In this manner, a high resolution, noise resistant, digital representation of each of capacitance value is obtained without the need for a high resolution A/D converter, or a high frequency oscillator circuit. One embodiment of the probe assembly employs a capacitor with two ground plates which provide symmetry to insure that accurate measurements are made thereby.

Two-phase flow induced parametric vibrations in structural systems

International Nuclear Information System (INIS)

Hara, Fumio

1980-01-01

This paper is divided into two parts concerning piping systems and a nuclear fuel pin system. The significant experimental results concerning the random vibration induced in an L-shaped pipe by air-water two-phase flow and the theoretical analysis of the vibration are described in the first part. It was clarified for the first time that the parametric excitation due to the periodic changes of system mass, centrifugal force and Coriolis force was the mechanism of exciting the vibration. Moreover, the experimental and theoretical analyses of the mechanism of exciting vibration by air-water two-phase flow in a straight, horizontal pipe were carried out, and the first natural frequency of the piping system was strongly related to the dominant frequency of void signals. The experimental results on the vibration of a nuclear fuel pin model in parallel air-water two-phase flow are reported in the latter part. The relations between vibrational strain variance and two-phase flow velocity or pressure fluctuation, and the frequency characteristics of vibrational strain variance were obtained. The theoretical analysis of the dynamic interaction between air-water two-phase flow and a fuel pin structure, and the vibrational instability of fuel pins in alternate air and water slugs or in large bubble flow are also reported. (Kako, I.)

Non-local two phase flow momentum transport in S BWR