An investigation of subchannel analysis models for single-phase and two-phase flow

Energy Technology Data Exchange (ETDEWEB)

Hwang, Dae Hyun

1996-01-01

The governing equations and lateral transport modelings of subchannel analysis code, which is the most widely used tool for the analysis of thermal hydraulics fields in reactor cores, have been thoroughly investigated in this study. The procedure for the derivation of subchannel integral balance equations from the local instantaneous phase equations was investigated by stages. The characteristics of governing equations according to the treatment of phase velocity were studies, and the equations based on the drift-flux equilibrium formulation have been derived. Turbulent mixing and void drift modeling, which affect considerably to the accuracy of subchannel analysis code, have been reviewed. In addition, some representative modelings of single-phase and two-phase turbulent mixing models have been introduced. (author). 5 tabs., 4 figs., 16 refs.

Natural circulation in single-phase and two-phase flow

International Nuclear Information System (INIS)

Cheung, F.B.; El-Genk, M.S.

1989-01-01

Natural circulation usually arises in a closed loop between a heat source and a heat sink were the fluid motion is driven by density difference. It may also occur in enclosures or cavities where the flow is induced primarily by temperature or concentration gradients within the fluid. The subject has recently received special attention by the heat transfer and nuclear reactor safety communities because of it importance to the areas of energy extraction, decay, heat removal in nuclear reactors, solar and geothermal heating, and cooling of electronic equipment. Although many new results and physical insights have been gained of the various natural circulation phenomena, a number of critical issues remain unresolved. These include, for example, transition from laminar to turbulent flow, buoyancy-induced turbulent flow modeling, change of flow regimes, flow field visualization, variable property effects, and flow instability. This symposium volume contains papers presented in the Natural Circulation in Single-Phase and Two-Phase Flow session at the 1989 Winter Annual Meeting of ASME, by authors from different countries including the United States, Japan, Canada, and Brazil. The papers deal with experimental and theoretical studies as well as state-of-the-art reviews, covering a broad spectrum of topics in natural circulation including: variable-conductance thermosyphons, microelectronic chip cooling, natural circulation in anisotropic porous media and in cavities, heat transfer in flat plat solar collectors, shutdown heat removal in fast reactors, cooling of light-water and heavy-water reactors. The breadth of papers contained in this volume clearly reflect the importance of the current interest in natural circulation as a means for passive cooling and heating

RELAP-7 Progress Report: A Mathematical Model for 1-D Compressible, Single-Phase Flow Through a Branching Junction

Energy Technology Data Exchange (ETDEWEB)

Berry, R. A. [Idaho National Lab. (INL), Idaho Falls, ID (United States)

2017-08-14

In the literature, the abundance of pipe network junction models, as well as inclusion of dissipative losses between connected pipes with loss coefficients, has been treated using the incompressible flow assumption of constant density. This approach is fundamentally, physically wrong for compressible flow with density change. This report introduces a mathematical modeling approach for general junctions in piping network systems for which the transient flows are compressible and single-phase. The junction could be as simple as a 1-pipe input and 1-pipe output with differing pipe cross-sectional areas for which a dissipative loss is necessary, or it could include an active component, between an inlet pipe and an outlet pipe, such as a pump or turbine. In this report, discussion will be limited to the former. A more general branching junction connecting an arbitrary number of pipes with transient, 1-D compressible single-phase flows is also presented. These models will be developed in a manner consistent with the use of a general equation of state like, for example, the recent Spline-Based Table Look-up method [1] for incorporating the IAPWS-95 formulation [2] to give accurate and efficient calculations for properties for water and steam with RELAP-7 [3].

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)

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

Suggestion of an average bidirectional flow tube for the measurement of single and two phase flow rate

International Nuclear Information System (INIS)

Yun, B.J.; Kang, K.H.; Euh, D.J.; Song, C.H.; Baek, W.P.

2005-01-01

Full text of publication follows: A new type instrumentation, average bidirectional flow tube, was suggested to apply to the single and two phase flow condition. Its working principle is similar to that of the Pitot tube. The pressure measured at the front of the flow tube is equal to the total pressure, while that measured at the rear tube is slightly less than static pressure of flow field due to the suction effect at the downstream. It gives an amplification effect of measured pressure difference at the flow tube. The proposed instrumentation has the characteristics that it could be applicable to low flow condition and measure bidirectional flow. It was tested in the air-water vertical and horizontal test sections which have 0.08 m inner diameter. The pressure difference across the average bidirectional flow tube, system pressure, average void fraction and injection phasic mass flow rates were measured on the measuring plane. Test was performed primarily in the single phase water and air flow condition to get the amplification factor k of the flow tube. The test was also performed in the air-water two phase flow condition and the covered flow regimes were bubbly, slug, churn turbulent flow in the vertical pipe and stratified flow in the horizontal pipe. In order to calculate the phasic and total mass flow rates from the measured differential pressure, Chexal drift-flux correlation and momentum exchange factor between the two phases were introduced. The test result shows that the suggested instrumentation with the measured void fraction, Chexal drift-flux correlation and Bosio and Malnes' momentum exchange model can predict the phasic mass flow rates within 15% error compared to the true values. A new momentum exchange model was also suggested and it gives up to 5% improvement of the measured mass flow rate compared to combination of Bosio and Malnes' momentum exchange model. (authors)

Modeling and numerical analysis of non-equilibrium two-phase flows

International Nuclear Information System (INIS)

Rascle, P.; El Amine, K.

1997-01-01

We are interested in the numerical approximation of two-fluid models of nonequilibrium two-phase flows described by six balance equations. We introduce an original splitting technique of the system of equations. This technique is derived in a way such that single phase Riemann solvers may be used: moreover, it allows a straightforward extension to various and detailed exchange source terms. The properties of the fluids are first approached by state equations of ideal gas type and then extended to real fluids. For the construction of numerical schemes , the hyperbolicity of the full system is not necessary. When based on suitable kinetic unwind schemes, the algorithm can compute flow regimes evolving from mixture to single phase flows and vice versa. The whole scheme preserves the physical features of all the variables which remain in the set of physical states. Several stiff numerical tests, such as phase separation and phase transition are displayed in order to highlight the efficiency of the proposed method. The document is a PhD thesis divided in 6 chapters and two annexes. They are entitled: 1. - Introduction (in French), 2. - Two-phase flow, modelling and hyperbolicity (in French), 3. - A numerical method using upwind schemes for the resolution of two-phase flows without exchange terms (in English), 4. - A numerical scheme for one-phase flow of real fluids (in English), 5. - An upwind numerical for non-equilibrium two-phase flows (in English), 6. - The treatment of boundary conditions (in English), A.1. The Perthame scheme (in English) and A.2. The Roe scheme (in English)

Critical flow rate in a single phase flow. Blocking concept

International Nuclear Information System (INIS)

Giot, Michel

1978-01-01

After referring to the phenomena accompanying the appearance of a critical flow rate in a nozzle and presenting equations governing single phase flows, the critical condition is defined. Several particular cases are then examined; the horizontal and vertical isentropic flow, Fanno's flow and Raleigh's and the isothermal flow. The entropy deviation is calculated on either side of a normal impact. To conclude, the link existing between the concepts of critical flow and the propagation rate of small perturbations is demonstrated. To do so, the method of perturbations, that of Prandtl and that of characteristic directions are applied in turn [fr

A combined single-multiphase flow formulation of the premixing phase using the level set method

International Nuclear Information System (INIS)

Leskovar, M.; Marn, J.

1999-01-01

The premixing phase of a steam explosion covers the interaction of the melt jet or droplets with the water prior to any steam explosion occurring. To get a better insight of the hydrodynamic processes during the premixing phase beside hot premixing experiments, where the water evaporation is significant, also cold isothermal premixing experiments are performed. The specialty of isothermal premixing experiments is that three phases are involved: the water, the air and the spheres phase, but only the spheres phase mixes with the other two phases whereas the water and air phases do not mix and remain separated by a free surface. Our idea therefore was to treat the isothermal premixing process with a combined single-multiphase flow model. In this combined model the water and air phase are treated as a single phase with discontinuous phase properties at the water air interface, whereas the spheres are treated as usually with a multiphase flow model, where the spheres represent the dispersed phase and the common water-air phase represents the continuous phase. The common water-air phase was described with the front capturing method based on the level set formulation. In the level set formulation, the boundary of two-fluid interfaces is modeled as the zero set of a smooth signed normal distance function defined on the entire physical domain. The boundary is then updated by solving a nonlinear equation of the Hamilton-Jacobi type on the whole domain. With this single-multiphase flow model the Queos isothermal premixing Q08 has been simulated. A numerical analysis using different treatments of the water-air interface (level set, high-resolution and upwind) has been performed for the incompressible and compressible case and the results were compared to experimental measurements.(author)

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

Single-phase flow and flow boiling of water in horizontal rectangular microchannels

OpenAIRE

Mirmanto

2013-01-01

This thesis was submitted for the degree of Doctor of Philosophy and awarded by Brunel University The current study is part of a long term experimental project devoted to investigating single-phase flow pressure drop and heat transfer, flow boiling pressure drop and heat transfer, flow boiling instability and flow visualization of de-ionized water flow in microchannels. The experimental facility was first designed and constructed by S. Gedupudi (2009) and in the present study; ...

Pressure drop and heat transfer of a mercury single-phase flow and an air-mercury two-phase flow in a helical tube under a strong magnetic field

International Nuclear Information System (INIS)

Takahashi, Minoru; Momozaki, Yoichi

2000-01-01

For the reduction of a large magneto-hydrodynamic (MHD) pressure drop of a liquid metal single-phase flow, a liquid metal two-phase flow cooling system has been proposed. As a fundamental study, MHD pressure drops and heat transfer characteristics of a mercury single-phase flow and an air-mercury two-phase flow were experimentally investigated. A strong transverse magnetic field relevant to the fusion reactor conditions was applied to the mercury single-phase flow and the air-mercury two-phase flow in a helically coiled tube that was inserted in the vertical bore of a solenoidal superconducting magnet. It was found that MHD pressure drops of a mercury single-phase flow in the helically coiled tube were nearly equal to those in a straight tube. The Nusselt number at an outside wall was higher than that at an inside wall both in the mercury single-phase flow in the absence and presence of a magnetic field. The Nusselt number of the mercury single-phase flow decreased, increased and again decreased with an increase in the magnetic flux density. MHD pressure drops did not decrease appreciably by injecting air into a mercury flow and changing the mercury flow into the air-mercury two-phase flow. Remarkable heat transfer enhancement did not appear by the air injection. The injection of air into the mercury flow enhanced heat transfer in the ranges of high mercury flow rate and low magnetic flux density, possibly due to the agitation effect of air bubbles. The air injection deteriorated heat transfer in the range of low mercury flow rates possibly because of the occupation of air near heating wall

Research on three-phase traffic flow modeling based on interaction range

Science.gov (United States)

Zeng, Jun-Wei; Yang, Xu-Gang; Qian, Yong-Sheng; Wei, Xu-Ting

2017-12-01

On the basis of the multiple velocity difference effect (MVDE) model and under short-range interaction, a new three-phase traffic flow model (S-MVDE) is proposed through careful consideration of the influence of the relationship between the speeds of the two adjacent cars on the running state of the rear car. The random slowing rule in the MVDE model is modified in order to emphasize the influence of vehicle interaction between two vehicles on the probability of vehicles’ deceleration. A single-lane model which without bottleneck structure under periodic boundary conditions is simulated, and it is proved that the traffic flow simulated by S-MVDE model will generate the synchronous flow of three-phase traffic theory. Under the open boundary, the model is expanded by adding an on-ramp, the congestion pattern caused by the bottleneck is simulated at different main road flow rates and on-ramp flow rates, which is compared with the traffic congestion pattern observed by Kerner et al. and it is found that the results are consistent with the congestion characteristics in the three-phase traffic flow theory.

Single- and two-phase flow modeling for coupled neutronics / thermal-hydraulics transient analysis of advanced sodium-cooled fast reactors

International Nuclear Information System (INIS)

Chenu, A.

2011-10-01

Nuclear power is nowadays in the front rank as regards helping to meet the growing worldwide energy demand while avoiding an excessive increase in greenhouse gas emissions. However, the operating nuclear power plants are mainly thermal-neutron reactors and, as such, can not be maintained on the basis of the currently identified uranium resources beyond one century at the present consumption rate. Sustainability of nuclear power thus involves closure of the fuel cycle through breeding. With a uranium-based fuel, breeding can only be achieved using a fast-neutron reactor. Sodium-cooled fast reactor (SFR) technology benefits from 400 reactor-years of accumulated experience and is thus a prime candidate for the implementation of so-called Generation-IV nuclear energy systems. In this context, the safety demonstration of SFRs remains a major Research and Development related issue. The current research aims at the development of a computational tool for the in-depth understanding of SFR core behaviour during accidental transients, particularly those including boiling of the coolant. An accurate modelling of the core physics during such transients requires the coupling between 3D neutron kinetics and thermal-hydraulics in the core, to account for the strong interactions between the two-phase coolant flow and power variations caused by the sodium void effect. The present study is specifically focused upon models for the representation of sodium two-phase flow. The extension of the thermal-hydraulics TRACE code, previously limited to the simulation of single-phase sodium flow, has been carried out through the implementation of equations-of-state and closure relations specific to sodium. The different correlations have then been implemented as options. From the validation study carried out, it has been possible to recommend a set of models which provide satisfactory results, while considering annular flow as the dominant regime up to dryout and a smooth breakdown of the

Two Phase Flow Split Model for Parallel Channels | Iloeje | Nigerian ...

African Journals Online (AJOL)

The model and code are capable of handling single and two phase flows, steady states and transients, up to ten parallel flow paths, simple and complicated geometries, including the boilers of fossil steam generators and nuclear power plants. A test calculation has been made with a simplified three-channel system ...

Mathematical modelling of two-phase flows

International Nuclear Information System (INIS)

Komen, E.M.J.; Stoop, P.M.

1992-11-01

A gradual shift from methods based on experimental correlations to methods based on mathematical models to study 2-phase flows can be observed. The latter can be used to predict dynamical behaviour of 2-phase flows. This report discusses various mathematical models for the description of 2-phase flows. An important application of these models can be found in thermal-hydraulic computer codes used for analysis of the thermal-hydraulic behaviour of water cooled nuclear power plants. (author). 17 refs., 7 figs., 6 tabs

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

An algebraic stress/flux model for two-phase turbulent flow

International Nuclear Information System (INIS)

Kumar, R.

1995-12-01

An algebraic stress model (ASM) for turbulent Reynolds stress and a flux model for turbulent heat flux are proposed for two-phase bubbly and slug flows. These mathematical models are derived from the two-phase transport equations for Reynolds stress and turbulent heat flux, and provide C μ , a turbulent constant which defines the level of eddy viscosity, as a function of the interfacial terms. These models also include the effect of heat transfer. When the interfacial drag terms and the interfacial momentum transfer terms are absent, the model reduces to a single-phase model used in the literature

Investigation for vertical, two-phase steam-water flow of three turbine models

International Nuclear Information System (INIS)

Silverman, S.; Goodrich, L.D.

1977-01-01

One of the basic quantities of interest during a loss-of-coolant experiment (LOCE) is the primary system mass flow rate. Presently, there are no transducers commercially available which continuously measure this parameter. Therefore, a transducer was designed at EG and G Idaho, Inc. which combines a drag-disc and turbine into a single unit. The basis for the design was that the drag-disc would measure momentum flux (rhoV 2 ), the turbine would measure velocity and the mass flow rate could then be calculated from the two quantities by assuming a flow profile. For two-phase flow, the outputs are approximately proportional to the desired parameter, but rather large errors can be expected under those assumptions. Preliminary evaluation of the experimental two- and single-phase calibration data has resulted in uncertainty estimates of +-8% of range for the turbine and +-20% of range for the drag-disc. In an effort to reduce the errors, further investigations were made to determine what the drag-disc and turbine really measure. In the present paper, three turbine models for vertical, two-phase, steam/water flow are investigated; the Aya Model, the Rouhani Model, and a volumetric flow model. Theoretical predictions are compared with experimental data for vertical, two-phase steam/water flow. For the purposes of the mass flow calculation, velocity profiles were assumed to be flat for the free-field condition. It is appreciated that this may not be true for all cases investigated, but for an initial inspection, flat profiles were assumed

Study of a new static mixer for two-phase and single-phase flows

International Nuclear Information System (INIS)

Foucrier, Michel

1996-01-01

The subject of this work is the study of OptimiX, a new static mixer, which is fully designed using an inverse method taking the final product features as input and based on the physical properties of the fluid to mix. The work began with the construction of an experimental loop which allowed us to qualify the mixer in two-phase and single-phase flow conditions. Next, a chemical method using a new test reaction and a micro-mixing model have been used to further characterise the mixer. This test reaction and the micro-mixing model have been developed by the 'Laboratoire des Sciences du Genie Chimique' of Nancy. The mixer OptimiX has proved to be an excellent device for both macro- and micro-mixing. The capability of this mixer to foster rapid reactions was also demonstrated. The well organised flow pattern of OptimiX, which results from its design, provides it with the unusual feature of being fully calculable. This work emphasizes the internal hydrodynamics of this mixer, justifies the universality of the design procedures, which validation is supported by the completed qualification work. (author) [fr

Multi-scale Modeling of Compressible Single-phase Flow in Porous Media using Molecular Simulation

KAUST Repository

Saad, Ahmed Mohamed

2016-01-01

potential model that accounts for the molecular quadrupole moment of fluids with non-spherical molecules such as CO2. The potential model was used to simulate the thermodynamic equilibrium properties for single-phase and two-phase systems using the canonical

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

Single and two-phase flow pressure drop for CANFLEX bundle

Energy Technology Data Exchange (ETDEWEB)

Park, Joo Hwan; Jun, Ji Su; Suk, Ho Chun [Korea Atomic Energy Research Institute, Taejon (Korea, Republic of); Dimmick, G R; Bullock, D E [Atomic Energy of Canada Limited, Ontario (Canada)

1999-12-31

Friction factor and two-phase flow frictional multiplier for a CANFLEX bundle are newly developed and presented in this paper. CANFLEX as a 43-element fuel bundle has been developed jointly by AECL/KAERI to provide greater operational flexibility for CANDU reactor operators and designers. Friction factor and two-phase flow frictional multiplier have been developed by using the experimental data of pressure drops obtained from two series of Freon-134a (R-134a) CHF tests with a string of simulated CANFLEX bundles in a single phase and a two-phase flow conditions. The friction factor for a CANFLEX bundle is found to be about 20% higher than that of Blasius for a smooth circular pipe. The pressure drop predicted by using the new correlations of friction factor and two-phase frictional multiplier are well agreed with the experimental pressure drop data of CANFLEX bundle within {+-} 5% error. 11 refs., 5 figs. (Author)

Single and two-phase flow pressure drop for CANFLEX bundle

Energy Technology Data Exchange (ETDEWEB)

Park, Joo Hwan; Jun, Ji Su; Suk, Ho Chun [Korea Atomic Energy Research Institute, Taejon (Korea, Republic of); Dimmick, G. R.; Bullock, D. E. [Atomic Energy of Canada Limited, Ontario (Canada)

1998-12-31

Friction factor and two-phase flow frictional multiplier for a CANFLEX bundle are newly developed and presented in this paper. CANFLEX as a 43-element fuel bundle has been developed jointly by AECL/KAERI to provide greater operational flexibility for CANDU reactor operators and designers. Friction factor and two-phase flow frictional multiplier have been developed by using the experimental data of pressure drops obtained from two series of Freon-134a (R-134a) CHF tests with a string of simulated CANFLEX bundles in a single phase and a two-phase flow conditions. The friction factor for a CANFLEX bundle is found to be about 20% higher than that of Blasius for a smooth circular pipe. The pressure drop predicted by using the new correlations of friction factor and two-phase frictional multiplier are well agreed with the experimental pressure drop data of CANFLEX bundle within {+-} 5% error. 11 refs., 5 figs. (Author)

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.

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

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

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

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

On the use of nuclear magnetic resonance to measure velocity and its fluctuations in single-phase and two-phase flows

International Nuclear Information System (INIS)

Jullien, Pierre

2013-01-01

This work deals with the use of NMR to measure velocity and its fluctuations in single-phase and two-phase flows. PGSE and imaging sequences have been used to determine the velocity distributions in upward turbulent pipe flows. NMR signals have been analysed in detail and the main artifacts have been characterized and suppressed. The measuring technique has been validated by comparison with a reference published data. A first comparison to 'homemade' hot-wire results in single-phase flow of water is presented and is very promising. Preliminary NMR results in two-phase flows emphasize the interest of NMR to benchmark velocity measurements in two-phase flows. Prospects of research have been identified, which will pave the way for the sequel of this research. (author) [fr

Single phase flow pressure drop and heat transfer in rectangular metallic microchannels

International Nuclear Information System (INIS)

Sahar, Amirah M.; Özdemir, Mehmed R.; Fayyadh, Ekhlas M.; Wissink, Jan; Mahmoud, Mohamed M.; Karayiannis, Tassos G.

2016-01-01

Numerical simulations were performed using Fluent 14.5 to investigate single phase flow and conjugate heat transfer in copper rectangular microchannels. Two different configurations were simulated: (1) single channel with hydraulic diameter of 0.561 mm and (2) multichannel configuration consisting of inlet and outlet manifolds and 25 channels with hydraulic diameter of 0.409 mm. In the single channel configuration, four numerical models were investigated namely, 2D thin-wall, 3D thin-wall (heated from the bottom), 3D thin-wall (three side heated) and 3D full conjugate models. In the multichannel configuration, only 3D full conjugate model was used. The simulation results of the single channel configuration were validated using experimental data of water as a test fluid while the results of the multichannel configuration were validated using experimental data of R134a refrigerant. In the multichannel configuration, flow distribution among the channels was also investigated. The 3D thin-wall model simulation was conducted at thermal boundary conditions similar to those assumed in the experimental data reduction (uniform heat flux) and showed excellent agreement with the experimental data. However, the results of the 3D full conjugate model demonstrated that there is a significant conjugate effect and the heat flux is not uniformly distributed along the channel resulting in significant deviation compared to the experimental data (more than 50%). Also, the results demonstrated that there is a significant difference between the 3D thin-wall and full conjugate models. The simulation of the multichannel configuration with an inlet manifold having gradual decrease in cross sectional area achieved very reasonable uniform flow distribution among the channels which will provide uniform heat transfer rates across the base of the microchannels.

Simulating single-phase and two-phase non-Newtonian fluid flow of a digital rock scanned at high resolution

Science.gov (United States)

Tembely, Moussa; Alsumaiti, Ali M.; Jouini, Mohamed S.; Rahimov, Khurshed; Dolatabadi, Ali

2017-11-01

Most of the digital rock physics (DRP) simulations focus on Newtonian fluids and overlook the detailed description of rock-fluid interaction. A better understanding of multiphase non-Newtonian fluid flow at pore-scale is crucial for optimizing enhanced oil recovery (EOR). The Darcy scale properties of reservoir rocks such as the capillary pressure curves and the relative permeability are controlled by the pore-scale behavior of the multiphase flow. In the present work, a volume of fluid (VOF) method coupled with an adaptive meshing technique is used to perform the pore-scale simulation on a 3D X-ray micro-tomography (CT) images of rock samples. The numerical model is based on the resolution of the Navier-Stokes equations along with a phase fraction equation incorporating the dynamics contact model. The simulations of a single phase flow for the absolute permeability showed a good agreement with the literature benchmark. Subsequently, the code is used to simulate a two-phase flow consisting of a polymer solution, displaying a shear-thinning power law viscosity. The simulations enable to access the impact of the consistency factor (K), the behavior index (n), along with the two contact angles (advancing and receding) on the relative permeability.

Two-phase wall function for modeling of turbulent boundary layer in subcooled boiling flow

International Nuclear Information System (INIS)

Bostjan Koncar; Borut Mavko; Yassin A Hassan

2005-01-01

Full text of publication follows: The heat transfer and phase-change mechanisms in the subcooled flow boiling are governed mainly by local multidimensional mechanisms near the heated wall, where bubbles are generated. The structure of such 'wall boiling flow' is inherently non-homogeneous and is further influenced by the two-phase flow turbulence, phase-change effects in the bulk, interfacial forces and bubble interactions (collisions, coalescence, break-up). In this work the effect of two-phase flow turbulence on the development of subcooled boiling flow is considered. Recently, the modeling of two-phase flow turbulence has been extensively investigated. A notable progress has been made towards deriving reliable models for description of turbulent behaviour of continuous (liquid) and dispersed phase (bubbles) in the bulk flow. However, there is a lack of investigation considering the modeling of two-phase flow boundary layer. In most Eulerian two-fluid models standard single-phase wall functions are used for description of turbulent boundary layer of continuous phase. That might be a good approximation at adiabatic flows, but their use for boundary layers with high concentration of dispersed phase is questionable. In this work, the turbulent boundary layer near the heated wall will be modeled with the so-called 'two-phase' wall function, which is based on the assumption of additional turbulence due to bubble-induced stirring in the boundary layer. In the two-phase turbulent boundary layer the wall function coefficients strongly depend on the void fraction. Moreover, in the turbulent boundary layer with nucleating bubbles, the bubble size variation also has a significant impact on the liquid phase. As a basis, the wall function of Troshko and Hassan (2001), developed for adiabatic bubbly flows will be used. The simulations will be performed by a general-purpose CFD code CFX-4.4 using additional models provided by authors. The results will be compared to the boiling

Factors that affect the calibration of turbines in single-phase flow

International Nuclear Information System (INIS)

Piper, T.C.

1977-05-01

Basic turbine operation in single-phase flow is related. Causes and relative magnitudes of retarding torque are given for two sizes of turbines when used for water flow measurement. An equation for slip caused by retarding torques is given. Evaluation of turbine slip behavior at the turbine low flow region shows that bearing retarding torques, change in flow patterns, or other effects can predominate in the relatively large changes in the calibration ''constant'' that occurs there. Fluid lubricity is singled out as an important fluid property in certain types of bearings and flow. Temperature induced changes in turbine size are shown to cause calibration changes if a turbine is used at a temperature significantly different than that at which it was calibrated

Factors that affect the calibration of turbines in single-phase flow

Energy Technology Data Exchange (ETDEWEB)

Piper, T. C.

1977-05-01

Basic turbine operation in single-phase flow is related. Causes and relative magnitudes of retarding torque are given for two sizes of turbines when used for water flow measurement. An equation for slip caused by retarding torques is given. Evaluation of turbine slip behavior at the turbine low flow region shows that bearing retarding torques, change in flow patterns, or other effects can predominate in the relatively large changes in the calibration ''constant'' that occurs there. Fluid lubricity is singled out as an important fluid property in certain types of bearings and flow. Temperature induced changes in turbine size are shown to cause calibration changes if a turbine is used at a temperature significantly different than that at which it was calibrated.

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.

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

Heat transfer and velocity characteristics of single- and two-phase flows in a subsonic model gun

International Nuclear Information System (INIS)

Bicen, A.F.; Khezzar, L.; Schmidt, M.; Whitelaw, J.H.

1989-01-01

Heat transfer and velocity measurements are reported for single- and two-phase flows in the wake of an in-bore projectile propelled by an inert gas at an initial gauge pressure of 8 bars to an exit velocity over 40 m/s in ∼ 33 ms. The results show that with the single phase the turbulent velocity boundary layers occupy over 20% of the barrel radius and that the wall heat transfer increases with distance from the breech and decreases with time during the shot. In the initial chamber, and later in the shot, the heat transfer results are close to those obtained from a convection correlation for a steady turbulent boundary layer, contrary to those at locations swept by the projectile, which are higher by up to 50% throughout the shot. The two-phase flow results show that 55-μm particles with loadings of 1.3% and 4% by volume initially lag the fluid and this lag increases with distance from the breech. Later in the shot the particles catch up and lead the decelerating fluid by an amount that is greater, with the higher particle loading and with a tendency for the particle velocity to increase around the edge of the boundary layer

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.

Microfluidic Manufacturing of Polymeric Nanoparticles: Comparing Flow Control of Multiscale Structure in Single-Phase Staggered Herringbone and Two-Phase Reactors.

Science.gov (United States)

Xu, Zheqi; Lu, Changhai; Riordon, Jason; Sinton, David; Moffitt, Matthew G

2016-12-06

We compare the microfluidic manufacturing of polycaprolactone-block-poly(ethylene oxide) (PCL-b-PEO) nanoparticles (NPs) in a single-phase staggered herringbone (SHB) mixer and in a two-phase gas-liquid segmented mixer. NPs generated from two different copolymer compositions in both reactors and at three different flow rates, along with NPs generated using a conventional bulk method, are compared with respect to morphologies, dimensions, and internal crystallinities. Our work, the first direct comparison between alternate microfluidic NP synthesis methods, shows three key findings: (i) NP morphologies and dimensions produced in the bulk are different from those produced in a microfluidic mixer, whereas NP crystallinities produced in the bulk and in the SHB mixer are similar; (ii) NP morphologies, dimensions, and crystallinities produced in the single-phase SHB and two-phase mixers at the lowest flow rate are similar; and (iii) NP morphologies, dimensions, and crystallinities change with flow rate in the two-phase mixer but not in the single-phase SHB mixer. These findings provide new insights into the relative roles of mixing and shear in the formation and flow-directed processing of polymeric NPs in microfluidics, informing future reactor designs for manufacturing NPs of low polydispersity and controlled multiscale structure and function.

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

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

A facility for the experimental investigation of single substance two phase flow

International Nuclear Information System (INIS)

Maeder, P.F.; Dickinson, D.A.; Nikitopoulos, D.E.; DiPippo, R.

1985-01-01

The paper describes a research facility dedicated to single-substance two-phase flow. The working fluid is dichlorotetrafluoroethane (or refrigerant R-114), allowing both operation at manageable pressures, temperatures and flowrates, and application of results to practical situations through similarity. Operation is in the blowdown mode. The control and data acquisition systems are fully automated and computer controlled. A range of flow conditions from predominantly liquid flow to high velocity, high void fraction choked flow can be attained

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

A novel deep reactive ion etched (DRIE) glass micro-model for two-phase flow experiments.

Science.gov (United States)

Karadimitriou, N K; Joekar-Niasar, V; Hassanizadeh, S M; Kleingeld, P J; Pyrak-Nolte, L J

2012-09-21

In the last few decades, micro-models have become popular experimental tools for two-phase flow studies. In this work, the design and fabrication of an innovative, elongated, glass-etched micro-model with dimensions of 5 × 35 mm(2) and constant depth of 43 microns is described. This is the first time that a micro-model with such depth and dimensions has been etched in glass by using a dry etching technique. The micro-model was visualized by a novel setup that allowed us to monitor and record the distribution of fluids throughout the length of the micro-model continuously. Quasi-static drainage experiments were conducted in order to obtain equilibrium data points that relate capillary pressure to phase saturation. By measuring the flow rate of water through the flow network for known pressure gradients, the intrinsic permeability of the micro-model's flow network was also calculated. The experimental results were used to calibrate a pore-network model and test its validity. Finally, we show that glass-etched micro-models can be valuable tools in single and/or multi-phase flow studies and their applications.

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.

Two-phase electro-hydrodynamic flow modeling by a conservative level set model.

Science.gov (United States)

Lin, Yuan

2013-03-01

The principles of electro-hydrodynamic (EHD) flow have been known for more than a century and have been adopted for various industrial applications, for example, fluid mixing and demixing. Analytical solutions of such EHD flow only exist in a limited number of scenarios, for example, predicting a small deformation of a single droplet in a uniform electric field. Numerical modeling of such phenomena can provide significant insights about EHDs multiphase flows. During the last decade, many numerical results have been reported to provide novel and useful tools of studying the multiphase EHD flow. Based on a conservative level set method, the proposed model is able to simulate large deformations of a droplet by a steady electric field, which is beyond the region of theoretic prediction. The model is validated for both leaky dielectrics and perfect dielectrics, and is found to be in excellent agreement with existing analytical solutions and numerical studies in the literature. Furthermore, simulations of the deformation of a water droplet in decyl alcohol in a steady electric field match better with published experimental data than the theoretical prediction for large deformations. Therefore the proposed model can serve as a practical and accurate tool for simulating two-phase EHD flow. © 2013 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim.

A theoretical model for measuring mass flowrate and quality of two phase flow by the noise of throttling set

International Nuclear Information System (INIS)

Tong Yunxian; Wang Wenran

1992-03-01

The mass flowrate and steam quality measuring of two phase flowrate is an essential issue in the tests of loss-of-coolant accident (LOCA). The spatial stochastic distribution of phase concentration would cause a differential pressure noise when two phase flow is crossing a throttling set. Under the assumption of that the variance of disperse phase concentration is proportional to its mean phase concentration and by using the separated flow model of two phase flow, it has demonstrated that the variance of noise of differential pressure square root is approximately proportional to the flowrate of disperse phase. Thus, a theoretical model for measuring mass flowrate and quality of two phase flow by noise measurement is developed. It indicates that there is a possibility to measure two phase flowrate and steam quality by using the simple theoretical model and a single throttling set

New theoretical model for two-phase flow discharged from stratified two-phase region through small break

International Nuclear Information System (INIS)

Yonomoto, Taisuke; Tasaka, Kanji

1988-01-01

A theoretical and experimental study was conducted to understand two-phase flow discharged from a stratified two-phase region through a small break. This problem is important for an analysis of a small break loss-of-coolant accident (LOCA) in a light water reactor (LWR). The present theoretical results show that a break quality is a function of h/h b , where h is the elevation difference between a bulk water level in the upstream region and break and b the suffix for entrainment initiation. This result is consistent with existing eperimental results in literature. An air-water experiment was also conducted changing a break orientation as an experimental parameter to develop and assess the model. Comparisons between the model and the experimental results show that the present model can satisfactorily predict the flow rate and the quality at the break without using any adjusting constant when liquid entrainment occurs in a stratified two-phase region. When gas entrainment occurs, the experimental data are correlated well by using a single empirical constant. (author)

A novel mechanical model for phase-separation in debris flows

Science.gov (United States)

Pudasaini, Shiva P.

2015-04-01

Understanding the physics of phase-separation between solid and fluid phases as a two-phase mass moves down slope is a long-standing challenge. Here, I propose a fundamentally new mechanism, called 'separation-flux', that leads to strong phase-separation in avalanche and debris flows. This new model extends the general two-phase debris flow model (Pudasaini, 2012) to include a separation-flux mechanism. The new flux separation mechanism is capable of describing and controlling the dynamically evolving phase-separation, segregation, and/or levee formation in a real two-phase, geometrically three-dimensional debris flow motion and deposition. These are often observed phenomena in natural debris flows and industrial processes that involve the transportation of particulate solid-fluid mixture material. The novel separation-flux model includes several dominant physical and mechanical aspects that result in strong phase-separation (segregation). These include pressure gradients, volume fractions of solid and fluid phases and their gradients, shear-rates, flow depth, material friction, viscosity, material densities, boundary structures, gravity and topographic constraints, grain shape, size, etc. Due to the inherent separation mechanism, as the mass moves down slope, more and more solid particles are brought to the front, resulting in a solid-rich and mechanically strong frontal surge head followed by a weak tail largely consisting of the viscous fluid. The primary frontal surge head followed by secondary surge is the consequence of the phase-separation. Such typical and dominant phase-separation phenomena are revealed here for the first time in real two-phase debris flow modeling and simulations. However, these phenomena may depend on the bulk material composition and the applied forces. Reference: Pudasaini, Shiva P. (2012): A general two-phase debris flow model. J. Geophys. Res., 117, F03010, doi: 10.1029/2011JF002186.

Numerical simulation for gas-liquid two-phase flow in pipe networks

International Nuclear Information System (INIS)

Li Xiaoyan; Kuang Bo; Zhou Guoliang; Xu Jijun

1998-01-01

The complex pipe network characters can not directly presented in single phase flow, gas-liquid two phase flow pressure drop and void rate change model. Apply fluid network theory and computer numerical simulation technology to phase flow pipe networks carried out simulate and compute. Simulate result shows that flow resistance distribution is non-linear in two phase pipe network

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.

Experimental and numerical study on single-phase flow characteristics of natural circulation system with heated narrow rectangular channel under rolling motion condition

International Nuclear Information System (INIS)

Yu, Shengzhi; Wang, Jianjun; Yan, Ming; Yan, Changqi; Cao, Xiaxin

2017-01-01

Highlights: • The phasic difference between flow rate and frictional pressure drop is negligible. • Effect mechanism of rolling motion on flow behaviors of NC is interpreted. • The startup model is proposed and verified. • Steady-state correlations are feasible to predict transient resistance. • The in-house code can simulate instantaneous flow behaviors of NC correctly. - Abstract: Effects of rolling motion on flow characteristics in a natural circulation system were investigated experimentally and numerically. The numerical results from validated code were mainly used to provide detailed information for the discussion and analysis of experimental results. The results indicate that under rolling motion condition, the phasic difference between flow rate and frictional pressure drop of narrow rectangular channel is negligible. Angular acceleration is the eigenvalue for the effects of rolling motion on flow rate under single-phase natural circulation condition. When angular acceleration is approximately equal, even though either the angle or the period of rolling motion is different, peak, trough and time-averaged values of flow rate are approximately equal. Under rolling motion and single-phase natural circulation conditions, the phenomenon that dimensionless time-averaged mass flow rate is smaller than that under steady state condition is controlled by the nonlinear relationship between mass flow rate and the resistance of loop. The factor also causes the result that the absolute difference of dimensionless flow rate between peak and steady state is smaller than that between trough and steady state. The startup model which is proposed in present paper can be used to predict the flow characteristics of single-phase natural circulation system at startup stage of rolling motion favorably. The self-developed code can simulate instantaneous flow characteristics of single-phase natural circulation system under rolling motion and steady state conditions

Mathematical model of two-phase flow in accelerator channel

Directory of Open Access Journals (Sweden)

О.Ф. Нікулін

2010-01-01

Full Text Available  The problem of  two-phase flow composed of energy-carrier phase (Newtonian liquid and solid fine-dispersed phase (particles in counter jet mill accelerator channel is considered. The mathematical model bases goes on the supposition that the phases interact with each other like independent substances by means of aerodynamics’ forces in conditions of adiabatic flow. The mathematical model in the form of system of differential equations of order 11 is represented. Derivations of equations by base physical principles for cross-section-averaged quantity are produced. The mathematical model can be used for estimation of any kinematic and thermodynamic flow characteristics for purposely parameters optimization problem solving and transfer functions determination, that take place in  counter jet mill accelerator channel design.

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

Similarity analysis and scaling criteria for LWRs under single-phase and two-phase natural circulation

International Nuclear Information System (INIS)

Ishii, M.; Kataoka, I.

1983-03-01

Scaling criteria for a natural circulation loop under single phase and two-phase flow conditions have been derived. For a single phase case the continuity, integral momentum, and energy equations in one-dimensional area average forms have been used. From this, the geometrical similarity groups, friction number, Richardson number, characteristic time constant ratio, Biot number, and heat source number are obtained. The Biot number involves the heat transfer coefficient which may cause some difficulties in simulating the turbulent flow regime. For a two-phase flow case, the similarity groups obtained from a perturbation analysis based on the one-dimensional drift-flux model have been used. The physical significance of the phase change number, subcooling number, drift-flux number, friction number are discussed and conditions imposed by these groups are evaluated. In the two-phase flow case, the critical heat flux is one of the most important transients which should be simulated in a scale model. The above results are applied to the LOFT facility in case of a natural circulation simulation. Some preliminary conclusions on the feasibility of the facility have been obtained

Similarity analysis and scaling criteria for LWRs under single-phase and two-phase natural circulation

Energy Technology Data Exchange (ETDEWEB)

Ishii, M.; Kataoka, I.

1983-03-01

Scaling criteria for a natural circulation loop under single phase and two-phase flow conditions have been derived. For a single phase case the continuity, integral momentum, and energy equations in one-dimensional area average forms have been used. From this, the geometrical similarity groups, friction number, Richardson number, characteristic time constant ratio, Biot number, and heat source number are obtained. The Biot number involves the heat transfer coefficient which may cause some difficulties in simulating the turbulent flow regime. For a two-phase flow case, the similarity groups obtained from a perturbation analysis based on the one-dimensional drift-flux model have been used. The physical significance of the phase change number, subcooling number, drift-flux number, friction number are discussed and conditions imposed by these groups are evaluated. In the two-phase flow case, the critical heat flux is one of the most important transients which should be simulated in a scale model. The above results are applied to the LOFT facility in case of a natural circulation simulation. Some preliminary conclusions on the feasibility of the facility have been obtained.

A review of solid-fluid selection options for optical-based measurements in single-phase liquid, two-phase liquid-liquid and multiphase solid-liquid flows

Science.gov (United States)

Wright, Stuart F.; Zadrazil, Ivan; Markides, Christos N.

2017-09-01

Experimental techniques based on optical measurement principles have experienced significant growth in recent decades. They are able to provide detailed information with high-spatiotemporal resolution on important scalar (e.g., temperature, concentration, and phase) and vector (e.g., velocity) fields in single-phase or multiphase flows, as well as interfacial characteristics in the latter, which has been instrumental to step-changes in our fundamental understanding of these flows, and the development and validation of advanced models with ever-improving predictive accuracy and reliability. Relevant techniques rely upon well-established optical methods such as direct photography, laser-induced fluorescence, laser Doppler velocimetry/phase Doppler anemometry, particle image/tracking velocimetry, and variants thereof. The accuracy of the resulting data depends on numerous factors including, importantly, the refractive indices of the solids and liquids used. The best results are obtained when the observational materials have closely matched refractive indices, including test-section walls, liquid phases, and any suspended particles. This paper reviews solid-liquid and solid-liquid-liquid refractive-index-matched systems employed in different fields, e.g., multiphase flows, turbomachinery, bio-fluid flows, with an emphasis on liquid-liquid systems. The refractive indices of various aqueous and organic phases found in the literature span the range 1.330-1.620 and 1.251-1.637, respectively, allowing the identification of appropriate combinations to match selected transparent or translucent plastics/polymers, glasses, or custom materials in single-phase liquid or multiphase liquid-liquid flow systems. In addition, the refractive indices of fluids can be further tuned with the use of additives, which also allows for the matching of important flow similarity parameters such as density and viscosity.

Analysis of flow induced valve operation and pressure wave propagation for single and two-phase flow conditions

International Nuclear Information System (INIS)

Nagel, H.

1986-01-01

The flow induced valve operation is calculated for single and two-phase flow conditions by the fluid dynamic computer code DYVRO and results are compared to experimental data. The analysis show that the operational behaviour of the valves is not only dependent on the condition of the induced flow, but also the pipe flow can cause a feedback as a result of the induced pressure waves. For the calculation of pressure wave propagation in pipes of which the operation of flow induced valves has a considerable influence it is therefore necessary to have a coupled analysis of the pressure wave propagation and the operational behaviour of the valves. The analyses of the fast transient transfer from steam to two-phase flow show a good agreement with experimental data. Hence even these very high loads on pipes resulting from such fluid dynamic transients can be calculated realistically. (orig.)

Validation of model predictions of pore-scale fluid distributions during two-phase flow

Science.gov (United States)

Bultreys, Tom; Lin, Qingyang; Gao, Ying; Raeini, Ali Q.; AlRatrout, Ahmed; Bijeljic, Branko; Blunt, Martin J.

2018-05-01

Pore-scale two-phase flow modeling is an important technology to study a rock's relative permeability behavior. To investigate if these models are predictive, the calculated pore-scale fluid distributions which determine the relative permeability need to be validated. In this work, we introduce a methodology to quantitatively compare models to experimental fluid distributions in flow experiments visualized with microcomputed tomography. First, we analyzed five repeated drainage-imbibition experiments on a single sample. In these experiments, the exact fluid distributions were not fully repeatable on a pore-by-pore basis, while the global properties of the fluid distribution were. Then two fractional flow experiments were used to validate a quasistatic pore network model. The model correctly predicted the fluid present in more than 75% of pores and throats in drainage and imbibition. To quantify what this means for the relevant global properties of the fluid distribution, we compare the main flow paths and the connectivity across the different pore sizes in the modeled and experimental fluid distributions. These essential topology characteristics matched well for drainage simulations, but not for imbibition. This suggests that the pore-filling rules in the network model we used need to be improved to make reliable predictions of imbibition. The presented analysis illustrates the potential of our methodology to systematically and robustly test two-phase flow models to aid in model development and calibration.

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.

Single- and two-phase flow characterization using optical fiber bragg gratings.

Science.gov (United States)

Baroncini, Virgínia H V; Martelli, Cicero; da Silva, Marco José; Morales, Rigoberto E M

2015-03-17

Single- and two-phase flow characterization using optical fiber Bragg gratings (FBGs) is presented. The sensor unit consists of the optical fiber Bragg grating positioned transversely to the flow and fixed in the pipe walls. The hydrodynamic pressure applied by the liquid or air/liquid flow to the optical fiber induces deformation that can be detected by the FBG. Given that the applied pressure is directly related to the mass flow, it is possible to establish a relationship using the grating resonance wavelength shift to determine the mass flow when the flow velocity is well known. For two phase flows of air and liquid, there is a significant change in the force applied to the fiber that accounts for the very distinct densities of these substances. As a consequence, the optical fiber deformation and the correspondent grating wavelength shift as a function of the flow will be very different for an air bubble or a liquid slug, allowing their detection as they flow through the pipe. A quasi-distributed sensing tool with 18 sensors evenly spread along the pipe is developed and characterized, making possible the characterization of the flow, as well as the tracking of the bubbles over a large section of the test bed. Results show good agreement with standard measurement methods and open up plenty of opportunities to both laboratory measurement tools and field applications.

Numerical modelling of single-phase flow in rough fractures with contacts

Science.gov (United States)

Olkiewicz, Piotr; Dabrowski, Marcin

2017-04-01

Fracture flow may dominate in rocks with low porosity and it can accompany both industrial and natural processes. Typical examples of such processes are natural flows in crystalline rocks and industrial flows in oil and gas production systems or hydraulic fracturing. Fracture flow provides an important mechanism for transporting mass and energy. The distribution of the apertures of fracture and contact area are the key parameters with regard to the fracture transmissivity. We use the method of correlated random fields [Mourzenko, 1996] to generate synthetic fracture geometry in 3D. The flow of an incompressible Newtonian viscous fluid in geological formation can be approximated by the Stokes, the Stokes-Brinkman or the Reynolds models. We use our own implementation of the finite element method based on MILAMIN [Dabrowski, 2008] to solve governing partial differential equation over domain. We compare the Stokes, the Stokes-Brinkamn and the Reynolds models for fracture flow based on systematic numerical simulations for a wide range of geometric parameters. Mismatch between the Reynolds and the Stokes models becomes significant with increasing fracture roughness or contact area. The Stokes-Brinkman model is more accurate than Reynolds models due to additional Laplacian term, which allows to fulfil no-slip boundary condition. We present condition when the Reynolds and the Stokes-Brinkman models are valid. In the last three decades many authors used the Reynolds equation for studying fracture flow because of its simplicity. We recommend using the Stokes-Brinkman model for fracture flow, which allows to fulfil no-slip boundary condition on asperities boundary and is more accurate for rough fractures than the Reynolds model.

Models for assessing the relative phase velocity in a two-phase flow. Status report

International Nuclear Information System (INIS)

Schaffrath, A.; Ringel, H.

2000-06-01

The knowledge of slip or drift flux in two phase flow is necessary for several technical processes (e.g. two phase pressure losses, heat and mass transfer in steam generators and condensers, dwell period in chemical reactors, moderation effectiveness of two phase coolant in BWR). In the following the most important models for two phase flow with different phase velocities (e.g. slip or drift models, analogy between pressure loss and steam quality, ε - ε models and models for the calculation of void distribution in reposing fluids) are classified, described and worked up for a further comparison with own experimental data. (orig.)

A double parameters measurement of steam-water two-phase flow with single orifice

International Nuclear Information System (INIS)

Zhong Shuoping; Tong Yunxian; Yu Meiying

1992-08-01

A double parameters measurement of steam-water two-phase flow with single orifice is described. An on-line measurement device based on micro-computer has been developed. The measured r.m.s error of steam quality is less than 6.5% and the measured relative r.m.s. error of mass flow rate is less than 9%

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.

Modeling of two-phase slug flow

International Nuclear Information System (INIS)

Fabre, J.; Line, A.

1992-01-01

When gas and liquid flow in a pipe, over a range of flow rates, a flow pattern results in which sequences of long bubbles, almost filling the pipe cross section, are successively followed by liquid slugs that may contain small bubbles. This flow pattern, usually called slug flow, is encountered in numerous practical situations, such as in the production of hydrocarbons in wells and their transportation in pipelines; the production of steam and water in geothermal power plants; the boiling and condensation in liquid-vapor systems of thermal power plants; emergency core cooling of nuclear reactors; heat and mass transfer between gas and liquid in chemical reactors. This paper provides a review of two phase slug flow modeling

A novel drag force coefficient model for gas–water two-phase flows under different flow patterns

Energy Technology Data Exchange (ETDEWEB)

Shang, Zhi, E-mail: shangzhi@tsinghua.org.cn

2015-07-15

Graphical abstract: - Highlights: • A novel drag force coefficient model was established. • This model realized to cover different flow patterns for CFD. • Numerical simulations were performed under wide range flow regimes. • Validations were carried out through comparisons to experiments. - Abstract: A novel drag force coefficient model has been developed to study gas–water two-phase flows. In this drag force coefficient model, the terminal velocities were calculated through the revised drift flux model. The revised drift flux is different from the traditional drift flux model because the natural curve movement of the bubble was revised through considering the centrifugal force. Owing to the revisions, the revised drift flux model was to extend to 3D. Therefore it is suitable for CFD applications. In the revised drift flux model, the different flow patterns of the gas–water two-phase flows were able to be considered. This model innovatively realizes the drag force being able to cover different flow patterns of gas–water two-phase flows on bubbly flow, slug flow, churn flow, annular flow and mist flow. Through the comparisons of the numerical simulations to the experiments in vertical upward and downward pipe flows, this model was validated.

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

Experimental study on the convective heat transfer enhancement in single-phase steam flow by a support grid

International Nuclear Information System (INIS)

Kim, Byoung Jae; Kim, Kihwan; Kim, Dong-Eok; Youn, Young-Jung; Park, Jong-Kuk; Moon, Sang-Ki; Song, Chul-Hwa

2014-01-01

Highlights: • The convective heat transfer enhancement by support grids is investigated. • Experiments were performed in a square array 2 × 2 rod bundle. • The enhancement was affected not only by the blockage ratio also by the Reynolds number. • For low Reynolds numbers, the enhancement depends on the Reynolds number (Re). • For high Reynolds numbers, the enhancement is nearly independent of Re. - Abstract: Single-phase flow occurs in the fuel rod bundle of a pressurized water reactor, during the normal operation period or at the early stage of the reflood phase in a loss-of-coolant accident scenario. In the former period, the flow is single-phase water flow, but in the latter case, the flow is single-phase steam flow. Support grids are required to maintain a proper geometry configuration of fuel rods within nuclear fuel assemblies. This study was conducted to elucidate the effects of support grids on the convective heat transfer in single-phase steam flow. Experiments were made in a square array 2 × 2 rod bundle. The four electrically-heating rods were maintained by support grids with mixing vanes creating a swirl flow. Two types of support grids were considered in this study. The two types are geometrically similar except the blockage ratio by different mixing vane angles. For all test runs, 2 kW power was supplied to each rod. The working fluid was superheated steam with Re = 2,301–39,594. The axial profile of the rod surface temperatures was measured, and the convective heat transfer enhancement by the presence of the support grids was examined. The peak heat transfer enhancement was a function of not only the blockage ratio but also the Reynolds number. Given the same blockage ratio, the heat transfer enhancement was sensitive to the Reynolds number in laminar flow, whereas it was nearly independent of the Reynolds number in turbulent flow

Regularized lattice Boltzmann model for immiscible two-phase flows with power-law rheology

Science.gov (United States)

Ba, Yan; Wang, Ningning; Liu, Haihu; Li, Qiang; He, Guoqiang

2018-03-01

In this work, a regularized lattice Boltzmann color-gradient model is developed for the simulation of immiscible two-phase flows with power-law rheology. This model is as simple as the Bhatnagar-Gross-Krook (BGK) color-gradient model except that an additional regularization step is introduced prior to the collision step. In the regularization step, the pseudo-inverse method is adopted as an alternative solution for the nonequilibrium part of the total distribution function, and it can be easily extended to other discrete velocity models no matter whether a forcing term is considered or not. The obtained expressions for the nonequilibrium part are merely related to macroscopic variables and velocity gradients that can be evaluated locally. Several numerical examples, including the single-phase and two-phase layered power-law fluid flows between two parallel plates, and the droplet deformation and breakup in a simple shear flow, are conducted to test the capability and accuracy of the proposed color-gradient model. Results show that the present model is more stable and accurate than the BGK color-gradient model for power-law fluids with a wide range of power-law indices. Compared to its multiple-relaxation-time counterpart, the present model can increase the computing efficiency by around 15%, while keeping the same accuracy and stability. Also, the present model is found to be capable of reasonably predicting the critical capillary number of droplet breakup.

Effects of rolling on characteristics of single-phase water flow in narrow rectangular ducts

International Nuclear Information System (INIS)

Xing Dianchuan; Yan Changqi; Sun Licheng; Xu Chao

2012-01-01

Highlights: ► Mass flow rate and friction pressure drop with different pressure head are compared. ► The effect of pressure head on flow fluctuation is considered theoretically. ► Time-mean and real-time friction pressure drop in different rolling motion are studied. ► Rolling motion influences the fluctuation of friction pressure drop in two aspects. ► New correlation for frictional coefficient in rolling motion is achieved. - Abstract: Experimental and theoretical studies of rolling effects on characteristics of single-phase water flow in narrow rectangular ducts are performed under ambient temperature and pressure. Two types of pressure head are supplied by elevate water tank and pump respectively. The results show that the frictional pressure drop under rolling condition fluctuates periodically, with its amplitude decaying as mean Reynolds number increase and the rolling amplitude decrease, while the amplitude is nearly invariable with rolling period. Rolling motion influences the fluctuation amplitude of frictional pressure drop in two aspects, on the one hand, rolling reduced periodical pulsing flow leads to the fluctuation of the frictional pressure drop, on the other hand, additional force acting on fluid near the wall due to the rolling motion makes local frictional resistance oscillate periodically. The mass flow rate oscillates periodically in rolling motion with the pressure head supplied by water tank, while its fluctuation is so weak that could be neglected for the case of the pressure head supplied by pump. An empirical correlation for the frictional coefficient under rolling condition is achieved, and the experimental data is well correlated. A mathematical model is also developed to study the effect of pressure head on mass flow rate fluctuation in rolling motion. The fluctuation amplitude of the mass flow rate decreases rapidly with a higher pressure head. Comparing with the vertical condition, rolling motion nearly has no effects on

Subsurface Flow Modeling in Single and Dual Continuum Anisotropic Porous Media using the Multipoint Flux Approximation Method

KAUST Repository

Negara, Ardiansyah

2015-05-01

Anisotropy of hydraulic properties of the subsurface geologic formations is an essential feature that has been established as a consequence of the different geologic processes that undergo during the longer geologic time scale. With respect to subsurface reservoirs, in many cases, anisotropy plays significant role in dictating the direction of flow that becomes no longer dependent only on driving forces like the pressure gradient and gravity but also on the principal directions of anisotropy. Therefore, there has been a great deal of motivation to consider anisotropy into the subsurface flow and transport models. In this dissertation, we present subsurface flow modeling in single and dual continuum anisotropic porous media, which include the single-phase groundwater flow coupled with the solute transport in anisotropic porous media, the two-phase flow with gravity effect in anisotropic porous media, and the natural gas flow in anisotropic shale reservoirs. We have employed the multipoint flux approximation (MPFA) method to handle anisotropy in the flow model. The MPFA method is designed to provide correct discretization of the flow equations for general orientation of the principal directions of the permeability tensor. The implementation of MPFA method is combined with the experimenting pressure field approach, a newly developed technique that enables the solution of the global problem breaks down into the solution of multitude of local problems. The numerical results of the study demonstrate the significant effects of anisotropy of the subsurface formations. For the single-phase groundwater flow coupled with the solute transport modeling in anisotropic porous media, the results shows the strong impact of anisotropy on the pressure field and the migration of the solute concentration. For the two-phase flow modeling with gravity effect in anisotropic porous media, it is observed that the buoyancy-driven flow, which emerges due to the density differences between the

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

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

Sampling device for withdrawing a representative sample from single and multi-phase flows

Science.gov (United States)

Apley, Walter J.; Cliff, William C.; Creer, James M.

1984-01-01

A fluid stream sampling device has been developed for the purpose of obtaining a representative sample from a single or multi-phase fluid flow. This objective is carried out by means of a probe which may be inserted into the fluid stream. Individual samples are withdrawn from the fluid flow by sampling ports with particular spacings, and the sampling parts are coupled to various analytical systems for characterization of the physical, thermal, and chemical properties of the fluid flow as a whole and also individually.

A New Appraoch to Modeling Immiscible Two-phase Flow in Porous Media

DEFF Research Database (Denmark)

Yuan, Hao; Shapiro, Alexander; Stenby, Erling Halfdan

In this work we present a systematic literature review regarding the macroscopic approaches to modeling immiscible two-phase flow in porous media, the formulation process of the incorporate PDE based on Film Model(viscous coupling), the calculation of saturation profile around the transition zone...... to modeling immiscible two-phase flow in porous media. The suggested approach to immiscible two-phase flow in porous media describes the dispersed mesoscopic fluids’ interfaces which are highly influenced by the injected interfacial energy and the local interfacial energy capacity. It reveals a new...... possibility of modeling two-phase flow through energy balance. The saturation profile generated through the suggested approach is different from those through other approaches....

Mathematical models for two-phase stratified pipe flow

Energy Technology Data Exchange (ETDEWEB)

Biberg, Dag

2005-06-01

The simultaneous transport of oil, gas and water in a single multiphase flow pipe line has for economical and practical reasons become common practice in the gas and oil fields operated by the oil industry. The optimal design and safe operation of these pipe lines require reliable estimates of liquid inventory, pressure drop and flow regime. Computer simulations of multiphase pipe flow have thus become an important design tool for field developments. Computer simulations yielding on-line monitoring and look ahead predictions are invaluable in day-to-day field management. Inaccurate predictions may have large consequences. The accuracy and reliability of multiphase pipe flow models are thus important issues. Simulating events in large pipelines or pipeline systems is relatively computer intensive. Pipe-lines carrying e.g. gas and liquefied gas (condensate) may cover distances of several hundred km in which transient phenomena may go on for months. The evaluation times associated with contemporary 3-D CFD models are thus not compatible with field applications. Multiphase flow lines are therefore normally simulated using specially dedicated 1-D models. The closure relations of multiphase pipe flow models are mainly based on lab data. The maximum pipe inner diameter, pressure and temperature in a multiphase pipe flow lab is limited to approximately 0.3 m, 90 bar and 60{sup o}C respectively. The corresponding field values are, however, much higher i.e.: 1 m, 1000 bar and 200{sup o}C respectively. Lab data does thus not cover the actual field conditions. Field predictions are consequently frequently based on model extrapolation. Applying field data or establishing more advanced labs will not solve this problem. It is in fact not practically possible to acquire sufficient data to cover all aspects of multiphase pipe flow. The parameter range involved is simply too large. Liquid levels and pressure drop in three-phase flow are e.g. determined by 13 dimensionless parameters

Criterion for traffic phases in single vehicle data and empirical test of a microscopic three-phase traffic theory

International Nuclear Information System (INIS)

Kerner, Boris S; Klenov, Sergey L; Hiller, Andreas

2006-01-01

Based on empirical and numerical microscopic analyses, the physical nature of a qualitatively different behaviour of the wide moving jam phase in comparison with the synchronized flow phase-microscopic traffic flow interruption within the wide moving jam phase-is found. A microscopic criterion for distinguishing the synchronized flow and wide moving jam phases in single vehicle data measured at a single freeway location is presented. Based on this criterion, empirical microscopic classification of different local congested traffic states is performed. Simulations made show that the microscopic criterion and macroscopic spatiotemporal objective criteria lead to the same identification of the synchronized flow and wide moving jam phases in congested traffic. Microscopic models in the context of three-phase traffic theory have been tested based on the microscopic criterion for the phases in congested traffic. It is found that microscopic three-phase traffic models can explain both microscopic and macroscopic empirical congested pattern features. It is obtained that microscopic frequency distributions for vehicle speed difference as well as fundamental diagrams and speed correlation functions can depend on the spatial co-ordinate considerably. It turns out that microscopic optimal velocity (OV) functions and time headway distributions are not necessarily qualitatively different, even if local congested traffic states are qualitatively different. The reason for this is that important spatiotemporal features of congested traffic patterns are lost in these as well as in many other macroscopic and microscopic traffic characteristics, which are widely used as the empirical basis for a test of traffic flow models, specifically, cellular automata traffic flow models

Moving Boudary Models for Dynamic Simulations of Two-phase Flows

DEFF Research Database (Denmark)

Jensen, Jakob Munch; Tummelscheit, H.

2002-01-01

. The Dymola Modelica translator can automatically reduce the DAE index and thus makes efficient simulation possible. Usually the flow entering a dry-expansion evaporator in a refrigeration system is two-phase, and there is thus no liquid region. The general MB model has a number of special cases where only...... model is used. The overall robustness and the simplicity of the MB model, makes it well suited for open loop as well as closed loop simulations of two-phase flows. Simulation results for an evaporator in a refrigeration system are shown. The open loop system is simulated both with the reduced MB...... but is less complex. The reduced MB-model is well suited for control purposes both for determining control parameters and for model based control strategies and examples of a controlled refrigeration system are shown. The general MB model divides the flow into three regions (liquid, two-phase and vapor...

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

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)

Modeling and numerical analysis of non-equilibrium two-phase flows; Modelisation et analyse numerique des ecoulements diphasiques en desequilibre

Energy Technology Data Exchange (ETDEWEB)

Rascle, P.; El Amine, K. [Electricite de France (EDF), Direction des Etudes et Recherches, 92 - Clamart (France)

1997-12-31

We are interested in the numerical approximation of two-fluid models of nonequilibrium two-phase flows described by six balance equations. We introduce an original splitting technique of the system of equations. This technique is derived in a way such that single phase Riemann solvers may be used: moreover, it allows a straightforward extension to various and detailed exchange source terms. The properties of the fluids are first approached by state equations of ideal gas type and then extended to real fluids. For the construction of numerical schemes , the hyperbolicity of the full system is not necessary. When based on suitable kinetic unwind schemes, the algorithm can compute flow regimes evolving from mixture to single phase flows and vice versa. The whole scheme preserves the physical features of all the variables which remain in the set of physical states. Several stiff numerical tests, such as phase separation and phase transition are displayed in order to highlight the efficiency of the proposed method. The document is a PhD thesis divided in 6 chapters and two annexes. They are entitled: 1. - Introduction (in French), 2. - Two-phase flow, modelling and hyperbolicity (in French), 3. - A numerical method using upwind schemes for the resolution of two-phase flows without exchange terms (in English), 4. - A numerical scheme for one-phase flow of real fluids (in English), 5. - An upwind numerical for non-equilibrium two-phase flows (in English), 6. - The treatment of boundary conditions (in English), A.1. The Perthame scheme (in English) and A.2. The Roe scheme (in English). 136 refs. This document represents a PhD thesis in the speciality Applied Mathematics presented par Khalid El Amine to the Universite Paris 6.

Drag reduction by polymer addition in single and two-phase gas-liquid flows in pipelines

Energy Technology Data Exchange (ETDEWEB)

Bizotto, Vanessa Cristina; Paes, Diogo Melo; Franca, Fernando de Almeida [Universidade Estadual de Campinas, SP (Brazil). Centro de Estudos de Petroleo. LabPetro]. E-mails: vanessa@cepetro.unicamp.br; diogopaes10@hotmail.com; Sabadini, Edvaldo [Universidade Estadual de Campinas, SP (Brazil). Inst. de Quimica]. E-mails: sabadini@iqm.unicamp.br; ffranca@fem.unicamp.br

2008-07-01

Turbulence mechanisms, as the eddies formation frequency and size, promote energy dissipation that appears as pressure drop in pipe flows. Adding minute amounts of polymers - ppm - of high molecular weight to the solution can lead to the reduction of the viscous dissipation. The formed macromolecules interact with the eddies, cause the eddies coherence breakdown, damp the energy transport and reduces the pressure drop. This phenomenon is known as the hydrodynamic drag reduction (DR, for short). Thus, for a given pipe flow rate there is decrease in pressure head, which is a desired operating strategy when transporting liquids. Studies on the hydrodynamic drag reduction in polymeric systems have been carried out in collaboration by the Chemistry Institute and the Petroleum Laboratory - LabPetro, UNICAMP. These studies have allowed microscopic approaches to the engineering scales, tackling the most usual processes - single phase flows, as well as gas-liquid two-phase flows in pipelines, which are quite common in the chemical and the petroleum industries. Tests conducted in the Chemistry Institute comprised over-the-bench experimentations made with a rotational double-gap type rheometer. These quick performed tests used small amount of polymers, and provided information on the additive concentration, the drag reduction and the solution mechanical stability along a turbulent shearing process. The results indicated that 17% is the limiting drag reduction achieved when a 2 ppm aqueous solution of polyacrylamide - PAM - was tested. These tests, besides giving preliminary estimations, are limited in terms of engineering application due to the low shearing rates applied by the viscometer. The tests performed at LabPetro comprised pressure drop measurements in actual pipe flows, both water single and air-water two-phase flows, using the previous knowledge acquired with the viscometer tests. In the former case, the Prandtl-von Karman map has been drawn to show the %DR in terms

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

RETRAN nonequilibrium two-phase flow model for operational transient analyses

International Nuclear Information System (INIS)

Paulsen, M.P.; Hughes, E.D.

1982-01-01

The field balance equations, flow-field models, and equation of state for a nonequilibrium two-phase flow model for RETRAN are given. The differential field balance model equations are: (1) conservation of mixture mass; (2) conservation of vapor mass; (3) balance of mixture momentum; (4) a dynamic-slip model for the velocity difference; and (5) conservation of mixture energy. The equation of state is formulated such that the liquid phase may be subcooled, saturated, or superheated. The vapor phase is constrained to be at the saturation state. The dynamic-slip model includes wall-to-phase and interphase momentum exchanges. A mechanistic vapor generation model is used to describe vapor production under bulk subcooling conditions. The speed of sound for the mixture under nonequilibrium conditions is obtained from the equation of state formulation. The steady-state and transient solution methods are described

An eddy viscosity model for flow in a tube bundle

International Nuclear Information System (INIS)

Soussan, D.; Grandotto, M.

1998-01-01

The work described in this paper is part of the development of GENEPI a 3-dimensional finite element code, designed for the thermalhydraulic analysis of steam generators. It focuses on the implementation of two-phase flow turbulence-induced viscosity in a tube bundle. The GENEPI code, as other industrial codes, uses the eddy viscosity concept introduced by Boussinesq for single phase flow. The concept assumes that the turbulent momentum transfer is similar to the viscous shear stresses. Eddy viscosity formulation is reasonably well known for single phase flows, especially in simple geometries (i.e., in smooth tube, around a single body, or behind a row of bars/tubes), but there exists very little information on it for two-phase flows. An analogy between single and two-phases is used to set up a model for eddy viscosity. The eddy viscosity model examined in this paper is used for a tube bundle geometry and, therefore, is extended to include anisotropy to the classic model. Each of the main flow directions (cross flow inline, cross flow staggered, and parallel flows) gives rise to a specific eddy viscosity formula. The results from a parametric study indicate that the eddy viscosity in the staggered flow is roughly 1.5 times as large as that for the inline cross flow, 60 times as large as that for the parallel flow, and 105 as large as that for the molecular viscosity. Then, the different terms are combined with each other to result in a global eddy viscosity model for a steam generator tube bundle flow. (author)

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

Three layer model analysis on two-phase critical flow through a converging nozzle

International Nuclear Information System (INIS)

Ochi, J.; Ayukawa, K.

1991-01-01

A three layer model is proposed for a two-phase critical flow through a converging nozzle in this paper. Most previous analyses of the two phase flow have been based on a homogeneous or a separated flow model as the conservation equations. These results were found to have large deviations from the actual measurements for two phase critical flows. The presented model is based on the assumption that a flow consists of three layers with a mixing region between gas and liquid phase layers. The effect of gas and liquid fraction occupied in the mixing layer was made clear from the numerical results. The measurements of the critical flow rate and the pressure profiles through a converging nozzle were made with air-water flow. The calculated results of these models are discussed in comparison with the experimental data for the flow rates and the pressure distributions under critical conditions

Numerical simulation of multi-dimensional two-phase flow based on flux vector splitting

Energy Technology Data Exchange (ETDEWEB)

Staedtke, H.; Franchello, G.; Worth, B. [Joint Research Centre - Ispra Establishment (Italy)

1995-09-01

This paper describes a new approach to the numerical simulation of transient, multidimensional two-phase flow. The development is based on a fully hyperbolic two-fluid model of two-phase flow using separated conservation equations for the two phases. Features of the new model include the existence of real eigenvalues, and a complete set of independent eigenvectors which can be expressed algebraically in terms of the major dependent flow parameters. This facilitates the application of numerical techniques specifically developed for high speed single-phase gas flows which combine signal propagation along characteristic lines with the conservation property with respect to mass, momentum and energy. Advantages of the new model for the numerical simulation of one- and two- dimensional two-phase flow are discussed.

Mass and heat transfer at the outer surface of helical coils under single and two phase flow

International Nuclear Information System (INIS)

Abdel-Aziz, M.H.; Nirdosh, I.; Sedahmed, G.H.

2016-01-01

Highlights: • The work aims to develop reactors which need rapid temperature control. • Mass and heat transfer at the outer surface of helical coils was studied experimentally. • The experiments were conducted under gas sparing, single and two phase flow. • Variables were helical tube diameter, physical properties, and gas and liquid velocity. • Results verification in terms of natural convection and surface renewal mechanism was explained. - Abstract: The mass transfer behavior of the outer surface of vertical helical coil was studied by the electrochemical technique under single phase flow, gas sparging and two phase flow. Variables studied were helical tube diameter, physical properties of the solution, solution velocity and superficial gas velocity. The mass transfer data were correlated by dimensionless equations. Mass transfer enhancement ratio in case of two phase flow ranged from 1.1 to 4.9 compared to single phase flow. Implication of the results for the design and operation of helical coil reactors used to conduct L–S exothermic diffusion controlled reactions which need rapid temperature control were outlined. In this case the inner coil surface will act as a cooler while the outer surface will act a reaction surface. Immobilized enzyme catalyzed biochemical reactions where heat sensitive materials may be involved represent an example for the reactions which can employ the helical coil reactor. Also the importance of the results in the design of and operation of diffusion controlled membrane processes which employ helical coil membrane was noted. In view of the analogy between heat and mass transfer the possibility of using the results in the design and operation of helical coil heat exchangers was highlighted.

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

interThermalPhaseChangeFoam—A framework for two-phase flow simulations with thermally driven phase change

Directory of Open Access Journals (Sweden)

Mahdi Nabil

2016-01-01

Full Text Available The volume-of-fluid (VOF approach is a mature technique for simulating two-phase flows. However, VOF simulation of phase-change heat transfer is still in its infancy. Multiple closure formulations have been proposed in the literature, each suited to different applications. While these have enabled significant research advances, few implementations are publicly available, actively maintained, or inter-operable. Here, a VOF solver is presented (interThermalPhaseChangeFoam, which incorporates an extensible framework for phase-change heat transfer modeling, enabling simulation of diverse phenomena in a single environment. The solver employs object oriented OpenFOAM library features, including Run-Time-Type-Identification to enable rapid implementation and run-time selection of phase change and surface tension force models. The solver is packaged with multiple phase change and surface tension closure models, adapted and refined from earlier studies. This code has previously been applied to study wavy film condensation, Taylor flow evaporation, nucleate boiling, and dropwise condensation. Tutorial cases are provided for simulation of horizontal film condensation, smooth and wavy falling film condensation, nucleate boiling, and bubble condensation. Validation and grid sensitivity studies, interfacial transport models, effects of spurious currents from surface tension models, effects of artificial heat transfer due to numerical factors, and parallel scaling performance are described in detail in the Supplemental Material (see Appendix A. By incorporating the framework and demonstration cases into a single environment, users can rapidly apply the solver to study phase-change processes of interest.

interThermalPhaseChangeFoam-A framework for two-phase flow simulations with thermally driven phase change

Science.gov (United States)

Nabil, Mahdi; Rattner, Alexander S.

The volume-of-fluid (VOF) approach is a mature technique for simulating two-phase flows. However, VOF simulation of phase-change heat transfer is still in its infancy. Multiple closure formulations have been proposed in the literature, each suited to different applications. While these have enabled significant research advances, few implementations are publicly available, actively maintained, or inter-operable. Here, a VOF solver is presented (interThermalPhaseChangeFoam), which incorporates an extensible framework for phase-change heat transfer modeling, enabling simulation of diverse phenomena in a single environment. The solver employs object oriented OpenFOAM library features, including Run-Time-Type-Identification to enable rapid implementation and run-time selection of phase change and surface tension force models. The solver is packaged with multiple phase change and surface tension closure models, adapted and refined from earlier studies. This code has previously been applied to study wavy film condensation, Taylor flow evaporation, nucleate boiling, and dropwise condensation. Tutorial cases are provided for simulation of horizontal film condensation, smooth and wavy falling film condensation, nucleate boiling, and bubble condensation. Validation and grid sensitivity studies, interfacial transport models, effects of spurious currents from surface tension models, effects of artificial heat transfer due to numerical factors, and parallel scaling performance are described in detail in the Supplemental Material (see Appendix A). By incorporating the framework and demonstration cases into a single environment, users can rapidly apply the solver to study phase-change processes of interest.

Experimental study on local resistance of two-phase flow through spacer grid with rod bundle

International Nuclear Information System (INIS)

Yan Chaoxing; Yan Changqi; Sun Licheng; Tian Qiwei

2015-01-01

The experimental study on local resistance of single-phase and two-phase flows through a spacer grid in a vertical channel with 3 × 3 rod bundle was carried out under the normal temperature and pressure. For the case of single-phase flow, the liquid Reynolds number covered the range of 290-18 007. For the case of two-phase flow, the ranges of gas and liquid superficial velocities were 0.013-3.763 m/s and 0.076-1.792 m/s, respectively. A correlation for predicting local resistance of single-phase flow was given based on experimental results. Eight classical two-phase viscosity formulae for homogeneous model were evaluated against the experimental data of two-phase flow. The results show that Dukler model predicts the experimental data well in the range of Re 1 < 9000 while McAdams correlation is the best one for Re 1 ≥ 9000. For all experimental data, Dukler model provides the best prediction with the mean relative error of 29.03%. A new correlation is fitted for the range of Re 1 < 9000 by considering mass quality, two- phase Reynolds number and liquid and gas densities, resulting in a good agreement with the experimental data. (authors)

Turbulence-induced bubble collision force modeling and validation in adiabatic two-phase flow using CFD

Energy Technology Data Exchange (ETDEWEB)

Sharma, Subash L., E-mail: sharma55@purdue.edu [School of Nuclear Engineering, Purdue University, West Lafayette, IN 47907-1290 (United States); Hibiki, Takashi; Ishii, Mamoru [School of Nuclear Engineering, Purdue University, West Lafayette, IN 47907-1290 (United States); Brooks, Caleb S. [Department of Nuclear, Plasma, and Radiological Engineering, University of Illinois, Urbana, IL 61801 (United States); Schlegel, Joshua P. [Nuclear Engineering Program, Missouri University of Science and Technology, Rolla, MO 65409 (United States); Liu, Yang [Nuclear Engineering Program, Department of Mechanical Engineering, Virginia Polytechnic Institute and State University, Blacksburg, VA 24061 (United States); Buchanan, John R. [Bechtel Marine Propulsion Corporation, Bettis Laboratory, West Mifflin, PA 15122 (United States)

2017-02-15

Highlights: • Void distribution in narrow rectangular channel with various non-uniform inlet conditions. • Modeling of void diffusion due to bubble collision force. • Validation of new modeling in adiabatic air–water two-phase flow in a narrow channel. - Abstract: The prediction capability of the two-fluid model for gas–liquid dispersed two-phase flow depends on the accuracy of the closure relations for the interfacial forces. In previous studies of two-phase flow Computational Fluid Dynamics (CFD), interfacial force models for a single isolated bubble has been extended to disperse two-phase flow assuming the effect in a swarm of bubbles is similar. Limited studies have been performed investigating the effect of the bubble concentration on the lateral phase distribution. Bubbles, while moving through the liquid phase, may undergo turbulence-driven random collision with neighboring bubbles without significant coalescence. The rate of these collisions depends upon the bubble approach velocity and bubble spacing. The bubble collision frequency is expected to be higher in locations with higher bubble concentrations, i.e., volume fraction. This turbulence-driven random collision causes the diffusion of the bubbles from high concentration to low concentration. Based on experimental observations, a phenomenological model has been developed for a “turbulence-induced bubble collision force” for use in the two-fluid model. For testing the validity of the model, two-phase flow data measured at Purdue University are utilized. The geometry is a 10 mm × 200 mm cross section channel. Experimentally, non-uniform inlet boundary conditions are applied with different sparger combinations to vary the volume fraction distribution across the wider dimension. Examining uniform and non-uniform inlet data allows for the influence of the volume fraction to be studied as a separate effect. The turbulence-induced bubble collision force has been implemented in ANSYS CFX. The

High-velocity two-phase flow two-dimensional modeling

International Nuclear Information System (INIS)

Mathes, R.; Alemany, A.; Thilbault, J.P.

1995-01-01

The two-phase flow in the nozzle of a LMMHD (liquid metal magnetohydrodynamic) converter has been studied numerically and experimentally. A two-dimensional model for two-phase flow has been developed including the viscous terms (dragging and turbulence) and the interfacial mass, momentum and energy transfer between the phases. The numerical results were obtained by a finite volume method based on the SIMPLE algorithm. They have been verified by an experimental facility using air-water as a simulation pair and a phase Doppler particle analyzer for velocity and droplet size measurement. The numerical simulation of a lithium-cesium high-temperature pair showed that a nearly homogeneous and isothermal expansion of the two phases is possible with small pressure losses and high kinetic efficiencies. In the throat region a careful profiling is necessary to reduce the inertial effects on the liquid velocity field

Identification of two-phase flow regimes by time-series modeling

International Nuclear Information System (INIS)

King, C.H.; Ouyang, M.S.; Pei, B.S.

1987-01-01

The identification of two-phase flow patterns in pipes or ducts is important to the design and operation of thermal-hydraulic systems, especially in the nuclear reactor cores of boiling water reactors or in the steam generators of pressurized water reactors. Basically, two-phase flow shows some fluctuating characteristics even at steady-state conditions. These fluctuating characteristics can be analyzed by statistical methods for obtaining flow signatures. There have been a number of experimental studies conducted that are concerned with the statistical properties of void fraction or pressure pulsation in two-phase flow. In this study, the authors propose a new technique of identifying the patterns of air-water two-phase flow in a vertical pipe. This technique is based on analyzing the statistic characteristics of the pressure signals of the test loop by time-series modeling

A residual-based a posteriori error estimator for single-phase Darcy flow in fractured porous media

KAUST Repository

Chen, Huangxin

2016-12-09

In this paper we develop an a posteriori error estimator for a mixed finite element method for single-phase Darcy flow in a two-dimensional fractured porous media. The discrete fracture model is applied to model the fractures by one-dimensional fractures in a two-dimensional domain. We consider Raviart–Thomas mixed finite element method for the approximation of the coupled Darcy flows in the fractures and the surrounding porous media. We derive a robust residual-based a posteriori error estimator for the problem with non-intersecting fractures. The reliability and efficiency of the a posteriori error estimator are established for the error measured in an energy norm. Numerical results verifying the robustness of the proposed a posteriori error estimator are given. Moreover, our numerical results indicate that the a posteriori error estimator also works well for the problem with intersecting fractures.

Scaling of two-phase flow transients using reduced pressure system and simulant fluid

International Nuclear Information System (INIS)

Kocamustafaogullari, G.; Ishii, M.

1987-01-01

Scaling criteria for a natural circulation loop under single-phase flow conditions are derived. Based on these criteria, practical applications for designing a scaled-down model are considered. Particular emphasis is placed on scaling a test model at reduced pressure levels compared to a prototype and on fluid-to-fluid scaling. The large number of similarty groups which are to be matched between modell and prototype makes the design of a scale model a challenging tasks. The present study demonstrates a new approach to this clasical problen using two-phase flow scaling parameters. It indicates that a real time scaling is not a practical solution and a scaled-down model should have an accelerated (shortened) time scale. An important result is the proposed new scaling methodology for simulating pressure transients. It is obtained by considerung the changes of the fluid property groups which appear within the two-phase similarity parameters and the single-phase to two-phase flow transition prameters. Sample calculations are performed for modeling two-phase flow transients of a high pressure water system by a low-pressure water system or a Freon system. It is shown that modeling is possible for both cases for simulation pressure transients. However, simulation of phase change transitions is not possible by a reduced pressure water system without distortion in either power or time. (orig.)

Developing two-phase flow modelling concepts for rock fractures

Energy Technology Data Exchange (ETDEWEB)

Keto, V. (Fortum Nuclear Services Oy, Espoo (Finland))

2010-01-15

The Finnish nuclear waste disposal company, Posiva Oy, is planning an underground repository for spent nuclear fuel to be constructed on the island of Olkiluoto on the south-west coast of Finland. One element of the site investigations conducted at Olkiluoto is the excavation of the underground rock characterisation facility (ONKALO) that will be extended to the final disposal depth (approximately -400 m). The bedrock around the excavated tunnel volume is fully saturated with groundwater, which water commonly contains a mixture of dissolved gases. These gases remain dissolved due to the high hydrostatic pressure. During tunnel excavation work the natural hydrostatic pressure field is disturbed and the water pressure will decrease close to the atmospheric pressure in the immediate vicinity of the tunnel. During this pressure drop two-phase flow conditions (combined flow of both water and gas) may develop in the vicinity of the underground opening, as the dissolved gas is exsoluted under the low pressure (the term exsolution refers here to release of the dissolved gas molecules from the water phase into a separate gas phase). This report steers towards concept development for numerical two-phase flow modeling for fractured rock. The focus is on the description of gas phase formation process under disturbed hydraulic conditions by exsolution of dissolved gases from groundwater, and on understanding the effects of a possibly formed gas phase on groundwater flow conditions in rock fractures. A mathematical model of three mutually coupled nonlinear partial differential equations for two-phase flow is presented and corresponding constitutional relationships are introduced and discussed. Illustrative numerical simulations are performed in a simplified setting using COMSOL Multiphysics 3.5a - software package. Shortcomings and conceptual problems are discussed. (orig.)

Developing two-phase flow modelling concepts for rock fractures

International Nuclear Information System (INIS)

Keto, V.

2010-01-01

The Finnish nuclear waste disposal company, Posiva Oy, is planning an underground repository for spent nuclear fuel to be constructed on the island of Olkiluoto on the south-west coast of Finland. One element of the site investigations conducted at Olkiluoto is the excavation of the underground rock characterisation facility (ONKALO) that will be extended to the final disposal depth (approximately -400 m). The bedrock around the excavated tunnel volume is fully saturated with groundwater, which water commonly contains a mixture of dissolved gases. These gases remain dissolved due to the high hydrostatic pressure. During tunnel excavation work the natural hydrostatic pressure field is disturbed and the water pressure will decrease close to the atmospheric pressure in the immediate vicinity of the tunnel. During this pressure drop two-phase flow conditions (combined flow of both water and gas) may develop in the vicinity of the underground opening, as the dissolved gas is exsoluted under the low pressure (the term exsolution refers here to release of the dissolved gas molecules from the water phase into a separate gas phase). This report steers towards concept development for numerical two-phase flow modeling for fractured rock. The focus is on the description of gas phase formation process under disturbed hydraulic conditions by exsolution of dissolved gases from groundwater, and on understanding the effects of a possibly formed gas phase on groundwater flow conditions in rock fractures. A mathematical model of three mutually coupled nonlinear partial differential equations for two-phase flow is presented and corresponding constitutional relationships are introduced and discussed. Illustrative numerical simulations are performed in a simplified setting using COMSOL Multiphysics 3.5a - software package. Shortcomings and conceptual problems are discussed. (orig.)

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.

Experimental investigations of single and two-phase flow in a heated rod bundle

International Nuclear Information System (INIS)

Barthel, Frank; Franz, Ronald; Hampel, Uwe; Technische Univ. Dresden

2013-01-01

An experimental facility for the study of boiling flows in a 3 x 3 rod bundle geometry was setup. The bundle resembles in essential geometrical parts the geometry in a pressurized water reactor fuel element. The facility is operated with a refrigerant fluid. Beside standard instrumentation for temperature, pressure and flow rate we employed particle image velocimetry for single phase flow studies, gamma ray densitometry for integral gas fraction measurement sand ultrafast X-ray tomography for the study of the void dynamics in the cross-section. Moreover extensive thermo-instrumentation allows axial rod surface temperature measurements for the central heated rod. First results will be discussed in this article. (orig.)

Contribution to the verification and the validation of an unsteady two-phase flow model

International Nuclear Information System (INIS)

Liu, Yujie

2013-01-01

This thesis contributes to the verification and the validation of the Baer-Nunziato (BN) model, to model water hammer phenomena in industrial piping systems. It consists of two parts, the first is to model water hammer flows with the BN model in Eulerian representation and the second is to extend this model to the ALE (Arbitrary Lagrangian Eulerian) formalism so as to take into account fluid-structure interaction (FSI). To model water hammer flows, closure laws of the BN model concerning the interfacial/source terms and the equations of state (EOS) were first studied. Then the whole system was simulated with a fractional step method including two steps, one for the resolution of the convective part, the other for the source terms. For the convective part, the Rusanov scheme was first checked, and some stability problems have been observed. Thus a more stable fractional step scheme has been proposed and verified. Regarding the source terms, four non-instantaneous relaxation schemes which represent the return to equilibrium of pressure, the transfers of momentum, heat and mass were successively applied. These schemes have been extended to 'generalized Stiffened Gas' EOS in order to represent phase-change. After regaining some typical phenomena associated with water hammer flows, the BN model was validated with the Simpson experiment, a classical water hammer test case, and the Canon experience, a rapid decompression of fluid in a high pressure duct. Moreover, the model was compared with two homogeneous models on these two experiments. Finally, an ALE version of the BN model was implemented, and verified on a case of wave propagation in a 'single' phase flow and a two-phase flow in a flexible pipe. The variation of wave propagation speed due to the coupling between the fluid and the structure has been well retrieved. The validation was performed on an experiment which examines the response of a pipe filled with water, subjected to a violent pressure peak (140 bar

Single-phase sodium pump model for LMFBR thermal-hydraulic analysis

International Nuclear Information System (INIS)

Madni, I.K.; Cazzoli, E.G.; Agrawal, A.K.

1979-01-01

A single-phase, homologous pump model has been developed for simulation of safety-related transients in LMFBR systems. Pump characteristics are modeled by homologous head and torque relations encompassing all regimes of operation. These relations were derived from independent model test results with a centrifugal pump of specific speed equal to 35 (SI units) or 1800 (gpm units), and are used to analyze the steady-state and transient behavior of sodium pumps in a number of LMFBR plants. Characteristic coefficients for the polynomials in all operational regimes are provided in a tabular form. The speed and flow dependence of head is included through solutions of the impeller and coolant dynamic equations. Results show the model to yield excellent agreement with experimental data in sodium for the FFTF prototype pump, and with vendor calculations for the CRBR pump. A sample pipe rupture calculation is also performed to demonstrate the necessity for modeling the complete pump characteristics

A state-of-the-art report on two-phase critical flow modelling

Energy Technology Data Exchange (ETDEWEB)

Jung, Jae Joon; Jang, Won Pyo; Kim, Dong Soo [Korea Atomic Energy Research Institute, Taejon (Korea, Republic of)

1993-09-01

This report reviews and analyses two-phase, critical flow models. The purposes of the report are (1) to make a knowledge base for the full understanding and best-estimate of two-phase, critical flow, (2) to analyse the model development trend and to derive the direction of further studies. A wide range of critical flow models are reviewed. Each model, in general, predicts critical flow well only within specified conditions. The critical flow models of best-estimate codes are special process model included in the hydrodynamic model. The results of calculations depend on the nodalization, discharge coefficient, and other user`s options. The following topics are recommended for continuing studies: improvement of two-fluid model, development of multidimensional model, data base setup and model error evaluation, and generalization of discharge coefficients. 24 figs., 5 tabs., 80 refs. (Author).

A state-of-the-art report on two-phase critical flow modelling

International Nuclear Information System (INIS)

Jung, Jae Joon; Jang, Won Pyo; Kim, Dong Soo

1993-09-01

This report reviews and analyses two-phase, critical flow models. The purposes of the report are (1) to make a knowledge base for the full understanding and best-estimate of two-phase, critical flow, (2) to analyse the model development trend and to derive the direction of further studies. A wide range of critical flow models are reviewed. Each model, in general, predicts critical flow well only within specified conditions. The critical flow models of best-estimate codes are special process model included in the hydrodynamic model. The results of calculations depend on the nodalization, discharge coefficient, and other user's options. The following topics are recommended for continuing studies: improvement of two-fluid model, development of multidimensional model, data base setup and model error evaluation, and generalization of discharge coefficients. 24 figs., 5 tabs., 80 refs. (Author)

A critical flow model for the Cathena thermalhydraulic code

International Nuclear Information System (INIS)

Popov, N.K.; Hanna, B.N.

1990-01-01

The calculation of critical flow rate, e.g., of choked flow through a break, is required for simulating a loss of coolant transient in a reactor or reactor-like experimental facility. A model was developed to calculate the flow rate through the break for given geometrical parameters near the break and fluid parameters upstream of the break for ordinary water, as well as heavy water, with or without non- condensible gases. This model has been incorporated in the CATHENA, one-dimensional, two-fluid thermalhydraulic code. In the CATHENA code a standard staggered-mesh, finite-difference representation is used to solve the thermalhydraulic equations. This model compares the fluid mixture velocity, calculated using the CATHENA momentum equations, with a critical velocity. When the mixture velocity is smaller than the critical velocity, the flow is assumed to be subcritical, and the model remains passive. When the fluid mixture velocity is higher than the critical velocity, the model sets the fluid mixture velocity equal to the critical velocity. In this paper the critical velocity at a link (momentum cell) is first estimated separately for single-phase liquid, two- phase, or single-phase gas flow condition at the upstream node (mass/energy cell). In all three regimes non-condensible gas can be present in the flow. For single-phase liquid flow, the critical velocity is estimated using a Bernoulli- type of equation, the pressure at the link is estimated by the pressure undershoot method

Modeling and measurement of interfacial area concentration in two-phase flow

International Nuclear Information System (INIS)

Paranjape, Sidharth; Ishii, Mamoru; Hibiki, Takashi

2010-01-01

This paper presents experimental and modeling approaches in characterizing interfacial structures in gas-liquid two-phase flow. For the modeling of the interfacial structure characterization, the interfacial area transport equation proposed earlier has been studied to provide a dynamic and mechanistic prediction tool for two-phase flow analysis. A state-of-the-art four-sensor conductivity probe technique has been developed to obtain detailed local interfacial structure information in a wide range of flow regimes spanning from bubbly to churn-turbulent flows. Newly obtained interfacial area data in 8 x 8 rod-bundle test section are also presented. This paper also reviews available models of the interfacial area sink and source terms and existing databases. The interfacial area transport equation has been benchmarked using condensation bubbly flow data.

A Finite-Volume approach for compressible single- and two-phase flows in flexible pipelines with fluid-structure interaction

Science.gov (United States)

Daude, F.; Galon, P.

2018-06-01

A Finite-Volume scheme for the numerical computations of compressible single- and two-phase flows in flexible pipelines is proposed based on an approximate Godunov-type approach. The spatial discretization is here obtained using the HLLC scheme. In addition, the numerical treatment of abrupt changes in area and network including several pipelines connected at junctions is also considered. The proposed approach is based on the integral form of the governing equations making it possible to tackle general equations of state. A coupled approach for the resolution of fluid-structure interaction of compressible fluid flowing in flexible pipes is considered. The structural problem is solved using Euler-Bernoulli beam finite elements. The present Finite-Volume method is applied to ideal gas and two-phase steam-water based on the Homogeneous Equilibrium Model (HEM) in conjunction with a tabulated equation of state in order to demonstrate its ability to tackle general equations of state. The extensive application of the scheme for both shock tube and other transient flow problems demonstrates its capability to resolve such problems accurately and robustly. Finally, the proposed 1-D fluid-structure interaction model appears to be computationally efficient.

Two-phase flow and pressure drop in T-junctions with horizontal run and vertical branch

International Nuclear Information System (INIS)

Katsaounis, A.

1987-01-01

Visual observations of single- and two-phase dividing flow and pressure drop measurements were performed in T-junctions with horizontal run and vertical branch. Both tees used were geometrically similar, in a scale of 1:4. The measurements were performed for plug/slug and stratified flow pattern regime in horizontal tube. Based on the single-phase form-resistance pressure drop correlation of Gardel a corresponded calculation model was developed for the two-phase flow verified by the own measurements. (orig.) [de

Hydrodynamic prediction of multidimensional single- and two-phase flow in rod arrays. Progress report, January 1-December 31, 1983

International Nuclear Information System (INIS)

Ebeling-Koning, D.B.; Robinson, J.T.; Todreas, N.E.

1984-01-01

The objective of this research is to develop comprehensive constitutive models for multidimensional two-phase flow in rod arrays. The constitutive parameters are the solid-fluid flow resistance and the gas-liquid interfacial momentum exchange force. This report covers work in four areas: (1) a correlation for flow resistance across banks of tubes which is independent of rod arrangement has been developed. The correlation was developed from data from three rod arrangements covering a Reynolds number range (based on superficial velocity) of 1 to 40,000; (2) complete pressure drop data for water flows in the laminar region in crossflow and 45 0 inclined rod arrays were taken; (3) the development of a model for the interfacial momentum exchange force in bubbly flows has been completed. This model has been validated against single bubble velocity data in inclined rod arrays. The model has been cast in a form suitable for implementation to two-fluid computer codes; and (4) rise velocities of bubbles in 0 0 , 45 0 , and 90 0 inclined rod arrays have been measured. This data should prove useful for the development of a bubble drag coefficient model for rod arrays

Numerical methods for limit problems in two-phase flow models

International Nuclear Information System (INIS)

Cordier, F.

2011-01-01

Numerical difficulties are encountered during the simulation of two-phase flows. Two issues are studied in this thesis: the simulation of phase transitions on one hand, and the simulation of both compressible and incompressible flows in the other hand. Un asymptotic study has shown that the loss of hyperbolicity of the bi fluid model was responsible for the difficulties encountered by the Roe scheme during the simulation of phase transitions. Robust and accurate polynomial schemes have thus been developed. To tackle the occasional lack of positivity of the solution, a numerical treatment based on adaptive diffusion was proposed and allowed to simulate with accuracy the test-cases of a boiling channel with creation of vapor and a tee-junction with separation of the phases. In a second part, an all-speed scheme for compressible and incompressible flows have been proposed. This pressure-based semi-implicit asymptotic preserving scheme is conservative, solves an elliptic equation on the pressure, and has been designed for general equations of state. The scheme was first developed for the full Euler equations and then extended to the Navier-Stokes equations. The good behaviour of the scheme in both compressible and incompressible regimes have been investigated. An extension of the scheme to the two-phase mixture model was implemented and demonstrated the ability of the scheme to simulate two-phase flows with phase change and a water-steam equation of state. (author) [fr

Numerical approach of multi-field two-phase flow models in the OVAP code

International Nuclear Information System (INIS)

Anela Kumbaro

2005-01-01

Full text of publication follows: A significant progress has been made in modeling the complexity of vapor-liquid two-phase flow. Different three-dimensional models exist in order to simulate the evolution of parameters which characterize a two-phase model. These models can be classified into various groups depending on the inter-field coupling. A hierarchy of increasing physical complexity can be defined. The simplest group corresponds to the homogeneous mixture models where no interactions are taken into account. Another group is constituted by the two-fluid models employing physically important interfacial forces between two-phases, liquid, and water. The last group is multi-field modeling where inter-field couplings can be taken into account at different degrees, such as the MUltiple Size Group modeling [2], the consideration of separate equations for the transport and generation of mass and momentum for each field under the assumption of the same energy for all the fields of the same phase, and a full multi-field two-phase model [1]. The numerical approach of the general three-dimensional two-phase flow is by complexity of the phenomena a very challenging task; the ideal numerical method should be at the same time simple in order to apply to any model, from equilibrium to multi-field model and conservative in order to respect the fundamental conservation physical laws. The approximate Riemann solvers have the good properties of conservation of mass, momentum and energy balance and have been extended successfully to two-fluid models [3]- [5]. But, the up-winding of the flux is based on the Eigen-decomposition of the two-phase flow model and the computation of the Eigen-structure of a multi-field model can be a high cost procedure. Our contribution will present a short review of the above two-phase models, and show numerical results obtained for some of them with an approximate Riemann solver and with lower-complexity alternative numerical methods that do not

Mechanical disequilibria in two-phase flow models: approaches by relaxation and by a reduced model

International Nuclear Information System (INIS)

Labois, M.

2008-10-01

This thesis deals with hyperbolic models for the simulation of compressible two-phase flows, to find alternatives to the classical bi-fluid model. We first establish a hierarchy of two-phase flow models, obtained according to equilibrium hypothesis between the physical variables of each phase. The use of Chapman-Enskog expansions enables us to link the different existing models to each other. Moreover, models that take into account small physical unbalances are obtained by means of expansion to the order one. The second part of this thesis focuses on the simulation of flows featuring velocity unbalances and pressure balances, in two different ways. First, a two-velocity two-pressure model is used, where non-instantaneous velocity and pressure relaxations are applied so that a balancing of these variables is obtained. A new one-velocity one-pressure dissipative model is then proposed, where the arising of second-order terms enables us to take into account unbalances between the phase velocities. We develop a numerical method based on a fractional step approach for this model. (author)

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

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.

Lattice Boltzmann model for three-phase viscoelastic fluid flow

Science.gov (United States)

Xie, Chiyu; Lei, Wenhai; Wang, Moran

2018-02-01

A lattice Boltzmann (LB) framework is developed for simulation of three-phase viscoelastic fluid flows in complex geometries. This model is based on a Rothman-Keller type model for immiscible multiphase flows which ensures mass conservation of each component in porous media even for a high density ratio. To account for the viscoelastic effects, the Maxwell constitutive relation is correctly introduced into the momentum equation, which leads to a modified lattice Boltzmann evolution equation for Maxwell fluids by removing the normal but excess viscous term. Our simulation tests indicate that this excess viscous term may induce significant errors. After three benchmark cases, the displacement processes of oil by dispersed polymer are studied as a typical example of three-phase viscoelastic fluid flow. The results show that increasing either the polymer intrinsic viscosity or the elastic modulus will enhance the oil recovery.

The use of a diffuse interface model to estimate effective transport properties for two-phase flows in porous media

International Nuclear Information System (INIS)

Fichot, Floriana; Duval, Fabiena; Garcia, Aureliena; Belloni, Julien; Quintard, Michel

2005-01-01

'interface' of finite width. Within the framework of the second gradient theory, diffuse interface models can be derived for both single-component and binary fluids. The latter model is often referred to as the Cahn-Hilliard or Phase-Field model. In this paper, the Cahn-Hilliard model is used to simulate immiscible two-phase flows in a representative geometry of a particle bed. The details of the modelling can be found in Duval [1]. The results are used to establish correlations for the two-phase pressure drop and the wall heat transfer coefficients as a function of void fraction and phase velocities. The results are compared with available analytical or experimental data. References: [1] F. Duval, F. Fichot, and M. Quintard. Diffuse interface methods for two-phase flow in porous media. In Trends in Numerical and Physical Modelling for Industrial Multiphase Flows, Cargese (France), September 2001; [2] F. Duval, F. Fichot, and M. Quintard. A local thermal non-equilibrium model for two-phase flows with phase change in porous media. International Journal of Heat and Mass Transfer, 47(3):613-639, 2004. (authors)

Dividing phases in two-phase flow and modeling of interfacial drag

Energy Technology Data Exchange (ETDEWEB)

Narumo, T.; Rajamaeki, M. [VTT Energy (Finland)

1997-07-01

Different models intended to describe one-dimensional two-phase flow are considered in this paper. The following models are introduced: conventional six-equation model, conventional model equipped with terms taking into account nonuniform transverse velocity distribution of the phases, several virtual mass models and a model in which the momentum equations have been derived by using the principles of Separation of the Flow According to Velocity (SFAV). The dynamics of the models have been tested by comparing their characteristic velocities to each other and against experimental data. The results show that the SFAV-model makes a hyperbolic system and predicts the propagation velocities of disturbances with the same order of accuracy as the best tested virtual mass models. Furthermore, the momentum interaction terms for the SFAV-model are considered. These consist of the wall friction terms and the interfacial friction term. The authors model wall friction with two independent terms describing the effect of each fluid on the wall separately. In the steady state, a relationship between the slip velocity and friction coefficients can be derived. Hence, the friction coefficients for the SFAV-model can be calculated from existing correlations, viz. from a drift-flux correlation and a wall friction correlation. The friction model was tested by searching steady-state distributions in a partial BWR fuel channel and comparing the relaxed values with the drift-flux correlation, which agreed very well with each other. In addition, response of the flow to a sine-wave disturbance in the water inlet flux was calculated as function of frequency. The results of the models differed from each other already with frequency of order 5 Hz, while the time constant for the relaxation, obtained from steady-state distribution calculation, would have implied significant differences appear not until with frequency of order 50 Hz.

Dividing phases in two-phase flow and modeling of interfacial drag

International Nuclear Information System (INIS)

Narumo, T.; Rajamaeki, M.

1997-01-01

Different models intended to describe one-dimensional two-phase flow are considered in this paper. The following models are introduced: conventional six-equation model, conventional model equipped with terms taking into account nonuniform transverse velocity distribution of the phases, several virtual mass models and a model in which the momentum equations have been derived by using the principles of Separation of the Flow According to Velocity (SFAV). The dynamics of the models have been tested by comparing their characteristic velocities to each other and against experimental data. The results show that the SFAV-model makes a hyperbolic system and predicts the propagation velocities of disturbances with the same order of accuracy as the best tested virtual mass models. Furthermore, the momentum interaction terms for the SFAV-model are considered. These consist of the wall friction terms and the interfacial friction term. The authors model wall friction with two independent terms describing the effect of each fluid on the wall separately. In the steady state, a relationship between the slip velocity and friction coefficients can be derived. Hence, the friction coefficients for the SFAV-model can be calculated from existing correlations, viz. from a drift-flux correlation and a wall friction correlation. The friction model was tested by searching steady-state distributions in a partial BWR fuel channel and comparing the relaxed values with the drift-flux correlation, which agreed very well with each other. In addition, response of the flow to a sine-wave disturbance in the water inlet flux was calculated as function of frequency. The results of the models differed from each other already with frequency of order 5 Hz, while the time constant for the relaxation, obtained from steady-state distribution calculation, would have implied significant differences appear not until with frequency of order 50 Hz

Development of Sodium Two Phase Flow Model for Kalimer Core Analysis

International Nuclear Information System (INIS)

Chang, W.P.; Hahn, Dohee

2002-01-01

An algorithm for sodium boiling is developed in order to extend the applicability of SSC-K, which is a main system analysis code for the KALIMER (Korea Advanced LIquid Metal Reactor) conceptual design. As the capability of the current SSC-K version is limited to simulation of only a single-phase sodium flow, its applicable range should not be enough to assess the fuel integrity under some of HCDA (Hypothetical Core Disruptive Accident) initiating events where sodium boiling is anticipated. The two-phase flow model similar to that used for the light water system is known to be no more effective directly to liquid metal reactors, because the phenomena observed between two reactor coolant systems are definitely different. The developing algorithm is based on a multiple-bubble slug ejection model, which allows a finite number of bubbles in a channel at any time. The present work is a continuous effort following the former study to confirm a qualitative acceptance on the model. Since the model has been applied only to the active fuel region in the former study, a part of its qualification seems to have already been demonstrated. For its application to the whole KALIMER core channel, however, the model needs to be examined the applicability to the fuel regions other than the active fuel. The present study primarily focuses on that point. In a result, although the model may be improved in a sense through the present study over the previous modeling, a clear limitation is also confirmed with the validity of the model. The further development, therefore, is required for this model to achieve its goal by resolving such limitations. (authors)

Single-stage three-phase AC to DC conversion with isolation and Bi-directional power flow

NARCIS (Netherlands)

Vermulst, B.J.D.; Duarte, J.L.; Wijnands, C.G.E.; Lomonova, E.A.

2014-01-01

An approach for three-phase AC to DC conversion is proposed, which consists of a single-stage while offering galvanic isolation, soft-switching, bi-directional power flow and a significant reduction of inductive and capacitive energy storage. Two elements enable this approach, namely a neutral

Phase-sensitive flow cytometer

Energy Technology Data Exchange (ETDEWEB)

Steinkamp, J.A.

1992-12-31

This report describes phase-sensitive flow cytometer (FCM) which provides additional FCM capability to use the fluorescence lifetime of one or more fluorochromes bound to single cells to provide additional information regarding the cells. The resulting fluorescence emission can be resolved into individual fluorescence signals if two fluorochromes are present or can be converted directly to a decay lifetime from a single fluorochrome. The excitation light for the fluorochromes is modulated to produce an amplitude modulated fluorescence pulse as the fluorochrome is excited in the FCM. The modulation signal also forms a reference signal that is phase-shifted a selected amount for subsequent mixing with the output modulated fluorescence intensity signal in phase-sensitive detection circuitry. The output from the phase-sensitive circuitry is then an individual resolved fluorochrome signal or a single fluorochrome decay lifetime, depending on the applied phase shifts.

Modeling and Stability Assessment of Single-Phase Grid Synchronization Techniques

DEFF Research Database (Denmark)

Golestan, Saeed; Guerrero, Josep M.; Vasquez, Juan

2018-01-01

(GSTs) is of vital importance. This task is most often based on obtaining a linear time-invariant (LTI) model for the GST and applying standard stability tests to it. Another option is modeling and dynamics/stability assessment of GSTs in the linear time-periodic (LTP) framework, which has received...... a very little attention. In this letter, the procedure of deriving the LTP model for single-phase GSTs is first demonstrated. The accuracy of the LTP model in predicting the GST dynamic behavior and stability is then evaluated and compared with that of the LTI one. Two well-known single-phase GSTs, i...

Comparison of Experimental and Numerical Two-Phase Flows in a Porous Micro-Model

Directory of Open Access Journals (Sweden)

Dustin Crandall

2009-12-01

Full Text Available Characterizing two-phase flow in porous media is important to provide estimates of sweep efficiency in enhanced oil recovery and storage estimates in potential geological CO2 sequestration repositories. To further the current understanding of two-phase flow in porous media a micro-model of interconnected channels was designed and fabricated using stereolithography to experimentally study gas-liquid flows. This flowcell was created with a wide variability of throat dimensions to represent naturally occurring porous media. Low flow rate experiments of immiscible two-phase drainage were performed within this cell. Additionally, a computational model for analyzing two-phase flows in the same flowcell was developed and used to simulate conditions not possible with our laboratory settings. The computational model was first tested for the identical conditions used in the experimental studies, and was shown to be in good agreement with the experimentally determined fractal dimension of the invading gas structure, time until breakthrough, and fluid saturation. The numerical model was then used to study two-phase air-water flows in flowcells with the same geometry and different gas-liquid-solid contact angles. The percent saturation of air and the motion of the fluids through the cell were found to vary with changes in these parameters. Finally, to simulate flows expected during geologic carbon sequestration, the fluid properties and interface conditions were set to model the flow of CO2 into a brine-saturated porous medium at representative subsurface conditions. The CO2 flows were shown to have larger gas saturations than the previous air into water studies. Thus the accuracy of the computational model was supported by the flowcell experiments, and the computational model extended the laboratory results to conditions not possible with the apparatus used in the experiments.

Three-Phase Flow and Capillarity in Porous Media

Energy Technology Data Exchange (ETDEWEB)

Helset, H.M.

1996-06-01

Multiphase flow occurs in most production operations on the reservoir rock of an oil field. Three-phase flow occurs in many oil recovery processes, such as gas or water alternating gas injection into water flooded reservoirs, thermal oil recovery and surfactant flooding. In this doctoral thesis, three-phase flow in porous media is modeled using the method of characteristics and simple wave theory. The general mathematical methods are used to model different flow processes: secondary migration of hydrocarbons, gravity segregated three-phase flow, and displacement experiments for relative permeabilities. The main conclusions from applying the methods to the migration process through a water-saturated carrier bed to a partially permeable seal are: (1) the three-phase solutions are generally dominated by the very low gas viscosity and gas density, (2) countercurrent flow and interference between the gas and oil mobilities can lead to trapping in situations where the permeability contrast between the seal and the carrier bed is insufficient to trap oil alone, (3) zones in which pressures exceed those from hydrostatic gradients (overpressured zones) exist even when the permeability contrast between the seal and the carrier bed is modest. The discussion of gravity segregated flow identifies stable displacement fronts, which usually occur only for a single value of the gas-water injection ratio. The new method developed for interpreting displacement experiments, with capillary pressure included, is based on the travelling wave solution of the flow equations. 185 refs., 69 figs., 5 tabs.

Computational fluid dynamics simulations of single-phase flow in a filter-press flow reactor having a stack of three cells

International Nuclear Information System (INIS)

Sandoval, Miguel A.; Fuentes, Rosalba; Walsh, Frank C.; Nava, José L.; Ponce de León, Carlos

2016-01-01

Highlights: • Computational fluid dynamic simulations in a filter-press stack of three cells. • The fluid velocity was different in each cell due to local turbulence. • The upper cell link pipe of the filter press cell acts as a fluid mixer. • The fluid behaviour tends towards a continuous mixing flow pattern. • Close agreement between simulations and experimental data was achieved. - Abstract: Computational fluid dynamics (CFD) simulations were carried out for single-phase flow in a pre-pilot filter press flow reactor with a stack of three cells. Velocity profiles and streamlines were obtained by solving the Reynolds-Averaged Navier-Stokes (RANS) equations with a standard k − ε turbulence model. The flow behaviour shows the appearance of jet flow at the entrance to each cell. At lengths from 12 to 15 cm along the cells channels, a plug flow pattern is developed at all mean linear flow rates studied here, 1.2 ≤ u ≤ 2.1 cm s −1 . The magnitude of the velocity profiles in each cell was different, due to the turbulence generated by the change of flow direction in the last fluid manifold. Residence time distribution (RTD) simulations indicated that the fluid behaviour tends towards a continuous mixing flow pattern, owing to flow at the output of each cell across the upper cell link pipe, which acts as a mixer. Close agreement between simulations and experimental RTD was obtained.

Mechanism of falling water limitation in two-phase counter flow through single hole vertical channel

International Nuclear Information System (INIS)

Sudo, Yukio; Ohnuki, Akira

1983-01-01

In the safety evaluation at the time of loss coolant accident, which is a credible accident in LWRs, recently main effort has been concentrated to the optimum evaluation calculation, and the grasp of vapor-liquid two-phase flow phenomena has become important. As one of the important phenomena, there is the limitation of falling water in two-phase counter flow through a vertical channel. This phenomenon is divided into the limitation of falling water stored in an upper plenum to a core through an upper core-supporting plate and a tie plate at the time of reflooding, and the limitation of falling emergency core-cooling water in downcomer channels at the time of reflooding in PWRs, under the presence of rising steam flow. In both cases, the evaluation of the quantity of falling water is important, because it contributes directly to core cooling. In this research, in order to clarify the mechanism of limitation of falling water in two-phase vertical counter flow, first, two-phase flow of air-water system through a single-hole vertical channel was taken up, and the effect of main parameters was experimentally studied. At the same time, the theoretical investigation was performed, and the comparison with the experimental results obtained so far was carried out. The different mechanisms for short and long channels gave the good results. (Kako, I.)

Computational fluid dynamics modeling of two-phase flow in a BWR fuel assembly

International Nuclear Information System (INIS)

Andrey Ioilev; Maskhud Samigulin; Vasily Ustinenko; Simon Lo; Adrian Tentner

2005-01-01

Full text of publication follows: The goal of this project is to develop an advanced Computational Fluid Dynamics (CFD) computer code (CFD-BWR) that allows the detailed analysis of the two-phase flow and heat transfer phenomena in a Boiling Water Reactor (BWR) fuel bundle under various operating conditions. This code will include more fundamental physical models than the current generation of sub-channel codes and advanced numerical algorithms for improved computational accuracy, robustness, and speed. It is highly desirable to understand the detailed two-phase flow phenomena inside a BWR fuel bundle. These phenomena include coolant phase changes and multiple flow regimes which directly influence the coolant interaction with fuel assembly and, ultimately, the reactor performance. Traditionally, the best analysis tools for the analysis of two-phase flow phenomena inside the BWR fuel assembly have been the sub-channel codes. However, the resolution of these codes is still too coarse for analyzing the detailed intra-assembly flow patterns, such as flow around a spacer element. Recent progress in Computational Fluid Dynamics (CFD), coupled with the rapidly increasing computational power of massively parallel computers, shows promising potential for the fine-mesh, detailed simulation of fuel assembly two-phase flow phenomena. However, the phenomenological models available in the commercial CFD programs are not as advanced as those currently being used in the sub-channel codes used in the nuclear industry. In particular, there are no models currently available which are able to reliably predict the nature of the flow regimes, and use the appropriate sub-models for those flow regimes. The CFD-BWR code is being developed as a customized module built on the foundation of the commercial CFD Code STAR-CD which provides general two-phase flow modeling capabilities. The paper describes the model development strategy which has been adopted by the development team for the

Description of turbulent velocity and temperature fields of single phase flow through tight rod bundles

International Nuclear Information System (INIS)

Monir, C.

1991-02-01

A two-dimensional procedure, VANTACY-II, describing the turbulent velocity and temperature fields for single phase flow in tight lattices is presented and validated. The flow is assumed to be steady, incrompressible and hydraulic and thermal fully developed. First, the state of art of turbulent momentum and heat transport in tight lattices is documented. It is shown that there is a necessity for experimental investigations in the field of turbulent heat transport. The presented new procedure is based on the turbulence model VELASCO-TUBS by NEELEN. The numerical solution of the balance equations is done by the finite element method code VANTACY by KAISER. The validation of the new procedure VANTACY-II is done by comparing the numerically calculated data for the velocity and temperature fields and for natural mixing with the experimental data of SEALE. The comparison shows a good agreement of experimental and numerically computed data. The observed differences can be mainly attributed to the model of the turbulent PRANDTL number used in the new procedure. (orig.) [de

Entropy analysis on non-equilibrium two-phase flow models

International Nuclear Information System (INIS)

Karwat, H.; Ruan, Y.Q.

1995-01-01

A method of entropy analysis according to the second law of thermodynamics is proposed for the assessment of a class of practical non-equilibrium two-phase flow models. Entropy conditions are derived directly from a local instantaneous formulation for an arbitrary control volume of a structural two-phase fluid, which are finally expressed in terms of the averaged thermodynamic independent variables and their time derivatives as well as the boundary conditions for the volume. On the basis of a widely used thermal-hydraulic system code it is demonstrated with practical examples that entropy production rates in control volumes can be numerically quantified by using the data from the output data files. Entropy analysis using the proposed method is useful in identifying some potential problems in two-phase flow models and predictions as well as in studying the effects of some free parameters in closure relationships

Entropy analysis on non-equilibrium two-phase flow models

Energy Technology Data Exchange (ETDEWEB)

Karwat, H.; Ruan, Y.Q. [Technische Universitaet Muenchen, Garching (Germany)

1995-09-01

A method of entropy analysis according to the second law of thermodynamics is proposed for the assessment of a class of practical non-equilibrium two-phase flow models. Entropy conditions are derived directly from a local instantaneous formulation for an arbitrary control volume of a structural two-phase fluid, which are finally expressed in terms of the averaged thermodynamic independent variables and their time derivatives as well as the boundary conditions for the volume. On the basis of a widely used thermal-hydraulic system code it is demonstrated with practical examples that entropy production rates in control volumes can be numerically quantified by using the data from the output data files. Entropy analysis using the proposed method is useful in identifying some potential problems in two-phase flow models and predictions as well as in studying the effects of some free parameters in closure relationships.

Study on hydrodynamic crisis of two-phase flow

International Nuclear Information System (INIS)

Nigmatulin, B.I.; Ivandaev, A.I.

1977-01-01

The phenomenon of hydrodynamic crisis (locking) of a two-phase flow is investigated. A model of a disperseannular flow with an effective monodisperse nucleus is used for describing the motion of a mixture under near-critical conditions. Main differential equations of a flow in a channel are given; in particular, the differential laws of variation of the effective diameters of drops in the nucleus as a result of mass exchange between the mixture components are singled out. Questions of concretization of the model are discussed. The conditions for the attainment of the maximum rate of flow of the gas through the channel are studied, as well as the effect of the flow prehistory on the formation of critical conditions in the outlet cross-section

Study on hydrodynamic crisis of two-phase flow

Energy Technology Data Exchange (ETDEWEB)

Nigmatulin, B I; Ivandaev, A I [Moskovskij Gosudarstvennyj Univ. (USSR). Nauchno-Issledovatel' skij Inst. Mekhaniki

1977-01-01

The phenomenon of hydrodynamic crisis (locking) of a two-phase flow is investigated. A model of a disperse annular flow with an effective monodisperse nucleus is used for describing the motion of a mixture under near-critical conditions. Main differential equations of a flow in a channel are given; in particular, the differential laws of variation of the effective diameters of drops in the nucleus as a result of mass exchange between the mixture components are singled out. Questions of concretization of the model are discussed. The conditions for the attainment of the maximum rate of flow of the gas through the channel are studied, as well as the effect of the flow prehistory on the formation of critical conditions in the outlet cross-section.

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

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

Modeling of bubble coalescence and disintegration in confined upward two-phase flow

International Nuclear Information System (INIS)

Sun Xiaodong; Kim, Seungjin; Ishii, Mamoru; Beus, Stephen G.

2004-01-01

This paper presents the modeling of bubble interaction mechanisms in the two-group interfacial area transport equation (IATE) for confined gas-liquid two-phase flow. The transport equation is applicable to bubbly, cap-turbulent, and churn-turbulent flow regimes. In the two-group IATE, bubbles are categorized into two groups: spherical/distorted bubbles as Group 1 and cap/slug/churn-turbulent bubbles as Group 2. Thus, two sets of equations are used to describe the generation and destruction rates of bubble number density, void fraction, and interfacial area concentration for the two groups of bubbles due to bubble expansion and compression, coalescence and disintegration, and phase change. Five major bubble interaction mechanisms are identified for the gas-liquid two-phase flow of interest, and are analytically modeled as the source/sink terms for the transport equation in the confined flow. These models include both intra-group and inter-group bubble interactions

Central upwind scheme for a compressible two-phase flow model.

Science.gov (United States)

Ahmed, Munshoor; Saleem, M Rehan; Zia, Saqib; Qamar, Shamsul

2015-01-01

In this article, a compressible two-phase reduced five-equation flow model is numerically investigated. The model is non-conservative and the governing equations consist of two equations describing the conservation of mass, one for overall momentum and one for total energy. The fifth equation is the energy equation for one of the two phases and it includes source term on the right-hand side which represents the energy exchange between two fluids in the form of mechanical and thermodynamical work. For the numerical approximation of the model a high resolution central upwind scheme is implemented. This is a non-oscillatory upwind biased finite volume scheme which does not require a Riemann solver at each time step. Few numerical case studies of two-phase flows are presented. For validation and comparison, the same model is also solved by using kinetic flux-vector splitting (KFVS) and staggered central schemes. It was found that central upwind scheme produces comparable results to the KFVS scheme.

Central upwind scheme for a compressible two-phase flow model.

Directory of Open Access Journals (Sweden)

Munshoor Ahmed

Full Text Available In this article, a compressible two-phase reduced five-equation flow model is numerically investigated. The model is non-conservative and the governing equations consist of two equations describing the conservation of mass, one for overall momentum and one for total energy. The fifth equation is the energy equation for one of the two phases and it includes source term on the right-hand side which represents the energy exchange between two fluids in the form of mechanical and thermodynamical work. For the numerical approximation of the model a high resolution central upwind scheme is implemented. This is a non-oscillatory upwind biased finite volume scheme which does not require a Riemann solver at each time step. Few numerical case studies of two-phase flows are presented. For validation and comparison, the same model is also solved by using kinetic flux-vector splitting (KFVS and staggered central schemes. It was found that central upwind scheme produces comparable results to the KFVS scheme.

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.35single bubble rising under buoyancy force in viscous fluids for a wide range of Eötvös and Morton numbers. Numerical results are compared with theoretical predictions and experimental results, and satisfactory agreement is shown.

Statistical descriptions of polydisperse turbulent two-phase flows

Energy Technology Data Exchange (ETDEWEB)

Minier, Jean-Pierre, E-mail: jean-pierre.minier@edf.fr

2016-12-15

Disperse two-phase flows are flows containing two non-miscible phases where one phase is present as a set of discrete elements dispersed in the second one. These discrete elements, or ‘particles’, can be droplets, bubbles or solid particles having different sizes. This situation encompasses a wide range of phenomena, from nano-particles and colloids sensitive to the molecular fluctuations of the carrier fluid to inertia particles transported by the large-scale motions of turbulent flows and, depending on the phenomenon studied, a broad spectrum of approaches have been developed. The aim of the present article is to analyze statistical models of particles in turbulent flows by addressing this issue as the extension of the classical formulations operating at a molecular or meso-molecular level of description. It has a three-fold purpose: (1) to bring out the thread of continuity between models for discrete particles in turbulent flows (above the hydrodynamical level of description) and classical mesoscopic formulations of statistical physics (below the hydrodynamical level); (2) to reveal the specific challenges met by statistical models in turbulence; (3) to establish a methodology for modeling particle dynamics in random media with non-zero space and time correlations. The presentation is therefore centered on organizing the different approaches, establishing links and clarifying physical foundations. The analysis of disperse two-phase flow models is developed by discussing: first, approaches of classical statistical physics; then, by considering models for single-phase turbulent flows; and, finally, by addressing current formulations for discrete particles in turbulent flows. This brings out that particle-based models do not cease to exist above the hydrodynamical level and offer great interest when combined with proper stochastic formulations to account for the lack of equilibrium distributions and scale separation. In the course of this study, general

Preliminary Two-Phase Terry Turbine Nozzle Models for RCIC Off-Design Operation Conditions

Energy Technology Data Exchange (ETDEWEB)

Zhao, Haihua [Idaho National Lab. (INL), Idaho Falls, ID (United States); O' Brien, James [Idaho National Lab. (INL), Idaho Falls, ID (United States)

2017-06-12

This report presents the effort to extend the single-phase analytical Terry turbine model to cover two-phase off-design conditions. The work includes: (1) adding well-established two-phase choking models – the Isentropic Homogenous Equilibrium Model (IHEM) and Moody’s model, and (2) theoretical development and implementation of a two-phase nozzle expansion model. The two choking models provide bounding cases for the two-phase choking mass flow rate. The new two-phase Terry turbine model uses the choking models to calculate the mass flow rate, the critical pressure at the nozzle throat, and steam quality. In the divergent stage, we only consider the vapor phase with a similar model for the single-phase case by assuming that the liquid phase would slip along the wall with a much slower speed and will not contribute the impulse on the rotor. We also modify the stagnation conditions according to two-phase choking conditions at the throat and the cross-section areas for steam flow at the nozzle throat and at the nozzle exit. The new two-phase Terry turbine model was benchmarked with the same steam nozzle test as for the single-phase model. Better agreement with the experimental data is observed than from the single-phase model. We also repeated the Terry turbine nozzle benchmark work against the Sandia CFD simulation results with the two-phase model for the pure steam inlet nozzle case. The RCIC start-up tests were simulated and compared with the single-phase model. Similar results are obtained. Finally, we designed a new RCIC system test case to simulate the self-regulated Terry turbine behavior observed in Fukushima accidents. In this test, a period inlet condition for the steam quality varying from 1 to 0 is applied. For the high quality inlet period, the RCIC system behaves just like the normal operation condition with a high pump injection flow rate and a nominal steam release rate through the turbine, with the net addition of water to the primary system; for

Compositional modeling of three-phase flow with gravity using higher-order finite element methods

KAUST Repository

Moortgat, Joachim

2011-05-11

A wide range of applications in subsurface flow involve water, a nonaqueous phase liquid (NAPL) or oil, and a gas phase, such as air or CO2. The numerical simulation of such processes is computationally challenging and requires accurate compositional modeling of three-phase flow in porous media. In this work, we simulate for the first time three-phase compositional flow using higher-order finite element methods. Gravity poses complications in modeling multiphase processes because it drives countercurrent flow among phases. To resolve this issue, we propose a new method for the upwinding of three-phase mobilities. Numerical examples, related to enhanced oil recovery and carbon sequestration, are presented to illustrate the capabilities of the proposed algorithm. We pay special attention to challenges associated with gravitational instabilities and take into account compressibility and various phase behavior effects, including swelling, viscosity changes, and vaporization. We find that the proposed higher-order method can capture sharp solution discontinuities, yielding accurate predictions of phase boundaries arising in computational three-phase flow. This work sets the stage for a broad extension of the higher-order methods for numerical simulation of three-phase flow for complex geometries and processes.

The structure of single-phase turbulent flows through closely spaced rod arrays

International Nuclear Information System (INIS)

Hooper, J.D.; Rehme, K.

1983-02-01

The axial and azimuthal turbulence intensity in the rod gap region has been shown, for developed single-phase turbulent flow through parallel rod arrays, to strongly increase with decreasing rod spacing. Two array geometries are reported, one constructed from a rectangular cross-section duct containing four rods and spaced at five p/d or w/d ratios. The second test section, constructed from six rods set in a regular square-pitch array, represented the interior flow region of a large array. The mean axial velocity, wall shear stress variation and axial pressure distribution were measured, together with hot-wire anemometer measurements of the Reynolds stresses. No significant non-zero secondary flow components were detected, using techniques capable of resolving secondary flow velocities to 1% of the local axial velocity. For the lowest p/d ratio of 1.036, cross-correlation measurements showed the presence of an energetic periodic azimuthal turbulent velocity component, correlated over a significant part of the flow area. The negligible contribution of secondary flows to the axial momentum balance, and the large azimuthal turbulent velocity component in the rod gap area, suggest a different mechanism than Reynolds stress gradient driven secondary flows for the turbulent transport process in the rod gap. (orig.) [de

A void fraction model for annular two-phase flow

Energy Technology Data Exchange (ETDEWEB)

Tandon, T.N.; Gupta, C.P.; Varma, H.K.

1985-01-01

An analytical model has been developed for predicting void fraction in two-phase annular flow. In the analysis, the Lockhart-Martinelli method has been used to calculate two-phase frictional pressure drop and von Karman's universal velocity profile is used to represent the velocity distribution in the annular liquid film. Void fractions predicted by the proposed model are generally in good agreement with a available experimental data. This model appears to be as good as Smith's correlation and better than the Wallis and Zivi correlations for computing void fraction.

A dual-porosity model for two-phase flow in deforming porous media

Science.gov (United States)

Shu, Zhengying

Only recently has one realized the importance of the coupling of fluid flow with rock matrix deformations for accurately modeling many problems in petroleum, civil, environmental, geological and mining engineering. In the oil industry, problems such as reservoir compaction, ground subsidence, borehole stability and sanding need to be simulated using a coupled approach to make more precise predictions than when each process is considered to be independent of the other. Due to complications associated with multiple physical processes and mathematical representation of a multiphase now system in deformable fractured reservoirs, very few references, if any, are available in the literature. In this dissertation, an approach, which is based on the dual-porosity concept and takes into account rock deformations, is presented to derive rigorously a set of coupled differential equations governing the behavior of fractured porous media and two-phase fluid flow. The finite difference numerical method, as an alternative method for finite element, is applied to discretize the governing equations both in time and space domains. Throughout the derived set of equations, the fluid pressures and saturations as well as the solid displacements are considered as the primary unknowns. The model is tested against the case of single-phase flow in a 1-D consolidation problem for which analytical solutions are available. An example of coupled two-phase fluid flow and rock deformations for a scenario of a one-dimensional, fractured porous medium is also discussed. The numerical model and simulator, RFIA (Rock Fluid InterAction), developed in this dissertation can be a powerful tool to solve difficult problems not only in petroleum engineering such as ground subsidence, borehole stability and sand control, but also in civil engineering such as groundwater flow through fractured bedrock and in environmental engineering such as waste deposit concerns in fractured and unconsolidated formations

Detailed pressure drop measurements in single-and two-phase adiabatic air-water turbulent flows in realistic BWR fuel assembly geometry with spacer grids

International Nuclear Information System (INIS)

Caraghiaur, Diana; Frid, Wiktor; Tillmark, Nils

2004-01-01

In recent years, advance numerical simulation tools based on CFD methods have been increasingly used in various multi-phase flow applications. One of these is two-phase flow in fuel assemblies of Boiling Water Reactors. The important and often missing aspect of this development is validation of CFD codes against proper experimental data. The purpose of the current paper is to present detailed pressure measurements over a spacer grid in low pressure adiabatic single- and bubbly two-phase flow, which will be used to further develop a CFD code for BWR fuel bundle analysis. The experiments have been carried out in a n asymmetric 24-rod sub-bundle, representing one quarter of a Westinghouse SVEA-96 nuclear reactor fuel assembly. Single-phase flow measurements have been performed at superficial velocities between 0.90-4.50 m/s and in the two-phase flow, which was simulated by air-water mixture, measurements have been performed at void fractions ranging from 4 to 12% and liquid superficial velocity of 4.50 m/s. In order to increase the number of measuring points, five pressure taps were drilled in one of the rods, which was easily moved vertically by a traversing system, covering most of the points in axial direction. Any of the rods in the bundle could be substitute by the pressure sensing rod and the measurements were made for five pressure taps facing-angles. A detailed pressure distribution comparison between single- and two-phase flows for different sub-channel positions and different flow conditions was performed over one of the spacers. In addition, single-phase pressure drop measurements in the upper part of the test section comprising two spacer grids have been carried out. (author)

Fractional single-phase-lagging heat conduction model for describing anomalous diffusion

Directory of Open Access Journals (Sweden)

T.N. Mishra

2016-03-01

Full Text Available The fractional single-phase-lagging (FSPL heat conduction model is obtained by combining scalar time fractional conservation equation to the single-phase-lagging (SPL heat conduction model. Based on the FSPL heat conduction model, anomalous diffusion within a finite thin film is investigated. The effect of different parameters on solution has been observed and studied the asymptotic behavior of the FSPL model. The analytical solution is obtained using Laplace transform method. The whole analysis is presented in dimensionless form. Numerical examples of particular interest have been studied and discussed in details.

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.

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.

A model of single and two-phase flow (critical or not) through cracks

International Nuclear Information System (INIS)

Seynhaeve, J.M.; Giot, M.; Granger, S.; Pages, D.

1995-07-01

The leaks through steam-generator cracks are the subject of research carried out in cooperation between EDF and UCL. A model to predict the mass flow rate with inlet subcooling has been developed, validated and published. The model takes into account the persistence of some metastable liquid in the crack. The present paper improves and extends the model, by making it applicable to all kinds of conditions prevailing in the S.G. tubes: not only subcooled water, but also saturated water, steam-water mixtures, saturated dry steam or superheated steam. Therefore, the flow at the crack inlet is analyzed and appropriate methods to initialize the numerical integration of the flow equations along the crack are proposed. The extensions of the model are still in the process of validation. However, a sensitivity analysis of its results has been made and is presented. (author)

An experimental and numerical study of developed single phase axial turbulent flow in a smooth rod bundle

International Nuclear Information System (INIS)

Hooper, J.D.

1977-01-01

A combined experimental and numerical model of a turbulent single phase coolant, flowing axially along the fuel pins of a nuclear reactor, was developed. The experimental rig represented two interconnected subchannels of a square array at a pitch/diameter ratio of 1.193. Air was the working fluid, and measurements were made of the mean radial velocity profiles, wall shear stress variation, turbulence velocity spectra and intensities. The numerically predicted wall shear distribution and mean velocity profiles, obtained using an empirical two-dimensional mixing length and eddy diffusivity concept to represent fluid turbulence, showed good agreement with the experimental results. (Author)

The questions of liquid metal two-phase flow modelling in the FBR core channels

International Nuclear Information System (INIS)

Martsiniouk, D.Ye.; Sorokin, A.P.

2000-01-01

The two-fluid model representation for calculations of two-phase flow characteristics in the FBR fuel pin bundles with liquid metal cooling is presented and analysed. Two conservation equations systems of the mass, momentum and energy have been written for each phase. Components accounted the mass-, momentum- and heat transfer throughout the interface occur in the macro-field equations after the averaging procedure realisation. The pattern map and correlations for two-fluid model in vertical liquid metal flows are presented. The description of processes interphase mass- and heat exchange and interphase friction is determined by the two-phase flow regime. The opportunity of the liquid metal two-phase flow regime definition is analysed. (author)

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)

On the use of nuclear magnetic resonance to characterize vertical two-phase bubbly flows

International Nuclear Information System (INIS)

Lemonnier, H.; Jullien, P.

2011-01-01

Research highlights: → We provide a complete theory of the PGSE measurement in single and two-phase flow. → Friction velocity can be directly determinated from measured velocity distributions. → Fast determination of moments shorten PGSE process with small loss of accuracy. → Turbulent diffusion measurements agree well with known trends and existing models. → We think NMR can be a tool to benchmark thermal anemometry in two-phase flow. - Abstract: Since the pioneering work of who showed that NMR can be used to measure accurately the mean liquid velocity and void fraction in two-phase pipe flow, it has been shown that NMR signal can also characterize the turbulent eddy diffusivity and velocity fluctuations. In this paper we provide an in depth validation of these statements together with a clarification of the nature of the mean velocity that is actually measured by NMR PFGSE sequence. The analysis shows that the velocity gradient at the wall is finely space-resolved and allows the determination of the friction velocity in single-phase flows. Next turbulent diffusion measurements in two-phase flows are presented, analyzed and compared to existing data and models. It is believed that NMR velocity measurement is sufficiently understood that it can be utilized to benchmark thermal anemometry in two-phase flows. Theoretical results presented in this paper also show how this can be undertaken.

Void fraction and flow regime determination by optical probe for boiling two-phase flow in a tube subchannel

International Nuclear Information System (INIS)

Cheng Huiping; Wu Hongtao; Ba Changxi; Yan Xiaoming; Huang Suyi

1995-12-01

In view of the need to determine void fraction and flow regime of vapor-liquid two-phase flow in the steam generator test model, domestic made optical probe was applied on a small-scale freon two-phase flow test rig. Optical probe signals were collected at a sampling rate up to 500 Hz and converted into digital form. Both the time signal, and the amplitude probability density function and FFT spectrum function calculated thereof were analysed in the time and frequency domains respectively. The threshold characterizing vapor or liquid contact with the probe tip was determined from the air-water two-phase flow pressure drop test results. Then, the boiling freon two-phase flow void fraction was determined by single threshold method, and compared with numerical heat transfer computation. Typical patterns which were revealed by the above-mentioned time signal and the functions were found corresponding to distinct flow regimes, as corroborated by visual observation. The experiment shows that the optical probe was a promising technique for two-phase flow void fraction measurement and flow regime identification (3 refs., 15 figs., 1 tab.)

A Schur complement method for compressible two-phase flow models

International Nuclear Information System (INIS)

Dao, Thu-Huyen; Ndjinga, Michael; Magoules, Frederic

2014-01-01

In this paper, we will report our recent efforts to apply a Schur complement method for nonlinear hyperbolic problems. We use the finite volume method and an implicit version of the Roe approximate Riemann solver. With the interface variable introduced in [4] in the context of single phase flows, we are able to simulate two-fluid models ([12]) with various schemes such as upwind, centered or Rusanov. Moreover, we introduce a scaling strategy to improve the condition number of both the interface system and the local systems. Numerical results for the isentropic two-fluid model and the compressible Navier-Stokes equations in various 2D and 3D configurations and various schemes show that our method is robust and efficient. The scaling strategy considerably reduces the number of GMRES iterations in both interface system and local system resolutions. Comparisons of performances with classical distributed computing with up to 218 processors are also reported. (authors)

Numerical simulation of the two-phase flows in a hydraulic coupling by solving VOF model

International Nuclear Information System (INIS)

Luo, Y; Zuo, Z G; Liu, S H; Fan, H G; Zhuge, W L

2013-01-01

The flow in a partially filled hydraulic coupling is essentially a gas-liquid two-phase flow, in which the distribution of two phases has significant influence on its characteristics. The interfaces between the air and the liquid, and the circulating flows inside the hydraulic coupling can be simulated by solving the VOF two-phase model. In this paper, PISO algorithm and RNG k–ε turbulence model were employed to simulate the phase distribution and the flow field in a hydraulic coupling with 80% liquid fill. The results indicate that the flow forms a circulating movement on the torus section with decreasing speed ratio. In the pump impeller, the air phase mostly accumulates on the suction side of the blades, while liquid on the pressure side; in turbine runner, air locates in the middle of the flow passage. Flow separations appear near the blades and the enclosing boundaries of the hydraulic coupling

An Eulerian two-phase model for steady sheet flow using large-eddy simulation methodology

Science.gov (United States)

Cheng, Zhen; Hsu, Tian-Jian; Chauchat, Julien

2018-01-01

A three-dimensional Eulerian two-phase flow model for sediment transport in sheet flow conditions is presented. To resolve turbulence and turbulence-sediment interactions, the large-eddy simulation approach is adopted. Specifically, a dynamic Smagorinsky closure is used for the subgrid fluid and sediment stresses, while the subgrid contribution to the drag force is included using a drift velocity model with a similar dynamic procedure. The contribution of sediment stresses due to intergranular interactions is modeled by the kinetic theory of granular flow at low to intermediate sediment concentration, while at high sediment concentration of enduring contact, a phenomenological closure for particle pressure and frictional viscosity is used. The model is validated with a comprehensive high-resolution dataset of unidirectional steady sheet flow (Revil-Baudard et al., 2015, Journal of Fluid Mechanics, 767, 1-30). At a particle Stokes number of about 10, simulation results indicate a reduced von Kármán coefficient of κ ≈ 0.215 obtained from the fluid velocity profile. A fluid turbulence kinetic energy budget analysis further indicates that the drag-induced turbulence dissipation rate is significant in the sheet flow layer, while in the dilute transport layer, the pressure work plays a similar role as the buoyancy dissipation, which is typically used in the single-phase stratified flow formulation. The present model also reproduces the sheet layer thickness and mobile bed roughness similar to measured data. However, the resulting mobile bed roughness is more than two times larger than that predicted by the empirical formulae. Further analysis suggests that through intermittent turbulent motions near the bed, the resolved sediment Reynolds stress plays a major role in the enhancement of mobile bed roughness. Our analysis on near-bed intermittency also suggests that the turbulent ejection motions are highly correlated with the upward sediment suspension flux, while

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

The difficult challenge of a two-phase CFD modelling for all flow regimes

International Nuclear Information System (INIS)

Bestion, D.

2014-01-01

Highlights: • The theoretical difficulties for modelling all flow regimes at CFD scale are identified. • The choice of the number of fields and of the time and space averaging or filtering are discussed and clarified. • Closure issues related to an all flow regime CFD model are listed and the main difficulties are identified. - Abstract: System thermalhydraulic codes model all two-phase flow regimes but they are limited to a macroscopic description. Two-phase CFD tools predict two-phase flow with a much finer space resolution but the current modelling capabilities are limited to dispersed bubbly or droplet flow and separate-phase flow. Much less experience exists on more complex flow regimes which combine the existence of dispersed fields with the presence of large interfaces such as a free surface or a film surface. A list of possible reactor issues which might benefit from an “all flow regime CFD model” is given. The first difficulty is to identify the various types of local flow configuration. It is shown that a 4-field model has much better capabilities than a two-fluid approach to identify most complex regimes. Then the choice between time averaging, space averaging, or even ensemble averaging is discussed. It is shown that only the RANS-2-fluid and a space-filtered 4-field model may be reasonably envisaged. The latter has the capabilities to identify all types of interfaces and should be privileged if a good accuracy is expected or if time fluctuations in intermittent flow have to be predicted while the former may be used when a high accuracy is not necessary and if time fluctuations in intermittent flow are not of interest. Finally the closure issue is presented including wall transfers, interfacial transfers, mass transfers between dispersed and continuous fields, and turbulent transfers. An important effort is required to model all interactions between sub-filter phenomena and the transfers from the sub-filter domain to the simulated domain. The

Virtual mass effects in two-phase flow. Topical report

International Nuclear Information System (INIS)

Cheng, L.Y.; Drew, D.A.; Lahey, R.T. Jr.

1978-03-01

The effect of virtual mass on phase separation during the acceleration of a two-phase mixture was studied. Virtual mass can be regarded as an induced inertia on the dispersed phase which is accelerating relative to the continuous phase, and it was found that the virtual mass acceleration is objective, implying an invariance with respect to reference frame. An objective form of the virtual acceleration was derived and required parameters were determined for limiting cases. Analyses determined that experiments on single bubble nozzle/diffuser flow cannot readily discriminate between various virtual mass acceleration models

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.

Numerical Simulations of Two-Phase Flow in a Self-Aerated Flotation Machine and Kinetics Modeling

KAUST Repository

Fayed, Hassan E.; Ragab, Saad

2015-01-01

A new boundary condition treatment has been devised for two-phase flow numerical simulations in a self-aerated minerals flotation machine and applied to a Wemco 0.8 m3 pilot cell. Airflow rate is not specified a priori but is predicted by the simulations as well as power consumption. Time-dependent simulations of two-phase flow in flotation machines are essential to understanding flow behavior and physics in self-aerated machines such as the Wemco machines. In this paper, simulations have been conducted for three different uniform bubble sizes (db = 0.5, 0.7 and 1.0 mm) to study the effects of bubble size on air holdup and hydrodynamics in Wemco pilot cells. Moreover, a computational fluid dynamics (CFD)-based flotation model has been developed to predict the pulp recovery rate of minerals from a flotation cell for different bubble sizes, different particle sizes and particle size distribution. The model uses a first-order rate equation, where models for probabilities of collision, adhesion and stabilization and collisions frequency estimated by Zaitchik-2010 model are used for the calculation of rate constant. Spatial distributions of dissipation rate and air volume fraction (also called void fraction) determined by the two-phase simulations are the input for the flotation kinetics model. The average pulp recovery rate has been calculated locally for different uniform bubble and particle diameters. The CFD-based flotation kinetics model is also used to predict pulp recovery rate in the presence of particle size distribution. Particle number density pdf and the data generated for single particle size are used to compute the recovery rate for a specific mean particle diameter. Our computational model gives a figure of merit for the recovery rate of a flotation machine, and as such can be used to assess incremental design improvements as well as design of new machines.

Numerical Simulations of Two-Phase Flow in a Self-Aerated Flotation Machine and Kinetics Modeling

Directory of Open Access Journals (Sweden)

Hassan Fayed

2015-03-01

Full Text Available A new boundary condition treatment has been devised for two-phase flow numerical simulations in a self-aerated minerals flotation machine and applied to a Wemco 0.8 m3 pilot cell. Airflow rate is not specified a priori but is predicted by the simulations as well as power consumption. Time-dependent simulations of two-phase flow in flotation machines are essential to understanding flow behavior and physics in self-aerated machines such as the Wemco machines. In this paper, simulations have been conducted for three different uniform bubble sizes (db = 0.5, 0.7 and 1.0 mm to study the effects of bubble size on air holdup and hydrodynamics in Wemco pilot cells. Moreover, a computational fluid dynamics (CFD-based flotation model has been developed to predict the pulp recovery rate of minerals from a flotation cell for different bubble sizes, different particle sizes and particle size distribution. The model uses a first-order rate equation, where models for probabilities of collision, adhesion and stabilization and collisions frequency estimated by Zaitchik-2010 model are used for the calculation of rate constant. Spatial distributions of dissipation rate and air volume fraction (also called void fraction determined by the two-phase simulations are the input for the flotation kinetics model. The average pulp recovery rate has been calculated locally for different uniform bubble and particle diameters. The CFD-based flotation kinetics model is also used to predict pulp recovery rate in the presence of particle size distribution. Particle number density pdf and the data generated for single particle size are used to compute the recovery rate for a specific mean particle diameter. Our computational model gives a figure of merit for the recovery rate of a flotation machine, and as such can be used to assess incremental design improvements as well as design of new machines.

Numerical Simulations of Two-Phase Flow in a Self-Aerated Flotation Machine and Kinetics Modeling

KAUST Repository

Fayed, Hassan E.

2015-03-30

A new boundary condition treatment has been devised for two-phase flow numerical simulations in a self-aerated minerals flotation machine and applied to a Wemco 0.8 m3 pilot cell. Airflow rate is not specified a priori but is predicted by the simulations as well as power consumption. Time-dependent simulations of two-phase flow in flotation machines are essential to understanding flow behavior and physics in self-aerated machines such as the Wemco machines. In this paper, simulations have been conducted for three different uniform bubble sizes (db = 0.5, 0.7 and 1.0 mm) to study the effects of bubble size on air holdup and hydrodynamics in Wemco pilot cells. Moreover, a computational fluid dynamics (CFD)-based flotation model has been developed to predict the pulp recovery rate of minerals from a flotation cell for different bubble sizes, different particle sizes and particle size distribution. The model uses a first-order rate equation, where models for probabilities of collision, adhesion and stabilization and collisions frequency estimated by Zaitchik-2010 model are used for the calculation of rate constant. Spatial distributions of dissipation rate and air volume fraction (also called void fraction) determined by the two-phase simulations are the input for the flotation kinetics model. The average pulp recovery rate has been calculated locally for different uniform bubble and particle diameters. The CFD-based flotation kinetics model is also used to predict pulp recovery rate in the presence of particle size distribution. Particle number density pdf and the data generated for single particle size are used to compute the recovery rate for a specific mean particle diameter. Our computational model gives a figure of merit for the recovery rate of a flotation machine, and as such can be used to assess incremental design improvements as well as design of new machines.

Stationary two-phase flow evaluation by the dynamic slip model

International Nuclear Information System (INIS)

Stevanovic, D.

1986-01-01

The equations which describe a dynamic slip model for stationary conditions are given in the paper. The basic solving procedure by the code DVOF4 is briefly described. The results are verified on the experiment FRIGG 313014. besides the void fraction and the vapor and liquid phase temperatures, the following parameters are plotted and explained: vapor phase generation rate, vapor and liquid phase velocities, slip between the phases, interfacial surface, friction drag between each phase and the wall, two-phase flow friction multiplier and pressure drop along the channel. (author)

Two-phase air-water stratified flow measurement using ultrasonic techniques

International Nuclear Information System (INIS)

Fan, Shiwei; Yan, Tinghu; Yeung, Hoi

2014-01-01

In this paper, a time resolved ultrasound system was developed for investigating two-phase air-water stratified flow. The hardware of the system includes a pulsed wave transducer, a pulser/receiver, and a digital oscilloscope. The time domain cross correlation method is used to calculate the velocity profile along ultrasonic beam. The system is able to provide velocities with spatial resolution of around 1mm and the temporal resolution of 200μs. Experiments were carried out on single phase water flow and two-phase air-water stratified flow. For single phase water flow, the flow rates from ultrasound system were compared with those from electromagnetic flow (EM) meter, which showed good agreement. Then, the experiments were conducted on two-phase air-water stratified flow and the results were given. Compared with liquid height measurement from conductance probe, it indicated that the measured velocities were explainable

Even distribution/dividing of single-phase fluids by symmetric bifurcation of flow channels

International Nuclear Information System (INIS)

Liu, Hong; Li, Peiwen

2013-01-01

Highlights: ► We addressed an issue of distributing a flow to a number of flow channels uniformly. ► The flow distribution is accomplished through bifurcation of channels. ► Some key parameters to the flow distribution uniformity have been identified. ► Flow uniformity was studied for several versions of flow distributor designs. ► A novel fluid packaging device of high efficiency was provided. -- Abstract: This study addresses a fundamental issue of distributing a single-phase fluid flow into a number of flow channels uniformly. A basic mechanism of flow distribution is accomplished through bifurcation of channels that symmetrically split one flow channel into two downstream channels. Applying the basic mechanism, cascades flow distributions are designed to split one flow into a large number of downstream flows uniformly. Some key parameters decisive to the flow distribution uniformity in such a system have been identified, and the flow distribution uniformity of air was studied for several versions of flow distributor designs using CFD analysis. The effect of the key parameters of the flow channel designs to the flow distribution uniformity was investigated. As an example of industrial application, a novel fluid packaging device of high efficiency was proposed and some CFD analysis results for the device were provided. The optimized flow distributor makes a very good uniform flow distribution which will significantly improve the efficiency of fluid packaging. The technology is expected to be of great significance to many industrial devices that require high uniformity of flow distribution

Contribution to the study of critical flow rates in a water-vapour two-phase flow

International Nuclear Information System (INIS)

Reocreux, Michel

1974-01-01

This research thesis aims at studying and analysing mechanisms involved in critical flows by adopting a theoretical and an experimental approach. After a recall of previous theoretical results and a discussion of their comparison with experimental results, the author outlines the main problems: the flow representation by a realistic model which takes all factors on which depend critical flow as well as many non critical flows into account, and the formulation of conditions to be met for a flow to be critical. Then, after a recall of the properties of critical single-phase flows, the author proposes an equation system. In the next part, he reports the development of an equation system for two-phase flows. The properties of the solutions of this system are studied to establish the general conditions required for a flow described by this system to be critical. These results are then applied to the equation system describing two-phase flows, and results are interpreted and discussed. In a second part, the author reports the experimental study by addressing experimental devices which could well produce the studied phenomenon, instrumentation and measurements, and the presentation and analysis of results [fr

Nonlinear Dynamic Model of Power Plants with Single-Phase Coolant Reactors

International Nuclear Information System (INIS)

Vollmer, H.

1968-12-01

The traditional way of developing dynamic models for a specific nuclear power plant and for specific purpose seems rather uneconomical, as much of the information often can not be utilized if the plant design or the required accuracy of the calculation is desired to be changed. It is therefore suggested that the model development may be made more systematic, general and flexible by - applying the 'box of bricks' system, where the main components of a nuclear power plant are treated separately and combined afterwards according to a given flow scheme, - a dynamic determination of the components which is as general as possible without taking into account those details which have a minor influence on the overall dynamics, - providing approximations of the more rigorous solution sufficient to meet the user s requirements on accuracy, - proper use of computers. A dynamic model for single-phase coolant reactor plants is established along these lines. By separation of the nonlinear and linear parts of the system, application of Laplace transformation and proper approximations, and the use of a hybrid computer it seems possible to determine the (nonlinear) dynamic behaviour of such a plant for perturbations which are not so large that phase changes of physical parameters occur, e. g. fuel does not melt. The model is applied to a steam cooled fast reactor power plant

Nonlinear Dynamic Model of Power Plants with Single-Phase Coolant Reactors

Energy Technology Data Exchange (ETDEWEB)

Vollmer, H

1968-12-15

The traditional way of developing dynamic models for a specific nuclear power plant and for specific purpose seems rather uneconomical, as much of the information often can not be utilized if the plant design or the required accuracy of the calculation is desired to be changed. It is therefore suggested that the model development may be made more systematic, general and flexible by - applying the 'box of bricks' system, where the main components of a nuclear power plant are treated separately and combined afterwards according to a given flow scheme, - a dynamic determination of the components which is as general as possible without taking into account those details which have a minor influence on the overall dynamics, - providing approximations of the more rigorous solution sufficient to meet the user s requirements on accuracy, - proper use of computers. A dynamic model for single-phase coolant reactor plants is established along these lines. By separation of the nonlinear and linear parts of the system, application of Laplace transformation and proper approximations, and the use of a hybrid computer it seems possible to determine the (nonlinear) dynamic behaviour of such a plant for perturbations which are not so large that phase changes of physical parameters occur, e. g. fuel does not melt. The model is applied to a steam cooled fast reactor power plant.

Flow visualization study of two-phase flow in a single bend outlet feeder pipe of a CANDU reactor

International Nuclear Information System (INIS)

Savalaxs, S.-A.; Lister, D.H.; Steward, F.R.

2005-01-01

In CANDU reactors, the feeder piping that is used to direct the high-temperature water coolant between the fuel channels and the steam generators is made of carbon steel. Since 1996, several CANDU stations have reported excessive corrosion of their outlet feeders. The first metre is particularity vulnerable because the piping there consists of single or double bends, which have relatively thin walls produced by the bending process. Early studies related the attack to the hydrodynamics of the coolant and verified that it was a type of flow-accelerated corrosion. In order to understand the hydrodynamics of the coolant in the outlet feeders by flow visualization, a full-scale transparent test section simulating the geometry and orientation of an outlet feeder bend with its upstream components was fabricated. The feeder consisted of a 54 mm diameter acrylic pipe with a 73 degree bend. This was connected to the upstream component with an acrylic simulation of a Grayloc flanged fitting. A test loop supplied room temperature water to the test section at flow rates up to 0.019 m3/s. Air could be injected into the water to give a mean volume fraction of up to 0.56. In this preliminary investigation, the size and velocity of air bubbles at different flow conditions and their distribution within the pipe bend were studied. Particular attention was paid to the flow pattern at the inside of the bend, where a CFD (computational fluid dynamics) code - Fluent 6.1-had failed to predict a liquid film in an earlier study. A high-speed digital video camera was used to determine the relation between bubble size and velocity. Such a relation should help to explain the discrepancy in the CFD modelling and provide the basis for accurate predictions of phase distribution in complex geometries at high flow rates. (authors)

Modeling and Control of a Single-Phase Marine Cooling System

DEFF Research Database (Denmark)

Hansen, Michael; Stoustrup, Jakob; Bendtsen, Jan Dimon

2013-01-01

This paper presents two model-based control design approaches for a single-phase marine cooling system. Models are derived from first principles and aim at describing significant system dynamics including nonlinearities and transport delays, while keeping the model complexity low. The two...

Interfacial shear modeling in two-phase annular flow

International Nuclear Information System (INIS)

Kumar, R.; Edwards, D.P.

1996-11-01

A new interfacial shear stress model called the law of the interface model, based on the law of the wall approach in turbulent flows, has been developed and locally applied in a fully developed, adiabatic, two-phase annular flow in a duct. Numerical results have been obtained using this model in conjunction with other models available in the literature that are required for the closure of the continuity and momentum equations. These results have been compared with droplet velocity data (using laser Doppler velocimetry and hot film anemometry), void fraction data (using gamma densitometry) and pressure drop data obtained in a R-134A refrigerant test facility. Droplet velocity results match the experimental data well, however, the prediction of the void fraction is less accurate. The poor prediction of void fraction, especially for the low void fraction cases, appears to be due to the lack of a good mechanistic model for entrainment

Interfacial shear modeling in two-phase annular flow

International Nuclear Information System (INIS)

Kumar, R.; Edwards, D.P.

1996-07-01

A new interfacial shear stress model called the law of the interface model, based on the law of the wall approach in turbulent flows, has been developed and locally applied in a fully developed, adiabatic, two-phase annular flow in a duct. Numerical results have been obtained using this model in conjunction with other models available in the literature that are required for the closure of the continuity and momentum equations. These results have been compared with droplet velocity data (using laser Doppler velocimetry and hot film anemometry), void fraction data (using gamma densitometry) and pressure drop data obtained in a R-134A refrigerant test facility. Droplet velocity results match the experimental data well, however, the prediction of the void fraction is less accurate. The poor prediction of void fraction, especially for the low void fraction cases, appears to be due to the lack of a good mechanistic model for entrainment

Numerical modeling of flow boiling instabilities using TRACE

International Nuclear Information System (INIS)

Kommer, Eric M.

2015-01-01

Highlights: • TRACE was used to realistically model boiling instabilities in single and parallel channel configurations. • Model parameters were chosen to exactly mimic other author’s work in order to provide for direct comparison of results. • Flow stability maps generated by the model show unstable flow at operating points similar to other authors. • The method of adjudicating when a flow is “unstable” is critical in this type of numerical study. - Abstract: Dynamic flow instabilities in two-phase systems are a vitally important area of study due to their effects on a great number of industrial applications, including heat exchangers in nuclear power plants. Several next generation nuclear reactor designs incorporate once through steam generators which will exhibit boiling flow instabilities if not properly designed or when operated outside design limits. A number of numerical thermal hydraulic codes attempt to model instabilities for initial design and for use in accident analysis. TRACE, the Nuclear Regulatory Commission’s newest thermal hydraulic code is used in this study to investigate flow instabilities in both single and dual parallel channel configurations. The model parameters are selected as to replicate other investigators’ experimental and numerical work in order to provide easy comparison. Particular attention is paid to the similarities between analysis using TRACE Version 5.0 and RELAP5/MOD3.3. Comparison of results is accomplished via flow stability maps non-dimensionalized via the phase change and subcooling numbers. Results of this study show that TRACE does indeed model two phase flow instabilities, with the transient response closely mimicking that seen in experimental studies. When compared to flow stability maps generated using RELAP, TRACE shows similar results with differences likely due to the somewhat qualitative criteria used by various authors to determine when the flow is truly unstable

A Simplified Micromechanical Modeling Approach to Predict the Tensile Flow Curve Behavior of Dual-Phase Steels

Science.gov (United States)

Nanda, Tarun; Kumar, B. Ravi; Singh, Vishal

2017-11-01

Micromechanical modeling is used to predict material's tensile flow curve behavior based on microstructural characteristics. This research develops a simplified micromechanical modeling approach for predicting flow curve behavior of dual-phase steels. The existing literature reports on two broad approaches for determining tensile flow curve of these steels. The modeling approach developed in this work attempts to overcome specific limitations of the existing two approaches. This approach combines dislocation-based strain-hardening method with rule of mixtures. In the first step of modeling, `dislocation-based strain-hardening method' was employed to predict tensile behavior of individual phases of ferrite and martensite. In the second step, the individual flow curves were combined using `rule of mixtures,' to obtain the composite dual-phase flow behavior. To check accuracy of proposed model, four distinct dual-phase microstructures comprising of different ferrite grain size, martensite fraction, and carbon content in martensite were processed by annealing experiments. The true stress-strain curves for various microstructures were predicted with the newly developed micromechanical model. The results of micromechanical model matched closely with those of actual tensile tests. Thus, this micromechanical modeling approach can be used to predict and optimize the tensile flow behavior of dual-phase steels.

Simultaneous thermal and optical imaging of two-phase flow in a micro-model.

Science.gov (United States)

Karadimitriou, N K; Nuske, P; Kleingeld, P J; Hassanizadeh, S M; Helmig, R

2014-07-21

In the study of non-equilibrium heat transfer in multiphase flow in porous media, parameters and constitutive relations, like heat transfer coefficients between phases, are unknown. In order to study the temperature development of a relatively hot invading immiscible non-wetting fluid and, ultimately, approximate heat transfer coefficients, a transparent micro-model is used as an artificial porous medium. In the last few decades, micro-models have become popular experimental tools for two-phase flow studies. In this work, the design of an innovative, elongated, PDMS (polydimethylsiloxane) micro-model with dimensions of 14.4 × 39 mm(2) and a constant depth of 100 microns is described. A novel setup for simultaneous thermal and optical imaging of flow through the micro-model is presented. This is the first time that a closed flow cell like a micro-model is used in simultaneous thermal and optical flow imaging. The micro-model is visualized by a novel setup that allowed us to monitor and record the distribution of fluids throughout the length of the micro-model continuously and also record the thermal signature of the fluids. Dynamic drainage and imbibition experiments were conducted in order to obtain information about the heat exchange between the phases. In this paper the setup as well as analysis and qualitative results are presented.

Accurate prediction of complex free surface flow around a high speed craft using a single-phase level set method

Science.gov (United States)

Broglia, Riccardo; Durante, Danilo

2017-11-01

This paper focuses on the analysis of a challenging free surface flow problem involving a surface vessel moving at high speeds, or planing. The investigation is performed using a general purpose high Reynolds free surface solver developed at CNR-INSEAN. The methodology is based on a second order finite volume discretization of the unsteady Reynolds-averaged Navier-Stokes equations (Di Mascio et al. in A second order Godunov—type scheme for naval hydrodynamics, Kluwer Academic/Plenum Publishers, Dordrecht, pp 253-261, 2001; Proceedings of 16th international offshore and polar engineering conference, San Francisco, CA, USA, 2006; J Mar Sci Technol 14:19-29, 2009); air/water interface dynamics is accurately modeled by a non standard level set approach (Di Mascio et al. in Comput Fluids 36(5):868-886, 2007a), known as the single-phase level set method. In this algorithm the governing equations are solved only in the water phase, whereas the numerical domain in the air phase is used for a suitable extension of the fluid dynamic variables. The level set function is used to track the free surface evolution; dynamic boundary conditions are enforced directly on the interface. This approach allows to accurately predict the evolution of the free surface even in the presence of violent breaking waves phenomena, maintaining the interface sharp, without any need to smear out the fluid properties across the two phases. This paper is aimed at the prediction of the complex free-surface flow field generated by a deep-V planing boat at medium and high Froude numbers (from 0.6 up to 1.2). In the present work, the planing hull is treated as a two-degree-of-freedom rigid object. Flow field is characterized by the presence of thin water sheets, several energetic breaking waves and plungings. The computational results include convergence of the trim angle, sinkage and resistance under grid refinement; high-quality experimental data are used for the purposes of validation, allowing to

Pressure drop and heat transfer characteristics for single-phase developing flow of water in rectangular microchannels

International Nuclear Information System (INIS)

Mirmanto; Kenning, D B R; Lewis, J S; Karayiannis, T G

2012-01-01

Experiments were conducted to investigate the pressure drop and heat transfer characteristics of single-phase flow of de-ionized water in single copper microchannels of hydraulic diameters 0.438 mm, 0.561 mm and 0.635 mm. The channel length was 62 mm. The experimental conditions covered a range of mass flux from 500 to 5000 kg/m 2 s in the laminar, transitional and low Reynolds number turbulent regimes. Pressure drop was measured for adiabatic flows with fluid inlet temperatures of 30°C, 60°C and 90°C. In the heat transfer tests, the heat flux ranged from 256 kW/m 2 to 519 kW/m 2 . Friction factors and Nusselt numbers determined from the measurements were higher than for fully-developed conditions, but in reasonable agreement with predictions made using published solutions for hydrodynamically and thermally developing flow. When entrance effects, experimental uncertainties, heat losses, inlet and exit losses, thermal boundary conditions and departure from laminar flow were considered, the results indicate that equations developed for flow and heat transfer in conventional size channels are applicable for water flows in microchannels of these sizes.

R 12 two-phase flow in throttle capillaries in critical flow conditions

International Nuclear Information System (INIS)

Petry, G.

1983-01-01

In this dissertation, the state of knowledge on two phase flow, its use and measurement processes are given from an extensive search of the literature. In the experimental part of the work, a continuously working experimental circuit was built up, by which single component two phase flow can be examined in critical flow conditions. Using the maintenance equations, a system of equations was produced, by which the content of steam flow, the content of steam volume and the slip between the phases at the end corssection of the capillary can be determined. The transfer of the experimental results into the Baker diagram shows that the experimental values lie in the region of mist, bubble and foam flow. (orig.) [de

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.

Cahn-Hilliard modeling of particles suspended in two-phase flows

NARCIS (Netherlands)

Choi, Y.J.; Anderson, P.D.

2012-01-01

In this paper, we present a model for the dynamics of particles suspended in two-phase flows by coupling the CahnHilliard theory with the extended finite element method (XFEM). In the CahnHilliard model the interface is considered to have a small but finite thickness, which circumvents explicit

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)

Governing equations for a seriated continuum: an unequal velocity model for two-phase flow

International Nuclear Information System (INIS)

Solbrig, C.W.; Hughes, E.D.

1975-05-01

The description of the flow of two-phase fluids is important in many engineering devices. Unexpected transient conditions which occur in these devices cannot, in general, be treated with single-component momentum equations. Instead, the use of momentum equations for each phase is necessary in order to describe the varied transient situations which can occur. These transient conditions can include phases moving in the opposite directions, such as steam moving upward and liquid moving downward, as well as phases moving in the same direction. The derivation of continuity and momentum equations for each phase and an overall energy equation for the mixture are presented. Terms describing interphase forces are described. A seriated (series of) continuum is distinguished from an interpenetrating medium by the representation of interphase friction with velocity differences in the former and velocity gradients in the latter. The seriated continuum also considers imbedded stationary solid surfaces such as occur in nuclear reactor cores. These stationary surfaces are taken into account with source terms. Sufficient constitutive equations are presented to form a complete set of equations. Methods are presented to show that all these coefficients are determinable from microscopic models and well known experimental results. Comparison of the present deviation with previous work is also given. The equations derived here may also be employed in certain multiphase, multicomponent flow applications. (U.S.)

Modelling a single phase voltage controlled rectifier using Laplace transforms

Science.gov (United States)

Kraft, L. Alan; Kankam, M. David

1992-01-01

The development of a 20 kHz, AC power system by NASA for large space projects has spurred a need to develop models for the equipment which will be used on these single phase systems. To date, models for the AC source (i.e., inverters) have been developed. It is the intent of this paper to develop a method to model the single phase voltage controlled rectifiers which will be attached to the AC power grid as an interface for connected loads. A modified version of EPRI's HARMFLO program is used as the shell for these models. The results obtained from the model developed in this paper are quite adequate for the analysis of problems such as voltage resonance. The unique technique presented in this paper uses the Laplace transforms to determine the harmonic content of the load current of the rectifier rather than a curve fitting technique. Laplace transforms yield the coefficient of the differential equations which model the line current to the rectifier directly.

Phase-field-based lattice Boltzmann modeling of large-density-ratio two-phase flows

Science.gov (United States)

Liang, Hong; Xu, Jiangrong; Chen, Jiangxing; Wang, Huili; Chai, Zhenhua; Shi, Baochang

2018-03-01

In this paper, we present a simple and accurate lattice Boltzmann (LB) model for immiscible two-phase flows, which is able to deal with large density contrasts. This model utilizes two LB equations, one of which is used to solve the conservative Allen-Cahn equation, and the other is adopted to solve the incompressible Navier-Stokes equations. A forcing distribution function is elaborately designed in the LB equation for the Navier-Stokes equations, which make it much simpler than the existing LB models. In addition, the proposed model can achieve superior numerical accuracy compared with previous Allen-Cahn type of LB models. Several benchmark two-phase problems, including static droplet, layered Poiseuille flow, and spinodal decomposition are simulated to validate the present LB model. It is found that the present model can achieve relatively small spurious velocity in the LB community, and the obtained numerical results also show good agreement with the analytical solutions or some available results. Lastly, we use the present model to investigate the droplet impact on a thin liquid film with a large density ratio of 1000 and the Reynolds number ranging from 20 to 500. The fascinating phenomena of droplet splashing is successfully reproduced by the present model and the numerically predicted spreading radius exhibits to obey the power law reported in the literature.

Modeling and numerical study of two phase flow

International Nuclear Information System (INIS)

Champmartin, A.

2011-01-01

This thesis describes the modelization and the simulation of two-phase systems composed of droplets moving in a gas. The two phases interact with each other and the type of model to consider directly depends on the type of simulations targeted. In the first part, the two phases are considered as fluid and are described using a mixture model with a drift relation (to be able to follow the relative velocity between the two phases and take into account two velocities), the two-phase flows are assumed at the equilibrium in temperature and pressure. This part of the manuscript consists of the derivation of the equations, writing a numerical scheme associated with this set of equations, a study of this scheme and simulations. A mathematical study of this model (hyperbolicity in a simplified framework, linear stability analysis of the system around a steady state) was conducted in a frame where the gas is assumed baro-tropic. The second part is devoted to the modelization of the effect of inelastic collisions on the particles when the time of the simulation is shorter and the droplets can no longer be seen as a fluid. We introduce a model of inelastic collisions for droplets in a spray, leading to a specific Boltzmann kernel. Then, we build caricatures of this kernel of BGK type, in which the behavior of the first moments of the solution of the Boltzmann equation (that is mass, momentum, directional temperatures, variance of the internal energy) are mimicked. The quality of these caricatures is tested numerically at the end. (author) [fr

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

Numerical modelling of two phase flow with hysteresis in heterogeneous porous media

Energy Technology Data Exchange (ETDEWEB)

Abreu, E. [Instituto Nacional de Matematica Pura e Aplicada (IMPA), Rio de Janeiro, RJ (Brazil); Furtado, F.; Pereira, F. [University of Wyoming, Laramie, WY (United States). Dept. of Mathematicsatics; Souza, G. [Universidade do Estado do Rio de Janeiro (UERJ), RJ (Brazil)

2008-07-01

Numerical simulators are necessary for the understanding of multiphase flow in porous media in order to optimize hydrocarbon recovery. In this work, the immiscible flow of two incompressible phases, a problem very common in waterflooding of petroleum reservoirs, is considered and numerical simulation techniques are presented. The system of equations which describe this type of flow form a coupled, highly nonlinear system of time-dependent partial differential equations (PDEs). The equation for the saturation of the invading fluid is a convection-dominated, degenerate parabolic PDE whose solutions typically exhibit sharp fronts (i.e., internal layers with strong gradients) and is very difficult to approximate numerically. It is well known that accurate modeling of convective and diffusive processes is one of the most daunting tasks in the numerical approximation of PDEs. Particularly difficult is the case where convection dominates diffusion. Specifically, we consider the injection problem for a model of two-phase (water/oil) flow in a core sample of porous rock, taking into account hysteresis effects in the relative permeability of the oil phase. (author)

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.

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

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

Theoretical and pragmatic modelling of governing equations for a two-phase flow in bubbly and annular flow regimes

International Nuclear Information System (INIS)

Bottoni, M.; Sengpiel, W.

1992-01-01

Starting from the rigorous formulation of the conservation equations for mass, momentum and enthalpy, derived for a two-phase flow by volume averaging microscopic balance equations over Eulerian control cells, the article discusses the formulation of the terms describing exchanges between the phases. Two flow regimes are taken into consideration, bubbly flow, applicable for small or medium void fractions, and annular flow, for large void fractions. When lack of knowledge of volume-averaged physical quantities make the rigorously formulated terms useless for computational purposes, modelling of these terms is discussed. 3 figs., 15 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

Electrical Capacitance Probe Characterization in Vertical Annular Two-Phase Flow

Directory of Open Access Journals (Sweden)

Grazia Monni

2013-01-01

Full Text Available The paper presents the experimental analysis and the characterization of an electrical capacitance probe (ECP that has been developed at the SIET Italian Company, for the measurement of two-phase flow parameters during the experimental simulation of nuclear accidents, as LOCA. The ECP is used to investigate a vertical air/water flow, characterized by void fraction higher than 95%, with mass flow rates ranging from 0.094 to 0.15 kg/s for air and from 0.002 to 0.021 kg/s for water, corresponding to an annular flow pattern. From the ECP signals, the electrode shape functions (i.e., the signals as a function of electrode distances in single- and two-phase flows are obtained. The dependence of the signal on the void fraction is derived and the liquid film thickness and the phase’s velocity are evaluated by means of rather simple models. The experimental analysis allows one to characterize the ECP, showing the advantages and the drawbacks of this technique for the two-phase flow characterization at high void fraction.

Two-phase critical flow models: a technical addendum to the CSNI state of the art report on critical flow modelling

International Nuclear Information System (INIS)

D'Auria, F.; Vigni, P.

1980-05-01

The purpose of this work was to obtain a comprehensive survey on the two-phase flow dynamics during accidental situations in nuclear reactors. About sixty theories regarding the two-phase flow calculation have been reviewed in this report with particular reference to their physical basis and assumptions; the aim is to control their applicability to nuclear safety problems. The main conclusions may be drawn as follows: the examined theories (perfect fluid, theories assuming thermodynamical equilibrium between liquid and vapor phases, non equilibrium models, etc.) are very different both for formulation and results; general validity of most theories is troublesome to check for the use of empirical coefficients. Moreover, according to the author's opinion, it is necessary to set up an organic program to obtain reliable experimental results in this field and to develop a model considering the whole blowdown transient

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

Multiphase reacting flows modelling and simulation

CERN Document Server

Marchisio, Daniele L

2007-01-01

The papers in this book describe the most widely applicable modeling approaches and are organized in six groups covering from fundamentals to relevant applications. In the first part, some fundamentals of multiphase turbulent reacting flows are covered. In particular the introduction focuses on basic notions of turbulence theory in single-phase and multi-phase systems as well as on the interaction between turbulence and chemistry. In the second part, models for the physical and chemical processes involved are discussed. Among other things, particular emphasis is given to turbulence modeling strategies for multiphase flows based on the kinetic theory for granular flows. Next, the different numerical methods based on Lagrangian and/or Eulerian schemes are presented. In particular the most popular numerical approaches of computational fluid dynamics codes are described (i.e., Direct Numerical Simulation, Large Eddy Simulation, and Reynolds-Averaged Navier-Stokes approach). The book will cover particle-based meth...

Challenges in modeling unstable two-phase flow experiments in porous micromodels

Science.gov (United States)

Meheust, Y.; Ferrari, A.; Jimenez-Martinez, J.; Le Borgne, T.; Lunati, I.

2014-12-01

The simulation of unstable invasion patterns in porous media flow is challenging since small perturbations tend to grow in time, so that slight differences in geometry or initial conditions potentially give rise to significantly different solutions. Here we present a detailed comparison of pore scale simulations and experiments of unstable primary drainage in porous micromodels. The porous medium consists of a Hele-Shaw cell containing cylindrical obstacles. Two experimental flow cells have been constructed by soft lithography, with different degrees of heterogeneity in the grain size distribution. To model two-phase flow at the pore scale, we solve Navier-Stokes equations for mass and momentum conservation in the discretized pore space and employ the Volume of Fluid (VOF) method to track the evolution of the interface. During drainage, if the defending fluid is the most viscous, viscous forces destabilize the interface, giving rise to the formation of preferential flow paths, in the form of a branched fingering structure. We test different numerical models (a 2D vertical integrated model and a full 3D model) and different initial conditions, studying their impact on the simulated spatial distributions of the fluid phases. Although due to the unstable nature of the invasion, small discrepancies between the experimental setup and the numerical model can result in different fluids patterns (see figure), simulations show a satisfactory agreement with the structures observed experimentally. To estimate the ability of the numerical approach to reproduce unstable displacement, we compare several quantities in both the statistical and deterministic sense. We demonstrate the impact of three main sources of uncertainty : i) the uncertainty on the pore space geometry, ii) the interface initialization and ii) three dimensional effects [1]. Simulations in weakly heterogeneous geometries are found to be more challenging because uncertainties on pore neck widths are on the same

Groundwater flow modelling under ice sheet conditions in Greenland (phase II)

International Nuclear Information System (INIS)

Jaquet, Olivier; Namar, Rabah; Siegel, Pascal; Jansson, Peter

2012-11-01

Within the framework of the GAP project, this second phase of geosphere modelling has enabled the development of an improved regional model that has led to a better representation of groundwater flow conditions likely to occur under ice sheet conditions. New data in relation to talik geometry and elevation, as well as to deformation zones were integrated in the geosphere model. In addition, more realistic hydraulic properties were considered for geosphere modelling; they were taken from the Laxemar site in Sweden. The geological medium with conductive deformation zones was modelled as a 3D continuum with stochastically hydraulic properties. Surface and basal glacial meltwater rates provided by a dynamic ice sheet model were assimilated into the groundwater flow model using mixed boundary conditions. The groundwater flow system is considered to be governed by infiltration of glacial meltwater in heterogeneous faulted crystalline rocks in the presence of permafrost and taliks. The characterisation of the permafrost-depth distribution was achieved using a coupled description of flow and heat transfer under steady state conditions. Using glaciological concepts and satellite data, an improved stochastic model was developed for the description at regional scale for the subglacial permafrost distribution in correlation with ice velocity and bed elevation data. Finally, the production of glacial meltwater by the ice sheet was traced for the determination of its depth and lateral extent. The major improvements are related to the type and handling of the subglacial boundary conditions. The use of meltwater rates provided by an ice sheet model applied as input to a mixed boundary condition enables to produce a more plausible flow field in the Eastern part of the domain, in comparison to previous modelling results (Jaquet et al. 2010). In addition, the integration of all potential taliks within the modelled domain provides a better characterisation of the likely groundwater

Groundwater flow modelling under ice sheet conditions in Greenland (phase II)

Energy Technology Data Exchange (ETDEWEB)

Jaquet, Olivier; Namar, Rabah; Siegel, Pascal [In2Earth Modelling Ltd, Lausanne (Switzerland); Jansson, Peter [Dept. of Physical Geography and Quaternary Geology, Stockholm Univ., Stockholm (Sweden)

2012-11-15

Within the framework of the GAP project, this second phase of geosphere modelling has enabled the development of an improved regional model that has led to a better representation of groundwater flow conditions likely to occur under ice sheet conditions. New data in relation to talik geometry and elevation, as well as to deformation zones were integrated in the geosphere model. In addition, more realistic hydraulic properties were considered for geosphere modelling; they were taken from the Laxemar site in Sweden. The geological medium with conductive deformation zones was modelled as a 3D continuum with stochastically hydraulic properties. Surface and basal glacial meltwater rates provided by a dynamic ice sheet model were assimilated into the groundwater flow model using mixed boundary conditions. The groundwater flow system is considered to be governed by infiltration of glacial meltwater in heterogeneous faulted crystalline rocks in the presence of permafrost and taliks. The characterisation of the permafrost-depth distribution was achieved using a coupled description of flow and heat transfer under steady state conditions. Using glaciological concepts and satellite data, an improved stochastic model was developed for the description at regional scale for the subglacial permafrost distribution in correlation with ice velocity and bed elevation data. Finally, the production of glacial meltwater by the ice sheet was traced for the determination of its depth and lateral extent. The major improvements are related to the type and handling of the subglacial boundary conditions. The use of meltwater rates provided by an ice sheet model applied as input to a mixed boundary condition enables to produce a more plausible flow field in the Eastern part of the domain, in comparison to previous modelling results (Jaquet et al. 2010). In addition, the integration of all potential taliks within the modelled domain provides a better characterisation of the likely groundwater

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.

Modeling of subcooled boiling in the vertical flow

International Nuclear Information System (INIS)

Koncar, B.; Mavko, B.

1999-01-01

A two-dimensional model of subcooled boiling in a vertical channel was developed. Its basic idea is that the vapor phase generation has a similar effect on the flow field as a hypothetical liquid phase generation. The bubble volume, generated due to evaporation process, was filled with liquid and included as a source term in the continuity equation for the liquid phase. Thus, the single-phase from of transport equations was preserved and bubbles were retained in the boundary layer near the heated surface. Time development of subcooled boiling was simulated and effects of governing physical mechanisms (evaporation, condensation, vapor-phase convection, vapor-phase diffusion) on the flow field and pressure drop were analyzed. The Results of the proposed two-dimensional model were compared with experimental data and RELAP5/MOD3.2 calculations. The presented model represents a contribution to the two-dimensional simulation of the subcooled boiling phenomenon.(author)

Seabrook simulator model upgrade: Implementation and validation of two-phase, nonequilibrium RCS and steam generator models

International Nuclear Information System (INIS)

Kao, S.

1990-01-01

A number of deficiencies in the original RCS and steam generator models on the Seabrook simulator were found to give unrealistic results under some off-normal and accident conditions. These deficiencies are attributed to the simplistic assumptions used in the original models, such as the homogeneous, equilibrium equations used in the pressurizer and steam generator models, and the single-phase flow model used in the RCS thermal-hydraulic model. To improve the fidelity of the simulator, efforts have been made to upgrade the RCS and steam generator models to include two-phase, nonequilibrium features. In the new RCS model, the following major assumptions are used to derive the finite difference form of the conservation equations: a donor-cell differencing scheme is adopted to allow flow reversal; a single pressure is used to evaluate properties; a single mass flow rate is assumed in each loop; enthalpy is assumed to vary linearly within each control volume; a homogeneous flow is assumed under two-phase conditions. The pressurizer is divided into a vapor region and a liquid region, each of which is represented by a set of mass and energy conservation equations. Interfacial mass and energy exchange mechanisms (condensation and flashing), thermal interactions between the vessel and fluids, and thermal nonequilibrium between the phases are included in the pressurizer model. The steam generator is divided into the vapor dome, riser, and downcomer regions. The assumptions applied are similar to those of the RCS and pressurizer models. A momentum model is incorporated to calculate the recirculation flow and simulate the downcomer level shrink/swell phenomenon. The new RCS and steam generator models are validated by comparing the simulator calculations against sister plant data and FSAR vendor analysis. The results show the new models give realistic and reliable calculations under off-normal and accident conditions

General correlations for pressure drop and heat transfer for single-phase turbulent flow in internally ribbed tubes

International Nuclear Information System (INIS)

Ravigururajan, T.S.; Bergles, A.E.

1985-01-01

General correlations for friction factors and heat transfer coefficients for single-phase turbulent flow in internally ribbed tubes are presented. Data from previous investigations are gathered for a wide range of tube parameters with e/d: 0.01 to 0.2; p/d: 0.1 to 7.0; α/90: 0.3 to 1.0, and flow parameters Re: 5000 to 250,000 and Pr: 0.66 to 37.6. The data were applied to a linear model to get normalized correlations that were then modified to fit tubes with extremely small parametric values. A shape function was included in the friction correlation to account for different rib profiles. The friction correlation predicts 96% of the data base to within +. 50% and 77% of the data base to within +. 20%. Corresponding figures for the heat transfer correlation are 99% and 69%. The present correlations are superior, for this extensive data base, to those presented by other investigators

The analysis of two-phase flow and heat transfer using a multidimensional, four field, two-fluid model

International Nuclear Information System (INIS)

Lahey, Richard T.; Drew, Donald A.

2001-01-01

This paper reviews the state-of-the-art in the prediction of multidimensional multiphase flow and heat transfer phenomena using a four field, two-fluid model. It is shown that accurate mechanistic computational fluid dynamic (CFD) predictions are possible for a wide variety of adiabatic and diabatic flows using this computational model. In particular, the model is able to predict the bubbly air/water upflow data of Serizawa (Serizawa, A., 1974. Fluid dynamic characteristics of two-phase flow. Ph.D. thesis, (Nuclear Engineering), Kyoto University, Japan), the downflow data of Wang et al. (Wang, S.K., Lee, S.J., Lahey Jr., R.T., Jones, O.C., 1987. 3-D turbulence structure and phase distribution measurements in bubbly two-phase flows. Int. J. Multiphase Flow 13 (3), 327-343), the isosceles triangle upflow data of Lopez de Bertodano et al. (Lopez de Bertodano, M., Lahey Jr., R.T., Jones, O.C., 1994b. Phase distribution in bubbly two-phase flow in vertical ducts. Int. J. Multiphase Flow 20 (5), 805-818), the heated annular R-113 subcooled boiling data of Velidandala, et al. (Velidandla, V., Pulta, S., Roy, P., Kaira, S.P., 1995. Velocity field in turbulent subcooled boiling flow. ASME Preprint HTD-314, 107-123) and the R-113 CHF data of Hino and Ueda (Hino, R., Ueda, T., 1985. Studies on heat transfer and flow characteristics in subcooled boiling-part 2, flow characteristics. Int. J. Multiphase Flow 11, 283-297). It can also predict external two-phase flows, such as those for spreading two-phase jets (Bonetto, F., Lahey Jr., R.T., 1993. An experimental study on air carryunder due to a plunging liquid jet. Int. J. Multiphase Flow 19 (2), 281-294) and multiphase flows around the hull of naval surface ships (Carrica, P.M., Bonetto, F., Drew, D.A., Lahey, R.T., 1999. A polydispersed model for bubbly two-phase flow around a surface ship. Int. J. Multiphase Flow 25 (2), 257-305)

Probabilistic physical characteristics of phase transitions at highway bottlenecks: incommensurability of three-phase and two-phase traffic-flow theories.

Science.gov (United States)

Kerner, Boris S; Klenov, Sergey L; Schreckenberg, Michael

2014-05-01

Physical features of induced phase transitions in a metastable free flow at an on-ramp bottleneck in three-phase and two-phase cellular automaton (CA) traffic-flow models have been revealed. It turns out that at given flow rates at the bottleneck, to induce a moving jam (F → J transition) in the metastable free flow through the application of a time-limited on-ramp inflow impulse, in both two-phase and three-phase CA models the same critical amplitude of the impulse is required. If a smaller impulse than this critical one is applied, neither F → J transition nor other phase transitions can occur in the two-phase CA model. We have found that in contrast with the two-phase CA model, in the three-phase CA model, if the same smaller impulse is applied, then a phase transition from free flow to synchronized flow (F → S transition) can be induced at the bottleneck. This explains why rather than the F → J transition, in the three-phase theory traffic breakdown at a highway bottleneck is governed by an F → S transition, as observed in real measured traffic data. None of two-phase traffic-flow theories incorporates an F → S transition in a metastable free flow at the bottleneck that is the main feature of the three-phase theory. On the one hand, this shows the incommensurability of three-phase and two-phase traffic-flow theories. On the other hand, this clarifies why none of the two-phase traffic-flow theories can explain the set of fundamental empirical features of traffic breakdown at highway bottlenecks.

Measurement of phase interaction in dispersed gas-particle two-phase flow by phase-doppler anemometry

Directory of Open Access Journals (Sweden)

Mergheni Ali Mohamed

2008-01-01

Full Text Available For simultaneous measurement of size and velocity distributions of continuous and dispersed phases in a two-phase flow a technique phase-Doppler anemometry was used. Spherical glass particles with a particle diameter range from 102 up to 212 µm were used. In this two-phase flow an experimental results are presented which indicate a significant influence of the solid particles on the flow characteristics. The height of influence of these effects depends on the local position in the jet. Near the nozzle exit high gas velocity gradients exist and therefore high turbulence production in the shear layer of the jet is observed. Here the turbulence intensity in the two-phase jet is decreased compared to the single-phase jet. In the developed zone the velocity gradient in the shear layer is lower and the turbulence intensity reduction is higher. .

Simulation of two-phase flows in vertical tubes with the CTFD code FLUBOX

International Nuclear Information System (INIS)

Graf, Udo; Papadimitriou, Pavlos

2007-01-01

The computational two-fluid dynamics (CTFD) code FLUBOX is developed at GRS for the multidimensional simulation of two-phase flows. The single-pressure two-fluid model is used as basis of the simulation. A basic mathematical property of the two-fluid model of FLUBOX is the hyperbolic character of the advection. The numerical solution methods of FLUBOX make explicit use of the hyperbolic structure of the coefficient matrices. The simulation of two-phase flow phenomena needs, apart from the conservation equations for each phase, an additional transport equation for the interfacial area concentration. The concentration of the interfacial area is one of the key parameters for the modeling of interfacial friction forces and interfacial transfer terms. A new transport equation for the interfacial area concentration is in development. It describes the dynamic change of the interfacial area concentration due to mass exchange and a force balance at the phase boundary. Results from FLUBOX calculations for different experiments of two-phase flows in vertical tubes are presented as part of the validation

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

Pressure drop and heat transfer of lithium single-phase flow under transverse magnetic field

International Nuclear Information System (INIS)

Takahashi, Minoru; Aritomi, Masanori; Inoue, Akira; Matsuzaki, Mitsuo

1996-01-01

Pressure drop and heat transfer characteristics of a lithium single-phase flow in a rectangular channel was investigated experimentally in the presence of a magnetic field. Friction loss coefficient under non-magnetic field and skin friction coefficient under magnetic field agreed well with the Blasius formula and a simple analytical expression, respectively. Nusselt number under non-magnetic field was slightly lower than the correlation by Hartnett and Irvine. Heat transfer was enhanced by increasing magnetic field above the Hartmann number of about 200. (author)

High Quality Model Predictive Control for Single Phase Grid Connected Photovoltaic Inverters

DEFF Research Database (Denmark)

Zangeneh Bighash, Esmaeil; Sadeghzadeh, Seyed Mohammad; Ebrahimzadeh, Esmaeil

2018-01-01

Single phase grid-connected inverters with LCL filter are widely used to connect the photovoltaic systems to the utility grid. Among the presented control schemes, predictive control methods are faster and more accurate but are more complex to implement. Recently, the model-predictive control...... algorithm for single-phase inverter has been presented, where the algorithm implementation is straightforward. In the proposed approach, all switching states are tested in each switching period to achieve the control objectives. However, since the number of the switching states in single-phase inverter...... is low, the inverter output current has a high total harmonic distortions. In order to reduce the total harmonic distortions of the injected current, this paper presents a high-quality model-predictive control for one of the newest structure of the grid connected photovoltaic inverter, i.e., HERIC...

Encyclopedia of two-phase heat transfer and flow III macro and micro flow boiling and numerical modeling fundamentals

CERN Document Server

2018-01-01

Set III of this encyclopedia is a new addition to the previous Sets I and II. It contains 26 invited chapters from international specialists on the topics of numerical modeling of two-phase flows and evaporation, fundamentals of evaporation and condensation in microchannels and macrochannels, development and testing of micro two-phase cooling systems for electronics, and various special topics (surface wetting effects, microfin tubes, two-phase flow vibration across tube bundles). The chapters are written both by renowned university researchers and by well-known engineers from leading corporate research laboratories. Numerous "must read" chapters cover the fundamentals of research and engineering practice on boiling, condensation and two-phase flows, two-phase heat transfer equipment, electronics cooling systems, case studies and so forth. Set III constitutes a "must have" reference together with Sets I and II for thermal engineering researchers and practitioners.

Design and development of drag-disc flow meter for measurement of transient two-phase flow

International Nuclear Information System (INIS)

Sreenivas Rao, G.; Kukreja, V.; Dolas, P.K.; Venkat Raj, V.

1989-01-01

Experiments have been carried out to test the suitability of drag-disc flowmeter for measuring two-phase flow. Calibration tests carried out under single-phase and two-phase flow conditions have confirmed the suitability of the drag-disc flowmeter. The experimental work and the results obtained are presented and discussed in the paper. (author). 6 figs

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.

A Novel Model Predictive Control for Single-Phase Grid-Connected Photovoltaic Inverters

DEFF Research Database (Denmark)

Zangeneh Bighash, Esmaeil; Sadeghzadeh, Seyed Mohammad; Ebrahimzadeh, Esmaeil

2017-01-01

Single-phase grid-connected inverters with LCL filter are widely used to connect photovoltaic systems to the utility grid. Among the existing control schemes, predictive control methods are faster and more accurate but also more complicated to implement. Recently, the Model Predictive Control (MPC......) algorithm for single-phase inverter has been presented, where the algorithm implementation is straightforward. In the MPC approach, all switching states are considered in each switching period to achieve the control objectives. However, since the number of switching states in single-phase inverters is small......, the inverter output current has a high Total Harmonic Distortions (THD). In order to reduce this, this paper presents an improved MPC for single-phase grid-connected inverters. In the proposed approach, the switching algorithm is changed and the number of the switching states is increased by means of virtual...

Experiments of steady state head and torque of centrifugal pumps in two-phase flow

International Nuclear Information System (INIS)

Minato, Akihiko; Tominaga, Kenji.

1988-01-01

Circulation pump behavior has large effect on coolant discharge flow rate in case of reactor pipe break. Experiment of two-phase pump performance was conducted as a joint study of Japanese BWR user utilities and makers. Two-phase head and torque of three centrifugal pumps in high temperature and high pressure (around 6 MPa) steam/water were measured. Head was decreased from single-phase characteristics when gas was mixed in liquid flow in condition with normal flow and normal rotation directions. When flow rate was large enough, two-phase head was about the same as single-phase one in reversal flow conditions. Two-phase head was smoothly increased as flowing steam volumetic concentration increased when flow rate was small and flow direction was reversal. Changes of torque with gas concentration were correspondent to those of head. This suggested that changes of interaction between flow and impellers due to phase slip effected on torque which caused head differences between single- and two-phase flows. Dependence of dimensionless head and torque of three test pumps on steam concentration were almost the same as each other. (author)

A study of critical two-phase flow models

International Nuclear Information System (INIS)

Siikonen, T.

1982-01-01

The existing computer codes use different boundary conditions in the calculation of critical two-phase flow. In the present study these boundary conditions are compared. It is shown that the boundary condition should be determined from the hydraulic model used in the computer code. The use of a correlation, which is not based on the hydraulic model used, leads often to bad results. Usually a good agreement with data is obtained in the calculation as far as the critical mass flux is concerned, but the agreement is not so good in the pressure profiles. The reason is suggested to be mainly in inadequate modeling of non-equilibrium effects. (orig.)

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

Deterministic sensitivity analysis of two-phase flow systems: forward and adjoint methods. Final report

International Nuclear Information System (INIS)

Cacuci, D.G.

1984-07-01

This report presents a self-contained mathematical formalism for deterministic sensitivity analysis of two-phase flow systems, a detailed application to sensitivity analysis of the homogeneous equilibrium model of two-phase flow, and a representative application to sensitivity analysis of a model (simulating pump-trip-type accidents in BWRs) where a transition between single phase and two phase occurs. The rigor and generality of this sensitivity analysis formalism stem from the use of Gateaux (G-) differentials. This report highlights the major aspects of deterministic (forward and adjoint) sensitivity analysis, including derivation of the forward sensitivity equations, derivation of sensitivity expressions in terms of adjoint functions, explicit construction of the adjoint system satisfied by these adjoint functions, determination of the characteristics of this adjoint system, and demonstration that these characteristics are the same as those of the original quasilinear two-phase flow equations. This proves that whenever the original two-phase flow problem is solvable, the adjoint system is also solvable and, in principle, the same numerical methods can be used to solve both the original and adjoint equations

Two-phase flow characteristic of inverted bubbly, slug and annular flow in post-critical heat flux region

International Nuclear Information System (INIS)

Ishii, M.; Denten, J.P.

1988-01-01

Inverted annular flow can be visualized as a liquid jet-like core surrounded by a vapor annulus. While many analytical and experimental studies of heat transfer in this regime have been performed, there is very little understanding of the basic hydrodynamics of the post-CHF flow field. However, a recent experimental study was done that was able to successfully investigate the effects of various steady-state inlet flow parameters on the post-CHF hydrodynamics of the film boiling of a single phase liquid jet. This study was carried out by means of a visual photographic analysis of an idealized single phase core inverted annular flow initial geometry (single phase liquid jet core surrounded by a coaxial annulus of gas). In order to extend this study, a subsequent flow visualization of an idealized two-phase core inverted annular flow geometry (two-phase central jet core, surrounded by a coaxial annulus of gas) was carried out. The objective of this second experimental study was to investigate the effect of steady-state inlet, pre-CHF two-phase jet core parameters on the hydrodynamics of the post-CHF flow field. In actual film boiling situations, two-phase flows with net positive qualities at the CHF point are encountered. Thus, the focus of the present experimental study was on the inverted bubbly, slug, and annular flow fields in the post dryout film boiling region. Observed post dryout hydrodynamic behavior is reported. A correlation for the axial extent of the transition flow pattern between inverted annular and dispersed droplet flow (the agitated regime) is developed. It is shown to depend strongly on inlet jet core parameters and jet void fraction at the dryout point. 45 refs., 9 figs., 4 tabs

Two-phase flow characteristic of inverted bubbly, slug, and annular flow in post-critical heat flux region

International Nuclear Information System (INIS)

Ishii, M.; Denten, J.P.

1989-01-01

Inverted annular flow can be visualized as a liquid jet-like core surrounded by a vapor annulus. While many analytical and experimental studies of heat transfer in this regime have been performed, there is very little understanding of the basic hydrodynamics of the post-critical heat flux (CHF) flow field. However, a recent experimental study was done that was able to successfully investigate the effects of various steady-state inlet flow parameters on the post-CHF hydrodynamics of the film boiling of a single phase liquid jet. This study was carried out by means of a visual photographic analysis of an idealized single phase core inverted annular flow initial geometry (single phase liquid jet core surrounded by a coaxial annulus of gas). In order to extend this study, a subsequent flow visualization of an idealized two-phase core inverted annular flow geometry (two-phase central jet core, surrounded by a coaxial annulus of gas) was carried out. The objective of this second experimental study was to investigate the effect of steady-state inlet, pre-CHF two-phase jet core parameters on the hydrodynamics of the post-CHF flow field. In actual film boiling situations, two-phase flows with net positive qualities at the CHF point are encountered. Thus, the focus of the present experimental study was on the inverted bubbly, slug, and annular flow fields in the post dryout film boiling region. Observed post dryout hydrodynamic behavior is reported. A correlation for the axial extent of the transition flow pattern between inverted annular and dispersed droplet flow (the agitated regime) is developed. It is shown to depend strongly on inlet jet core parameters and jet void fraction at the dryout point

A gas kinetic scheme for the Baer–Nunziato two-phase flow model

International Nuclear Information System (INIS)

Pan, Liang; Zhao, Guiping; Tian, Baolin; Wang, Shuanghu

2012-01-01

Numerical methods for the Baer–Nunziato (BN) two-phase flow model have attracted much attention in recent years. In this paper, we present a new gas kinetic scheme for the BN two-phase flow model containing non-conservative terms in the framework of finite volume method. In the view of microscopic aspect, a generalized Bhatnagar–Gross–Krook (BGK) model which matches with the BN model is constructed. Based on the integral solution of the generalized BGK model, we construct the distribution functions at the cell interface. Then numerical fluxes can be obtained by taking moments of the distribution functions, and non-conservative terms are explicitly introduced into the construction of numerical fluxes. In this method, not only the complex iterative process of exact solutions is avoided, but also the non-conservative terms included in the equation can be handled well.

Modeling of two-phase flow with thermal and mechanical non-equilibrium

International Nuclear Information System (INIS)

Houdayer, G.; Pinet, B.; Le Coq, G.; Reocreux, M.; Rousseau, J.C.

1977-01-01

To improve two-phase flow modeling by taking into account thermal and mechanical non-equilibrium a joint effort on analytical experiment and physical modeling has been undertaken. A model describing thermal non-equilibrium effects is first presented. A correlation of mass transfer has been developed using steam water critical flow tests. This model has been used to predict in a satisfactory manner blowdown tests. It has been incorporated in CLYSTERE system code. To take into account mechanical non-equilibrium, a six equations model is written. To get information on the momentum transfers special nitrogen-water tests have been undertaken. The first results of these studies are presented

Two-phase wall friction model for the trace computer code

International Nuclear Information System (INIS)

Wang Weidong

2005-01-01

The wall drag model in the TRAC/RELAP5 Advanced Computational Engine computer code (TRACE) has certain known deficiencies. For example, in an annular flow regime, the code predicts an unphysical high liquid velocity compared to the experimental data. To address those deficiencies, a new wall frictional drag package has been developed and implemented in the TRACE code to model the wall drag for two-phase flow system code. The modeled flow regimes are (1) annular/mist, (2) bubbly/slug, and (3) bubbly/slug with wall nucleation. The new models use void fraction (instead of flow quality) as the correlating variable to minimize the calculation oscillation. In addition, the models allow for transitions between the three regimes. The annular/mist regime is subdivided into three separate regimes for pure annular flow, annular flow with entrainment, and film breakdown. For adiabatic two-phase bubbly/slug flows, the vapor phase primarily exists outside of the boundary layer, and the wall shear uses single-phase liquid velocity for friction calculation. The vapor phase wall friction drag is set to zero for bubbly/slug flows. For bubbly/slug flows with wall nucleation, the bubbles are presented within the hydrodynamic boundary layer, and the two-phase wall friction drag is significantly higher with a pronounced mass flux effect. An empirical correlation has been studied and applied to account for nucleate boiling. Verification and validation tests have been performed, and the test results showed a significant code improvement. (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.

Direct numerical simulation of reactor two-phase flows enabled by high-performance computing

Energy Technology Data Exchange (ETDEWEB)

Fang, Jun; Cambareri, Joseph J.; Brown, Cameron S.; Feng, Jinyong; Gouws, Andre; Li, Mengnan; Bolotnov, Igor A.

2018-04-01

Nuclear reactor two-phase flows remain a great engineering challenge, where the high-resolution two-phase flow database which can inform practical model development is still sparse due to the extreme reactor operation conditions and measurement difficulties. Owing to the rapid growth of computing power, the direct numerical simulation (DNS) is enjoying a renewed interest in investigating the related flow problems. A combination between DNS and an interface tracking method can provide a unique opportunity to study two-phase flows based on first principles calculations. More importantly, state-of-the-art high-performance computing (HPC) facilities are helping unlock this great potential. This paper reviews the recent research progress of two-phase flow DNS related to reactor applications. The progress in large-scale bubbly flow DNS has been focused not only on the sheer size of those simulations in terms of resolved Reynolds number, but also on the associated advanced modeling and analysis techniques. Specifically, the current areas of active research include modeling of sub-cooled boiling, bubble coalescence, as well as the advanced post-processing toolkit for bubbly flow simulations in reactor geometries. A novel bubble tracking method has been developed to track the evolution of bubbles in two-phase bubbly flow. Also, spectral analysis of DNS database in different geometries has been performed to investigate the modulation of the energy spectrum slope due to bubble-induced turbulence. In addition, the single-and two-phase analysis results are presented for turbulent flows within the pressurized water reactor (PWR) core geometries. The related simulations are possible to carry out only with the world leading HPC platforms. These simulations are allowing more complex turbulence model development and validation for use in 3D multiphase computational fluid dynamics (M-CFD) codes.

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.

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

Two-phase flow in the cooling circuit of a cryogenic rocket engine

Science.gov (United States)

Preclik, D.

1992-07-01

Transient two-phase flow was investigated for the hydrogen cooling circuit of the HM7 rocket engine. The nuclear reactor code ATHLET/THESEUS was adapted to cryogenics and applied to both principal and prototype experiments for validation and simulation purposes. The cooling circuit two-phase flow simulation focused on the hydrogen prechilling and pump transient phase prior to ignition. Both a single- and a multichannel model were designed and employed for a valve leakage flow, a nominal prechilling flow, and a prechilling with a subsequent pump-transient flow. The latter case was performed in order to evaluate the difference between a nominal and a delayed turbo-pump start-up. It was found that an extension of the nominal prechilling sequence in the order of 1 second is sufficient to finally provide for liquid injection conditions of hydrogen which, as commonly known, is undesirable for smooth ignition and engine starting transients.

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

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

Computational simulation of flow and heat transfer in single-phase natural circulation loops

International Nuclear Information System (INIS)

Pinheiro, Larissa Cunha

2017-01-01

Passive decay heat removal systems based on natural circulation are essential assets for the new Gen III+ nuclear power reactors and nuclear spent fuel pools. The aim of the present work is to study both laminar and turbulent flow and heat transfer in single-phase natural circulation systems through computational fluid dynamics simulations. The working fluid is considered to be incompressible with constant properties. In the way, the Boussinesq Natural Convection Hypothesis was applied. The model chosen for the turbulence closure problem was the k -- εThe commercial computational fluid dynamics code ANSYS CFX 15.0 was used to obtain the numerical solution of the governing equations. Two single-phase natural circulation circuits were studied, a 2D toroidal loop and a 3D rectangular loop, both with the same boundary conditions of: prescribed heat flux at the heater and fixed wall temperature at the cooler. The validation and verification was performed with the numerical data provided by DESRAYAUD et al. [1] and the experimental data provided by MISALE et al. [2] and KUMAR et al. [3]. An excellent agreement between the Reynolds number (Re) and the modified Grashof number (Gr_m), independently of Prandtl Pr number was observed. However, the convergence interval was observed to be variable with Pr, thus indicating that Pr is a stability governing parameter for natural circulation. Multiple steady states was obtained for Pr = 0,7. Finally, the effect of inclination was studied for the 3D circuit, both in-plane and out-of-plane inclinations were verified for the steady state laminar regime. As a conclusion, the Re for the out-of-plane inclination was in perfect agreement with the correlation found for the zero inclination system, while for the in-plane inclined system the results differ from that of the corresponding vertical loop. (author)

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

KAUST Repository

Shao, Sihong; Qian, Tiezheng

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

Numerical approximation of a binary fluid-surfactant phase field model of two-phase incompressible flow

KAUST Repository

Zhu, Guangpu

2018-04-17

In this paper, we consider the numerical approximation of a binary fluid-surfactant phase field model of two-phase incompressible flow. The nonlinearly coupled model consists of two Cahn-Hilliard type equations and incompressible Navier-Stokes equations. Using the Invariant Energy Quadratization (IEQ) approach, the governing system is transformed into an equivalent form, which allows the nonlinear potentials to be treated efficiently and semi-explicitly. we construct a first and a second-order time marching schemes, which are extremely efficient and easy-to-implement, for the transformed governing system. At each time step, the schemes involve solving a sequence of linear elliptic equations, and computations of phase variables, velocity and pressure are totally decoupled. We further establish a rigorous proof of unconditional energy stability for the semi-implicit schemes. Numerical results in both two and three dimensions are obtained, which demonstrate that the proposed schemes are accurate, efficient and unconditionally energy stable. Using our schemes, we investigate the effect of surfactants on droplet deformation and collision under a shear flow. The increase of surfactant concentration can enhance droplet deformation and inhibit droplet coalescence.

Numerical approximation of a binary fluid-surfactant phase field model of two-phase incompressible flow

KAUST Repository

Zhu, Guangpu; Kou, Jisheng; Sun, Shuyu; Yao, Jun; Li, Aifen

2018-01-01

In this paper, we consider the numerical approximation of a binary fluid-surfactant phase field model of two-phase incompressible flow. The nonlinearly coupled model consists of two Cahn-Hilliard type equations and incompressible Navier-Stokes equations. Using the Invariant Energy Quadratization (IEQ) approach, the governing system is transformed into an equivalent form, which allows the nonlinear potentials to be treated efficiently and semi-explicitly. we construct a first and a second-order time marching schemes, which are extremely efficient and easy-to-implement, for the transformed governing system. At each time step, the schemes involve solving a sequence of linear elliptic equations, and computations of phase variables, velocity and pressure are totally decoupled. We further establish a rigorous proof of unconditional energy stability for the semi-implicit schemes. Numerical results in both two and three dimensions are obtained, which demonstrate that the proposed schemes are accurate, efficient and unconditionally energy stable. Using our schemes, we investigate the effect of surfactants on droplet deformation and collision under a shear flow. The increase of surfactant concentration can enhance droplet deformation and inhibit droplet coalescence.

A modified variable physical properties model, for analyzing nanofluids flow and heat transfer over nonlinearly stretching sheet

Directory of Open Access Journals (Sweden)

Pooria Akbarzadeh

2017-07-01

Full Text Available In this paper, the problem of laminar nanofluid flow which results from the nonlinear stretching of a flat sheet is investigated numerically. In this paper, a modified variable physical properties model for analyzing nanofluids flow and heat transfer is introduced. In this model, the effective viscosity, density, and thermal conductivity of the solid-liquid mixture (nanofluids which are commonly utilized in the homogenous single-phase model, are locally combined with the prevalent single-phase model. A numerical similarity solution is considered which depends on the local Prandtl number, local Brownian motion number, local Lewis number, and local thermophoresis number. The results are compared to the prevalent single-phase model. This comparison depicts that the prevalent single-phase model has a considerable deviation for predicting the behavior of nanofluids flow especially in dimensionless temperature and nanoparticle volume fraction. In addition the effect of the governing parameters such as Prandtl number, the Brownian motion number, the thermophoresis parameter, the Lewis number, and etc. on the velocity, temperature, and volume fraction distribution and the dimensionless heat and mass transfer rates are examined.

Advancement and Application of Multi-Phase CFD Modeling to High Speed Supercavitating Flows

Science.gov (United States)

2013-08-13

October 2008 - December 2013 4. TITLE AND SUBTITLE Advancement and Application of Multi-Phase CFD Modeling to High Speed Supercavitating Flows...influence cavity hysteresis behavior. These observations are used to guide improved supercavitating -vehicle analyses including numerical predictions...experiments, and modeling 15. SUBJECT TERMS supercavitation , computational fluid dynamics, multiphase flow 16. SECURITY CLASSIFICATION OF: a

A method to couple HEM and HRM two-phase flow models

Energy Technology Data Exchange (ETDEWEB)

Herard, J.M.; Hurisse, O. [Elect France, Div Rech and Dev, Dept Mecan Fluides Energies and Environm, F-78401 Chatou (France); Hurisse, O. [Univ Aix Marseille 1, Ctr Math and Informat, Lab Anal Topol and Probabil, CNRS, UMR 6632, F-13453 Marseille 13 (France); Ambroso, A. [CEA Saclay, DEN, DM2S, SFME, LETR, 91 - Gif sur Yvette (France)

2009-04-15

We present a method for the unsteady coupling of two distinct two-phase flow models (namely the Homogeneous Relaxation Model, and the Homogeneous Equilibrium Model) through a thin interface. The basic approach relies on recent works devoted to the interfacial coupling of CFD models, and thus requires to introduce an interface model. Many numerical test cases enable to investigate the stability of the coupling method. (authors)

A method to couple HEM and HRM two-phase flow models

International Nuclear Information System (INIS)

Herard, J.M.; Hurisse, O.; Hurisse, O.; Ambroso, A.

2009-01-01

We present a method for the unsteady coupling of two distinct two-phase flow models (namely the Homogeneous Relaxation Model, and the Homogeneous Equilibrium Model) through a thin interface. The basic approach relies on recent works devoted to the interfacial coupling of CFD models, and thus requires to introduce an interface model. Many numerical test cases enable to investigate the stability of the coupling method. (authors)

Assessment of TRAC-PF1/MOD3 Mark-22 assembly model using SRL ''A'' tank single-assembly flow experiments

International Nuclear Information System (INIS)

Fischer, S.R.; Lam, K.; Lin, J.C.

1991-01-01

This paper summarizes the results of an assessment of our TRAC-PF1/MOD3 Mark-22 prototype fuel assembly model against single-assembly data obtained from the ''A'' Tank single-assembly tests that were performed at the Savannah River Laboratory. We felt the data characterize prototypic assembly behavior over a range of air-water flow conditions of interest for loss-of-coolant accident (LOCA) calculations. This study was part of a benchmarking effort performed to evaluate and validate a multiple-assembly, full-plant model that is being developed by Los Alamos National Laboratory to study various aspects of the Savannah River plant operating conditions, including LOCA transients, using TRAC-PF1/MOD3 Version 1.10. The results of this benchmarking effort demonstrate that TRAC-PF1/MOD3 is capable pf calculating plenum conditions and assembly flows during conditions thought to be typical of the Emergency Cooling System (ECS) phase of a LOCA. 10 refs., 12 fig

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

Traffic flow behavior at a single-lane urban roundabout

Science.gov (United States)

Lakouari, N.; Oubram, O.; Ez-Zahraouy, H.; Cisneros-Villalobos, L.; Velásquez-Aguilar, J. G.

In this paper, we propose a stochastic cellular automata model to study the traffic behavior at a single-lane roundabout. Vehicles can enter the interior lane or exit from it via N intersecting lane, the boundary conditions are stochastic. The traffic is controlled by a self-organized scheme. It has turned out that depending on the rules of insertion to the roundabout, five distinct traffic phases can appear, namely, free flow, congestion, maximum current, jammed and gridlock. The transition between the free flow and the gridlock is forbidden. The density profiles are used to study the traffic pattern at the interior lane of the roundabout. In order to quantify the interactions between vehicles in the interior lane of the roundabout, the velocity correlation coefficient (VCC) is also studied. Besides, the spatiotemporal diagrams corresponding to the entry/exit lanes are derived numerically. Furthermore, we have investigated the effect of displaying signal (PIn), as the PIn decreases, the maximum current increases at the expense of the free flow and the jamming phase. Finally, we have investigated the effect of the braking probability P on the interior lane of the roundabout. We have found that the increase of P raises the spontaneous jam formation on the ring. Thus, enlarges the maximum current and the jamming phase while the free flow phase decreases.

Development and calibration of instruments for measurements in transient two-phase flow

International Nuclear Information System (INIS)

Banerjee, S.; Heidrick, T.R.

1981-01-01

For validation and development of theoretical models for transient two-phase flow, it is necessary to measure local and cross-sectionally averaged thermalhydraulic parameters. Of these parameters, void fraction and mass velocity are the most difficult to measure. In this paper, we present our recent work on various techniques for determining these quantities. The possibility of determining flow regime by using fast neutron transmission is discussed. The development of a miniaturized electrical resistivity probe for measuring local void fraction is described, together with calibrations obtained by integrating the void fraction profile and comparing the cross-sectionally averaged void fraction with direct measurements using two quick closing valves. Results on the calibration of combinations of full-flow turbine meters, Pitot tube rakes and gamma densitometers for measuring cross-sectionally averaged mass velocity in steady steam-water flow are presented. The results are interpreted with a simple model using single-phase calibration factors for the Pitot tube rakes and turbine meters. Calibration experiments were also done in transient steam-water flows and interpretation of the results with the steady state models is also discussed

Studies on shock phenomena in two-phase flow, (4). Characteristics in channel flow consisting of bubbly mixture and liquid in series

Energy Technology Data Exchange (ETDEWEB)

Akagawa, Koji; Fujii, Terushige; Ito, Yutaka; Hiraki, Sei

1982-04-01

The research carried out so far was related to the case in which the mean void ratio in a pipe distributed almost invariably in axial direction. However, in actual piping system, the distribution of void ratio sometimes changes in axial direction such as evaporating tubes. In this study, in order to clarify the basic characteristics of shock phenomena in a piping system in which the density of two-phase flow changes in axial direction, experiment was carried out on air and water two-component bubbly flow, in which single phase was in upstream, and two-phase flow with constant void ratio in axial direction was in downstream. Also, the theoretical study on the phenomena was performed. The experimental setup and experimental method, the result of the waveform of pressure response, the behavior of pressure waves at the interface of two-phase flow and single phase flow, the qualitative analysis of the waveform of pressure response, and the analysis of pressure rise are reported. By the sudden closure of a valve, the pressure in two-phase flow rose by the initial potential surge, thereafter stepped pressure rise was observed. This phenomenon can be explained by the reflection of pressure waves at the interface of two-phase flow and single phase flow.

Analysis of forced convective transient boiling by homogeneous model of two-phase flow

International Nuclear Information System (INIS)

Kataoka, Isao

1985-01-01

Transient forced convective boiling is of practical importance in relation to the accident analysis of nuclear reactor etc. For large length-to-diameter ratio, the transient boiling characteristics are predicted by transient two-phase flow calculations. Based on homogeneous model of two-phase flow, the transient forced convective boiling for power and flow transients are analysed. Analytical expressions of various parameters of transient two-phase flow have been obtained for several simple cases of power and flow transients. Based on these results, heat flux, velocity and time at transient CHF condition are predicted analytically for step and exponential power increases, and step, exponential and linear velocity decreases. The effects of various parameters on heat flux, velocity and time at transient CHF condition have been clarified. Numerical approach combined with analytical method is proposed for more complicated cases. Solution method for pressure transient are also described. (author)

Computational fluid dynamics modeling of two-phase flow in a BWR fuel assembly. Final CRADA Report

International Nuclear Information System (INIS)

Tentner, A.

2009-01-01

A direct numerical simulation capability for two-phase flows with heat transfer in complex geometries can considerably reduce the hardware development cycle, facilitate the optimization and reduce the costs of testing of various industrial facilities, such as nuclear power plants, steam generators, steam condensers, liquid cooling systems, heat exchangers, distillers, and boilers. Specifically, the phenomena occurring in a two-phase coolant flow in a BWR (Boiling Water Reactor) fuel assembly include coolant phase changes and multiple flow regimes which directly influence the coolant interaction with fuel assembly and, ultimately, the reactor performance. Traditionally, the best analysis tools for this purpose of two-phase flow phenomena inside the BWR fuel assembly have been the sub-channel codes. However, the resolution of these codes is too coarse for analyzing the detailed intra-assembly flow patterns, such as flow around a spacer element. Advanced CFD (Computational Fluid Dynamics) codes provide a potential for detailed 3D simulations of coolant flow inside a fuel assembly, including flow around a spacer element using more fundamental physical models of flow regimes and phase interactions than sub-channel codes. Such models can extend the code applicability to a wider range of situations, which is highly important for increasing the efficiency and to prevent accidents.

Simulation of the phenomenon of single-phase and two-phase natural circulation

International Nuclear Information System (INIS)

Castrillo, Lazara Silveira

1998-02-01

Natural convection phenomenon is often used to remove the residual heat from the surfaces of bodies where the heat is generated e.g. during accidents or transients of nuclear power plants. Experimental study of natural circulation can be done in small scale experimental circuits and the results can be extrapolated for larger operational facilities. The numerical analysis of transients can be carried out by using large computational codes that simulate the thermohydraulic behavior in such facilities. The computational code RELAP5/MOD2, (Reactor Excursion and Leak Analysis Program) was developed by U.S. Nuclear Regulatory Commissions's. Division of Reactor Safety Research with the objective of analysis of transients and postulated accidents in the light water reactor (LWR) systems, including small and large ruptures with loss of coolant accidents (LOCA's). The results obtained by the simulation of single-phase and two-phase natural circulation, using the RELAP5/MOD2, are presented in this work. The study was carried out using the experimental circuit built at the 'Departamento de Engenharia Quimica da Escola Politecnica da Universidade de Sao Paulo'. In the circuit, two experiments were carried out with different conditions of power and mass flow, obtaining a single-phase regime with a level of power of 4706 W and flow of 5.10 -5 m 3 /s (3 l/min) and a two-phase regime with a level of power of 6536 W and secondary flow 2,33.10 -5 m 3 /s (1,4 l/min). The study allowed tio evaluate the capacity of the code for representing such phenomena as well as comparing the transients obtained theoretically with the experimental results. The comparative analysis shows that the code represents fairly well the single-phase transient, but the results for two-phase transients, starting from the nodalization and calibration used for the case single-phase transient, did not reproduce faithfully some experimental results. (author)

Magnetic resonance velocity imaging of liquid and gas two-phase flow in packed beds.

Science.gov (United States)

Sankey, M H; Holland, D J; Sederman, A J; Gladden, L F

2009-02-01

Single-phase liquid flow in porous media such as bead packs and model fixed bed reactors has been well studied by MRI. To some extent this early work represents the necessary preliminary research to address the more challenging problem of two-phase flow of gas and liquid within these systems. In this paper, we present images of both the gas and liquid velocities during stable liquid-gas flow of water and SF(6) within a packing of 5mm spheres contained within columns of diameter 40 and 27 mm; images being acquired using (1)H and (19)F observation for the water and SF(6), respectively. Liquid and gas flow rates calculated from the velocity images are in agreement with macroscopic flow rate measurements to within 7% and 5%, respectively. In addition to the information obtained directly from these images, the ability to measure liquid and gas flow fields within the same sample environment will enable us to explore the validity of assumptions used in numerical modelling of two-phase flows.

Local two-phase modeling of the water-steam flows occurring in steam generators

International Nuclear Information System (INIS)

Denefle, Romain

2013-01-01

The present study is related to the need of modeling the two-phase flows occurring in a steam generator (liquid at inlet and vapour at outlet). The choice is made to investigate a hybrid modeling of the flow, considering the gas phase as two separated fields, each one being modeled with different closure laws. In so doing, the small and spherical bubbles are modeled through a dispersed approach within the two-fluid model, and the distorted bubbles are simulated with an interface locating method. The main outcome is about the implementation, the verification and the validation of the model dedicated to the large and distorted bubbles, as well as the coupling of the two approaches for the gas, allowing the presentation of demonstration calculations using the so-called hybrid approach. (author)

BACCHUS-3D/SP. A computer programme for the three-dimensional description of sodium single-phase flow in bundle geometry

International Nuclear Information System (INIS)

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

1983-07-01

The computer programme BACCHUS implemented at KfK includes a steady-state version, a two-dimensional and a three-dimensional transient single-phase flow version describing the thermal-hydraulic behaviour of the coolant (sodium or water) in bundle geometry under nominal or accident conditions. All versions are coupled with a pin model describing the temperature distribution in fuel (or electrical heaters) and cladding. The report describes the programme from the viewpoints of the geometrical model, the mathematical foundations and the numerical treatment of the basic equations. Although emphasis is put on the three-dimensional version, the two-dimensional and the steady state versions are also documented in self-consistent sections. (orig.) [de

Validation of Friction Models in MARS-MultiD Module with Two-Phase Cross Flow Experiment

International Nuclear Information System (INIS)

Choi, Chi-Jin; Yang, Jin-Hwa; Cho, Hyoung-Kyu; Park, Goon-Cher; Euh, Dong-Jin

2015-01-01

In the downcomer of Advanced Power Reactor 1400 (APR1400) which has direct vessel injection (DVI) lines as an emergency core cooling system, multidimensional two-phase flow may occur due to the Loss-of-Coolant-Accident (LOCA). The accurate prediction about that is high relevance to evaluation of the integrity of the reactor core. For this reason, Yang performed an experiment that was to investigate the two-dimensional film flow which simulated the two-phase cross flow in the upper downcomer, and obtained the local liquid film velocity and thickness data. From these data, it could be possible to validate the multidimensional modules of system analysis codes. In this study, MARS-MultiD was used to simulate the Yang's experiment, and obtained the local variables. Then, the friction models used in MARS-MultiD were validated by comparing the two-phase flow experimental results with the calculated local variables. In this study, the two-phase cross flow experiment was modeled by the MARS-MultiD. Compared with the experimental results, the calculated results by the code properly presented mass conservation which could be known from the relation between the liquid film velocity and thickness at the same flow rate. The magnitude and direction of the liquid film, however, did not follow well with experimental results. According to the results of Case-2, wall friction should be increased, and interfacial friction should be decreased in MARS-MultiD. These results show that it is needed to modify the friction models in the MARS-MultiD to simulate the two-phase cross 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)

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)

CFD modelling of polydispersed bubbly two-phase flow around an obstacle

International Nuclear Information System (INIS)

Krepper, Eckhard; Beyer, Matthias; Frank, Thomas; Lucas, Dirk; Prasser, Horst-Michael

2009-01-01

A population balance model (the Inhomogeneous MUSIG model) has recently been developed in close cooperation between ANSYS-CFX and Forschungszentrum Dresden-Rossendorf and implemented into the CFD-Code CFX [Krepper, E., Lucas, D., Prasser, H.-M, 2005. On the modelling of bubbly flow in vertical pipes. Nucl. Eng. Des. 235, 597-611; Frank, T., Zwart, P.J., Shi, J.-M., Krepper, E., Rohde, U., 2005. Inhomogeneous MUSIG Model-a population balance approach for polydispersed bubbly flows, International Conference 'Nuclear Energy for New Europe 2005', Bled, Slovenia, September 5-8, 2005; Krepper, E., Beyer, M., Frank, Th., Lucas, D., Prasser, H.-M., 2007. Application of a population balance approach for polydispersed bubbly flows, 6th Int. Conf. on Multiphase Flow Leipzig 2007, (paper 378)]. The current paper presents a brief description of the model principles. The capabilities of this model are discussed via the example of a bubbly flow around a half-moon shaped obstacle arranged in a 200 mm pipe. In applying the Inhomogeneous MUSIG approach, a deeper understanding of the flow structures is possible and the model allows effects of polydispersion to be investigated. For the complex flow around the obstacle, the general structure of the flow was well reproduced in the simulations. This test case demonstrates the complicated interplay between size dependent bubble migration and the effects of bubble coalescence and breakup on real flows. The closure models that characterize the bubble forces responsible for the simulation of bubble migration show agreement with the experimental observations. However, clear deviations occur for bubble coalescence and fragmentation. The models applied here, which describe bubble fragmentation and coalescence could be proved as a weakness in the validity of numerous CFD analyses of vertical upward two-phase pipe flow. Further work on this topic is under way.

Experimental study of flow monitoring instruments in air-water, two-phase downflow

International Nuclear Information System (INIS)

Sheppard, J.D.; Hayes, P.H.; Wynn, M.C.

1976-01-01

The performance of a turbine meter, target flow meter (drag disk), and a gamma densitometer was studied in air-water, two-phase vertical downflow. Air and water were metered into an 0.0889-m-ID (3.5-in.) piping system; air flows ranged from 0.007 to 0.3 m 3 /sec (16 to 500 scfm) and water flows ranged from 0.0006 to 0.03 m 3 /sec (10 to 500 gpm). The study included effects of flow rate, quality, flow regime, and flow dispersion on the mean and fluctuating components of the instrument signals. Wire screen flow dispersers located at the inlet to the test section had a significant effect on the readings of the drag disk and gamma densitometer, but had little effect on the turbine. Further, when flow dispersers were used, mass flow rates determined from the three instrument readings and a two-velocity, slip flow model showed good agreement with actual mass flow rate over a three-fold range in quality; mass flows determined with the drag disk and densitometer readings assuming homogeneous flow were nearly as accurate. However, when mass flows were calculated using the turbine and densitometer or turbine and drag disk readings assuming homogeneous flow, results were scattered and relatively inaccurate compared to the actual mass flows. Turbine meter data were used with a two-velocity turbine model and continuity relationships for each phase to determine the void fraction and mean phase velocities in the test section. The void fraction was compared with single beam gamma densitometer results and fluid momentum calculated from a two-velocity model was compared with drag disk readings

A Navier-Stokes/Cahn-Hilliard model for the simulation of three phase immiscible incompressible flows

International Nuclear Information System (INIS)

Celine Lapuerta; Bruno Piar; Franck Boyer; Philippe Angot; Michel Quintard

2005-01-01

This paper presents a Navier-Stokes/Cahn-Hilliard model designed for incompressible flows of three immiscible phases, characterized by different surface tensions and without phase change. This physical context is relevant to study the late phase of a hypothetical severe accident in a nuclear pressurized water reactor. Thanks to a suitable choice of a free energy and a particular form of the Cahn-Hilliard equation, the evolution of the three phases is described by only two order parameters. Moreover, this model allows the simulation of purely two phase flows as a limiting case: no artificial apparition of the third phase occurs if this later is physically absent which contrasts with others models of the literature. We examine the spreading of a liquid lens at the interface between two stratified phases. We present results showing that the method gives correct contact angles and pressure jumps, at equilibrium. (authors)

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

Experimental study and theoretical modelling of two-phase flow in a converging diverging nozzle

International Nuclear Information System (INIS)

Selmer-Olsen, Stale

1991-01-01

A theoretical and experimental study of high quality two-phase flows in converging-diverging nozzles is presented. The main objectives are the prediction of critical (choked) flow rates and the evolution of characteristic parameters towards the nozzle outlet. First, a thorough analysis of available models shows the importance of a correct modelling of the mechanical and thermal interactions between the gas and liquid phases. As a second step, a purely dispersed flow model is considered. The solution algorithm which is utilized describes accurately the critical (choked) flow conditions as well as the topology of the solutions. The dispersed flow model accounts for effects on the gas flow rate of the upstream and the downstream pressures, the liquid flow rate and the nozzle geometry. The pressure profile along the nozzle and the location of the critical cross-section are also well predicted. The flow is shown to switch from critical to sub-critical when the liquid flow rate is increased, all other control parameters at the inlet and the outlet maintained. This new finding is interpreted as a result of the possible location of the critical cross-section anywhere in the diverging part of the nozzle. Moreover, the experiments show that the critical (choked) gas flow rate depends on the inlet configuration of gas/liquid. In the third step, a careful analysis of the data is used as a basis for proposing a new dispersed-annular flow model. This model accounts for the liquid flowing both as a liquid film and as entrained droplets in the core, non-developed flow is accounted for as well as flow separation in the diffuser. Finally, advanced local measuring techniques of pressure, film thickness and film velocity have been developed in the course of the work. In particular film thickness measurements allowed the development of the flow structure to be understood. (author) [fr

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

Introduction of phase field model to crimson code for two-phase flow simulation using unstructured mesh

International Nuclear Information System (INIS)

Tanaka, Nobuatsu; Maseguchi, Ryo; Ogawara, Takuya

2008-01-01

This study is concerned with improvement of numerical code called CRIMSON (Civa RefIned Multiphase SimulatiON), which has been developed to evaluate multi-phase flow behaviors based on the recent CFD (computational fluid dynamics) technologies. The CRIMSON employs a finite-volume method combined with the high order interpolation scheme, CIVA (cubic-interpolation with area/volume coordinates). The CRIMSON solves gas-liquid two phases by a unified scheme of CUP (combined unified procedure). The conventional CIVA method has two problems of interface blurring in long-term calculation and non-conservativeness. In this study, the problems were solved by introducing the ideas of the level set method and the phase field method. We verified out method by applying it to some popular benchmark problems of single bubble rising and collapse of water column problems. (author)

Behavior of pumps conveying two-phase liquid flow

International Nuclear Information System (INIS)

Grison, Pierre; Lauro, J.-F.

1979-01-01

Determination of the two-phase flow (critical or otherwise) through a pump is an essential requirement for complete description of a loss of primary coolant accident in a PWR plant. Theoretical and experimental research at Electricite de France on this subject is described and problems associated with the introduction of a two-phase fluid (with mass transfer) are discussed, with an attempt to single out new phenomena involved and establish their effect on pump behavior. A complementary experimental investigation is described and the results of tests at pressures and temperatures up to 120 bars and 320 0 C respectively are compared with the theoretical model data [fr

Behavior of pumps conveying two-phase liquid flow

Energy Technology Data Exchange (ETDEWEB)

Grison, P; Lauro, J F [Electricite de France, 78 - Chatou. Direction des Etudes et Recherches

1979-01-01

Determination of the two-phase flow (critical or otherwise) through a pump is an essential requirement for complete description of a loss of primary coolant accident in a PWR plant. Theoretical and experimental research at Electricite de France on this subject is described and problems associated with the introduction of a two-phase fluid (with mass transfer) are discussed, with an attempt to single out new phenomena involved and establish their effect on pump behavior. A complementary experimental investigation is described and the results of tests at pressures and temperatures up to 120 bars and 320/sup 0/C respectively are compared with the theoretical model data.

A numerical method for a model of two-phase flow in a coupled free flow and porous media system

KAUST Repository

Chen, Jie; Sun, Shuyu; Wang, Xiaoping

2014-01-01

In this article, we study two-phase fluid flow in coupled free flow and porous media regions. The model consists of coupled Cahn-Hilliard and Navier-Stokes equations in the free fluid region and the two-phase Darcy law in the porous medium region. We propose a Robin-Robin domain decomposition method 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. Numerical examples are presented to illustrate the effectiveness of this method. © 2014 Elsevier Inc.

A numerical method for a model of two-phase flow in a coupled free flow and porous media system

KAUST Repository

Chen, Jie

2014-07-01

In this article, we study two-phase fluid flow in coupled free flow and porous media regions. The model consists of coupled Cahn-Hilliard and Navier-Stokes equations in the free fluid region and the two-phase Darcy law in the porous medium region. We propose a Robin-Robin domain decomposition method 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. Numerical examples are presented to illustrate the effectiveness of this method. © 2014 Elsevier Inc.

Equations of motion for two-phase flow in a pin bundle of a nuclear reactor

International Nuclear Information System (INIS)

Chawla, T.C.; Ishii, M.

1978-01-01

By performing Eulerian area averaging over a channel area of the local continuity, momentum, and energy equations for single phase turbulent flow and assuming each phase in two-phase flows to be continuum but coupled by the appropriate 'jump' conditions at the interface, the corresponding axial macroscopic balances for two-fluid model in a pin bundle are obtained. To determine the crossflow, a momentum equation in transverse (to the gap between the pins) direction is obtained for each phase by carrying out Eulerian segment averaging of the local momentum equation, where the segment is taken parallel to the gap. By considering the mixture as a whole, a diffusion model based on drift-flux velocity is formulated. In the axial direction it is expressed in terms of three mixture conservation equations of mass, momentum, and energy with one additional continuity equation for the vapor phase. For the determination of crossflow, transverse momentum equation for a mixture is obtained. It is considered that the previous formulation of the two-phase flow based on the 'slip' flow model and the integral subchannel balances using finite control volumes is inadequate in that the model is heuristic and, a priori, assumes the order of magnitude of the terms, also the model is incomplete and incorrect when applied to two-phase mixtures in thermal non-equilibrium such as during accidental depressurization of a water cooled reactor. The governing equations presented are shown to be a very formal and sound physical basis and are indispensable for physically correct methods of analyzing two-phase flows in a pin bundle. (author)

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

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

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

TOPFLOW-experiments, model development and validation for the qualification of CFD-odes for two-phase flows. Final report

International Nuclear Information System (INIS)

Lucas, D.; Beyer, M.; Banowski, M.; Seidel, T.; Krepper, E.; Liao, Y.; Apanasevich, P.; Gauss, F.; Ma, T.

2016-12-01

This report summarizes the main results obtained in frame of the project. The aim of the project was the qualification of CFD-methods for two-phase flows with phase transfer relevant for nuclear safety research. To reach this aim CFD-grade experimental data are required. Such data can be obtained at the TOPFLOW facility because of the combination of experiments in scales and at parameters which are relevant for nuclear safety research with innovative measuring techniques. The experimental part of this project comprises investigations on flows in vertical pipes using the ultrafast X-ray tomography, on flows with and without phase transfer in a special test basin and on counter-current flow limitation in a model of a PWR hot leg. These experiments are only briefly presented in this report since detailed documentations are given in separated reports for all of these 3 experimental series. One important results of the activities devoted on CFD qualification is the establishment of the baseline model concept and the definition of the baseline model for poly-disperse bubbly flows. This is an important contribution to improve the predictive capabilities of CFD-models basing on the two- or multi-fluid approach. On the other hand, the innovative Generalized Two-Phase Flow concept (GENTOP) aims on an extension of the range of applicability of CFD-methods. In many relevant flow situations different morphologies of the phases or different flow pattern occur simultaneously in one flow domain. In addition transitions between these morphologies may occur. The GENTOP-concept for the first time a framework was established which allows the simulation of such flow situations in a consistent manner. Other activities of the project aim on special model developments to improve the simulation capabilities for flows with phase transfer.

Models for turbulent flows with variable density and combustion