Characteristics of two-phase flows in large diameter channels

Energy Technology Data Exchange (ETDEWEB)

Schlegel, J.P., E-mail: schlegelj@mst.edu [Department of Mining and Nuclear Engineering, Missouri University of Science and Technology, 301 W 14th St., Rolla, MO 65401 (United States); Hibiki, T.; Ishii, M. [School of Nuclear Engineering, Purdue University, 400 Central Dr., West Lafayette, IN 47907 (United States)

2016-12-15

Two-phase flows in large diameter channels have a great deal of importance in a wide variety of industrial applications. Nuclear systems, petroleum refineries, and chemical processes make extensive use of larger systems. Flows in such channels have very different properties from flows in smaller channels which are typically used in experimental research. In this paper, the various differences between flows in large and small channels are highlighted using the results of previous experimental and analytical research. This review is followed by a review of recent experiments in and model development for flows in large diameter channels performed by the authors. The topics of these research efforts range from void fraction and interfacial area concentration measurement to flow regime identification and modeling, drift-flux modeling for high void fraction conditions, and evaluation of interfacial area transport models for large diameter channels.

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

International Nuclear Information System (INIS)

Davis, C.B.

1998-01-01

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

Comparison of swirling strengths derived from two- and three-dimensional velocity fields in channel flow

Science.gov (United States)

Chen, Huai; Li, Danxun; Bai, Ruonan; Wang, Xingkui

2018-05-01

Swirling strength is an effective vortex indicator in wall turbulence, and it can be determined based on either two-dimensional (2D) or three-dimensional (3D) velocity fields, written as λci2D and λci3D, respectively. A comparison between λci2D and λci3D has been made in this paper in sliced XY, YZ, and XZ planes by using 3D DNS data of channel flow. The magnitude of λci2D in three orthogonal planes differs in the inner region, but the difference tends to diminish in the outer flow. The magnitude of λci3D exceeds each λci2D, and the square of λci3D is greater than the summation of squares of three λci2D. Extraction with λci2D in XY, YZ, and XZ planes yields different population densities and vortex sizes, i.e., in XZ plane, the vortices display the largest population density and the smallest size, and in XY and YZ planes the vortices are similar in size but fewer vortices are extracted in the XY plane in the inner layer. Vortex size increases inversely with the threshold used for growing the vortex region from background turbulence. When identical thresholds are used, the λci3D approach leads to a slightly smaller population density and a greater vortex radius than the λci2D approach. A threshold of 0.8 for the λci3D approach is approximately equivalent to a threshold of 1.5 for the λci2D approach.

On the large-scale structure and spectral dynamics of two-dimensional turbulence in a periodic channel

NARCIS (Netherlands)

Kramer, W.; Clercx, H.J.H.; Heijst, van G.J.F.

2008-01-01

This paper reports on a numerical study of forced two-dimensional turbulence in a periodic channel with flat no-slip walls. Since corners or curved domain boundaries, met in the standard rectangular, square or circular geometries, are absent in this geometry, the (statistical) analysis of the flow

A coupled three dimensional model of vanadium redox flow battery for flow field designs

International Nuclear Information System (INIS)

Yin, Cong; Gao, Yan; Guo, Shaoyun; Tang, Hao

2014-01-01

A 3D (three-dimensional) model of VRB (vanadium redox flow battery) with interdigitated flow channel design is proposed. Two different stack inlet designs, single-inlet and multi-inlet, are structured in the model to study the distributions of fluid pressure, electric potential, current density and overpotential during operation of VRB cell. Electrolyte flow rate and stack channel dimension are proved to be the critical factors affecting flow distribution and cell performance. The model developed in this paper can be employed to optimize both VRB stack design and system operation conditions. Further improvements of the model concerning current density and electrode properties are also suggested in the paper. - Highlights: • A coupled three-dimensional model of vanadium redox flow cell is proposed. • Interdigitated flow channels with two different manifold designs are simulated. • Manifold structure affects uniformity of distribution patterns significantly. • Increased electrolyte flow rate improves cell performance for both designs. • Decreased channel size and enlarged land width enhance cell voltage

Analysis of Poiseuille Flow Property in Two-Dimensional Mi-cro Channels of Microfluidic Pneumatic Micro-Valve

Science.gov (United States)

Yang, Shaohua; Long, Wei; Chen, Yajun

2018-03-01

In this paper, the control mechanism and mathematical description of the microfluidic flow in the microfluidic process of the PDMS membrane type pneumatic micro-valve were studied. The velocity and pressure variation law of the velocity field inside micro valve was analyzed by numerical simulation method. The influence of the two kinds of inlet drive modes on the working effect and the pressure flow characteristics of the pneumatic micro-valve was studied. The structure of the elastic solid valve diaphragm under the dual action of the airway and the liquid channel was analyzed. Deformation and stress distribution. The results show that the gas flow in the gas flow channel under the diaphragm by the vacuum part of the role of the formation of a suction gas vortex, pressure-driven mode was easier under the diaphragm to produce a strong gas vortex, resulting in internal and external pressure to promote diaphragm cut-off liquid channel; In the pressure pneumatic mode, the stress at both ends of the diaphragm was smaller, the membrane was not easy to tear failure.

Flow Patterns and Thermal Drag in a One-Dimensional Inviscid Channel with Heating or Cooling

Institute of Scientific and Technical Information of China (English)

无

1993-01-01

In this paper investigations on the flow patterns and the thermal drag phenomenon in one -dimensional inviscid channel flow with heating or cooling are described and discussed:expressions of flow rate ratio and thermal drag coefficient for different flow patterns and its physical mechanism are presented.

Numerical study on flow rate limitation of open capillary channel flow through a wedge

Directory of Open Access Journals (Sweden)

Ting-Ting Zhang

2016-04-01

Full Text Available The flow characteristics of slender-column flow in wedge-shaped channel under microgravity condition are investigated in this work. The one-dimensional theoretical model is applied to predict the critical flow rate and surface contour of stable flow. However, the one-dimensional model overestimates the critical flow rate for not considering the extra pressure loss. Then, we develop a three-dimensional simulation method with OpenFOAM, a computational fluid dynamics tool, to simulate various phenomena in wedge channels with different lengths. The numerical results are verified with the capillary channel flow experimental data on the International Space Station. We find that the three-dimensional simulation perfectly predicts the critical flow rates and surface contours under various flow conditions. Meanwhile, the general behaviors in subcritical, critical, and supercritical flow are studied in three-dimensional simulation considering variations of flow rate and open channel length. The numerical techniques for three-dimensional simulation is validated for a wide range of configurations and is hopeful to provide valuable guidance for capillary channel flow experiment and efficient liquid management in space.

Two-dimensional turbulent flows on a bounded domain

NARCIS (Netherlands)

Kramer, W.

2006-01-01

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

Flow channeling in a single fracture as a two-dimensional strongly heterogeneous permeable medium

International Nuclear Information System (INIS)

Tsang, Y.W.; Tsang, C.F.

1990-01-01

Recent interest in the evaluation of contaminant transport in bedrock aquifers and in the performance assessment of geologic nuclear waste repositories has motivated many studies of fluid flow and tracer transport in fractured rocks. Until recently, numerical modeling of fluid flow in the fractured medium commonly makes the assumption that each fracture may be idealized as a pair of parallel plates separated by a constant distance which represents the aperture of the fracture. More recent theoretical work has taken into account that the aperture in a real rock fracture in fact takes on a range of values. Evidence that flow in fractures tends to coalesce in preferred paths has been found in the field. Current studies of flow channeling in a fracture as a result of the variable apertures may also be applicable to flow and transport in a strongly heterogenous porous medium. This report includes the methodology used to study the flow channelling and tracer transport in a single fracture consisting of variable apertures. Relevant parameters that control flow channeling are then identified and the relationship of results to the general problem of flow in a heterogenous porous medium are discussed

Two-phase flow boiling pressure drop in small channels

International Nuclear Information System (INIS)

Sardeshpande, Madhavi V.; Shastri, Parikshit; Ranade, Vivek V.

2016-01-01

Highlights: • Study of typical 19 mm steam generator tube has been undertaken in detail. • Study of two phase flow boiling pressure drop, flow instability and identification of flow regimes using pressure fluctuations is the main focus of present work. • Effect of heat and mass flux on pressure drop and void fraction was studied. • Flow regimes identified from pressure fluctuations data using FFT plots. • Homogeneous model predicted pressure drop well in agreement. - Abstract: Two-phase flow boiling in small channels finds a variety of applications in power and process industries. Heat transfer, boiling flow regimes, flow instabilities, pressure drop and dry out are some of the key issues related to two-phase flow boiling in channels. In this work, the focus is on pressure drop in two-phase flow boiling in tubes of 19 mm diameter. These tubes are typically used in steam generators. Relatively limited experimental database is available on 19 mm ID tube. Therefore, in the present work, the experimental set-up is designed for studying flow boiling in 19 mm ID tube in such a way that any of the different flow regimes occurring in a steam generator tube (from pre-heating of sub-cooled water to dry-out) can be investigated by varying inlet conditions. The reported results cover a reasonable range of heat and mass flux conditions such as 9–27 kW/m 2 and 2.9–5.9 kg/m 2 s respectively. In this paper, various existing correlations are assessed against experimental data for the pressure drop in a single, vertical channel during flow boiling of water at near-atmospheric pressure. A special feature of these experiments is that time-dependent pressures are measured at four locations along the channel. The steady-state pressure drop is estimated and the identification of boiling flow regimes is done with transient characteristics using time series analysis. Experimental data and corresponding results are compared with the reported correlations. The results will be

Linear predictions of supercritical flow instability in two parallel channels

International Nuclear Information System (INIS)

Shah, M.

2008-01-01

A steady state linear code that can predict thermo-hydraulic instability boundaries in a two parallel channel system under supercritical conditions has been developed. Linear and non-linear solutions of the instability boundary in a two parallel channel system are also compared. The effect of gravity on the instability boundary in a two parallel channel system, by changing the orientation of the system flow from horizontal flow to vertical up-flow and vertical down-flow has been analyzed. Vertical up-flow is found to be more unstable than horizontal flow and vertical down flow is found to be the most unstable configuration. The type of instability present in each flow-orientation of a parallel channel system has been checked and the density wave oscillation type is observed in horizontal flow and vertical up-flow, while the static type of instability is observed in a vertical down-flow for the cases studied here. The parameters affecting the instability boundary, such as the heating power, inlet temperature, inlet and outlet K-factors are varied to assess their effects. This study is important for the design of future Generation IV nuclear reactors in which supercritical light water is proposed as the primary coolant. (author)

Sub-nanometre channels embedded in two-dimensional materials

KAUST Repository

Han, Yimo

2017-12-04

Two-dimensional (2D) materials are among the most promising candidates for next-generation electronics due to their atomic thinness, allowing for flexible transparent electronics and ultimate length scaling1. Thus far, atomically thin p–n junctions2,3,4,5,6,7,8, metal–semiconductor contacts9,10,11, and metal–insulator barriers12,13,14 have been demonstrated. Although 2D materials achieve the thinnest possible devices, precise nanoscale control over the lateral dimensions is also necessary. Here, we report the direct synthesis of sub-nanometre-wide one-dimensional (1D) MoS2 channels embedded within WSe2 monolayers, using a dislocation-catalysed approach. The 1D channels have edges free of misfit dislocations and dangling bonds, forming a coherent interface with the embedding 2D matrix. Periodic dislocation arrays produce 2D superlattices of coherent MoS2 1D channels in WSe2. Using molecular dynamics simulations, we have identified other combinations of 2D materials where 1D channels can also be formed. The electronic band structure of these 1D channels offers the promise of carrier confinement in a direct-gap material and the charge separation needed to access the ultimate length scales necessary for future electronic applications.

Three-dimensional numerical simulations of turbulent cavitating flow in a rectangular channel

Science.gov (United States)

Iben, Uwe; Makhnov, Andrei; Schmidt, Alexander

2018-05-01

Cavitation is a phenomenon of formation of bubbles (cavities) in liquid as a result of pressure drop. Cavitation plays an important role in a wide range of applications. For example, cavitation is one of the key problems of design and manufacturing of pumps, hydraulic turbines, ship's propellers, etc. Special attention is paid to cavitation erosion and to performance degradation of hydraulic devices (noise, fluctuations of the mass flow rate, etc.) caused by the formation of a two-phase system with an increased compressibility. Therefore, development of a model to predict cavitation inception and collapse of cavities in high-speed turbulent flows is an important fundamental and applied task. To test the algorithm three-dimensional simulations of turbulent flow of a cavitating liquid in a rectangular channel have been conducted. The obtained results demonstrate the efficiency and robustness of the formulated model and the algorithm.

A three-dimensional model of PEM fuel cells with serpentine flow channels

International Nuclear Information System (INIS)

Nguyen, P.T.; Berning, T.; Bang, M.; Djilali, N.

2003-01-01

A three-dimensional computational model of PEM fuel cell with serpentine flow field channels is presented in this paper. This model presents a comprehensive account for all important transport phenomena in fuel cell such as heat transfer, mass transfer, electrode kinetics, and potential fields in the membrane and gas diffusion layers. A new approach of solving for the potential losses across the cell was also developed in this model. The dependency of local current density on oxygen concentration and activation overpotential is fully addressed in this model. The computational domain consists of serpentine gas flow channels, porous gas diffusion layers, catalyst layers, and a membrane. Results obtained from this model are in good agreement with experimental results. (author)

Measurements of local two-phase flow parameters in a boiling flow channel

International Nuclear Information System (INIS)

Yun, Byong Jo; Park, Goon-CherI; Chung, Moon Ki; Song, Chul Hwa

1998-01-01

Local two-phase flow parameters were measured lo investigate the internal flow structures of steam-water boiling flow in an annulus channel. Two kinds of measuring methods for local two-phase flow parameters were investigated. These are a two-conductivity probe for local vapor parameters and a Pitot cube for local liquid parameters. Using these probes, the local distribution of phasic velocities, interfacial area concentration (IAC) and void fraction is measured. In this study, the maximum local void fraction in subcooled boiling condition is observed around the heating rod and the local void fraction is smoothly decreased from the surface of a heating rod to the channel center without any wall void peaking, which was observed in air-water experiments. The distributions of local IAC and bubble frequency coincide with those of local void fraction for a given area-averaged void fraction. (author)

Cross-equatorial flow through an abyssal channel under the complete Coriolis force: Two-dimensional solutions

Science.gov (United States)

Stewart, A. L.; Dellar, P. J.

The component of the Coriolis force due to the locally horizontal component of the Earth's rotation vector is commonly neglected, under the so-called traditional approximation. We investigate the role of this "non-traditional" component of the Coriolis force in cross-equatorial flow of abyssal ocean currents. We focus on the Antarctic Bottom Water (AABW), which crosses from the southern to the northern hemisphere through the Ceara abyssal plain in the western Atlantic ocean. The bathymetry in this region resembles a northwestward channel, connecting the Brazil Basin in the south to the Guyana Basin in the north. South of the equator, the AABW leans against the western continental rise, consistent with a northward flow in approximate geostrophic balance. The AABW then crosses to the other side of the abyssal channel as it crosses the equator, and flows into the northern hemisphere leaning towards the east against the Mid-Atlantic Ridge. The non-traditional component of the Coriolis force is strongest close to the equator. The traditional component vanishes at the equator, being proportional to the locally vertical component of the Earth's rotation vector. The weak stratification of the abyssal ocean, and subsequent small internal deformation radius, defines a relatively short characteristic horizontal lengthscale that tends to make non-traditional effects more prominent. Additionally, the steep gradients of the channel bathymetry induce large vertical velocities, which are linked to zonal accelerations by the non-traditional components of the Coriolis force. We therefore expect non-traditional effects to play a substantial role in cross-equatorial transport of the AABW. We present asymptotic steady solutions for non-traditional shallow water flow through an idealised abyssal channel, oriented at an oblique angle to the equator. The current enters from the south, leaning up against the western side of the channel in approximate geostrophic balance, and crosses the

Finite-time barriers to front propagation in two-dimensional fluid flows

Science.gov (United States)

Mahoney, John R.; Mitchell, Kevin A.

2015-08-01

Recent theoretical and experimental investigations have demonstrated the role of certain invariant manifolds, termed burning invariant manifolds (BIMs), as one-way dynamical barriers to reaction fronts propagating within a flowing fluid. These barriers form one-dimensional curves in a two-dimensional fluid flow. In prior studies, the fluid velocity field was required to be either time-independent or time-periodic. In the present study, we develop an approach to identify prominent one-way barriers based only on fluid velocity data over a finite time interval, which may have arbitrary time-dependence. We call such a barrier a burning Lagrangian coherent structure (bLCS) in analogy to Lagrangian coherent structures (LCSs) commonly used in passive advection. Our approach is based on the variational formulation of LCSs using curves of stationary "Lagrangian shear," introduced by Farazmand et al. [Physica D 278-279, 44 (2014)] in the context of passive advection. We numerically validate our technique by demonstrating that the bLCS closely tracks the BIM for a time-independent, double-vortex channel flow with an opposing "wind."

Sub-Nanometer Channels Embedded in Two-Dimensional Materials

KAUST Repository

Han, Yimo

2017-07-31

Two-dimensional (2D) materials are among the most promising candidates for next-generation electronics due to their atomic thinness, allowing for flexible transparent electronics and ultimate length scaling1. Thus far, atomically-thin p-n junctions2-7, metal-semiconductor contacts8-10, and metal-insulator barriers11-13 have been demonstrated. While 2D materials achieve the thinnest possible devices, precise nanoscale control over the lateral dimensions are also necessary. Although external one-dimensional (1D) carbon nanotubes14 can be used to locally gate 2D materials, this adds a non-trivial third dimension, complicating device integration and flexibility. Here, we report the direct synthesis of sub-nanometer 1D MoS2 channels embedded within WSe2 monolayers, using a dislocation-catalyzed approach. The 1D channels have edges free of misfit dislocations and dangling bonds, forming a coherent interface with the embedding 2D matrix. Periodic dislocation arrays produce 2D superlattices of coherent MoS2 1D channels in WSe2. Molecular dynamics (MD) simulations have identified other combinations of 2D materials that could form 1D channels. Density function theory (DFT) calculation predicts these 1D channels display type II band alignment needed for carrier confinement and charge separation to access the ultimate length scales necessary for future electronic applications.

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

Science.gov (United States)

Helman, James L.; Hesselink, Lambertus

1990-01-01

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

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

Science.gov (United States)

Bhateja, Ashish; Khakhar, Devang V.

2018-06-01

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

Effect of aspect ratio on relationship between flow resistance and flow regime of two-phase flow in rectangular channel

International Nuclear Information System (INIS)

Yan Chaoxing; Yan Changqi; Sun Licheng; Xing Dianchuan; Wang Yang

2013-01-01

On the basis of visual observation, the effects of aspect ratio on relationship between flow resistance and flow regime were investigated experimentally for two-phase flow in three rectangular channels with the same cross-section width of 43 mm and different heights of 1.41, 3 and 10 mm, respectively. According to the criteria in terms of restriction factor C o , the former two channels belong to narrow channel, whereas the last one is conventional channel. The experimental results show that the two-phase pressure drops in rectangular channel with different aspect ratios have different variation trends with the increase of the gas velocity. For the 10 mm channel, the gravitational pressure drop makes the major percentage of total pressure drop at low gas velocity while the frictional pressure drop is dominant for the 1.41 mm and 3 mm channels. With the increase of the gas flow rate, the frictional pressure drop contributes more to total pressure drop. The range of churn flow can be distinguished from its pressure drop characteristic in 10 mm channel. (authors)

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

International Nuclear Information System (INIS)

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

2014-01-01

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

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)

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

International Nuclear Information System (INIS)

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

1992-01-01

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

Two-phase flow instabilities in a vertical annular channel

Energy Technology Data Exchange (ETDEWEB)

Babelli, I.; Nair, S.; Ishii, M. [Purdue Univ., West Lafayette, IN (United States)

1995-09-01

An experimental test facility was built to study two-phase flow instabilities in vertical annular channel with emphasis on downward flow under low pressure and low flow conditions. The specific geometry of the test section is similar to the fuel-target sub-channel of the Savannah River Site (SRS) Mark 22 fuel assembly. Critical Heat Flux (CHF) was observed following flow excursion and flow reversal in the test section. Density wave instability was not recorded in this series of experimental runs. The results of this experimental study show that flow excursion is the dominant instability mode under low flow, low pressure, and down flow conditions. The onset of instability data are plotted on the subcooling-Zuber (phase change) numbers stability plane.

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

International Nuclear Information System (INIS)

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

2008-01-01

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

Three dimensional computation of turbulent flow in meandering channels

Energy Technology Data Exchange (ETDEWEB)

Van Thinh Nguyen

2000-07-01

In this study a finite element calculation procedure together with two-equation turbulent model k-{epsilon} and mixing length are applied to the problem of simulating 3D turbulent flow in closed and open meandering channels. Near the wall a special approach is applied in order to overcome the weakness of the standard k-{epsilon} in the viscous sub-layer. A specialized shape function is used in the special near wall elements to capture accurately the strong variations of the mean flow variables in the viscosity-affected near wall region. Based on the analogy of water and air flows, a few characteristics of hydraulic problems can be examined in aerodynamic models, respectively. To study the relationships between an aerodynamic and a hydraulic model many experiments have been carried out by Federal Waterway Engineering and Research Institute of Karlsruhe, Germany. In order to test and examine the results of these physical models, an appropriated numerical model is necessary. The numerical mean will capture the limitations of the experimental setup. The similarity and the difference between an aerodynamic and a hydraulic model will be found out by the results of numerical computations and will be depicted in this study. Despite the presence of similarities between the flow in closed channels and the flow in open channels, it should be stated that the presence of a free surface in the open channel introduces serious complications to three dimensional computation. A new unknown, which represents the position of nodes on this free surface, is introduced. A special approach is required for solving this unknown. A procedure surface tracking is applied to the free surface boundary like a moving boundary. Grid nodes on the free surface are free to move in such a way that they belong to the spines, which are the generator lines to define the allowed motion of the nodes on the free surface. (orig.) [German] Die numerische Simulation ist heute ein wichtiges Hilfsmittel fuer die

Non-Darcy behavior of two-phase channel flow.

Science.gov (United States)

Xu, Xianmin; Wang, Xiaoping

2014-08-01

We study the macroscopic behavior of two-phase flow in porous media from a phase-field model. A dissipation law is first derived from the phase-field model by homogenization. For simple channel geometry in pore scale, the scaling relation of the averaged dissipation rate with the velocity of the two-phase flow can be explicitly obtained from the model which then gives the force-velocity relation. It is shown that, for the homogeneous channel surface, Dacry's law is still valid with a significantly modified permeability including the contribution from the contact line slip. For the chemically patterned surfaces, the dissipation rate has a non-Darcy linear scaling with the velocity, which is related to a depinning force for the patterned surface. Our result offers a theoretical understanding on the prior observation of non-Darcy behavior for the multiphase flow in either simulations or experiments.

One-dimensional model of steady, compressible channel flow with mass, momentum, and energy addition

International Nuclear Information System (INIS)

Johnston, S.C.

1976-09-01

A one-dimensional model of steady, compressible channel flow with mass, momentum and energy addition is discussed. An exact solution to the governing equations was found and from it a similarity parameter relating dimensionless mass, momentum and energy addition identified. This similarity parameter is used to make two flows having different dimensionless mass, momentum and energy additions equivalent. Application of the similarity parameter to the LASL Intense Neutron Source experiment and the Sandia simulation of that experiment results in an expression relating the dimensionless mass addition of combustible gas required in the Sandia experiment to dimensionless energy addition in the LASL experiment. Results of the analysis indicate that the Sandia experiment can realistically simulate the energy addition in the LASL Intense Neutron Source experiment

Experimental study on downward two-phase flow in narrow rectangular channel

Energy Technology Data Exchange (ETDEWEB)

Kim, T.H.; Jeong, J.H. [Pusan National Univ., Busan (Korea, Republic of)

2014-07-01

Adiabatic vertical two-phase flow of air and water through narrow rectangular channels was investigated. This study involved the observation of flow using a high speed camera and flow regimes were determined by image processing program using a MATLAB. The flows regimes in channel with downward flow are similar to those found by previous studies with upward flow. The flow regimes in downward flow at low liquid velocity are different from the previous studies in upward flow. The flow regimes can be classified into bubbly, cap-bubbly, slug and churn flow. (author)

Transient two-dimensional flow in porous media

International Nuclear Information System (INIS)

Sharpe, L. Jr.

1979-01-01

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

Soap film flows: Statistics of two-dimensional turbulence

International Nuclear Information System (INIS)

Vorobieff, P.; Rivera, M.; Ecke, R.E.

1999-01-01

Soap film flows provide a very convenient laboratory model for studies of two-dimensional (2-D) hydrodynamics including turbulence. For a gravity-driven soap film channel with a grid of equally spaced cylinders inserted in the flow, we have measured the simultaneous velocity and thickness fields in the irregular flow downstream from the cylinders. The velocity field is determined by a modified digital particle image velocimetry method and the thickness from the light scattered by the particles in the film. From these measurements, we compute the decay of mean energy, enstrophy, and thickness fluctuations with downstream distance, and the structure functions of velocity, vorticity, thickness fluctuation, and vorticity flux. From these quantities we determine the microscale Reynolds number of the flow R λ ∼100 and the integral and dissipation scales of 2D turbulence. We also obtain quantitative measures of the degree to which our flow can be considered incompressible and isotropic as a function of downstream distance. We find coarsening of characteristic spatial scales, qualitative correspondence of the decay of energy and enstrophy with the Batchelor model, scaling of energy in k space consistent with the k -3 spectrum of the Kraichnan endash Batchelor enstrophy-scaling picture, and power-law scalings of the structure functions of velocity, vorticity, vorticity flux, and thickness. These results are compared with models of 2-D turbulence and with numerical simulations. copyright 1999 American Institute of Physics

Proton transport in a membrane protein channel: two-dimensional infrared spectrum modeling.

NARCIS (Netherlands)

Liang, C.; Knoester, J.; Jansen, T.L.Th.A.

2012-01-01

We model the two-dimensional infrared (2DIR) spectrum of a proton channel to investigate its applicability as a spectroscopy tool to study the proton transport process in biological systems. Proton transport processes in proton channels are involved in numerous fundamental biochemical reactions.

Experimental and numerical characterization of the water flow in spacer-filled channels of spiral-wound membranes

KAUST Repository

Bucs, Szilard; Valladares Linares, Rodrigo; Marston, Jeremy O.; Radu, Andrea I.; Vrouwenvelder, Johannes S.; Picioreanu, Cristian

2015-01-01

Micro-scale flow distribution in spacer-filled flow channels of spiral-wound membrane modules was determined with a particle image velocimetry system (PIV), aiming to elucidate the flow behaviour in spacer-filled flow channels. Two-dimensional water

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

International Nuclear Information System (INIS)

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

1996-01-01

Particle image velocimetry (PIV) is a nonintrusive measurement technique that can be used to study the structure of various fluid flows. PIV is used to measure the time-varying, full-field velocity data of a particle-seeded flow field within either a two-dimensional plane or three-dimensional volume. PIV is a very efficient measurement technique since it can obtain both qualitative and quantitative spatial information about the flow field being studied. The quantitative spatial velocity information can be further processed into information of flow parameters such as vorticity and turbulence over extended areas. The objective of this study was to apply recent advances and improvements in the PIV flow measurement technique to the full-field, nonintrusive analysis of a three-dimensional, two-phase fluid flow system in such a manner that both components of the two-phase system could be experimentally quantified

Analysis of two-phase flow and boiling heat transfer in inclined channel of core-catcher

International Nuclear Information System (INIS)

Tahara, M.; Suzuki, Y.; Abe, N.; Kurita, T.; Hamazaki, R.; Kojima, Y.

2008-01-01

Passive Corium Cooling System (CCS) provides a function of ex-vessel debris cooling and molten core stabilization during a severe accident. CCS features inclined cooling channels arranged axi-symmetrically below the core-catcher basin. In order to estimate the coolability of the inclined cooling channel, it is indispensable to identify the flow pattern of the two-phase flow in the cooling channel. Several former studies for the two-phase flow pattern in the inclined channel are referred. Taitel and Dukler (1976) developed a prediction method of the flow pattern transition in horizontal and near horizontal tubes. Barnea et al. (1980) showed the flow pattern map of upward flow with 10 degrees inclination. Sakaguti et al. (1996) observed the two-phase flow patterns in the horizontal pipe connected with slightly upward pipe, in which the flow pattern in the pipe with a bending part was expressed by the combination of a basic flow pattern and some auxiliary flow patterns. Then we investigated these studies In order to identify the flow patterns observed in the inclined cooling channel of CCS. Furthermore we experimentally observed the flow patterns in the inclined cooling channel with various inlet conditions. As a result of the investigation and observation, typical flow patterns in the inclined cooling channel were identified. Two typical flow patterns were observed depending on the steam flow rate, one of which is 'elongated bubble 'flow, and the other is 'churn with collapsing backward and upward slug 'flow The flow and heat transfer in the inclined channel of CCS is analyzed by using a two-phase analysis code employing two-fluid model in which the constitutive equations for the two-phase flow in inclined channels are incorporated. That is, drift flux parameter for each of the elongated bubble flow, and the churn with collapsing backward and upward slug flow are incorporated to the two-phase analysis code, which are based on the rising velocity of the long bubble in

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.

Numerical Simulation of Turbulent Half-corrugated Channel Flow by Hydrophilic and Hydrophobic Surfaces

Directory of Open Access Journals (Sweden)

M. R. Rastan

2018-03-01

Full Text Available In the first part of the present study, a two dimensional half-corrugated channel flow is simulated at Reynolds number of 104, in no-slip condition (hydrophilic surfaces( using various low Reynolds turbulence models as well as standard k-ε model; and an appropriate turbulence model (k-ω 1998 model( is proposed. Then, in order to evaluate the proposed solution method in simulation of flow adjacent to hydrophobic surfaces, turbulent flow is simulated in simple channel and the results are compared with the literature. Finally, two dimensional half-corrugated channel flow at Reynolds number of 104 is simulated again in vicinity of hydrophobic surfaces for varoius slip lengths. The results show that this method is capable of drag reduction in such a way that an increase of 200 μm in slip length leads to a massive drag reduction up to 38%. In addition, to access a significant drag reduction in turbulent flows, the non-dimensionalized slip length should be larger than the minimum.

Experimental investigation on flow patterns of gas-liquid two-phase upward flow through packed channel with spheres

International Nuclear Information System (INIS)

Zhang Nan; Sun Zhongning; Zhao Zhongnan

2011-01-01

Experiments of visualized two-phase upward flow were conducted in the packed channel, which filled with 3, 5, 8 mm in diameter of glass sphere respectively. The gas superficial velocity ranges from 0.005 to 1.172 m/s. The liquid superficial velocity ranges from 0.004 to 0.093 m/s. Four representative flow patterns were observed as bubbly flow, cluster flow, liquid-pulse flow and churn-pulse flow, and corresponding flow pattern maps were also presented. It is found that the pulse flow region is dominant. The comparisons of flow pattern map between packed channel and non-packed channel show that the bubbly flow region in packed channel is narrower than that of non-packed channel due to the packing. The comparisons of flow pattern maps for three different packing sizes show that the cluster flow region expands with the increase of the packing diameter. In the low liquid superficial velocity, the cluster flow directly changes to churn-pulse flow in the packed channel with 8 mm packing. (authors)

Three dimensional numerical investigation of flow mixing in curved tubes and mass transfer in T-channel junction

International Nuclear Information System (INIS)

Pandey, Pradeep; Nayak, A.K.; Vijayan, P.K.

2014-01-01

Three dimensional flow patterns appearing in geometries such as curved pipes and T-channel junctions have important applications and are attractive for research. Unlike the flow in a straight tube, fluid motion in a curved tube is not parallel to the axis of bend, owing to the presence of centrifugal effects. It is characterized by a secondary flow in a cross-sectional plane normal to the main flow. Consequently, secondary flow separation near the inner wall is observed in the developing region. The strength of the secondary flow is greatly influenced by the curvature ratio and in turn, a non-dimensional parameter called the Dean Number. Secondary flow increases flow resistance, resulting in a larger pressure drop along the bend. The location of the maximum axial velocity gets shifted towards the outer wall. Flow in a T-channel junction is also a configuration of great significance. The simulations of the present work show that flow at low Reynolds numbers (Re ≤ 115) is steady and symmetric. For low Reynolds numbers, flow in the downstream channel remains highly segregated about the centerline. The appearance of vortices in the T-channel junction does little to redistribute concentration when flow remains symmetric. With increasing Reynolds number, transition takes place towards asymmetric flow. The incoming flow field gets redistributed at the center-plane and the dividing streamline becomes increasingly distorted. The flow field is characterized by thin elongated fluid interfaces across which momentum diffusion takes place. Flow at higher Reynolds numbers (Re ≥ 250) becomes unsteady in which unstable stagnation stream traces move periodically leftward and rightward at top and bottom walls. Trajectories of mass-less particles show greater dwelling in the junction as compared to those of finite mass particle. The numerical simulation is carried out in the present work using ANUPRAVAHA, a general purpose CFD solver developed at IIT Kanpur in collaboration with

Numerical study on the two-dimensional flows of plasma and ionizing gas using trial particles

International Nuclear Information System (INIS)

Brushlinskij, K.V.; Kozlov, A.N.; Morozov, A.I.

1985-01-01

Two-dimensional flows of plasma and ionized αs in a channel between two coaxial electrodes are considered in the MHD-model with account of Hall effect. Stationary solutions of the problem on the flow are obtained either analytically in approximation of a ''smooth'' channel - for ideal conducting plasma, or numerically using the methos of establishment - in the ge-neral case of finite conductivity. A method of further numerical analysis of some peculiarities of flow is suggested in the paper. It is based on studying dynamics of single ''test'' particles in fields of the main MHD plasma flow. Trajectory of the test ion is calculated with account for interaction forces with earlier determined electromagentic field and friction responsible for Coulomb collisions with particles of the background flow. The calculations display trajectories of test particles with different masses, initial positions and initial rates. They are shown to be dose to current lines of background medium in plasma of finite conductivity, that testified to the virtue of effectiveness of the MHD-model. In case of ideal conductivity trajectories of test and background particles can noticeably differ from one another. Stabilization effects of motion of particles accidentally knocked out from the flow and separation of pariticles of different mass by electromao.netic forces are considered

An experimental study of two-phase flow instability on two parallel channel with low steam quality

International Nuclear Information System (INIS)

Jiang Shengyao; Wu shaorong; Bo Jinhai; Yao Meisheng; Han Bing; Zhang Youjie

1988-01-01

An experimental result of two-phase flow instability on two parallel channel natural circulation with low steam quality is presented. The comparison of instability in the single channel and that in parallel channel is given. The effect of unequal inlet resistance coefficient and unequal power on the parallel channel instability is described and the behaviour of instability with equal exit steam quality in the two channel is investigated

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

Flow model for open-channel reach or network

Science.gov (United States)

Schaffranek, R.W.

1987-01-01

Formulation of a one-dimensional model for simulating unsteady flow in a single open-channel reach or in a network of interconnected channels is presented. The model is both general and flexible in that it can be used to simulate a wide range of flow conditions for various channel configurations. It is based on a four-point (box), implicit, finite-difference approximation of the governing nonlinear flow equations with user-definable weighting coefficients to permit varying the solution scheme from box-centered to fully forward. Unique transformation equations are formulated that permit correlation of the unknowns at the extremities of the channels, thereby reducing coefficient matrix and execution time requirements. Discharges and water-surface elevations computed at intermediate locations within a channel are determined following solution of the transformation equations. The matrix of transformation and boundary-condition equations is solved by Gauss elimination using maximum pivot strategy. Two diverse applications of the model are presented to illustrate its broad utility. (USGS)

Two-phase flow regimes for counter-current air-water flows in narrow rectangular channels

International Nuclear Information System (INIS)

Kim, Byong Joo; Sohn, Byung Hu; Jeong, Si Young

2001-01-01

A study of counter-current two-phase flow in narrow rectangular channels has been performed. Two-phase flow regimes were experimentally investigated in a 760 mm long and 100 mm wide test section with 2.0 and 5.0 mm gap widths. The resulting flow regime maps were compared with the existing transition criteria. The experimental data and the transition criteria of the models showed relatively good agreement. However, the discrepancies between the experimental data and the model predictions of the flow regime transition became pronounced as the gap width increased. As the gap width increased the transition gas superficial velocities increased. The critical void fraction for the bubbly-to-slug transition was observed to be about 0.25. The two-phase distribution parameter for the slug flow was larger for the narrower channel. The uncertainties in the distribution parameter could lead to a disagreement in slug-to-churn transition between the experimental findings and the transition criteria. For the transition from churn to annular flow the effect of liquid superficial velocity was found to be insignificant

Pressure drop of magnetohydrodynamic two-phase annular flow in rectangular channel

International Nuclear Information System (INIS)

Kumamaru, Hiroshige; Fujiwara, Yoshiki; Ogita, Kenji

1999-01-01

Numerical calculations have been performed on magnetohydrodynamic (MHD) two-phase annular flow in a rectangular channel with a small aspect ratio, i.e.a small ratio of the channel side perpendicular to the applied magnetic field and the side parallel to the field. Results of the present calculation agree nearly with Inoue et al.'s experimental results in the region of large liquid Reynolds numbers and large Hartmann numbers. Calculation results also show that the pressure drop ratio, i.e. the ratio of pressure drop of two-phase flow to that of single-phase flow under the same liquid flow rate and applied magnetic field, becomes lower than ∼0.02 for conditions of a fusion reactor plant. (author)

Experimental study for flow regime of downward air-water two-phase flow in a vertical narrow rectangular channel

Energy Technology Data Exchange (ETDEWEB)

Kim, T. H.; Yun, B. J.; Jeong, J. H. [Pusan National University, Geunjeong-gu, Busan (Korea, Republic of)

2015-05-15

Studies were mostly about flow in upward flow in medium size circular tube. Although there are great differences between upward and downward flow, studies on vertical upward flow are much more active than those on vertical downward flow in a channel. In addition, due to the increase of surface forces and friction pressure drop, the pattern of gas-liquid two-phase flow bounded to the gap of inside the rectangular channel is different from that in a tube. The downward flow in a rectangular channel is universally applicable to cool the plate type nuclear fuel in research reactor. The sub-channel of the plate type nuclear fuel is designed with a few millimeters. Downward air-water two-phase flow in vertical rectangular channel was experimentally observed. The depth, width, and length of the rectangular channel is 2.35 mm, 66.7 mm, and 780 mm, respectively. The test section consists of transparent acrylic plates confined within a stainless steel frame. The flow patterns of the downward flow in high liquid velocity appeared to be similar to those observed in previous studies with upward flow. In downward flow, the transition lines for bubbly-slug and slug-churn flow shift to left in the flow regime map constructed with abscissa of the superficial gas velocity and ordinate of the superficial liquid velocity. The flow patterns observed with downward flow at low liquid velocity are different from those with upward flow.

Integral transformation of the Navier-Stokes equations for laminar flow in channels of arbitrary two-dimensional geometry

International Nuclear Information System (INIS)

Perez Guerrero, Jesus Salvador

1995-01-01

Laminar developing flow in channels of arbitrary geometry was studied by solving the Navier-Stokes equations in the stream function-only formulation through the Generalized Integral Transform Technique (GITT). The stream function is expanded in an infinite system based on eigenfunctions obtained by considering solely the diffusive terms of the original formulation. The Navier-Stokes equations are transformed into an infinite system of ordinary differential equations, by using the transformation and inversion formulae. For computational purposes, the infinite series is truncated, according to an automatic error control procedure. The ordinary differential is solved through well-established scientific subroutines from widely available mathematical libraries. The classical problem of developing flow between parallel-plates is analysed first, as for both uniform and irrotational inlet conditions. The effect of truncating the duct length in the accuracy of the obtained solution is studied. A convergence analysis of the results obtained by the GITT is performed and compared with results obtained by finite difference and finite element methods, for different values of Reynolds number. The problem of flow over a backward-facing step then follows. Comparisons with experimental results in the literature indicate an excellent agreement. The numerical co-validation was established for a test case, and perfect agreement is reached against results considered as benchmarks in the recent literature. The results were shown to be physically more reasonable than others obtained by purely numerical methods, in particular for situations where three-dimensional effects are identified. Finally, a test problem for an irregular by shoped duct was studied and compared against results found in the literature, with good agreement and excellent convergence rates for the stream function field along the whole channel, for different values of Reynolds number. (author)

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.

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

International Nuclear Information System (INIS)

Naymik, T.G.

1978-01-01

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

Two dimensional MHD flows between porous boundaries

International Nuclear Information System (INIS)

Gratton, F.T.

1994-01-01

Similarity solutions of dissipative MHD equations representing conducting fluids injected through porous walls and flowing out in both directions from the center of the channel, are studied as a function of four non dimensional parameters, Reynolds number R e , magnetic Reynolds number R m , Alfvenic Mach number, M A , and pressure gradient coefficient, C. The effluence is restrained by an external magnetic field normal to the walls. When R m m >>1, the solution may model a collision of plasmas of astrophysical interest. In this case the magnetic field lines help to drive the outflow acting jointly with the pressure gradient. The law for C as a function of the other parameters is given for several asymptotic limits. (author). 3 refs, 6 figs

Two-dimensional heat flow apparatus

Science.gov (United States)

McDougall, Patrick; Ayars, Eric

2014-06-01

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

Analysis of Two Phase Natural Circulation Flow in the Cooling Channel of the PECS

Energy Technology Data Exchange (ETDEWEB)

Park, R. J; Ha, K. S; Rhee, B. W; Kim, H. Y [Korea Atomic Energy Research Institute, Daejeon (Korea, Republic of)

2014-10-15

Decay heat and sensible heat of the relocated and spread corium are removed by the natural circulation flow at the bottom and side wall of the core catcher and the top water cooling of the corium. The coolant in the inclined channel absorbs the decay heat and sensible heat transferred from the corium through the structure of the core catcher body and flows up to the pool as a two phase mixture. On the other hand, some of the pool water will flow into the inlet of the downcomer piping, and will flow into the inclined cooling channel of the core catcher by gravity. As shown in Fig. 1, the engineered cooling channel is designed to provide effective long-term cooling and stabilization of the corium mixture in the core catcher body while facilitating steam venting in the PECS. To maintain the integrity of the ex-vessel core catcher, however, it is necessary that the coolant be sufficiently circulated along the inclined cooling channel to avoid CHF (Critical Heat Flux) on the heating surface of the cooling channel. For this reason, a verification experiment on the cooling capability of the EU-APR1400 core catcher has been performed in the CE (Cooling Experiment)-PECS facility at KAERI. Preliminary simulations of two-phase natural circulation in the CE-PECS were performed to predict two-phase flow characteristics and to determine the natural circulation mass flow rate in the flow channel. In this study, simulations of two-phase natural circulation in a real core catcher of the PECS have been performed to determine the natural circulation mass flow rate in the flow channel using the RELAP5/MOD3 computer code.

Micromachined two dimensional resistor arrays for determination of gas parameters

NARCIS (Netherlands)

van Baar, J.J.J.; Verwey, Willem B.; Dijkstra, Mindert; Dijkstra, Marcel; Wiegerink, Remco J.; Lammerink, Theodorus S.J.; Krijnen, Gijsbertus J.M.; Elwenspoek, Michael Curt

A resistive sensor array is presented for two dimensional temperature distribution measurements in a micromachined flow channel. This allows simultaneous measurement of flow velocity and fluid parameters, like thermal conductivity, diffusion coefficient and viscosity. More general advantages of

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

Science.gov (United States)

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

2017-09-01

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

Kinematics and statistics of dense, slow granular flow through vertical channels

Science.gov (United States)

Ananda, K. S.; Moka, Sudheshna; Nott, Prabhu R.

We have investigated the flow of dry granular materials through vertical channels in the regime of dense slow flow using video imaging of the particles adjacent to a transparent wall. Using an image processing technique based on particle tracking velocimetry, the video movies were analysed to obtain the velocities of individual particles. Experiments were conducted in two- and three-dimensional channels. In the latter, glass beads and mustard seeds were used as model granular materials, and their translational velocities were measured. In the former, aluminium disks with a dark diametral stripe were used and their translational velocities and spin were measured. Experiments in the three-dimensional channels were conducted for a range of the channel width W, and for smooth and rough sidewalls. As in earlier studies, we find that shearing takes place predominantly in thin layers adjacent to the walls, while the rest of the material appears to move as a plug. However, there are large velocity fluctuations even in the plug, where the macroscopic deformation rate is negligibly small. The thickness of the shear layer, scaled by the particle diameter dp, increases weakly with W/dp. The experimental data for the velocity field are in good agreement with the Cosserat plasticity model proposed recently. We also measured the mean spin of the particles in the two-dimensional channel, and its deviation from half the vorticity. There is a clear, measurable deviation, which too is in qualitative agreement with the Cosserat plasticity model. The statistics of particle velocity and spin fluctuations in the two-dimensional channel were analysed by determining their probability distribution function, and their spatial and temporal correlation. They were all found to be broadly similar to previous observations for three-dimensional channels, but some differences are evident. The spatial correlation of the velocity fluctuations are much stronger in the two-dimensional channel, implying

Theoretical investigation of flow regime for boiling water two-phase flow in horizontal rectangular narrow channels

International Nuclear Information System (INIS)

Zhang Chunwei; Qiu Suizheng; Yan Mingyu; Wang Bulei; Nie Changhua

2005-01-01

The flow regime transition criteria for the boiling water two-phase flow in horizontal rectangular narrow channels (1 x 20 mm, 2 x 20 mm) were theoretically explored. The discernible flow patterns were bubble, intermittent slug, churn, annular and steam-water separation flow. By using two-fluid model, equations of conservation of momentum were established for the two-phase flow. New flow-regime criteria were obtained and agreed well with the experiment data. (authors)

Biased and flow driven Brownian motion in periodic channels

Science.gov (United States)

Martens, S.; Straube, A.; Schmid, G.; Schimansky-Geier, L.; Hänggi, P.

2012-02-01

In this talk we will present an expansion of the common Fick-Jacobs approximation to hydrodynamically as well as by external forces driven Brownian transport in two-dimensional channels exhibiting smoothly varying periodic cross-section. We employ an asymptotic analysis to the components of the flow field and to stationary probability density for finding the particles within the channel in a geometric parameter. We demonstrate that the problem of biased Brownian dynamics in a confined 2D geometry can be replaced by Brownian motion in an effective periodic one-dimensional potential ψ(x) which takes the external bias, the change of the local channel width, and the flow velocity component in longitudinal direction into account. In addition, we study the influence of the external force magnitude, respectively, the pressure drop of the fluid on the particle transport quantities like the averaged velocity and the effective diffusion coefficient. The critical ratio between the external force and pressure drop where the average velocity equals zero is identified and the dependence of the latter on the channel geometry is derived. Analytic findings are confirmed by numerical simulations of the particle dynamics in a reflection symmetric sinusoidal channel.

Flow of two stratified fluids in an open channel with addition of fluids along the channel length

International Nuclear Information System (INIS)

Gardner, G.C.

1980-01-01

It is shown that two stably stratified fluids flowing in an open channel have two critical flow conditions. The one at higher flowrates is equivalent to the choked flow condition of a single fluid over a broad-crested weir, when the Froude number is unity. The lower critical condition imposes restrictions, which define the system if fluids are added progressively along the channel length and the flowrates increase from low to high values. However, if the flowrate does not become sufficiently large to pass through the lower critical condition, this condition will then define a form of choking, which again determines the system. It is shown that an important special case, with the proportional flowrates of the two fluids kept constant, has an analytical solution in which the relative depths of the fluids is a constant along the channel. Other systems must be solved numerically. (orig.)

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

Three-dimensional features on oscillating microbubbles streaming flows

Science.gov (United States)

Rossi, Massimiliano; Marin, Alvaro G.; Wang, Cheng; Hilgenfeldt, Sascha; Kähler, Christian J.

2013-11-01

Ultrasound-driven oscillating micro-bubbles have been used as active actuators in microfluidic devices to perform manifold tasks such as mixing, sorting and manipulation of microparticles. A common configuration consists in side-bubbles, created by trapping air pockets in blind channels perpendicular to the main channel direction. This configuration results in bubbles with a semi-cylindrical shape that creates a streaming flow generally considered quasi two-dimensional. However, recent experiments performed with three-dimensional velocimetry methods have shown how microparticles can present significant three-dimensional trajectories, especially in regions close to the bubble interface. Several reasons will be discussed such as boundary effects of the bottom/top wall, deformation of the bubble interface leading to more complex vibrational modes, or bubble-particle interactions. In the present investigation, precise measurements of particle trajectories close to the bubble interface will be performed by means of 3D Astigmatic Particle Tracking Velocimetry. The results will allow us to characterize quantitatively the three-dimensional features of the streaming flow and to estimate its implications in practical applications as particle trapping, sorting or mixing.

Three-Dimensional, Numerical Investigation of Flow and Heat Transfer in Rectangular Channels Subject to Partial Blockage

KAUST Repository

Salama, Amgad; El-Amin, Mohamed; Sun, Shuyu

2014-01-01

Numerical simulation of flow and heat transfer in two adjacent channels is conducted with one of the channels partially blocked. This system simulates typical channels of a material testing reactor. The blockage is assumed due to the buckling of one of the channel plates inward along its width. The blockage ratio considered in this work is defined as the ratio between the cross-sectional area of the blocked and the unblocked channel. In this work, we consider a blockage ratio of approximately 40%. However, the blockage is different along the width of the channel, ranging from 0% at the end of the channel to 90% in the middle. The channel walls are sandwiching volumetric heat sources that vary spatially as chopped cosine functions. Interesting patterns are highlighted and investigated. The reduction in the flow area of one channel results in the flow redistributing among the two channels according to the changes in their hydraulic conductivities. The results of the numerical simulations show that the maximum wall temperature in the blocked channel is well below the boiling temperature at the operating pressure.

Three-Dimensional, Numerical Investigation of Flow and Heat Transfer in Rectangular Channels Subject to Partial Blockage

KAUST Repository

Salama, Amgad

2014-08-25

Numerical simulation of flow and heat transfer in two adjacent channels is conducted with one of the channels partially blocked. This system simulates typical channels of a material testing reactor. The blockage is assumed due to the buckling of one of the channel plates inward along its width. The blockage ratio considered in this work is defined as the ratio between the cross-sectional area of the blocked and the unblocked channel. In this work, we consider a blockage ratio of approximately 40%. However, the blockage is different along the width of the channel, ranging from 0% at the end of the channel to 90% in the middle. The channel walls are sandwiching volumetric heat sources that vary spatially as chopped cosine functions. Interesting patterns are highlighted and investigated. The reduction in the flow area of one channel results in the flow redistributing among the two channels according to the changes in their hydraulic conductivities. The results of the numerical simulations show that the maximum wall temperature in the blocked channel is well below the boiling temperature at the operating pressure.

A mathematical model of the flow and bed topography in curved channels

NARCIS (Netherlands)

Olesen, K.W.

1985-01-01

A two-dimensional horizontal mathematical model of the flow and bed topography in alluvial channel bends is presented. The applicability of the model is restricted to channels of which the width-depth ratio is large, the Froude number is small, bed load is dominant and grain sorting effects are

Development of One Dimensional Hyperbolic Coupled Solver for Two-Phase Flows

International Nuclear Information System (INIS)

Kim, Eoi Jin; Kim, Jong Tae; Jeong, Jae June

2008-08-01

The purpose of this study is a code development for one dimensional two-phase two-fluid flows. In this study, the computations of two-phase flow were performed by using the Roe scheme which is one of the upwind schemes. The upwind scheme is widely used in the computational fluid dynamics because it can capture discontinuities clearly such as a shock. And this scheme is applicable to multi-phase flows by the extension methods which were developed by Toumi, Stadtke, etc. In this study, the extended Roe upwind scheme by Toumi for two-phase flow was implemented in the one-dimensional code. The scheme was applied to a shock tube problem and a water faucet problem. This numerical method seems efficient for non oscillating solutions of two phase flow problems, and also capable for capturing discontinuities

Development of One Dimensional Hyperbolic Coupled Solver for Two-Phase Flows

Energy Technology Data Exchange (ETDEWEB)

Kim, Eoi Jin; Kim, Jong Tae; Jeong, Jae June

2008-08-15

The purpose of this study is a code development for one dimensional two-phase two-fluid flows. In this study, the computations of two-phase flow were performed by using the Roe scheme which is one of the upwind schemes. The upwind scheme is widely used in the computational fluid dynamics because it can capture discontinuities clearly such as a shock. And this scheme is applicable to multi-phase flows by the extension methods which were developed by Toumi, Stadtke, etc. In this study, the extended Roe upwind scheme by Toumi for two-phase flow was implemented in the one-dimensional code. The scheme was applied to a shock tube problem and a water faucet problem. This numerical method seems efficient for non oscillating solutions of two phase flow problems, and also capable for capturing discontinuities.

Flow around turbulence promoters in parallel channel, (2)

International Nuclear Information System (INIS)

Shiina, Yasuaki

1983-01-01

Effects of walls on shedding vortex in developed channel flow were investigated putting a cylinder at the center of channels or on a wall for the value of w/d from 2 to 4. Results were compared with the uniform flow result. When a cylinder was put at the center of the channels, non-dimensional frequency plotted against Reynolds number agreed with the uniform flow result at low Reynolds number. However, it increased rapidly with Reynolds number, then it lay considerably above the uniform flow results at high Reynolds number. When a cylinder was put on a wall, non-dimensional frequency was considerably lower than the uniform flow result in the cases of w/d = 3 and 4. In the case of w/d = 2, however, frequency was higher than the uniform flow result at high Reynolds number. In all cases in the present study, the transition Reynolds number increased with decrease in the value of w/d. These results indicate that the increase in shedding frequency was due to the shift in velocity distribution from Poiseuille parabora in the wake region, which obviously increased with Reynolds number and with decrease in channel width. (author)

Predicting transition in two- and three-dimensional separated flows

International Nuclear Information System (INIS)

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

2008-01-01

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

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

Energy Technology Data Exchange (ETDEWEB)

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

2015-12-01

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

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

International Nuclear Information System (INIS)

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

2015-01-01

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

Three dimensional peristaltic flow of hyperbolic tangent fluid in non-uniform channel having flexible walls

Directory of Open Access Journals (Sweden)

M. Ali Abbas

2016-03-01

Full Text Available In this present analysis, three dimensional peristaltic flow of hyperbolic tangent fluid in a non-uniform channel has been investigated. We have considered that the pressure is uniform over the whole cross section and the interial effects have been neglected. For this purpose we consider laminar flow under the assumptions of long wavelength (λ→∞ and creeping flow (Re→0 approximations. The attained highly nonlinear equations are solved with the help of Homotopy perturbation method. The influence of various physical parameters of interest is demonstrated graphically for wall tension, mass characterization, damping nature of the wall, wall rigidity, wall elastance, aspect ratio and the Weissenberg number. In this present investigation we found that the magnitude of the velocity is maximum in the center of the channel whereas it is minimum near the walls. Stream lines are also drawn to discuss the trapping mechanism for all the physical parameters. Comparison has also been presented between Newtonian and non-Newtonian fluid.

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

Science.gov (United States)

Park, S.; Iversen, J. D.

1984-01-01

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

Three-dimensional investigations of the threading regime in a microfluidic flow-focusing channel

Science.gov (United States)

Gowda, Krishne; Brouzet, Christophe; Lefranc, Thibault; Soderberg, L. Daniel; Lundell, Fredrik

2017-11-01

We study the flow dynamics of the threading regime in a microfluidic flow-focusing channel through 3D numerical simulations and experiments. Making strong filaments from cellulose nano-fibrils (CNF) could potentially steer to new high-performance bio-based composites competing with conventional glass fibre composites. CNF filaments can be obtained through hydrodynamic alignment of dispersed CNF by using the concept of flow-focusing. The aligned structure is locked by diffusion of ions resulting in a dispersion-gel transition. Flow-focusing typically refers to a microfluidic channel system where the core fluid is focused by the two sheath fluids, thereby creating an extensional flow at the intersection. In this study, threading regime corresponds to an extensional flow field generated by the water sheath fluid stretching the dispersed CNF core fluid and leading to formation of long threads. The experimental measurements are performed using optical coherence tomography (OCT) and 3D numerical simulations with OpenFOAM. The prime focus is laid on the 3D characteristics of thread formation such as wetting length of core fluid, shape, aspect ratio of the thread and velocity flow-field in the microfluidic channel.

Counter-current gas-liquid two-phase flow in a narrow rectangular channel

International Nuclear Information System (INIS)

Sohn, Byung Hu; Kim, Byong Joo

2000-01-01

A study of counter-current two-phase flow in a narrow rectangular channel has been performed. Two-phase flow patterns and void fractions were experimentally studied in a 760 mm long and 100 mm wide test section with 3.0 mm gap. The resulting data have been compared to previous transition criteria and empirical correlations. The comparison of experimental data to the transition criteria developed by Taitel and Barnea showed good agreement for the bubbly-to-slug transition. For the criteria of Mishima and Ishii to be applicable to the slug to churn transition, a new model seems to be needed for the accurate prediction of the distribution parameter for the counter-current flow in narrow rectangular channels. For the churn-to-annular transition the model of Taitel and Barnea was found to be close to the experimental data. However the model should be improved in conjunction with the channel geometry to accurately predict the counter-current flow limitation and flow transition. It was verified the distribution parameter was well-correlated by the drift-flux model. The distribution parameter for the present study was found to be about 1.2 for all flow regimes except 1.0 for an annular flow. (author)

Accretion torques due to three-dimensional channelled flows in magnetic cataclysmic variables

International Nuclear Information System (INIS)

Campbell, C.G.

1986-01-01

Angular momentum transfer due to three-dimensional magnetically channelled accretion flows in cataclysmic binaries is considered. The white dwarf experiences a torque due to the twist in that part of its magnetic field which interacts with the accretion stream. The channelling process can also enhance angular momentum exchange between the stream and the orbit by increasing the gravitational torques. The components of the accretion torque are calculated for an arbitrary static magnetic orientation of the white dwarf, and their variation with orientation is presented. For high inclinations of the accreting pole to the orbital plane the component of the accretion torque parallel to this plane can be comparable to its perpendicular component. It is shown that the parallel component of the torque is still significant relative to the perpendicular component if material links to the white dwarf's magnetic field well away from the L 1 region. (author)

Analytical and numerical solution of three-dimensional channel flow in presence of a sinusoidal fluid injection and a chemical reaction

Directory of Open Access Journals (Sweden)

Sahin Ahmed

2015-06-01

Full Text Available Modeling of three-dimensional channel flow in a chemically-reacting fluid between two long vertical parallel flat plates in the presence of a transverse magnetic field is presented. The stationary plate is subjected to a transverse sinusoidal injection velocity distribution while the uniformly moving plate is subjected to a constant suction and slip boundary conditions. Due to this type of injection velocity, the flow becomes three dimensional. Comparisons with previously published work are performed and the results are found to be in excellent agreement. An increase in the permeability/magnetic parameter is found to escalate the velocity near the plate in motion. Growing Reynolds number or magnetic parameter enhances the x-component and reduces the z-component of the skin-friction at the wall at rest. The acquired knowledge in our study can be used by designers to control MHD flow as suitable for certain applications which include laminar magneto-aerodynamics, materials processing and MHD propulsion thermo-fluid dynamics.

Statistical mechanics of two-dimensional and geophysical flows

International Nuclear Information System (INIS)

Bouchet, Freddy; Venaille, Antoine

2012-01-01

The theoretical study of the self-organization of two-dimensional and geophysical turbulent flows is addressed based on statistical mechanics methods. This review is a self-contained presentation of classical and recent works on this subject; from the statistical mechanics basis of the theory up to applications to Jupiter’s troposphere and ocean vortices and jets. Emphasize has been placed on examples with available analytical treatment in order to favor better understanding of the physics and dynamics. After a brief presentation of the 2D Euler and quasi-geostrophic equations, the specificity of two-dimensional and geophysical turbulence is emphasized. The equilibrium microcanonical measure is built from the Liouville theorem. Important statistical mechanics concepts (large deviations and mean field approach) and thermodynamic concepts (ensemble inequivalence and negative heat capacity) are briefly explained and described. On this theoretical basis, we predict the output of the long time evolution of complex turbulent flows as statistical equilibria. This is applied to make quantitative models of two-dimensional turbulence, the Great Red Spot and other Jovian vortices, ocean jets like the Gulf-Stream, and ocean vortices. A detailed comparison between these statistical equilibria and real flow observations is provided. We also present recent results for non-equilibrium situations, for the studies of either the relaxation towards equilibrium or non-equilibrium steady states. In this last case, forces and dissipation are in a statistical balance; fluxes of conserved quantity characterize the system and microcanonical or other equilibrium measures no longer describe the system.

Effect of the Presence of Semi-circular Cylinders on Heat Transfer From Heat Sources Placed in Two Dimensional Channel

Directory of Open Access Journals (Sweden)

Ahmed W. Mustava

2013-04-01

Full Text Available The effect of a semi-circular cylinders in a two dimensional channel on heat transfer by forced convection from two heat sources with a constant temperature has been studied numerically. Each channel contains two heat sources; one on the upper surface of the channel and the other on the lower surface of the channel. There is semi-circular cylinder under the source in upper surface and there is semi-circular cylinder above the source in lower surface. The location of the second heat source with its semi-cylinder has been changed and keeps the first source with its semi- cylinder at the same location. The flow and temperature field are studied numerically with different values of Reynolds numbers and for different spacing between the centers of the semi-cylinders. The laminar flow field is analyzed numerically by solving the steady forms of the two-dimensional incompressible Navier- Stokes and energy equations.  The Cartesian velocity components and pressure on a collocated (non-staggered grid are used as dependent variables in the momentum equations, which discretized by finite volume method, body fitted coordinates are used to represent the complex channel geometry accurately, and grid generation technique based on elliptic partial differential equations is employed. SIMPLE algorithm is used to adjust the velocity field to satisfy the conservation of mass.  The range of Reynolds number is (Re= 100 – 800 and the range of the spacing between the semi-cylinders is(1-4 and the Prandtl number is 0.7.The results showed that increasing the spacing between the semi-cylinders increases the average of Nusselt number of the first heat source for all Reynolds numbers. As well as the results show that the best case among the cases studied to enhance the heat transfer is when the second heat source and its semi-cylinder located on at the distance (S=1.5 from the first half of the cylinder and the Reynolds number is greater than (Re ≥ 400 because of the

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

International Nuclear Information System (INIS)

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

2007-01-01

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

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

Energy Technology Data Exchange (ETDEWEB)

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

2016-12-15

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

Topology optimization of Channel flow problems

DEFF Research Database (Denmark)

Gersborg-Hansen, Allan; Sigmund, Ole; Haber, R. B.

2005-01-01

function which measures either some local aspect of the velocity field or a global quantity, such as the rate of energy dissipation. We use the finite element method to model the flow, and we solve the optimization problem with a gradient-based math-programming algorithm that is driven by analytical......This paper describes a topology design method for simple two-dimensional flow problems. We consider steady, incompressible laminar viscous flows at low to moderate Reynolds numbers. This makes the flow problem non-linear and hence a non-trivial extension of the work of [Borrvall&Petersson 2002......]. Further, the inclusion of inertia effects significantly alters the physics, enabling solutions of new classes of optimization problems, such as velocity--driven switches, that are not addressed by the earlier method. Specifically, we determine optimal layouts of channel flows that extremize a cost...

Couple stress fluid flow in a rotating channel with peristalsis

Science.gov (United States)

Abd elmaboud, Y.; Abdelsalam, Sara I.; Mekheimer, Kh. S.

2018-04-01

This article describes a new model for obtaining closed-form semi-analytical solutions of peristaltic flow induced by sinusoidal wave trains propagating with constant speed on the walls of a two-dimensional rotating infinite channel. The channel rotates with a constant angular speed about the z - axis and is filled with couple stress fluid. The governing equations of the channel deformation and the flow rate inside the channel are derived using the lubrication theory approach. The resulting equations are solved, using the homotopy perturbation method (HPM), for exact solutions to the longitudinal velocity distribution, pressure gradient, flow rate due to secondary velocity, and pressure rise per wavelength. The effect of various values of physical parameters, such as, Taylor's number and couple stress parameter, together with some interesting features of peristaltic flow are discussed through graphs. The trapping phenomenon is investigated for different values of parameters under consideration. It is shown that Taylor's number and the couple stress parameter have an increasing effect on the longitudinal velocity distribution till half of the channel, on the flow rate due to secondary velocity, and on the number of closed streamlines circulating the bolus.

Modeling water droplet condensation and evaporation in DNS of turbulent channel flow

Energy Technology Data Exchange (ETDEWEB)

Russo, E; Kuerten, J G M; Geld, C W M van der [Department of Mechanical Engineering, Eindhoven University of Technology, PO Box 513, 5600 MB Eindhoven (Netherlands); Geurts, B J, E-mail: e.russo@tue.nl [Faculty EEMCS, University of Twente, PO Box 217, 7500 AE Enschede (Netherlands)

2011-12-22

In this paper a point particle model for two-way coupling in water droplet-laden incompressible turbulent flow of air is proposed. The model is based on conservation laws and semi-empirical correlations. It has been implemented and tested in a DNS code based for turbulent channel flow with an Eulerian-Lagrangian approach. The two-way coupling is investigated in terms of the effects of mass and heat transfer on the droplets distributions along the channel wall-normal direction and by comparison of the droplet temperature statistics with respect to the case without evaporation and condensation. A remarkable conclusion is that the presence of evaporating and condensing droplets results in an increase in the non-dimensional heat transfer coefficient of the channel flow represented by the Nusselt number.

Lumped parameter modeling of a two-phase thermal-hydraulic channel with interface tracking

International Nuclear Information System (INIS)

Jo, J.H.; Kaufman, J.M.; Ruger, C.J.; Stein, S.

1978-01-01

A nonhomogenous, thermal nonequilibrium model for one-dimensional two-phase flow in a heated channel has been formulated in lumped parameter form. The channel is divided into a variable number of flow regimes separated by moving interfaces. The model can be used to predict the behavior of a LWR core and both primary and secondary sides of a steam generator under transient conditions. (author)

Research on one-dimensional two-phase flow

International Nuclear Information System (INIS)

Adachi, Hiromichi

1988-10-01

In Part I the fundamental form of the hydrodynamic basic equations for a one-dimensional two-phase flow (two-fluid model) is described. Discussions are concentrated on the treatment of phase change inertial force terms in the equations of motion and the author's equations of motion which have a remarkable uniqueness on the following three points. (1) To express force balance of unit mass two-phase fluid instead of that of unit volume two-phase fluid. (2) To pick up the unit existing mass and the unit flowing mass as the unit mass of two-phase fluid. (3) To apply the kinetic energy principle instead of the momentum low in the evaluation of steady inertial force term. In these three, the item (1) is for excluding a part of momentum change or kinetic energy change due to mass change of the examined part of fluid, which is independent of force. The item (2) is not to introduce a phenomenological physical model into the evaluation of phase change inertial force term. And the item (3) is for correctly applying the momentum law taking into account the difference of representative velocities between the main flow fluid (vapor phase or liquid phase) and the phase change part of fluid. In Part II, characteristics of various kinds of high speed two-phase flow are clarified theoretically by the basic equations derived. It is demonstrated that the steam-water two-phase critical flow with violent flashing and the airwater two-phase critical flow without phase change can be described with fundamentally the same basic equations. Furthermore, by comparing the experimental data from the two-phase critical discharge test and the theoretical prediction, the two-phase discharge coefficient, C D , for large sharp-edged orifice is determined as the value which is not affected by the experimental facility characteristics, etc. (author)

Electro-optic tunable multi-channel filter in two-dimensional ferroelectric photonic crystals

International Nuclear Information System (INIS)

Fu, Yulan; Zhang, Jiaxiang; Hu, Xiaoyong; Gong, Qihuang

2010-01-01

An electro-optic tunable multi-channel filter is presented, which is based on a two-dimensional ferroelectric photonic crystal made of barium titanate. The filtering properties of the photonic crystal filter can be tuned by an applied voltage or by adjusting the structural parameters. The channel shifts about 30 nm under excitation of an applied voltage of 54.8 V. The influences of the structural disorders caused by the perturbations in the radius or the position of air holes on the filtering properties are also analyzed

Dynamics of vortex interactions in two-dimensional flows

DEFF Research Database (Denmark)

Juul Rasmussen, J.; Nielsen, A.H.; Naulin, V.

2002-01-01

The dynamics and interaction of like-signed vortex structures in two dimensional flows are investigated by means of direct numerical solutions of the two-dimensional Navier-Stokes equations. Two vortices with distributed vorticity merge when their distance relative to their radius, d/R-0l. is below...... a critical value, a(c). Using the Weiss-field, a(c) is estimated for vortex patches. Introducing an effective radius for vortices with distributed vorticity, we find that 3.3 ... is effectively producing small scale structures and the relation to the enstrophy "cascade" in developed 2D turbulence is discussed. The influence of finite viscosity on the merging is also investigated. Additionally, we examine vortex interactions on a finite domain, and discuss the results in connection...

THERMOSS: a thermohydraulic model of flow stagnation in a horizontal fuel channel

International Nuclear Information System (INIS)

Gulshani, P.; Caplan, M.Z.; Spinks, N.J.

1984-01-01

Following a postulated inlet-side small break in the CANDU reactor, emergency coolant is injected to refull the horizontal fuel channels and remove the decay heat. As part of the accident analysis, the effects of loss of forced circulation during the accident are predicted. A break size exists for which, at the end of pump rundown, the break force balances the natural circulation force and the channel flow is reduced to near zero. The subcooled, stagnant channel condition is referred to as the standing-start condition. Subsequently, the channel coolant boils and stratifies. Eventually the steam flow from the channel heats up the endfitting to the saturation temperature and reaches the vertical feeder. The resulting buoyancy-induced flow then refills the channel. One dimensional, two-fluid conservation equations are solved in closed form to predict the duration of stagnation. In this calculation the channel water level is an important intermediate variable because it determines the amount of steam production

Noise control of subsonic cavity flows using plasma actuated receptive channels

International Nuclear Information System (INIS)

Gupta, Arnob Das; Roy, Subrata

2014-01-01

We introduce a passive receptive rectangular channel at the trailing edge of an open rectangular cavity to reduce the acoustic tones generated due to coherent shear layer impingement. The channel is numerically tested at Mach 0.3 using an unsteady three-dimensional large eddy simulation. Results show reduction in pressure fluctuations in the cavity due to which sound pressure levels are suppressed. Two linear dielectric barrier discharge plasma actuators are placed inside the channel to enhance the flow through it. Specifically, acoustic suppression of 7 dB was obtained for Mach 0.3 flow with the plasma actuated channel. Also, the drag coefficient for the cavity reduced by over three folds for the channel and over eight folds for the plasma actuated channel. Such a channel can be useful in noise and drag reduction for various applications, including weapons bay, landing gear and branched piping systems. (fast track communication)

VOF modelling of gas–liquid flow in PEM water electrolysis cell micro-channels

DEFF Research Database (Denmark)

Lafmejani, Saeed Sadeghi; Olesen, Anders Christian; Kær, Søren Knudsen

2017-01-01

In this study, the gaseliquid flow through an interdigitated anode flow field of a PEM water electrolysis cell (PEMEC) is analysed using a three-dimensional, transient, computational fluid dynamics (CFD) model. To account for two-phase flow, the volume of fluid (VOF) method in ANSYS Fluent 17...... of the channel. The model is capable of revealing effect of different bubble shapes/lengths in the outgoing channel. Shape and the sequence of the bubbles affect the water flow distribution in the ATL. The model presented in this work is the first step in the development of a comprehensive CFD model...

Derivation of a correlation for Drag coefficient in two-dimensional bounded supercavitating flows, using artificial neural networks

Energy Technology Data Exchange (ETDEWEB)

Shafaghat, R.; Hosseinalipour, S.M.; Derakhshani, S.M.E. [Iran University of Science and Technology, Department of Mechanical Engineering, Tehran (Iran)

2010-07-15

Artificial neural networks (ANNs) are used as a new approach for the determination of the relations between drag coefficient and Cavitation Number with cavity geometry in supercavitating flows which have been most widely used in the hydrodynamics researches. Also the result of the ANNs as a cost function potentially will be used in an optimization algorithm. Instead of complex differential equations and limited experimental data, faster and simpler solutions were obtained using equations derived from the ANN model. For training of the ANN the numerical results are used that are obtained from a boundary element method (BEM). At this problem, a two-dimensional supercavitation potential inviscid flow pasts a symmetric two-dimensional cavitator, which is placed perpendicular to the flow in a channel of infinite width and immediately a cavity is formed behind the cavitator. It was found that the coefficient of multiple determination (R{sup 2}-value) between the actual and ANN predicted data is equal to about 0.9998 for the drag coefficient and Cavitation number. As seen from the obtained results, the calculated cavity geometry for all drag coefficients and Cavitation Numbers are obviously within acceptable limits. (orig.)

Experimental and numerical characterization of the water flow in spacer-filled channels of spiral-wound membranes

KAUST Repository

Bucs, Szilard

2015-09-25

Micro-scale flow distribution in spacer-filled flow channels of spiral-wound membrane modules was determined with a particle image velocimetry system (PIV), aiming to elucidate the flow behaviour in spacer-filled flow channels. Two-dimensional water velocity fields were measured in a flow cell (representing the feed spacer-filled flow channel of a spiral wound reverse osmosis membrane module without permeate production) at several planes throughout the channel height. At linear flow velocities (volumetric flow rate per cross-section of the flow channel considering the channel porosity, also described as crossflow velocities) used in practice (0.074 and 0.163 m∙s-1) the recorded flow was laminar with only slight unsteadiness in the upper velocity limit. At higher linear flow velocity (0.3 m∙s-1) the flow was observed to be unsteady and with recirculation zones. Measurements made at different locations in the flow cell exhibited very similar flow patterns within all feed spacer mesh elements, thus revealing the same hydrodynamic conditions along the length of the flow channel. Three-dimensional (3-D) computational fluid dynamics simulations were performed using the same geometries and flow parameters as the experiments, based on steady laminar flow assumption. The numerical results were in good agreement (0.85-0.95 Bray-Curtis similarity) with the measured flow fields at linear velocities of 0.074 and 0.163 m∙s-1, thus supporting the use of model-based studies in the optimization of feed spacer geometries and operational conditions of spiral wound membrane systems.

Experimental and numerical characterization of the water flow in spacer-filled channels of spiral-wound membranes.

Science.gov (United States)

Bucs, Szilard S; Linares, Rodrigo Valladares; Marston, Jeremy O; Radu, Andrea I; Vrouwenvelder, Johannes S; Picioreanu, Cristian

2015-12-15

Micro-scale flow distribution in spacer-filled flow channels of spiral-wound membrane modules was determined with a particle image velocimetry system (PIV), aiming to elucidate the flow behaviour in spacer-filled flow channels. Two-dimensional water velocity fields were measured in a flow cell (representing the feed spacer-filled flow channel of a spiral wound reverse osmosis membrane module without permeate production) at several planes throughout the channel height. At linear flow velocities (volumetric flow rate per cross-section of the flow channel considering the channel porosity, also described as crossflow velocities) used in practice (0.074 and 0.163 m·s(-1)) the recorded flow was laminar with only slight unsteadiness in the upper velocity limit. At higher linear flow velocity (0.3 m·s(-1)) the flow was observed to be unsteady and with recirculation zones. Measurements made at different locations in the flow cell exhibited very similar flow patterns within all feed spacer mesh elements, thus revealing the same hydrodynamic conditions along the length of the flow channel. Three-dimensional (3-D) computational fluid dynamics simulations were performed using the same geometries and flow parameters as the experiments, based on steady laminar flow assumption. The numerical results were in good agreement (0.85-0.95 Bray-Curtis similarity) with the measured flow fields at linear velocities of 0.074 and 0.163 m·s(-1), thus supporting the use of model-based studies in the optimization of feed spacer geometries and operational conditions of spiral wound membrane systems. Copyright © 2015 Elsevier Ltd. All rights reserved.

Laser--Doppler anemometry technique applied to two-phase dispersed flows in a rectangular channel

International Nuclear Information System (INIS)

Lee, S.L.; Srinivasan, J.

1979-01-01

A new optical technique using Laser--Doppler anemometry has been applied to the local measurement of turbulent upward flow of a dilute water droplet--air two-phase dispersion in a vertical rectangular channel. Individually examined were over 20,000 droplet signals coming from each of a total of ten transversely placed measuring points, the closest of which to the channel wall was 250 μ away from the wall. Two flows of different patterns due to different imposed flow conditions were investigated, one with and the other without a liquid film formed on the channel wall. Reported are the size and number density distribution and the axial and lateral velocity distributions for the droplets as well as the axial and lateral velocity distributions for the air

Advanced numerical methods for three dimensional two-phase flow calculations

Energy Technology Data Exchange (ETDEWEB)

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

1997-07-01

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

Advanced numerical methods for three dimensional two-phase flow calculations

International Nuclear Information System (INIS)

Toumi, I.; Caruge, D.

1997-01-01

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

3-Dimensional numerical study of cooling performance of a heat sink with air-water flow through mini-channel

Science.gov (United States)

Majumder, Sambit; Majumder, Abhik; Bhaumik, Swapan

2016-07-01

The present microelectronics market demands devices with high power dissipation capabilities having enhanced cooling per unit area. The drive for miniaturizing the devices to even micro level dimensions is shooting up the applied heat flux on such devices, resulting in complexity in heat transfer and cooling management. In this paper, a method of CPU processor cooling is introduced where active and passive cooling techniques are incorporated simultaneously. A heat sink consisting of fins is designed, where water flows internally through the mini-channel fins and air flows externally. Three dimensional numerical simulations are performed for large set of Reynolds number in laminar region using finite volume method for both developing flows. The dimensions of mini-channel fins are varied for several aspect ratios such as 1, 1.33, 2 and 4. Constant temperature (T) boundary condition is applied at heat sink base. Channel fluid temperature, pressure drop are analyzed to obtain best cooling option in the present study. It has been observed that as the aspect ratio of the channel decreases Nusselt number decreases while pressure drop increases. However, Nusselt number increases with increase in Reynolds number.

Tunable double-channel filter based on two-dimensional ferroelectric photonic crystals

International Nuclear Information System (INIS)

Jiang, Ping; Ding, Chengyuan; Hu, Xiaoyong; Gong, Qihuang

2007-01-01

A tunable double-channel filter is presented, which is based on a two-dimensional nonlinear ferroelectric photonic crystal made of cerium doped barium titanate. The filtering properties of the photonic crystal filter can be tuned by adjusting the defect structure or by a pump light. The influences of the structure disorders caused by the perturbations in the radius or the position of air holes on the filtering properties are also analyzed

Tunable double-channel filter based on two-dimensional ferroelectric photonic crystals

Energy Technology Data Exchange (ETDEWEB)

Jiang, Ping [State Key Laboratory for Mesoscopic Physics, Department of Physics, Peking University, Beijing 100871 (China); Ding, Chengyuan [State Key Laboratory for Mesoscopic Physics, Department of Physics, Peking University, Beijing 100871 (China); Hu, Xiaoyong [State Key Laboratory for Mesoscopic Physics, Department of Physics, Peking University, Beijing 100871 (China)]. E-mail: xiaoyonghu@pku.edu.cn; Gong, Qihuang [State Key Laboratory for Mesoscopic Physics, Department of Physics, Peking University, Beijing 100871 (China)]. E-mail: qhgong@pku.edu.cn

2007-04-02

A tunable double-channel filter is presented, which is based on a two-dimensional nonlinear ferroelectric photonic crystal made of cerium doped barium titanate. The filtering properties of the photonic crystal filter can be tuned by adjusting the defect structure or by a pump light. The influences of the structure disorders caused by the perturbations in the radius or the position of air holes on the filtering properties are also analyzed.

Investigating the Performance of One- and Two-dimensional Flood Models in a Channelized River Network: A Case Study of the Obion River System

Science.gov (United States)

Kalyanapu, A. J.; Dullo, T. T.; Thornton, J. C.; Auld, L. A.

2015-12-01

Obion River, is located in the northwestern Tennessee region, and discharges into the Mississippi River. In the past, the river system was largely channelized for agricultural purposes that resulted in increased erosion, loss of wildlife habitat and downstream flood risks. These impacts are now being slowly reversed mainly due to wetland restoration. The river system is characterized by a large network of "loops" around the main channels that hold water either from excess flows or due to flow diversions. Without data on each individual channel, levee, canal, or pond it is not known where the water flows from or to. In some segments along the river, the natural channel has been altered and rerouted by the farmers for their irrigation purposes. Satellite imagery can aid in identifying these features, but its spatial coverage is temporally sparse. All the alterations that have been done to the watershed make it difficult to develop hydraulic models, which could predict flooding and droughts. This is especially true when building one-dimensional (1D) hydraulic models compared to two-dimensional (2D) models, as the former cannot adequately simulate lateral flows in the floodplain and in complex terrains. The objective of this study therefore is to study the performance of 1D and 2D flood models in this complex river system, evaluate the limitations of 1D models and highlight the advantages of 2D models. The study presents the application of HEC-RAS and HEC-2D models developed by the Hydrologic Engineering Center (HEC), a division of the US Army Corps of Engineers. The broader impacts of this study is the development of best practices for developing flood models in channelized river systems and in agricultural watersheds.

Numerical prediction of augmented turbulent heat transfer in an annular fuel channel with repeated two-dimensional square ribs

International Nuclear Information System (INIS)

Takase, K.

1996-01-01

The square-ribbed fuel rod for high temperature gas-cooled reactors was designed and developed so as to enhance the turbulent heat transfer in comparison with the previous standard fuel rod. The turbulent heat transfer characteristics in an annular fuel channel with repeated two-dimensional square ribs were analysed numerically on a fully developed incompressible flow using the k-ε turbulence model and the two-dimensional axisymmetrical coordinate system. Numerical analyses were carried out under the conditions of Reynolds numbers from 3000 to 20000 and ratios of square-rib pitch to height of 10, 20 and 40 respectively. The predictions of the heat transfer coefficients agreed well within an error of 10% for the square-rib pitch to height ratio of 10, 20% for 20 and 25% for 40 respectively, with the heat transfer empirical correlations obtained from the experimental data due to the simulated square-ribbed fuel rods. Therefore it was found that the effect of heat transfer augmentation due to the square ribs could be predicted by the present numerical simulations and the mechanism could be explained by the change in the turbulence kinematic energy distribution along the flow direction. (orig.)

Two-dimensional hydrodynamic modeling to quantify effects of peak-flow management on channel morphology and salmon-spawning habitat in the Cedar River, Washington

Science.gov (United States)

Czuba, Christiana; Czuba, Jonathan A.; Gendaszek, Andrew S.; Magirl, Christopher S.

2010-01-01

The Cedar River in Washington State originates on the western slope of the Cascade Range and provides the City of Seattle with most of its drinking water, while also supporting a productive salmon habitat. Water-resource managers require detailed information on how best to manage high-flow releases from Chester Morse Lake, a large reservoir on the Cedar River, during periods of heavy precipitation to minimize flooding, while mitigating negative effects on fish populations. Instream flow-management practices include provisions for adaptive management to promote and maintain healthy aquatic habitat in the river system. The current study is designed to understand the linkages between peak flow characteristics, geomorphic processes, riverine habitat, and biological responses. Specifically, two-dimensional hydrodynamic modeling is used to simulate and quantify the effects of the peak-flow magnitude, duration, and frequency on the channel morphology and salmon-spawning habitat. Two study reaches, representative of the typical geomorphic and ecologic characteristics of the Cedar River, were selected for the modeling. Detailed bathymetric data, collected with a real-time kinematic global positioning system and an acoustic Doppler current profiler, were combined with a LiDAR-derived digital elevation model in the overbank area to develop a computational mesh. The model is used to simulate water velocity, benthic shear stress, flood inundation, and morphologic changes in the gravel-bedded river under the current and alternative flood-release strategies. Simulations of morphologic change and salmon-redd scour by floods of differing magnitude and duration enable water-resource managers to incorporate model simulation results into adaptive management of peak flows in the Cedar River. PDF version of a presentation on hydrodynamic modelling in the Cedar River in Washington state. Presented at the American Geophysical Union Fall Meeting 2010.

Approximate characteristics for one-dimensional two-phase flows

International Nuclear Information System (INIS)

Sarayloo, A.; Peddleson, J.

1985-01-01

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

Numerical prediction of turbulent heat transfer augmentation in an annular fuel channel with two-dimensional square ribs

International Nuclear Information System (INIS)

Takase, Kazuyuki

1996-01-01

The square-ribbed fuel rod for high temperature gas-cooled reactors was developed in order to enhance the turbulent heat transfer in comparison with the standard fuel rod. To evaluate the heat transfer performance of the square-ribbed fuel rod, the turbulent heat transfer coefficients in an annular fuel channel with repeated two-dimensional square ribs were analyzed numerically on a fully developed incompressible flow using the k - ε turbulence model and the two-dimensional axisymmetrical coordinate system. Numerical analyses were carried out for a range of Reynolds numbers from 3000 to 20000 and ratios of square-rib pitch to height of 10, 20 and 40, respectively. The predicted values of the heat transfer coefficients agreed within an error of 10% for the square-rib pitch to height ratio of 10, 20% for 20 and 25% for 40, respectively, with the heat transfer empirical correlations obtained from the experimental data. It was concluded by the present study that the effect of the heat transfer augmentation by square ribs could be predicted sufficiently by the present numerical simulations and also a part of its mechanism could be explained by means of the change in the turbulence kinematic energy distribution along the flow direction. (author)

An experimental investigation of flow instability between two heated parallel channels with supercritical water

Energy Technology Data Exchange (ETDEWEB)

Xi, Xi; Xiao, Zejun, E-mail: fabulous_2012@sina.com; Yan, Xiao; Li, Yongliang; Huang, Yanping

2014-10-15

Highlights: • Flow instability experiment between two heated channels with supercritical water is carried out. • Two kinds of out of phase flow instability are found and instability boundaries under different working conditions are obtained. • Dynamics characteristics of flow instability are analyzed. - Abstract: Super critical water reactor (SCWR) is the generation IV nuclear reactor in the world. Under normal operation, water enters SCWR from cold leg with a temperature of 280 °C and then leaves the core with a temperature of 500 °C. Due to the sharp change of temperature, there is a huge density change in the core, which could result in potential flow instability and the safety of reactor would be threatened consequently. So it is necessary to carry out relevant investigation in this field. An experimental investigation which concerns with out of phase flow instability between two heated parallel channels with supercritical water has been carried out in this paper. Due to two INCONEL 625 pipes with a thickness of 6.5 mm are adopted, more experimental results are attained. To find out the influence of axial power shape on the onset of flow instability, each heated channel is divided into two sections and the heating power of each section can be controlled separately. Finally the instability boundaries are obtained under different inlet temperatures, axial power shapes, total inlet mass flow rates and system pressures. The dynamics characteristics of out of phase oscillation are also analyzed.

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

International Nuclear Information System (INIS)

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

1997-01-01

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

Flow Structure and Channel Morphology at a Confluent-Meander Bend

Science.gov (United States)

Riley, J. D.; Rhoads, B. L.

2009-12-01

Flow structure and channel morphology in meander bends have been well documented. Channel curvature subjects flow through a bend to centrifugal acceleration, inducing a counterbalancing pressure-gradient force that initiates secondary circulation. Transverse variations in boundary shear stress and bedload transport parallel cross-stream movement of high velocity flow and determine spatial patterns of erosion along the outer bank and deposition along the inner bank. Laboratory experiments and numerical modeling of confluent-meander bends, a junction planform that develops when a tributary joins a meandering river along the outer bank of a bend, suggest that flow and channel morphology in such bends deviate from typical patterns. The purpose of this study is to examine three-dimensional (3-D) flow structure and channel morphology at a natural confluent-meander bend. Field data were collected in southeastern Illinois where Big Muddy Creek joins the Little Wabash River near a local maximum of curvature along an elongated meander loop. Measurements of 3-D velocity components were obtained with an acoustic Doppler current profiler (ADCP) for two flow events with differing momentum ratios. Channel bathymetry was also resolved from the four-beam depths of the ADCP. Analysis of velocity data reveals a distinct shear layer flanked by dual helical cells within the bend immediately downstream of the confluence. Flow from the tributary confines flow from the main channel along the inner part of the channel cross section, displacing the thalweg inward, limiting the downstream extent of the point bar, protecting the outer bank from erosion and enabling bar-building along this bank. Overall, this pattern of flow and channel morphology is quite different from typical patterns in meander bends, but is consistent with a conceptual model derived from laboratory experiments and numerical modeling.

Analytical study of narrow channel flow for a spallation target system design

Energy Technology Data Exchange (ETDEWEB)

Islam, Md. Shafiqul; Monde, Masanori [Saga Univ., Saga (Japan); Terada, Atsuhiko; Kinoshita, Hidetaka; Hino, Ryutaro [Japan Atomic Energy Research Inst., Tokai, Ibaraki (Japan). Tokai Research Establishment

2001-07-01

Heat transfer and pressure drop characteristics under fully developed turbulent water flow condition were analyzed over a two-dimensional narrow rectangular channel whose height is H=1.2 mm. The channel configuration and water flow condition simulate forced convection cooling of a spallation target system components design such as a solid target and a proton beam window. The high-Reynolds number form of the standard k - {epsilon} and RNG k - {epsilon} models employing wall functions for the Reynolds number (Re) range of 7,000 to 22,000 were used in the analyses. As for heat transfer characteristics of a smooth channel, the Nusselt number obtained by the standard k - {epsilon} model agreed very well with the Dittus-Boelter correlation. No significant differences in friction factors for the smooth channel were observed for these two models, which agreed well with the Blasius correlation. However, the standard k - {epsilon} model could not predict friction factors well for the rib-roughened channel. (author)

Analytical study of narrow channel flow for a spallation target system design

International Nuclear Information System (INIS)

Islam, Md. Shafiqul; Monde, Masanori; Terada, Atsuhiko; Kinoshita, Hidetaka; Hino, Ryutaro

2001-07-01

Heat transfer and pressure drop characteristics under fully developed turbulent water flow condition were analyzed over a two-dimensional narrow rectangular channel whose height is H=1.2 mm. The channel configuration and water flow condition simulate forced convection cooling of a spallation target system components design such as a solid target and a proton beam window. The high-Reynolds number form of the standard k - ε and RNG k - ε models employing wall functions for the Reynolds number (Re) range of 7,000 to 22,000 were used in the analyses. As for heat transfer characteristics of a smooth channel, the Nusselt number obtained by the standard k - ε model agreed very well with the Dittus-Boelter correlation. No significant differences in friction factors for the smooth channel were observed for these two models, which agreed well with the Blasius correlation. However, the standard k - ε model could not predict friction factors well for the rib-roughened channel. (author)

Decaying quasi-two-dimensional viscous flow on a square domain

DEFF Research Database (Denmark)

Konijnenberg, J.A. van de; Flor, J.B.; Heijst, G.J.F. van

1998-01-01

A comparison is made between experimental, numerical and analytical results for the two-dimensional flow on a square domain. The experiments concern the flow at the interface of a two-layer stratified fluid, evoked by either stirring the fluid with a rake, or by injecting additional fluid...... at the interface. Two numerical simulations were performed with initial conditions and boundary conditions that correspond approximately with those met in the experiments. The analytical results concern the calculation of the lowest modes of a decaying Stokes flow on a square domain. At late times...... relationship between vorticity and stream function in the experiments and the simulations. (C) 1998 American Institute of Physics....

Dimensional reduction of a general advection–diffusion equation in 2D channels

Science.gov (United States)

Kalinay, Pavol; Slanina, František

2018-06-01

Diffusion of point-like particles in a two-dimensional channel of varying width is studied. The particles are driven by an arbitrary space dependent force. We construct a general recurrence procedure mapping the corresponding two-dimensional advection-diffusion equation onto the longitudinal coordinate x. Unlike the previous specific cases, the presented procedure enables us to find the one-dimensional description of the confined diffusion even for non-conservative (vortex) forces, e.g. caused by flowing solvent dragging the particles. We show that the result is again the generalized Fick–Jacobs equation. Despite of non existing scalar potential in the case of vortex forces, the effective one-dimensional scalar potential, as well as the corresponding quasi-equilibrium and the effective diffusion coefficient can be always found.

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

Two-dimensional analysis of motion artifacts, including flow effects

International Nuclear Information System (INIS)

Litt, A.M.; Brody, A.S.; Spangler, R.A.; Scott, P.D.

1990-01-01

The effects of motion on magnetic resonance images have been theoretically analyzed for the case of a point-like object in simple harmonic motion and for other one-dimensional trajectories. The authors of this paper extend this analysis to a generalized two-dimensional magnetization with an arbitrary motion trajectory. The authors provide specific solutions for the clinically relevant cases of the cross-sections of cylindrical objects in the body, such as the aorta, which has a roughly one-dimensional, simple harmonic motion during respiration. By extending the solution to include inhomogeneous magnetizations, the authors present a model which allows the effects of motion artifacts and flow artifacts to be analyzed simultaneously

Two dimensional modelling of flood flows and suspended sediment transport: the case of Brenta River

Science.gov (United States)

D'Alpaos, L.; Martini, P.; Carniello, L.

2003-04-01

The paper deals with numerical modelling of flood waves and suspended sediment in plain river basins. The two dimensional depth integrated momentum and continuity equations, modified to take into account of the bottom irregularities that strongly affect the hydrodynamic and the continuity in partially dry areas (for example, during the first stages of a plain flooding and in tidal flows), are solved with a standard Galerkin finite element method using a semi-implicit numerical scheme and considering the role both of the small channel network and the regulation dispositive on the flooding wave propagation. Transport of suspended sediment and bed evolution are coupled with the flood propagation through the convection-dispersion equation and the Exner's equation. Results of a real case study are presented in which the effects of extreme flood of Brenta River (Italy) are examinated. The flooded areas (urban and rural areas) are identified and a mitigation solution based on a diversion channel flowing into Venice Lagoon is proposed. We show that this solution strongly reduces the flood risk in the downstream areas and can provide an important sediment source to the Venice Lagoon. Finally, preliminary results of the sediment dispersion in the Venice Lagoon are presented.

Effect of flow field on open channel flow properties using numerical investigation and experimental comparison

Energy Technology Data Exchange (ETDEWEB)

Khazaee, I. [Department of Mechanical Engineering, Torbat-e-jam branch, Islamic Azad University, Torbat-e-jam (Iran, Islamic Republic of); Mohammadiun, M. [Department of Mechanical Engineering, Shahrood branch, Islamic Azad University, Shahrood (Iran, Islamic Republic of)

2012-07-01

In this paper a complete three-dimensional and two phase CFD model for flow distribution in an open channel investigated. The finite volume method (FVM) with a dynamic Sub grid-scale was carried out for seven cases of different aspect ratios, different inclination angles or slopes and convergence-divergence condition. The volume of fluid (VOF) method was used to allow the free-surface to deform freely with the underlying turbulence. The discharge through open channel flow is often evaluated by velocity-area integration method from the measurement of velocity at discrete locations in the measuring section. The variation of velocity along horizontal and vertical directions is thus very important to decide the location of the sensors. The aspect ratio of the channel, slope of the channel and divergence- convergence of the channel have investigated and the results show that the depth of water at the end of the channel is higher at AR=0.8 against the AR=0.4 and AR=1.2. Also it is clear that by increasing the inclination angle or slope of the channel in case1, case4 and case5 the depth of the water increases. Also it is clear that the outlet mass flow rate is at a minimum value at a range of inclination angle of the channel.

Field Investigation of Flow Structure and Channel Morphology at Confluent-Meander Bends

Science.gov (United States)

Riley, J. D.; Rhoads, B. L.

2007-12-01

The movement of water and sediment through drainage networks is inevitably influenced by the convergence of streams and rivers at channel confluences. These focal components of fluvial systems produce a complex hydrodynamic environment, where rapid changes in flow structure and sediment transport occur to accommodate the merging of separate channel flows. The inherent geometric and hydraulic change at confluences also initiates the development of distinct geomorphic features, reflected in the bedform and shape of the channel. An underlying assumption of previous experimental and theoretical models of confluence dynamics has been that converging streams have straight channels with angular configurations. This generalized conceptualization was necessary to establish confluence planform as symmetrical or asymmetrical and to describe subsequent flow structure and geomorphic features at confluences. However, natural channels, particularly those of meandering rivers, curve and bend. This property and observation of channel curvature at natural junctions have led to the hypothesis that natural stream and river confluences tend to occur on the concave outer bank of meander bends. The resulting confluence planform, referred to as a confluent-meander bend, was observed over a century ago but has received little scientific attention. This paper examines preliminary data on three-dimensional flow structure and channel morphology at two natural confluent-meander bends of varying size and with differing tributary entrance locations. The large river confluence of the Vermilion River and Wabash River in west central Indiana and the comparatively small junction of the Little Wabash River and Big Muddy Creek in southeastern Illinois are the location of study sites for field investigation. Measurements of time-averaged three-dimensional velocity components were obtained at these confluences with an acoustic Doppler current profiler for flow events with differing momentum ratios. Bed

Multi-dimensional two-fluid flow computation. An overview

International Nuclear Information System (INIS)

Carver, M.B.

1992-01-01

This paper discusses a repertoire of three-dimensional computer programs developed to perform critical analysis of single-phase, two-phase and multi-fluid flow in reactor components. The basic numerical approach to solving the governing equations common to all the codes is presented and the additional constitutive relationships required for closure are discussed. Particular applications are presented for a number of computer codes. (author). 12 refs

Waves in separated two-phase flow

International Nuclear Information System (INIS)

Pols, R.M.

1998-06-01

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

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

International Nuclear Information System (INIS)

Mendes de Moura, L.F.

1989-04-01

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

Two-dimensional magnetohydrodynamic equilibria with flow and studies of equilibria fluctuations

International Nuclear Information System (INIS)

Agim, Y.Z.

1989-08-01

A set of reduced ideal MHD equations is derived to investigate equilibria of plasmas with mass flow in general two-dimensional geometry. These equations provide a means of investigating the effects of flow on self-consistent equilibria in a number of new two-dimensional configurations such as helically symmetric configurations with helical axis, which are relevant to stellarators, as well as axisymmetric configurations. It is found that as in the axisymmetric case, general two-dimensional flow equilibria are governed by a second-order quasi-linear partial differential equation for a magnetic flux function, which is coupled to a Bernoulli-type equation for the density. The equation for the magnetic flux function becomes hyperbolic at certain critical flow speeds which follow from its characteristic equation. When the equation is hyperbolic, shock phenomena may exist. As a particular example, unidirectional flow along the lines of symmetry is considered. In this case, the equation mentioned above is always elliptic. An exact solution for the case of helically symmetric unidirectional flow is found and studied to determine flow effects on the magnetic topology. In second part of this thesis, magnetic fluctuations due to the thermally excited MHD waves are investigated using fluid and kinetic models to describe stable, uniform, compressible plasma in the range above the drift wave frequency and below the ion cyclotron frequency. It is shown that the fluid model with resistivity yields spectral densities which are roughly Lorentzian, exhibit equipartition with no apparent cutoff in wavenumber space and a Bohm-type diffusion coefficient. Under certain conditions, the ensuing transport may be comparable to classical values. For a phenomenological cutoff imposed on the spectrum, the typical fluctuating-to-equilibrium magnetic field ratio is found to be of the order of 10 -10

Simulation of the solidification in a channel of a water-cooled glass flow

Directory of Open Access Journals (Sweden)

G. E. Ovando Chacon

2014-12-01

Full Text Available A computer simulation study of a laminar steady-state glass flow that exits from a channel cooled with water is reported. The simulations are carried out in a two-dimensional, Cartesian channel with a backward-facing step for three different angles of the step and different glass outflow velocities. We studied the interaction of the fluid dynamics, phase change and thermal behavior of the glass flow due to the heat that transfers to the cooling water through the wall of the channel. The temperature, streamline, phase change and pressure fields are obtained and analyzed for the glass flow. Moreover, the temperature increments of the cooling water are characterized. It is shown that, by reducing the glass outflow velocity, the solidification is enhanced; meanwhile, an increase of the step angle also improves the solidification of the glass flow.

Investigation of one-dimensional interfacial area transport for vertical upward air–water two-phase flow in an annular channel at elevated pressures

International Nuclear Information System (INIS)

Ozar, B.; Brooks, C.S.; Euh, D.J.; Hibiki, T.; Ishii, M.

2013-01-01

Highlights: • Interfacial area transport equation (IATE) for a rectangular duct is modified for an annulus. • IATE predicts interfacial area transport in bubbly-to-churn flow. • Scalability of IATE to elevated pressure conditions is validated. • Detailed 1D interfacial area transport data are presented. • Detailed interfacial area transport mechanisms are discussed. -- Abstract: The interfacial area transport of vertical, upward, air–water two-phase flows in an annular channel has been investigated at different system pressures. The inner and outer diameters of the annular channel were 19.1 mm and 38.1 mm, respectively. Twenty three inlet flow conditions were selected, which covered bubbly, cap-bubbly, and churn-turbulent flows. These flow conditions also overlapped with twelve conditions of a previous study for comparison. The local flow parameters, such as void fractions, interfacial area concentrations (IAC), and bubble interface velocities, were measured at nine radial positions for the three axial locations and converted into area-averaged parameters. The axial evolutions of local flow structure were interpreted in terms of bubble coalescence, breakup, expansion of the gas-phase due to pressure drop and system pressure. An assessment of interfacial area transport equation (IATE) was made and compared with the experimental data. A discussion of the comparison between model prediction and the experimental results were made

Investigation of one-dimensional interfacial area transport for vertical upward air–water two-phase flow in an annular channel at elevated pressures

Energy Technology Data Exchange (ETDEWEB)

Ozar, B., E-mail: ozar@fauske.com [School of Nuclear Engineering, Purdue University, 400 Central Drive, West Lafayette, IN 47907-2017 (United States); Brooks, C.S. [School of Nuclear Engineering, Purdue University, 400 Central Drive, West Lafayette, IN 47907-2017 (United States); Euh, D.J. [Korea Atomic Energy Research Institute, 150 Deokjin, Yuseong, Daejeon 305-353 (Korea, Republic of); Hibiki, T.; Ishii, M. [School of Nuclear Engineering, Purdue University, 400 Central Drive, West Lafayette, IN 47907-2017 (United States)

2013-10-15

Highlights: • Interfacial area transport equation (IATE) for a rectangular duct is modified for an annulus. • IATE predicts interfacial area transport in bubbly-to-churn flow. • Scalability of IATE to elevated pressure conditions is validated. • Detailed 1D interfacial area transport data are presented. • Detailed interfacial area transport mechanisms are discussed. -- Abstract: The interfacial area transport of vertical, upward, air–water two-phase flows in an annular channel has been investigated at different system pressures. The inner and outer diameters of the annular channel were 19.1 mm and 38.1 mm, respectively. Twenty three inlet flow conditions were selected, which covered bubbly, cap-bubbly, and churn-turbulent flows. These flow conditions also overlapped with twelve conditions of a previous study for comparison. The local flow parameters, such as void fractions, interfacial area concentrations (IAC), and bubble interface velocities, were measured at nine radial positions for the three axial locations and converted into area-averaged parameters. The axial evolutions of local flow structure were interpreted in terms of bubble coalescence, breakup, expansion of the gas-phase due to pressure drop and system pressure. An assessment of interfacial area transport equation (IATE) was made and compared with the experimental data. A discussion of the comparison between model prediction and the experimental results were made.

Hermite-Pade approximation approach to hydromagnetic flows in convergent-divergent channels

International Nuclear Information System (INIS)

Makinde, O.D.

2005-10-01

The problem of two-dimensional, steady, nonlinear flow of an incompressible conducting viscous fluid in convergent-divergent channels under the influence of an externally applied homogeneous magnetic field is studied using a special type of Hermite-Pade approximation approach. This semi-numerical scheme offers some advantages over solutions obtained by using traditional methods such as finite differences, spectral method, shooting method, etc. It reveals the analytical structure of the solution function and the important properties of overall flow structure including velocity field, flow reversal control and bifurcations are discussed. (author)

Two-dimensional electron flow in pulsed power transmission lines and plasma opening switches

International Nuclear Information System (INIS)

Church, B.W.; Longcope, D.W.; Ng, C.K.; Sudan, R.N.

1991-01-01

The operation of magnetically insulated transmission lines (MITL) and the interruption of current in a plasma opening switch (POS) are determined by the physics of the electrons emitted by the cathode surface. A mathematical model describes the self-consistent two-dimensional flow of an electron fluid. A finite element code, FERUS, has been developed to solve the two equations which describe Poisson's and Ampere's law in two dimensions. The solutions from this code are obtained for parameters where the electron orbits are considerably modified by the self-magnetic field of the current. Next, the self-insulated electron flow in a MITL with a step change in cross-section is studied using a conventional two-dimensional fully electromagnetic particle-in-cell code, MASK. The equations governing two-dimensional quasi-static electron flow are solved numerically by a third technique which is suitable for predicting current interruption in a POS. The object of the study is to determine the critical load impedance, Z CL , required for current interruption for a given applied voltage, cathode voltage and plasma length. (author). 9 refs, 5 figs

Characterization of coherent structures in three-dimensional turbulent flows using the finite-size Lyapunov exponent

International Nuclear Information System (INIS)

Bettencourt, João H; López, Cristóbal; Hernández-García, Emilio

2013-01-01

In this paper, we use the finite-size Lyapunov exponent (FSLE) to characterize Lagrangian coherent structures in three-dimensional (3D) turbulent flows. Lagrangian coherent structures act as the organizers of transport in fluid flows and are crucial to understand their stirring and mixing properties. Generalized maxima (ridges) of the FSLE fields are used to locate these coherent structures. 3D FSLE fields are calculated in two phenomenologically distinct turbulent flows: a wall-bounded flow (channel flow) and a regional oceanic flow obtained by the numerical solution of the primitive equations where two-dimensional (2D) turbulence dominates. In the channel flow, autocorrelations of the FSLE field show that the structure is substantially different from the near wall to the mid-channel region and relates well to the more widely studied Eulerian coherent structure of the turbulent channel flow. The ridges of the FSLE field have complex shapes due to the 3D character of the turbulent fluctuations. In the oceanic flow, strong horizontal stirring is present and the flow regime is similar to that of 2D turbulence where the domain is populated by coherent eddies that interact strongly. This in turn results in the presence of high FSLE lines throughout the domain leading to strong non-local mixing. The ridges of the FSLE field are quasi-vertical surfaces, indicating that the horizontal dynamics dominates the flow. Indeed, due to rotation and stratification, vertical motions in the ocean are much less intense than horizontal ones. This suppression is absent in the channel flow, as the 3D character of the FSLE ridges shows. This article is part of a special issue of Journal of Physics A: Mathematical and Theoretical devoted to ‘Lyapunov analysis: from dynamical systems theory to applications’. (paper)

Effects of Parallel Channel Interactions on Two-Phase Flow Split in ...

African Journals Online (AJOL)

The tests would aid the development of a realistic transient computer model for tracking the distribution of two-phase flows into the multiple parallel channels of a Nuclear Reactor, during Loss of Coolant Accidents (LOCA), and were performed at the General Electric Nuclear Energy Division Laboratory, California. The test ...

Three-dimensional microbubble streaming flows

Science.gov (United States)

Rallabandi, Bhargav; Marin, Alvaro; Rossi, Massimiliano; Kaehler, Christian; Hilgenfeldt, Sascha

2014-11-01

Streaming due to acoustically excited bubbles has been used successfully for applications such as size-sorting, trapping and focusing of particles, as well as fluid mixing. Many of these applications involve the precise control of particle trajectories, typically achieved using cylindrical bubbles, which establish planar flows. Using astigmatic particle tracking velocimetry (APTV), we show that, while this two-dimensional picture is a useful description of the flow over short times, a systematic three-dimensional flow structure is evident over long time scales. We demonstrate that this long-time three-dimensional fluid motion can be understood through asymptotic theory, superimposing secondary axial flows (induced by boundary conditions at the device walls) onto the two-dimensional description. This leads to a general framework that describes three-dimensional flows in confined microstreaming systems, guiding the design of applications that profit from minimizing or maximizing these effects.

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

International Nuclear Information System (INIS)

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

1976-11-01

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

Investigating flow patterns in a channel with complex obstacles using the lattice Boltzmann method

Energy Technology Data Exchange (ETDEWEB)

Yojina, Jiraporn; Ngamsaad, Waipot; Nuttavut, Narin; Triampo, Darapond; Lenbury, Yongwimon; Sriyab, Somchai; Triampo, Wannapong [Faculty of Science, Mahidol University, Bangkok (Thailand); Kanthang, Paisan [Rajamangala University of Technology, Bangkok (Thailand)

2010-10-15

In this work, mesoscopic modeling via a computational lattice Boltzmann method (LBM) is used to investigate the flow pattern phenomena and the physical properties of the flow field around one and two square obstacles inside a two-dimensional channel with a fixed blockage ratio,{beta} =14 , centered inside a 2D channel, for a range of Reynolds numbers (Re) from 1 to 300. The simulation results show that flow patterns can initially exhibit laminar flow at low Re and then make a transition to periodic, unsteady, and, finally, turbulent flow as the Re get higher. Streamlines and velocity profiles and a vortex shedding pattern are observed. The Strouhal numbers are calculated to characterize the shedding frequency and flow dynamics. The effect of the layouts or configurations of the obstacles are also investigated, and the possible connection between the mixing process and the appropriate design of a chemical mixing system is discussed

On the occurrence of burnout downstream of a flow obstacle in boiling two-phase upward flow within a vertical annular channel

International Nuclear Information System (INIS)

Mori, Shoji; Tominaga, Akira; Fukano, Tohru

2004-01-01

If a flow obstruction such as a spacer is set in a boiling two-phase flow within an annular channel, the inner tube of which is used as a heater, the temperature on the surface of the heater tube is severely affected by the existence of the spacer. In some case the spacer has a cooling effect, and in the other case it causes the dryout of the cooling liquid film on the heating surface resulting in the burnout of the tube. But the burnout mechanism near the spacer is not still clear. In the present paper we discus the influence of the flow obstacle on the occurrence of burnout downstream of the flow obstacle in boiling two-phase upward flow within a vertical annular channel. (author)

The effect of the flow direction inside the header on two-phase flow distribution in parallel vertical channels

International Nuclear Information System (INIS)

Marchitto, A.; Fossa, M.; Guglielmini, G.

2012-01-01

Uniform fluid distribution is essential for efficient operation of chemical-processing equipment such as contactors, reactors, mixers, burners and in most refrigeration equipment, where two phases are acting together. To obtain optimum distribution, proper consideration must be given to flow behaviour in the distributor, flow conditions upstream and downstream of the distributor, and the distribution requirements (fluid or phase) of the equipment. Even though the principles of single phase distribution have been well developed for more than three decades, they are frequently not taken in the right account by equipment designers when a mixture is present, and a significant fraction of process equipment consequently suffers from maldistribution. The experimental investigation presented in this paper is aimed at understanding the main mechanisms which drive the flow distribution inside a two-phase horizontal header in order to design improved distributors and to optimise the flow distribution inside compact heat exchanger. Experimentation was devoted to establish the influence of the inlet conditions and of the channel/distributor geometry on the phase/mass distribution into parallel vertical channels. The study is carried out with air–water mixtures and it is based on the measurement of component flow rates in individual channels and on pressure drops across the distributor. The effects of the operating conditions, the header geometry and the inlet port nozzle were investigated in the ranges of liquid and gas superficial velocities of 0.2–1.2 and 1.5–16.5 m/s, respectively. In order to control the main flow direction inside the header, different fitting devices were tested; the insertion of a co-axial, multi-hole distributor inside the header has confirmed the possibility of greatly improving the liquid and gas flow distribution by the proper selection of position, diameter and number of the flow openings between the supplying distributor and the system of

Three-dimensional multi-relaxation-time lattice Boltzmann front-tracking method for two-phase flow

International Nuclear Information System (INIS)

Xie Hai-Qiong; Zeng Zhong; Zhang Liang-Qi

2016-01-01

We developed a three-dimensional multi-relaxation-time lattice Boltzmann method for incompressible and immiscible two-phase flow by coupling with a front-tracking technique. The flow field was simulated by using an Eulerian grid, an adaptive unstructured triangular Lagrangian grid was applied to track explicitly the motion of the two-fluid interface, and an indicator function was introduced to update accurately the fluid properties. The surface tension was computed directly on a triangular Lagrangian grid, and then the surface tension was distributed to the background Eulerian grid. Three benchmarks of two-phase flow, including the Laplace law for a stationary drop, the oscillation of a three-dimensional ellipsoidal drop, and the drop deformation in a shear flow, were simulated to validate the present model. (paper)

Flow and sediment transport across oblique channels

DEFF Research Database (Denmark)

Hjelmager Jensen, Jacob; Madsen, Erik Østergaard; Fredsøe, Jørgen

1998-01-01

A 3D numerical investigation of flow across channels aligned obliquely to the main flow direction has been conducted. The applied numerical model solves the Reynolds-averaged Navier-Stokes equations using the k-ε model for turbulence closure on a curvilinear grid. Three momentum equations...... are solved, but the computational domain is 2D due to a uniformity along the channel alignment. Two important flow features arise when the flow crosses the channel: (i) the flow will be refracted in the direction of the channel alignment. This may be described by a depth-averaged model. (ii) due to shear...

A thermal analysis computer programme package for the estimation of KANUPP coolant channel flows and outlet header temperature distribution

International Nuclear Information System (INIS)

Siddiqui, M.S.

1992-06-01

COFTAN is a computer code for actual estimation of flows and temperatures in the coolant channels of a pressure tube heavy water reactor. The code is being used for Candu type reactor with coolant flowing 208 channels. The simulation model first performs the detailed calculation of flux and power distribution based on two groups diffusion theory treatment on a three dimensional mesh and then channel powers, resulting from the summation of eleven bundle powers in each of the 208 channels, are employed to make actual estimation of coolant flows using channel powers and channel outlet temperature monitored by digital computers. The code by using the design flows in individual channels and applying a correction factor based on control room monitored flows in eight selected channels, can also provide a reserve computational tool of estimating individual channel outlet temperatures, thus providing an alternate arrangements for checking Rads performance. 42 figs. (Orig./A.B.)

Contribution of large-scale coherent structures towards the cross flow in two interconnected channels

International Nuclear Information System (INIS)

Mahmood, A.; Rohde, M.; Hagen, T.H.J.J. van der; Mudde, R.F.

2009-01-01

Single phase cross flow through a gap region joining two vertical channels has been investigated experimentally for Reynolds numbers, based on the channels hydraulic diameter, ranging from 850 to 21000. The flow field in the gap region is investigated by 2D-PIV and the inter channel mass transfer is quantified by the tracer injection method. Experiments carried out for variable gap heights and shape show the existence of a street of large-scale counter rotating vortices on either side of the channel-gap interface, resulting from the mean velocity gradient in the gap and the main channel region. The appearance of the coherent vortices is subject to a threshold associated with the difference between the maximum and the minimum average stream wise velocities in the channel and the gap region, respectively. The auto power spectral density of the cross velocity component in the gap region exhibits a slope of -3 in the inertial range, indicating the 2D nature of these vortices. The presence of the large-scale vortices enhances the mass transfer through the gap region by approximately 63% of the mass transferred by turbulent mixing alone. The inter-channel mass transfer, due to cross flow, is found to be dependent not only on the large-scale vortices characteristics, but also on the gap geometry. (author)

On the occurrence of burnout downstream of the flow obstacle in boiling two-phase upward flow within a vertical annular channel

International Nuclear Information System (INIS)

Mori, Shoji; Fukano, Tohru

2003-01-01

If a flow obstruction such as a spacer is set in a boiling two-phase flow within an annular channel, the inner tube of which is used as a heater, the temperature on the surface of the heater tube is severely affected by the existence of the spacer. In some cases the spacer has a cooling effect, and in the other case it causes the dryout of the cooling liquid film on the heating surface resulting in the burnout of the tube. But the thermo-fluid dynamic mechanism to cause burnout near the spacer is not still clear. In the present paper we discuss the influence of the flow obstacle on the occurrence of burnout downstream of the flow obstacle in boiling two-phase upward flow within a vertical annular channel. (author)

Three-dimensional (3-D) video systems: bi-channel or single-channel optics?

Science.gov (United States)

van Bergen, P; Kunert, W; Buess, G F

1999-11-01

This paper presents the results of a comparison between two different three-dimensional (3-D) video systems, one with single-channel optics, the other with bi-channel optics. The latter integrates two lens systems, each transferring one half of the stereoscopic image; the former uses only one lens system, similar to a two-dimensional (2-D) endoscope, which transfers the complete stereoscopic picture. In our training centre for minimally invasive surgery, surgeons were involved in basic and advanced laparoscopic courses using both a 2-D system and the two 3-D video systems. They completed analog scale questionnaires in order to record a subjective impression of the relative convenience of operating in 2-D and 3-D vision, and to identify perceived deficiencies in the 3-D system. As an objective test, different experimental tasks were developed, in order to measure performance times and to count pre-defined errors made while using the two 3-D video systems and the 2-D system. Using the bi-channel optical system, the surgeon has a heightened spatial perception, and can work faster and more safely than with a single-channel system. However, single-channel optics allow the use of an angulated endoscope, and the free rotation of the optics relative to the camera, which is necessary for some operative applications.

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

Experimental Investigation on Zonal Structure in Drag-Reducing Channel Flow with Surfactant Additives

Directory of Open Access Journals (Sweden)

Masaaki Motozawa

2011-01-01

Full Text Available The spatial structure of a drag-reducing channel flow with surfactant additives in a two-dimensional channel was investigated experimentally. We carried out detailed measurements of the instantaneous velocity in the streamwise wall-normal plane and streamwise spanwise plane by using particle image velocimetry (PIV. The surfactant used in this experiment is a kind of cationic surfactant CTAC. The weight concentrations of the CTAC solution were 25 and 40 ppm on the flow. We considered the effects of Reynolds number ranging from 10000 to 25000 and the weight concentration of CTAC. The results of this paper showed that in the drag-reducing flow, there appeared an area where the root mean square of streamwise velocity fluctuation and the vorticity fluctuation sharply decreased. This indicated that two layers with different turbulent structure coexisted on the boundary of this area. Moreover, these layers had characteristic flow structures, as confirmed by observation of the instantaneous vorticity fluctuation map.

Complete energetic description of hydrokinetic turbine impact on flow channel dynamics

Science.gov (United States)

Brasseale, E.; Kawase, M.

2016-02-01

Energy budget analysis on tidal channels quantifies and demarcates the impacts of marine renewables on environmental fluid dynamics. Energy budget analysis assumes the change in total kinetic energy within a volume of fluid can be described by the work done by each force acting on the flow. In a numerically simulated channel, the balance between energy change and work done has been validated up to 5% error.The forces doing work on the flow include pressure, turbulent dissipation, and stress from the estuary floor. If hydrokinetic turbines are installed in an estuarine channel to convert tidal energy into usable power, the dynamics of the channel change. Turbines provide additional pressure work against the flow of the channel which will slow the current and lessen turbulent dissipation and bottom stress. These losses may negatively impact estuarine circulation, seafloor scour, and stratification.The environmental effects of turbine deployment have been quantified using a three dimensional, Reynolds-averaged, Navier-Stokes model of an idealized flow channel situated between the ocean and a large estuarine basin. The channel is five kilometers wide, twenty kilometers long and fifty meters deep, and resolved to a grid size of 10 meters by 10 meters by 1 meter. Tidal currents are simulated by an initial difference in sea surface height across the channel of 160 centimeters from the channel entrance to the channel exit. This creates a pressure gradient which drives flow through the channel. Tidal power turbines are represented as disks that force the channel in proportion to the strength of the current. Three tidal turbines twenty meters in diameters have been included in the model to simulate the impacts of a pilot scale test deployment.This study is the first to appreciate the energetic impact of marine renewables in a three dimensional model through the energy equation's constituent terms. This study provides groundwork for understanding and predicting the

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.

Rarefied gas flows through a curved channel: Application of a diffusion-type equation

Science.gov (United States)

Aoki, Kazuo; Takata, Shigeru; Tatsumi, Eri; Yoshida, Hiroaki

2010-11-01

Rarefied gas flows through a curved two-dimensional channel, caused by a pressure or a temperature gradient, are investigated numerically by using a macroscopic equation of convection-diffusion type. The equation, which was derived systematically from the Bhatnagar-Gross-Krook model of the Boltzmann equation and diffuse-reflection boundary condition in a previous paper [K. Aoki et al., "A diffusion model for rarefied flows in curved channels," Multiscale Model. Simul. 6, 1281 (2008)], is valid irrespective of the degree of gas rarefaction when the channel width is much shorter than the scale of variations of physical quantities and curvature along the channel. Attention is also paid to a variant of the Knudsen compressor that can produce a pressure raise by the effect of the change of channel curvature and periodic temperature distributions without any help of moving parts. In the process of analysis, the macroscopic equation is (partially) extended to the case of the ellipsoidal-statistical model of the Boltzmann equation.

Steady state flow analysis of two-phase natural circulation in multiple parallel channel loop

International Nuclear Information System (INIS)

Bhusare, V.H.; Bagul, R.K.; Joshi, J.B.; Nayak, A.K.; Kannan, Umasankari; Pilkhwal, D.S.; Vijayan, P.K.

2016-01-01

Highlights: • Liquid circulation velocity increases with increasing superficial gas velocity. • Total two-phase pressure drop decreases with increasing superficial gas velocity. • Channels with larger driving force have maximum circulation velocities. • Good agreement between experimental and model predictions. - Abstract: In this work, steady state flow analysis has been carried out experimentally in order to estimate the liquid circulation velocities and two-phase pressure drop in air–water multichannel circulating loop. Experiments were performed in 15 channel circulating loop. Single phase and two-phase pressure drops in the channels have been measured experimentally and have been compared with theoretical model of Joshi et al. (1990). Experimental measurements show good agreement with model.

Numerical simulations of heat transfer in an annular fuel channel with three-dimensional spacer ribs set up periodically under a fully developed turbulent flow

International Nuclear Information System (INIS)

Takase, Kazuyuki; Akino, Norio

1996-06-01

Thermal-hydraulic characteristics of an annular fuel channel with spacer ribs for high temperature gas-cooled reactors were analyzed numerically by three-dimensional heat transfer computations under a fully developed turbulent flow. The two-equations κ-ε turbulence model was applied to the present turbulent analysis. In particular, the κ-ε turbulence model constants and the turbulent Prandtl number were improved from the previous standard values proposed by Jones and Launder in order to obtain heat transfer predictions with higher accuracy. Consequently, heat transfer coefficients and friction factors in the spacer-ribbed fuel channel were predicted with sufficient accuracy in the range of Reynolds number exceeding 3000. It was clarified quantitatively from the present study that main mechanism for the heat transfer augmentation in the spacer-ribbed fuel channel was combined effects of the turbulence promoter effect by the spacer ribs and the velocity acceleration effect by a reduction in the channel cross-section. (author)

Structure of two-phase adiabatic flow in air sparging regime in vertical cylindrical channel with water

Directory of Open Access Journals (Sweden)

V. I. Solonin

2014-01-01

Full Text Available The article presents a research of two-phase adiabatic flow in air sparging regime in vertical cylindrical channel filled with water. A purpose of the work is to obtain experimental data for further analysis of a character of the moving phases. Research activities used the optic methods PIV (Particle Image Visualization because of their noninvasiveness to obtain data without disturbing effect on the flow. A laser sheet illuminated the fluorescence particles, which were admixed in water along the channel length. A digital camera recorded their motion for a certain time interval that allowed building the velocity vector fields. As a result, gas phase velocity components typical for a steady area of the channel and their relations for various intensity of volume air rate were obtained. A character of motion both for an air bubble and for its surrounding liquid has been conducted. The most probable direction of phases moving in the channel under sparging regime is obtained by building the statistic scalar fields. The use of image processing enabled an analysis of the initial area of the air inlet into liquid. A characteristic curve of the bubbles offset from the axis for various intensity of volume gas rate and channel diameter is defined. A character of moving phases is obtained by building the statistic scalar fields. The values of vertical components of liquid velocity in the inlet part of channel are calculated. Using the obtained data of the gas phase velocities a true void fraction was calculated. It was compared with the values of void fraction, calculated according to the liquid level change in the channel. Obtained velocities were compared with those of the other researchers, and a small difference in their values was explained by experimental conditions. The article is one of the works to research the two-phase flows with no disturbing effect on them. Obtained data allow us to understand a character of moving the two-phase flows in

Computation of two-dimensional isothermal flow in shell-and-tube heat exchangers

International Nuclear Information System (INIS)

Carlucci, L.N.; Galpin, P.F.; Brown, J.D.; Frisina, V.

1983-07-01

A computational procedure is outlined whereby two-dimensional isothermal shell-side flow distributions can be calculated for tube bundles having arbitrary boundaries and flow blocking devices, such as sealing strips, defined in arbitrary locations. The procedure is described in some detail and several computed results are presented to illustrate the robustness and generality of the method

Mixed finite element simulations in two-dimensional groundwater flow problems

International Nuclear Information System (INIS)

Kimura, Hideo

1989-01-01

A computer code of groundwater flow in two-dimensional porous media based on the mixed finite element method was developed for accurate approximations of Darcy velocities in safety evaluation of radioactive waste disposal. The mixed finite element procedure solves for both the Darcy velocities and pressure heads simultaneously in the Darcy equation and continuity equation. Numerical results of a single well pumping at a constant rate in a uniform flow field showed that the mixed finite element method gives more accurate Darcy velocities nearly 50 % on average error than standard finite element method. (author)

Numerical study on channel size effect for proton exchange membrane fuel cell with serpentine flow field

International Nuclear Information System (INIS)

Wang Xiaodong; Yan Weimon; Duan Yuanyuan; Weng Fangbor; Jung Guobin; Lee Chiyuan

2010-01-01

This work numerically investigates the effect of the channel size on the cell performance of proton exchange membrane (PEM) fuel cells with serpentine flow fields using a three-dimensional, two-phase model. The local current densities in the PEM, oxygen mass flow rates and liquid water concentrations at the interface of the cathode gas diffusion layer and catalyst layer were analyzed to understand the channel size effect. The predictions show that smaller channel sizes enhance liquid water removal and increase oxygen transport to the porous layers, which improve cell performance. Additionally, smaller channel sizes also provide more uniform current density distributions in the cell. However, as the channel size decreases, the total pressure drops across the cell increases, which leads to more pump work. With taking into account the pressure losses, the optimal cell performance occurs for a cell with a flow channel cross-sectional area of 0.535 x 0.535 mm 2 .

Separation prediction in two dimensional boundary layer flows using artificial neural networks

International Nuclear Information System (INIS)

Sabetghadam, F.; Ghomi, H.A.

2003-01-01

In this article, the ability of artificial neural networks in prediction of separation in steady two dimensional boundary layer flows is studied. Data for network training is extracted from numerical solution of an ODE obtained from Von Karman integral equation with approximate one parameter Pohlhousen velocity profile. As an appropriate neural network, a two layer radial basis generalized regression artificial neural network is used. The results shows good agreements between the overall behavior of the flow fields predicted by the artificial neural network and the actual flow fields for some cases. The method easily can be extended to unsteady separation and turbulent as well as compressible boundary layer flows. (author)

Slip-line field analysis of metal flow during two dimensional forging

International Nuclear Information System (INIS)

Fenton, R.G.; Khataan, H.A.

1981-01-01

A method of computation and a computer software package were developed for solving problems of two dimensional plastic flow between symmetrical dies of any specified shape. The load required to initiate plastic flow, the stress and velocity distributions in the plastic region of the metal, and the pressure distribution acting on the die are determined. The method can be used to solve any symmetrical plane strain flow problem regardless of the complexity of the die. The accurate solution obtained by this efficient method can provide valuable help to forging die designers. (Author) [pt

A New Void Fraction Measurement Method for Gas-Liquid Two-Phase Flow in Small Channels

Directory of Open Access Journals (Sweden)

Huajun Li

2016-01-01

Full Text Available Based on a laser diode, a 12 × 6 photodiode array sensor, and machine learning techniques, a new void fraction measurement method for gas-liquid two-phase flow in small channels is proposed. To overcome the influence of flow pattern on the void fraction measurement, the flow pattern of the two-phase flow is firstly identified by Fisher Discriminant Analysis (FDA. Then, according to the identification result, a relevant void fraction measurement model which is developed by Support Vector Machine (SVM is selected to implement the void fraction measurement. A void fraction measurement system for the two-phase flow is developed and experiments are carried out in four different small channels. Four typical flow patterns (including bubble flow, slug flow, stratified flow and annular flow are investigated. The experimental results show that the development of the measurement system is successful. The proposed void fraction measurement method is effective and the void fraction measurement accuracy is satisfactory. Compared with the conventional laser measurement systems using standard laser sources, the developed measurement system has the advantages of low cost and simple structure. Compared with the conventional void fraction measurement methods, the proposed method overcomes the influence of flow pattern on the void fraction measurement. This work also provides a good example of using low-cost laser diode as a competent replacement of the expensive standard laser source and hence implementing the parameter measurement of gas-liquid two-phase flow. The research results can be a useful reference for other researchers’ works.

A New Void Fraction Measurement Method for Gas-Liquid Two-Phase Flow in Small Channels.

Science.gov (United States)

Li, Huajun; Ji, Haifeng; Huang, Zhiyao; Wang, Baoliang; Li, Haiqing; Wu, Guohua

2016-01-27

Based on a laser diode, a 12 × 6 photodiode array sensor, and machine learning techniques, a new void fraction measurement method for gas-liquid two-phase flow in small channels is proposed. To overcome the influence of flow pattern on the void fraction measurement, the flow pattern of the two-phase flow is firstly identified by Fisher Discriminant Analysis (FDA). Then, according to the identification result, a relevant void fraction measurement model which is developed by Support Vector Machine (SVM) is selected to implement the void fraction measurement. A void fraction measurement system for the two-phase flow is developed and experiments are carried out in four different small channels. Four typical flow patterns (including bubble flow, slug flow, stratified flow and annular flow) are investigated. The experimental results show that the development of the measurement system is successful. The proposed void fraction measurement method is effective and the void fraction measurement accuracy is satisfactory. Compared with the conventional laser measurement systems using standard laser sources, the developed measurement system has the advantages of low cost and simple structure. Compared with the conventional void fraction measurement methods, the proposed method overcomes the influence of flow pattern on the void fraction measurement. This work also provides a good example of using low-cost laser diode as a competent replacement of the expensive standard laser source and hence implementing the parameter measurement of gas-liquid two-phase flow. The research results can be a useful reference for other researchers' works.

An experimental study of gravity-driven countercurrent two-phase flow in horizontal and inclined channels

International Nuclear Information System (INIS)

Lillibridge, K.H.; Ghiaasiaan, S.M.; Abdel-Khalik, S.I.

1994-01-01

Countercurrent two-phase flow in horizontal and inclined channels, connecting a sealed liquid-filled reservoir to the atmosphere, is experimentally studied. This type of gravity-driven countercurrent two-phase flow can occur during the operation of passive safety coolant injection systems of advanced reactors. It can also occur in the pressurizer surge line of pressurized water reactors during severe accidents when the hot leg becomes voided. Four distinct flow regimes are identified: (a) stratified countercurrent, which mainly occurs when the channel is horizontal; (b) intermittent stratified-slug; (c) oscillating, which occurs when the angle of inclination is ≥30 deg; and (d) annular countercurrent. The characteristics of each regime and their sensitivity to important geometric parameters are examined. The superficial velocities in the stratified countercurrent and oscillating regimes are empirically correlated

CANDU channel flow verification

International Nuclear Information System (INIS)

Mazalu, N.; Negut, Gh.

1997-01-01

The purpose of this evaluation was to obtain accurate information on each channel flow that enables us to assess precisely the level of reactor thermal power and, for reasons of safety, to establish which channel is boiling. In order to assess the channel flow parameters, computer simulations were done with the NUCIRC code and the results were checked by measurements. The complete channel flow measurements were made in the zero power cold condition. In hot conditions there were made flow measurements using the Shut Down System 1 (SDS 1) flow devices from 0.1 % F.P. up to 100 % F.P. The NUCIRC prediction for CANDU channel flows and the measurements by Ultrasonic Flow Meter at zero power cold conditions and SDS 1 flow channel measurements at different reactor power levels showed an acceptable agreement. The 100 % F.P. average errors for channel flow of R, shows that suitable NUCIRC flow assessment can be made. So, it can be done a fair prediction of the reactor power distribution. NUCIRC can predict accurately the onset of boiling and helps to warn at the possible power instabilities at high powers or it can detect the flow blockages. The thermal hydraulic analyst has in NUCIRC a suitable tool to do accurate predictions for the thermal hydraulic parameters for different steady state power levels which subsequently leads to an optimal CANDU reactor operation. (authors)

Three-dimensional numerical modeling of turbulent single-phase and two-phase flow in curved pipes

International Nuclear Information System (INIS)

Xin, R.C.; Dong, Z.F.; Ebadian, M.A.

1996-01-01

In this study, three-dimensional single-phase and two-phase flows in curved pipes have been investigated numerically. Two different pipe configurations were computed. When the results of the single-phase flow simulation were compared with the experimental data, a fairly good agreement was achieved. A flow-developing process has been suggested in single-phase flow, in which the turbulence is stronger near the outer tube wall than near the inner tube wall. For two-phase flow, the Eulerian multiphase model was used to simulate the phase distribution of a three-dimensional gas-liquid bubble flow in curved pipe. The RNG/κ-ε turbulence model was used to determine the turbulence field. An inlet gas void fraction of 5 percent was simulated. The gas phase effects on the liquid phase flow velocity have been examined by comparing the results of single-phase flow and two-phase flow. The findings show that for the downward flow in the U bend, the gas concentrates at the inner portion of the cross section at φ = π/18 - π/6 in most cases. The results of the phase distribution simulation are compared to experimental observations qualitatively and topologically

GIS-based two-dimensional numerical simulation of rainfall-induced debris flow

Directory of Open Access Journals (Sweden)

C. Wang

2008-02-01

Full Text Available This paper aims to present a useful numerical method to simulate the propagation and deposition of debris flow across the three dimensional complex terrain. A depth-averaged two-dimensional numerical model is developed, in which the debris and water mixture is assumed to be continuous, incompressible, unsteady flow. The model is based on the continuity equations and Navier-Stokes equations. Raster grid networks of digital elevation model in GIS provide a uniform grid system to describe complex topography. As the raster grid can be used as the finite difference mesh, the continuity and momentum equations are solved numerically using the finite difference method. The numerical model is applied to simulate the rainfall-induced debris flow occurred in 20 July 2003, in Minamata City of southern Kyushu, Japan. The simulation reproduces the propagation and deposition and the results are in good agreement with the field investigation. The synthesis of numerical method and GIS makes possible the solution of debris flow over a realistic terrain, and can be used to estimate the flow range, and to define potentially hazardous areas for homes and road section.

To the elementary theory of critical (maximum) flow rate of two-phase mixture in channels with various sections

International Nuclear Information System (INIS)

Nigmatulin, B.I.; Soplenkov, K.I.

1978-01-01

On the basis of the concepts of two-phase dispersive flow with various structures (bubble, vapour-drop etc) in the framework of the two-speed and two-temperature one-dimension stationary model of the current with provision for phase transitions the conditions, under which a critical (maximum) flow rate of two-phase mixture is achieved during its outflowing from a channel with the pre-set geometry, have been determined. It is shown, that for the choosen set of two-phase flow equations with the known parameters of deceleration and structure one of the critical conditions is satisfied: either solution of the set of equations corresponding a critical flow rate is a special one, i.e. passes through a special point locating between minimum and outlet channel sections where the carrying phase velocity approaches the value of decelerated sound speed in the mixture or the determinator of the initial set of equations equals zero for the outlet channel sections, i.e. gradients of the main flow parameters tend to +-infinity in this section, and carrying phase velocity also approaches the value of the decelerated sound velocity in the mixture

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

International Nuclear Information System (INIS)

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

2004-01-01

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

Two-dimensional numerical simulation of flow around three-stranded rope

Science.gov (United States)

Wang, Xinxin; Wan, Rong; Huang, Liuyi; Zhao, Fenfang; Sun, Peng

2016-08-01

Three-stranded rope is widely used in fishing gear and mooring system. Results of numerical simulation are presented for flow around a three-stranded rope in uniform flow. The simulation was carried out to study the hydrodynamic characteristics of pressure and velocity fields of steady incompressible laminar and turbulent wakes behind a three-stranded rope. A three-cylinder configuration and single circular cylinder configuration are used to model the three-stranded rope in the two-dimensional simulation. The governing equations, Navier-Stokes equations, are solved by using two-dimensional finite volume method. The turbulence flow is simulated using Standard κ-ɛ model and Shear-Stress Transport κ-ω (SST) model. The drag of the three-cylinder model and single cylinder model is calculated for different Reynolds numbers by using control volume analysis method. The pressure coefficient is also calculated for the turbulent model and laminar model based on the control surface method. From the comparison of the drag coefficient and the pressure of the single cylinder and three-cylinder models, it is found that the drag coefficients of the three-cylinder model are generally 1.3-1.5 times those of the single circular cylinder for different Reynolds numbers. Comparing the numerical results with water tank test data, the results of the three-cylinder model are closer to the experiment results than the single cylinder model results.

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.

Two-phase fluid flow measurements in small diameter channels using real-time neutron radiography

International Nuclear Information System (INIS)

Carlisle, B.S.; Johns, R.C.; Hassan, Y.A.

2004-01-01

A series of real-time, neutron radiography, experiments are ongoing at the Texas A and M Nuclear Science Center Reactor (NSCR). These tests determine the resolving capabilities for radiographic imaging of two phase water and air flow regimes through small diameter flow channels. Though both film and video radiographic imaging is available, the real-time video imaging was selected to capture the dynamic flow patterns with results that continue to improve. (author)

On the large-scale structure and spectral dynamics of two-dimensional turbulence in a periodic channel

NARCIS (Netherlands)

Kramer, W.; Clercx, H.J.H.; van Heijst, G.J.F.

2008-01-01

This paper reports on a numerical study of forced two-dimensional turbulence in a periodic channel with flat no-slip walls. Since corners or curved domain boundaries, which are met in the standard rectangular, square, or circular geometries, are absent in this geometry, the (statistical) analysis of

HYTRAN: hydraulic transient code for investigating channel flow stability

International Nuclear Information System (INIS)

Kao, H.S.; Cardwell, W.R.; Morgan, C.D.

1976-01-01

HYTRAN is an analytical program used to investigate the possibility of hydraulic oscillations occurring in a reactor flow channel. The single channel studied is ordinarily the hot channel in the reactor core, which is parallel to other channels and is assumed to share a constant pressure drop with other channels. Since the channel of highest thermal state is studied, provision is made for two-phase flow that can cause a flow instability in the channel. HYTRAN uses the CHATA(1) program to establish a steady-state condition. A heat flux perturbation is then imposed on the channel, and the flow transient is calculated as a function of time

Studies on diversion cross-flow between two parallel channels communicating by a lateral slot. II

International Nuclear Information System (INIS)

Tapucu, A.; Merilo, M.

1977-01-01

The axial pressure variations of two parallel channels with single phase flows communicating by a long lateral slot have been studied experimentally. Using mass and momentum conservation principles, the axial pressure variations have been derived in terms of two parameters ksub(d) and ksub(r), for donor and recipient channels, respectively. These parameters include the combined effect of fluid transferred from donor to recipient channel, and drag force brought on by the connection gap, and are functions of the velocities and slot geometry parameters. A pressure difference oscillation between channels along the slot has been detected which is sinusoidal with wave lengths which seem to be a function of the gap clearance. (Auth.)

Two-phase flow instability and bifurcation analysis of inclined multiple uniformly heated channels - 15107

International Nuclear Information System (INIS)

Mishra, A.M.; Paul, S.; Singh, S.; Panday, V.

2015-01-01

In this paper the two-phase flow instability analysis of multiple heated channels with various inclinations is studied. In addition, the bifurcation analysis is also carried out to capture the nonlinear dynamics of the system and to identify the regions in parameter space for which subcritical and supercritical bifurcations exist. In order to carry out the analysis, the system is mathematically represented by nonlinear Partial Differential Equation (PDE) for mass, momentum and energy in single as well as two-phase region. Then converted into Ordinary Differential Equation (ODE) using weighted residual method. Also, coupling equation is being used under the assumption that pressure drop in each channel is the same and the total mass flow rate is equal to sum of the individual mass flow rates. The homogeneous equilibrium model is used for the analysis. Stability Map is obtained in terms of phase change number (Npch) and Subcooling Number (Nsb) by solving a set of nonlinear, coupled algebraic equations obtained at equilibrium using Newton Raphson Method. MATLAB Code is verified by comparing it with results obtained by Matcont (Open source software) under same parametric values. Numerical simulations of the time-dependent, nonlinear ODEs are carried out for selected points in the operating parameter space to obtain the actual damped and growing oscillations in the channel inlet flow velocity which confirms the stability region across the stability map. Generalized Hopf (GH) points are observed for different inclinations, they are also points for subcritical and supercritical bifurcations. (authors)

Stochastic and collisional diffusion in two-dimensional periodic flows

International Nuclear Information System (INIS)

Doxas, I.; Horton, W.; Berk, H.L.

1990-05-01

The global effective diffusion coefficient D* for a two-dimensional system of convective rolls with a time dependent perturbation added, is calculated. The perturbation produces a background diffusion coefficient D, which is calculated analytically using the Menlikov-Arnold integral. This intrinsic diffusion coefficient is then enhanced by the unperturbed flow, to produce the global effective diffusion coefficient D*, which we can calculate theoretically for a certain range of parameters. The theoretical value agrees well with numerical simulations. 23 refs., 4 figs

Effect on two-phase flow frictional pressure drop characteristic in narrow rectangular channel at fluctuant condition

International Nuclear Information System (INIS)

Li Changwei; Cao Xiaxin; Sun Licheng; Jin Guangyuan

2013-01-01

Based on the data of two-phase flow in narrow rectangular channel, the influence of the two-phase flow friction characteristic under the different fluctuant states was analyzed. Through analyzing the experimental data, it is shown that the fluctuant amplitude of the friction pressure drop is affected slightly by the fluctuant period in narrow rectangular channel, but the frequency of the friction pressure drop fluctuation is changed. However, the change of fluctuant period is of little effect on the average frictional pressure drop. Comparing the φ l 2 (φ g 2 )-X variation curves at static condition with the ones at fluctuant condition, using the L-M method, it's found that the two phase frictional pressure drop in the narrow rectangular channel under the fluctuant state can be calculated by the φ l 2 (φ g 2 )-X variation curve at static condition. (authors)

Local gas- and liquid-phase measurements for air-water two-phase flows in a rectangular channel

International Nuclear Information System (INIS)

Zhou, X.; Sun, X.; Williams, M.; Fu, Y.; Liu, Y.

2014-01-01

Local gas- and liquid-phase measurements of various gas-liquid two-phase flows, including bubbly, cap-bubbly, slug, and churn-turbulent flows, were performed in an acrylic vertical channel with a rectangular cross section of 30 mm x 10 mm and height of 3.0 m. All the measurements were carried out at three measurement elevations along the flow channel, with z/D h = 9, 72, and 136, respectively, to study the flow development. The gas-phase velocity, void fraction, and bubble number frequency were measured using a double-sensor conductivity probe. A high-speed imaging system was utilized to perform the flow regime visualization and to provide additional quantitative information of the two-phase flow structure. An image processing scheme was developed to obtain the gas-phase velocity, void fraction, Sauter mean diameter, bubble number density, and interfacial area concentration. The liquid-phase velocity and turbulence measurements were conducted using a particle image velocimetry-planar laser-induced fluorescence (PIV-PLIF) system, which enables whole-field and high-resolution data acquisition. An optical phase separation method, which uses fluorescent particles and optical filtration technique, is adopted to extract the velocity information of the liquid phase. An image pre-processing scheme is imposed on the raw PIV images acquired to remove noises due to the presence of bubble residuals and optically distorted particles in the images captured by the PIV-PLIF system. Due to the better light access and less bubble distortion in the narrow rectangular channel, the PIV-PLIF system were able to perform reasonably well in flows of even higher void fractions as compared to the situations with circular pipe test sections. The flow conditions being studied covered various flow regime transitions, void fractions, and liquid-phase flow Reynolds numbers. The obtained experimental data can also be used to validate two-phase CFD results. (author)

A general theory of two- and three-dimensional rotational flow in subsonic and transonic turbomachines

Science.gov (United States)

Wu, Chung-Hua

1993-01-01

This report represents a general theory applicable to axial, radial, and mixed flow turbomachines operating at subsonic and supersonic speeds with a finite number of blades of finite thickness. References reflect the evolution of computational methods used, from the inception of the theory in the 50's to the high-speed computer era of the 90's. Two kinds of relative stream surfaces, S(sub 1) and S(sub 2), are introduced for the purpose of obtaining a three-dimensional flow solution through the combination of two-dimensional flow solutions. Nonorthogonal curvilinear coordinates are used for the governing equations. Methods of computing transonic flow along S(sub 1) and S(sub 2) stream surfaces are given for special cases as well as for fully three-dimensional transonic flows. Procedures pertaining to the direct solutions and inverse solutions are presented. Information on shock wave locations and shapes needed for computations are discussed. Experimental data from a Deutsche Forschungs- und Versuchsanstalt fur Luft- und Raumfahrt e.V. (DFVLR) rotor and from a Chinese Academy of Sciences (CAS) transonic compressor rotor are compared with the computed flow properties.

Two-dimensional steady unsaturated flow through embedded elliptical layers

Science.gov (United States)

Bakker, Mark; Nieber, John L.

2004-12-01

New analytic element solutions are presented for unsaturated, two-dimensional steady flow in vertical planes that include nonoverlapping impermeable elliptical layers and elliptical inhomogeneities. The hydraulic conductivity, which is represented by an exponential function of the pressure head, differs between the inside and outside of an elliptical inhomogeneity; both the saturated hydraulic conductivity and water retention parameters are allowed to differ between the inside and outside. The Richards equation is transformed, through the Kirchhoff transformation and a second standard transformation, into the modified Helmholtz equation. Analytic element solutions are obtained through separation of variables in elliptical coordinates. The resulting equations for the Kirchhoff potential consist of infinite sums of products of exponentials and modified Mathieu functions. In practical applications the series are truncated but still fulfill the differential equation exactly; boundary conditions are met approximately but up to machine accuracy, provided that enough terms are used. The pressure head, saturation, and flow may be computed analytically at any point in the vadose zone. Examples are given of the shadowing effect of an impermeable elliptical layer in a uniform flow field and funnel-type flow between two elliptical inhomogeneities. The presented solutions may be applied to study transport processes in vadose zones containing many impermeable elliptical layers or elliptical inhomogeneities.

Two-phase flow models

International Nuclear Information System (INIS)

Delaje, Dzh.

1984-01-01

General hypothesis used to simplify the equations, describing two-phase flows, are considered. Two-component and one-component models of two-phase flow, as well as Zuber and Findlay model for actual volumetric steam content, and Wallis model, describing the given phase rates, are presented. The conclusion is made, that the two-component model, in which values averaged in time are included, is applicable for the solving of three-dimensional tasks for unsteady two-phase flow. At the same time, using the two-component model, including values, averaged in space only one-dimensional tasks for unsteady two-phase flow can be solved

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

Energy Technology Data Exchange (ETDEWEB)

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

1997-06-01

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

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

International Nuclear Information System (INIS)

Goldberg, L.F.

1990-08-01

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

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

International Nuclear Information System (INIS)

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

1977-08-01

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

Obtaining of Analytical Relations for Hydraulic Parameters of Channels With Two Phase Flow Using Open CFD Toolbox

Science.gov (United States)

Varseev, E.

2017-11-01

The present work is dedicated to verification of numerical model in standard solver of open-source CFD code OpenFOAM for two-phase flow simulation and to determination of so-called “baseline” model parameters. Investigation of heterogeneous coolant flow parameters, which leads to abnormal friction increase of channel in two-phase adiabatic “water-gas” flows with low void fractions, presented.

Convective heat transfer from rough surfaces with two-dimensional ribs - transitional and laminar flow

International Nuclear Information System (INIS)

Dalle Donne, M.; Meyer, L.

1978-01-01

Measurements of friction factor and heat transfer coefficients for two rods of 18.9 mm 0.D. with two-dimensional roughness, each in two different outer smooth tubes have been performed in turbulent and laminar flow. The turbulent flow results indicate that the flow was not thermally fully established, the isothermal data however agree reasonably well with our previously obtained general correlation. Laminar flow results can be correlated best when the Reynolds and Greatz numbers are evaluated at the temperature average between the temperature of the inner rod surface and of the outer smooth surface of the annulus, the average being weighted over the two surfaces. (orig.) [de

Transient response in granular quasi-two-dimensional bounded heap flow.

Science.gov (United States)

Xiao, Hongyi; Ottino, Julio M; Lueptow, Richard M; Umbanhowar, Paul B

2017-10-01

We study the transition between steady flows of noncohesive granular materials in quasi-two-dimensional bounded heaps by suddenly changing the feed rate. In both experiments and simulations, the primary feature of the transition is a wedge of flowing particles that propagates downstream over the rising free surface with a wedge front velocity inversely proportional to the square root of time. An additional longer duration transient process continues after the wedge front reaches the downstream wall. The entire transition is well modeled as a moving boundary problem with a diffusionlike equation derived from local mass balance and a local relation between the flux and the surface slope.

Screened Raman response in two-dimensional d(x2-y2)-wave superconductors: Relative intensities in different symmetry channels

DEFF Research Database (Denmark)

Wenger, F.; Käll, M.

1997-01-01

We analyze the Raman-scattering response in a two-dimensional d(x2-y2)-wave superconductor and point out a strong suppression of relative intensity in the screened A(1g) channel compared to the B-1g channel for a generic tight-binding model. This is in contrast with the observed behavior in high...

Wall-cooling-induced mixed-convection flow recirculation in a vertical square-array multi-rod channel

International Nuclear Information System (INIS)

Luangdilok, W.; Todreas, N.E.

1989-01-01

This work investigated the structure of penetrative flow recirculation and associated flow conditions in a multi-rod channel induced by interassembly heat transfer that causes cooling through channel walls. Three investigation approaches, experimental, numerical, and analytical were employed in a complimentary fashion. Physical experiments involved water flow visualization and temperature measurement in a 4x4 rod square channel. Numerical experiments involved 3-dimensional simulations of water and sodium flow in a 2x2-rod channels. An approximate reverse flow model including Prandtl number effect was developed. A correlating equation based on the model and experiments was verified for water to correctly predict the trend of the 4x4-rod experimental penetration depth data. (orig.)

Synchrotron 4-dimensional imaging of two-phase flow through porous media.

Science.gov (United States)

Kim, F H; Penumadu, D; Patel, P; Xiao, X; Garboczi, E J; Moylan, S P; Donmez, M A

2016-01-01

Near real-time visualization of complex two-phase flow in a porous medium was demonstrated with dynamic 4-dimensional (4D) (3D + time) imaging at the 2-BM beam line of the Advanced Photon Source (APS) at Argonne National Laboratory. Advancing fluid fronts through tortuous flow paths and their interactions with sand grains were clearly captured, and formations of air bubbles and capillary bridges were visualized. The intense X-ray photon flux of the synchrotron facility made 4D imaging possible, capturing the dynamic evolution of both solid and fluid phases. Computed Tomography (CT) scans were collected every 12 s with a pixel size of 3.25 µm. The experiment was carried out to improve understanding of the physics associated with two-phase flow. The results provide a source of validation data for numerical simulation codes such as Lattice-Boltzmann, which are used to model multi-phase flow through porous media.

A two-phase flow regime map for a MAPLE-type nuclear research reactor fuel channel: Effect of hexagonal finned bundle

International Nuclear Information System (INIS)

Harvel, G.D.; Chang, J.S.

1997-01-01

A two-phase flow regime map is developed experimentally and theoretically for a vertical hexagonal flow channel with and without a 36-finned rod hexagonal bundle. This type of flow channel is of interest to MAPLE-type nuclear research reactors. The flow regime maps are determined by visual observations and observation of waveforms shown by a capacitance-type void fraction meter. The experimental results show that the inclusion of the finned hexagonal bundle shifts the flow regime transition boundaries toward higher water flow rates. Existing flow regime maps based on pipe flow require slight modifications when applied to the hexagonal flow channel with and without a MAPLE-type finned hexagonal bundle. The proposed theoretical model agrees well with experimental results

Study of gas-water flow in horizontal rectangular channels

Science.gov (United States)

Chinnov, E. A.; Ron'shin, F. V.; Kabov, O. A.

2015-09-01

The two-phase flow in the narrow short horizontal rectangular channels 1 millimeter in height was studied experimentally. The features of formation of the two-phase flow were studied in detail. It is shown that with an increase in the channel width, the region of the churn and bubble regimes increases, compressing the area of the jet flow. The areas of the annular and stratified flow patterns vary insignificantly.

Comparative study of the two-fluid momentum equations for multi-dimensional bubbly flows: Modification of Reynolds stress

Energy Technology Data Exchange (ETDEWEB)

Lee, Seung Jun; Park, Ik Kyu; Yoon, Han Young [Thermal-Hydraulic Safety Research Division, Korea Atomic Energy Research Institute, Daejeon (Korea, Republic of); Jae, Byoung [School of Mechanical Engineering, Chungnam National University, Daejeon (Korea, Republic of)

2017-01-15

Two-fluid equations are widely used to obtain averaged behaviors of two-phase flows. This study addresses a problem that may arise when the two-fluid equations are used for multi-dimensional bubbly flows. If steady drag is the only accounted force for the interfacial momentum transfer, the disperse-phase velocity would be the same as the continuous-phase velocity when the flow is fully developed without gravity. However, existing momentum equations may show unphysical results in estimating the relative velocity of the disperse phase against the continuous-phase. First, we examine two types of existing momentum equations. One is the standard two-fluid momentum equation in which the disperse-phase is treated as a continuum. The other is the averaged momentum equation derived from a solid/ fluid particle motion. We show that the existing equations are not proper for multi-dimensional bubbly flows. To resolve the problem mentioned above, we modify the form of the Reynolds stress terms in the averaged momentum equation based on the solid/fluid particle motion. The proposed equation shows physically correct results for both multi-dimensional laminar and turbulent flows.

Inception of supraglacial channelization under turbulent flow conditions

Science.gov (United States)

Mantelli, E.; Camporeale, C.; Ridolfi, L.

2013-12-01

Glacier surfaces exhibit an amazing variety of meltwater-induced morphologies, ranging from small scale ripples and dunes on the bed of supraglacial channels to meandering patterns, till to large scale drainage networks. Even though the structure and geometry of these morphologies play a key role in the glacier melting processes, the physical-based modeling of such spatial patterns have attracted less attention than englacial and subglacial channels. In order to partially fill this gap, our work concerns the large scale channelization occurring on the ice slopes and focuses on the role of turbulence on the wavelength selection processes during the channelization inception. In a recent study[1], two of us showed that the morphological instability induced by a laminar film flowing over an ice bed is characterized by transversal length scales of order of centimeters. Being these scales much smaller than the spacing observed in the channelization of supraglacial drainage networks (that are of order of meters) and considering that the water films flowing on glaciers can exhibit Reynolds numbers larger than 104, we investigated the role of turbulence in the inception of channelization. The flow-field is modeled by means of two-dimensional shallow water equations, where Reynolds stresses are also considered. In the depth-averaged heat balance equation an incoming heat flux from air is assumed and forced convection heat exchange with the wall is taken into account, in addition to convection and diffusion in the liquid. The temperature profile in the ice is finally coupled to the liquid through Stefan equation. We then perform a linear stability analysis and, under the assumption of small Stefan number, we solve the differential eigenvalue problem analytically. As main outcome of such an analysis, the morphological instability of the ice-water interface is detected and investigated in a wide range of the independent parameters: longitudinal and transversal wavenumbers

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

International Nuclear Information System (INIS)

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

1997-01-01

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

Channel flow structure measurements using particle image velocimetry

International Nuclear Information System (INIS)

Norazizi Mohamed; Noraeini Mokhtar; Aziz Ibrahim; Ramli Abu Hassan

1996-01-01

Two different flow structures in a laboratory channel were examined using a flow visualization technique, known as Particle Image Velocimetry (PIV). The first channel flow structure was that of a steady flow over a horizontal channel bottom. Photographs of particle displacements were taken in the boundary layer in a plane parallel to the flow. These photographs were analyzed to give simultaneous measurements of two components of the velocity at hundreds of points in the plane. Averaging these photographs gave the velocity profile a few millimeters from the bottom of the channel to the water surface. The results gave good agreement with the known boundary layer theory. This technique is extended to the study of the structure under a progressive wave in the channel. A wavelength of the propagating wave is divided into sections by photographing it continously for a number of frames. Each frame is analyzed and a velocity field under this wave at various phase points were produced with their respective directions. The results show that velocity vectors in a plane under the wave could be achieved instantaneously and in good agreement with the small amplitude wave theory

Numerical investigation of flow instability in parallel channels with supercritical water

International Nuclear Information System (INIS)

Shitsi, Edward; Debrah, Seth Kofi; Agbodemegbe, Vincent Yao; Ampomah-Amoako, Emmanuel

2017-01-01

Highlights: •Supercritical flow instability in parallel channels is investigated. •Flow dynamics and heat transfer characteristics are analyzed. •Mass flow rate, pressure, heating power, and axial power shape have significant effects on flow instability. •Numerical results are validated with experimental results. -- Abstract: SCWR is one of the selected Gen IV reactors purposely for electricity generation in the near future. It is a promising technology with higher efficiency compared to current LWRs but without the challenges of heat transfer and its associated flow instability. Supercritical flow instability is mainly caused by sharp change in the coolant properties around the pseudo-critical point of the working fluid and research into this phenomenon is needed to address concerns of flow instability at supercritical pressures. Flow instability in parallel channels at supercritical pressures is investigated in this paper using a three dimensional (3D) numerical tool (STAR-CCM+). The dynamics characteristics such as amplitude and period of out-of-phase inlet mass flow oscillation at the heated channel inlet, and heat transfer characteristic such as maximum outlet temperature of the heated channel outlet temperature oscillation are discussed. Influences of system parameters such as axial power shape, pressure, mass flow rate, and gravity are discussed based on the obtained mass flow and temperature oscillations. The results show that the system parameters have significant effect on the amplitude of the mass flow oscillation and maximum temperature of the heated outlet temperature oscillation but have little effect on the period of the mass flow oscillation. The amplitude of mass flow oscillation and maximum temperature of the heated channel outlet temperature oscillation increase with heating power. The numerical results when compared to experiment data show that the 3D numerical tool (STAR-CCM+) could capture dynamics and heat transfer characteristics of

Two dimensional modelling of flood flows and suspended sedimenttransport: the case of the Brenta River, Veneto (Italy

Directory of Open Access Journals (Sweden)

P. Martini

2004-01-01

Full Text Available The paper presents a numerical model for the simulation of flood waves and suspended sediment transport in a lowland river basin of North Eastern Italy. The two dimensional depth integrated momentum and continuity equations are modified to take into account the bottom irregularities that strongly affect the hydrodynamics in partially dry areas, as for example, in the first stages of an inundation process or in tidal flow. The set of equations are solved with a standard Galerkin finite element method using a semi-implicit numerical scheme where the effects of both the small channel network and the regulation devices on the flood wave propagation are accounted for. Transport of suspended sediment and bed evolution are coupled with the hydrodynamics using an appropriate form of the advection-dispersion equation and Exner's equation. Applications to a case study are presented in which the effects of extreme flooding on the Brenta River (Italy are examined. Urban and rural flood risk areas are identified and the effects of a alleviating action based on a diversion channel flowing into Venice Lagoon are simulated. The results show that this solution strongly reduces the flood risk in the downstream areas and can provide an important source of sediment for the Venice Lagoon. Finally, preliminary results of the sediment dispersion due to currents and waves in the Venice Lagoon are presented.

Two dimensional modelling of flood flows and suspended sedimenttransport: the case of the Brenta River, Veneto (Italy)

Science.gov (United States)

Martini, P.; Carniello, L.; Avanzi, C.

2004-03-01

The paper presents a numerical model for the simulation of flood waves and suspended sediment transport in a lowland river basin of North Eastern Italy. The two dimensional depth integrated momentum and continuity equations are modified to take into account the bottom irregularities that strongly affect the hydrodynamics in partially dry areas, as for example, in the first stages of an inundation process or in tidal flow. The set of equations are solved with a standard Galerkin finite element method using a semi-implicit numerical scheme where the effects of both the small channel network and the regulation devices on the flood wave propagation are accounted for. Transport of suspended sediment and bed evolution are coupled with the hydrodynamics using an appropriate form of the advection-dispersion equation and Exner's equation. Applications to a case study are presented in which the effects of extreme flooding on the Brenta River (Italy) are examined. Urban and rural flood risk areas are identified and the effects of a alleviating action based on a diversion channel flowing into Venice Lagoon are simulated. The results show that this solution strongly reduces the flood risk in the downstream areas and can provide an important source of sediment for the Venice Lagoon. Finally, preliminary results of the sediment dispersion due to currents and waves in the Venice Lagoon are presented.

A soap film shock tube to study two-dimensional compressible flows

Energy Technology Data Exchange (ETDEWEB)

Wen, C.Y.; Chen, Y.M.; Chang-Jian, S.K. [Dept. of Mechanical Engineering, Da-Yeh University Chang-Hwa (Taiwan)

2001-07-01

A new experimental approach to the study of the two-dimensional compressible flow phenomena is presented. In this technique, a variety of compressible flows were generated by bursting plane vertical soap films. An aureole and a ''shock wave'' preceding the rim of the expanding hole were clearly observed using traditional high-speed flash photography and a fast line-scan charge coupled device (CCD) camera. The moving shock wave images obtained from the line-scan CCD camera were similar to the x-t diagrams in gas dynamics. The moving shock waves cause thickness jumps and induce supersonic flows. Photographs of the supersonic flows over a cylinder and a wedge are presented. The results suggest clearly the feasibility of the ''soap film shock tube''. (orig.)

Improving flow distribution in influent channels using computational fluid dynamics.

Science.gov (United States)

Park, No-Suk; Yoon, Sukmin; Jeong, Woochang; Lee, Seungjae

2016-10-01

Although the flow distribution in an influent channel where the inflow is split into each treatment process in a wastewater treatment plant greatly affects the efficiency of the process, and a weir is the typical structure for the flow distribution, to the authors' knowledge, there is a paucity of research on the flow distribution in an open channel with a weir. In this study, the influent channel of a real-scale wastewater treatment plant was used, installing a suppressed rectangular weir that has a horizontal crest to cross the full channel width. The flow distribution in the influent channel was analyzed using a validated computational fluid dynamics model to investigate (1) the comparison of single-phase and two-phase simulation, (2) the improved procedure of the prototype channel, and (3) the effect of the inflow rate on flow distribution. The results show that two-phase simulation is more reliable due to the description of the free-surface fluctuations. It should first be considered for improving flow distribution to prevent a short-circuit flow, and the difference in the kinetic energy with the inflow rate makes flow distribution trends different. The authors believe that this case study is helpful for improving flow distribution in an influent channel.

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

International Nuclear Information System (INIS)

Neves Conti, T. das.

1983-01-01

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

Numerical simulation of secondary flow in bubbly turbulent flow in sub-channel

International Nuclear Information System (INIS)

Ikeno, Tsutomu; Kataoka, Isao

2009-01-01

Secondary flow in bubbly turbulent flow in sub-channel was simulated by using an algebraic turbulence stress model. The mass, momentum, turbulence energy and bubble diffusion equations were used as fundamental equation. The basis for these equations was the two-fluid model: the equation of liquid phase was picked up from the equation system theoretically derived for the gas-liquid two-fluid turbulent flow. The fundamental equation was transformed onto a generalized coordinate system fitted to the computational domain in sub-channel. It was discretized for the SIMPLE algorism using the finite-volume method. The shape of sub-channel causes a distortion of the computational mesh, and orthogonal nature of the mesh is sometimes broken. An iterative method to satisfy a requirement for the contra-variant velocity was introduced to represent accurate symmetric boundary condition. Two-phase flow at a steady state was simulated for different magnitude of secondary flow and void fraction. The secondary flow enhanced the momentum transport in sub-channel and accelerated the liquid phase in the rod gap. This effect was slightly mitigated when the void fraction increased. The acceleration can contribute to effective cooling in the rod gap. The numerical result implied a phenomenon of industrial interest. This suggested that experimental approach is necessary to validate the numerical model and to identify the phenomenon. (author)

Isotachophoresis system having larger-diameter channels flowing into channels with reduced diameter and with selectable counter-flow

Energy Technology Data Exchange (ETDEWEB)

Mariella, Jr., Raymond P.

2018-03-06

An isotachophoresis system for separating a sample containing particles into discrete packets including a flow channel, the flow channel having a large diameter section and a small diameter section; a negative electrode operably connected to the flow channel; a positive electrode operably connected to the flow channel; a leading carrier fluid in the flow channel; a trailing carrier fluid in the flow channel; and a control for separating the particles in the sample into discrete packets using the leading carrier fluid, the trailing carrier fluid, the large diameter section, and the small diameter section.

Numerical study of mixed convection heat transfer enhancement in a channel with active flow modulation

Science.gov (United States)

Billah, Md. Mamun; Khan, Md Imran; Rahman, Mohammed Mizanur; Alam, Muntasir; Saha, Sumon; Hasan, Mohammad Nasim

2017-06-01

A numerical study of steady two dimensional mixed convention heat transfer phenomena in a rectangular channel with active flow modulation is carried out in this investigation. The flow in the channel is modulated via a rotating cylinder placed at the center of the channel. In this study the top wall of the channel is subjected to an isothermal low temperature while a discrete isoflux heater is positioned on the lower wall. The fluid flow under investigation is assumed to have a Prandtl number of 0.71 while the Reynolds No. and the Grashof No. are varied in wide range for four different situations such as: i) plain channel with no cylinder, ii) channel with stationary cylinder, iii) channel with clockwise rotating cylinder and iv) channel with counter clockwise rotating cylinder. The results obtained in this study are presented in terms of the distribution of streamlines, isotherms in the channel while the heat transfer process from the heat source is evaluated in terms of the local Nusselt number, average Nusselt number. The outcomes of this study also indicate that the results are strongly dependent on the type of configuration and direction of rotation of the cylinder and that the average Nusselt number value rises with an increase in Reynolds and Grashof numbers but the correlation between these parameters at higher values of Reynolds and Grashof numbers becomes weak.

Noninteracting beams of ballistic two-dimensional electrons

International Nuclear Information System (INIS)

Spector, J.; Stormer, H.L.; Baldwin, K.W.; Pfeiffer, L.N.; West, K.W.

1991-01-01

We demonstrate that two beams of two-dimensional ballistic electrons in a GaAs-AlGaAs heterostructure can penetrate each other with negligible mutual interaction analogous to the penetration of two optical beams. This allows electrical signal channels to intersect in the same plane with negligible crosstalk between the channels

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

International Nuclear Information System (INIS)

Chawla, T.C.; Ishii, M.

1980-01-01

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

Flow Patterns in an Open Channel Confluence with Increasingly Dominant Tributary Inflow

Directory of Open Access Journals (Sweden)

Laurent Schindfessel

2015-08-01

Full Text Available Despite the ratio of incoming discharges being recognized as a key parameter in open-channel confluence hydrodynamics, little is known about the flow patterns when the tributary provides more than 90% of the total discharge. This paper offers a systematic study of flow features when the tributary becomes increasingly dominant in a 90° confluence with a fixed concordant bed. Large-eddy simulations are used to investigate the three-dimensional complex flow patterns for three different discharge ratios. It is found that the tributary flow impinges on the opposing bank when the tributary flow becomes sufficiently dominant, causing a recirculating eddy in the upstream channel of the confluence, which induces significant changes in the incoming velocity distribution. Moreover, it results in stronger helicoidal cells in the downstream channel, along with zones of upwelling flow. In turn, the changed flow patterns also influence the mixing layer and the flow recovery. Finally, intermittent events of stronger upwelling flow are discerned. Improved understanding of flow patterns at confluences where the tributary is dominant is applicable to both engineering and earth sciences.

Present status of numerical analysis on transient two-phase flow

International Nuclear Information System (INIS)

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

1987-01-01

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

Novel target design algorithm for two-dimensional optical storage (TwoDOS)

NARCIS (Netherlands)

Huang, Li; Chong, T.C.; Vijaya Kumar, B.V.K.; Kobori, H.

2004-01-01

In this paper we introduce the Hankel transform based channel model of Two-Dimensional Optical Storage (TwoDOS) system. Based on this model, the two-dimensional (2D) minimum mean-square error (MMSE) equalizer has been derived and applied to some simple but common cases. The performance of the 2D

MMOSS-I: a CANDU multiple-channel thermosyphoning flow stability model

Energy Technology Data Exchange (ETDEWEB)

Gulshani, P [Atomic Energy of Canada Ltd., Mississauga, ON (Canada); Huynh, H [Hydro-Quebec, Montreal, PQ (Canada)

1996-12-31

This paper presents a multiple-channel flow stability model, dubbed MMOSS, developed to predict the conditions for the onset of flow oscillations in a CANDU-type multiple-channel heat transport system under thermosyphoning conditions. The model generalizes that developed previously to account for the effects of any channel flow reversal. Two-phase thermosyphoning conditions are predicted by thermalhydraulic codes for some postulated accident scenarios in CANDU. Two-phase thermosyphoning experiments in the multiple-channel RD-14M facility have indicated that pass-to-pass out-of-phase oscillations in the loop conditions caused the flow in some of the heated channels to undergo sustained reversal in direction. This channel flow reversal had significant effects on the channel and loop conditions. It is, therefore, important to understand the nature of the oscillations and be able to predict the conditions for the onset of the oscillations or for stable flow in RD-14M and the reactor. For stable flow conditions, oscillation-induced channel flow reversal is not expected. MMOSS was developed for a figure-of-eight system with any number of channels. The system characteristic equation was derived from a linearization of the conservation equations. In this paper, the MMOSS characteristic equation is solved for a system of N identical channel assemblies. The resulting model is called MMOSS-I. This simplification provides valuable physical insight and reasonably accurate results. MMOSS-I and a previously-developed steady-state model THERMOSYPHON are used to predict thermosyphoning flow stability maps for RD-14M and the Gentilly 2 reactor. (author). 11 refs., 7 figs.

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

CERN Document Server

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

1999-01-01

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

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

Science.gov (United States)

Na, Yun Whan

2008-07-01

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

Three-Dimensional Flow Generated by a Partially Penetrating Well in a Two-Aquifer System

Science.gov (United States)

Sepulveda, N.

2007-12-01

An analytical solution is presented for three-dimensional (3D) flow in a confined aquifer and the overlying storative semiconfining layer and unconfined aquifer. The equation describing flow caused by a partially penetrating production well is solved analytically to provide a method to accurately determine the hydraulic parameters in the confined aquifer, semiconfining layer, and unconfined aquifer from aquifer-test data. Previous solutions for a partially penetrating well did not account for 3D flow or storativity in the semiconfining unit. The 3D and two- dimensional (2D) flow solutions in the semiconfining layer are compared for various hydraulic conductivity ratios between the aquifer and the semiconfining layer. Analysis of the drawdown data from an aquifer test in central Florida showed that the 3D solution in the semiconfining layer provides a more unique identification of the hydraulic parameters than the 2D solution. The analytical solution could be used to analyze, with higher accuracy, the effect that pumping water from the lower aquifer in a two-aquifer system has on wetlands.

Advanced porous electrodes with flow channels for vanadium redox flow battery

Science.gov (United States)

Bhattarai, Arjun; Wai, Nyunt; Schweiss, Ruediger; Whitehead, Adam; Lim, Tuti M.; Hng, Huey Hoon

2017-02-01

Improving the overall energy efficiency by reducing pumping power and improving flow distribution of electrolyte, is a major challenge for developers of flow batteries. The use of suitable channels can improve flow distribution through the electrodes and reduce flow resistance, hence reducing the energy consumption of the pumps. Although several studies of vanadium redox flow battery have proposed the use of bipolar plates with flow channels, similar to fuel cell designs, this paper presents the use of flow channels in the porous electrode as an alternative approach. Four types of electrodes with channels: rectangular open channel, interdigitated open cut channel, interdigitated circular poked channel and cross poked circular channels, are studied and compared with a conventional electrode without channels. Our study shows that interdigitated open channels can improve the overall energy efficiency up to 2.7% due to improvement in flow distribution and pump power reduction while interdigitated poked channel can improve up to 2.5% due to improvement in flow distribution.

Determination of the onset nonlinearity hydrodynamic characteristics at two-phase flow in parallel vertical channels

International Nuclear Information System (INIS)

Jovic, V.; Afgan, N.; Jovic, L.; Spasojevic, D.

1993-01-01

The paper presents results of the experimental and theoretical analyses of linear and nonlinear characteristics of adiabatic two-phase water-air flow in vertical parallel channels. Regime character changes and linear to nonlinear dynamic characteristics transfer conditions were defined. (author)

Numerical approximation of flow in a symmetric channel with vibrating walls

Directory of Open Access Journals (Sweden)

Sváček P.

2010-12-01

Full Text Available In this paper the numerical solution of two dimensional fluid-structure interaction problem is addressed. The fluid motion is modelled by the incompressible unsteady Navier-Stokes equations. The spatial discretization by stabilized finite element method is used. The motion of the computational domain is treated with the aid of Arbitrary Lagrangian Eulerian (ALE method. The time-space problem is solved with the aid of multigrid method. The method is applied onto a problem of interaction of channel flow with moving walls, which models the air flow in the glottal region of the human vocal tract. The pressure boundary conditions and the effects of the isotropic and anisotropic mesh refinement are discussed. The numerical results are presented.

Experimental investigation of flow over two-dimensional multiple hill models.

Science.gov (United States)

Li, Qing'an; Maeda, Takao; Kamada, Yasunari; Yamada, Keisuke

2017-12-31

The aim of this study is to investigate the flow field characteristics in ABL (Atmospheric Boundary Layer) flow over multiple hills and valleys in two-dimensional models under neutral conditions. Active turbulence grids and boundary layer generation frame were used to simulate the natural winds in wind tunnel experiments. As a result, the mean wind velocity, the velocity vector diagram and turbulence intensity around the hills were investigated by using a PIV (Particle Image Velocimetry) system. From the measurement results, it was known that the average velocity was increased along the upstream slope of upside hill, and then separated at the top of the hills, the acceleration region of U/U ref >1 was generated at the downstream of the hill. Meanwhile, a large clockwise circulation flow was generated between the two hill models. Moreover, the turbulence intensity showed small value in the circulation flow regions. Compared to 1H model, the turbulence intensity in the mainstream direction showed larger value than that in the vertical direction. This paper provided a better understanding of the wind energy distribution on the terrain for proper selection of suitable sites for installing wind farms in the ABL. Copyright © 2017 Elsevier B.V. All rights reserved.

Numerical Simulation of Gas-Solid Two-Phase Flow for Four-Channels Pulverized Swirling Burner

Directory of Open Access Journals (Sweden)

Defu LI

2013-05-01

Full Text Available This article presents a mathematical model of cold gas-solid two-phase flow which is based on the cement rotary kiln in service. By altering the parameters of air supply system of four- channels pulverized burner, investigations are taken of that motion trajectory and particle distributions in the very turbulent field. The results show that motion trail of most particles in rotary kiln is a combination process of gradual diffusion and slow sedimentation; increasing internal flow velocity would aggravate coal particles to diffuse; external flow velocity should be controlled in a reasonable range.

Effects of stream discharge, alluvial depth and bar amplitude on hyporheic flow in pool-riffle channels

Science.gov (United States)

Daniele Tonina; John M. Buffington

2011-01-01

Hyporheic flow results from the interaction between streamflow and channel morphology and is an important component of stream ecosystems because it enhances water and solute exchange between the river and its bed. Hyporheic flow in pool-riffle channels is particularly complex because of three-dimensional topography that spans a range of partially to fully submerged...

Kolmogorov flow in two dimensional strongly coupled dusty plasma

Energy Technology Data Exchange (ETDEWEB)

Gupta, Akanksha; Ganesh, R., E-mail: ganesh@ipr.res.in; Joy, Ashwin [Institute for Plasma Research, Bhat Gandhinagar, Gujarat 382 428 (India)

2014-07-15

Undriven, incompressible Kolmogorov flow in two dimensional doubly periodic strongly coupled dusty plasma is modelled using generalised hydrodynamics, both in linear and nonlinear regime. A complete stability diagram is obtained for low Reynolds numbers R and for a range of viscoelastic relaxation time τ{sub m} [0 < τ{sub m} < 10]. For the system size considered, using a linear stability analysis, similar to Navier Stokes fluid (τ{sub m} = 0), it is found that for Reynolds number beyond a critical R, say R{sub c}, the Kolmogorov flow becomes unstable. Importantly, it is found that R{sub c} is strongly reduced for increasing values of τ{sub m}. A critical τ{sub m}{sup c} is found above which Kolmogorov flow is unconditionally unstable and becomes independent of Reynolds number. For R < R{sub c}, the neutral stability regime found in Navier Stokes fluid (τ{sub m} = 0) is now found to be a damped regime in viscoelastic fluids, thus changing the fundamental nature of transition of Kolmogorov flow as function of Reynolds number R. A new parallelized nonlinear pseudo spectral code has been developed and is benchmarked against eigen values for Kolmogorov flow obtained from linear analysis. Nonlinear states obtained from the pseudo spectral code exhibit cyclicity and pattern formation in vorticity and viscoelastic oscillations in energy.

HANARO core channel flow-rate measurement

Energy Technology Data Exchange (ETDEWEB)

Kim, Heon Il; Chae, Hee Tae; Im, Don Soon; Kim, Seon Duk [Korea Atomic Energy Research Institute, Taejon (Korea, Republic of)

1996-06-01

HANARO core consists of 23 hexagonal flow tubes and 16 cylindrical flow tubes. To get the core flow distribution, we used 6 flow-rate measuring dummy fuel assemblies (instrumented dummy fuel assemblies). The differential pressures were measured and converted to flow-rates using the predetermined relationship between AP and flow-rate for each instrumented dummy fuel assemblies. The flow-rate for the cylindrical flow channels shows +-7% relative errors and that for the hexagonal flow channels shows +-3.5% relative errors. Generally the flow-rates of outer core channels show smaller values compared to those of inner core. The channels near to the core inlet pipe and outlet pipes also show somewhat lower flow-rates. For the lower flow channels, the thermal margin was checked by considering complete linear power histories. From the experimental results, the gap flow-rate was estimated to be 49.4 kg/s (cf. design flow of 50 kg/s). 15 tabs., 9 figs., 10 refs. (Author) .new.

Dynamics of the free surface of stratified two-phase flows in channels with rectangular cross-sections

International Nuclear Information System (INIS)

Vallee, Christophe

2012-01-01

Stratified two-phase flows were investigated at different test facilities with horizontal test sections in order to provide an experimental database for the development and validation of computational fluid dynamics (CFD) codes. These channels were designed with rectangular cross-sections to enable optimal observation conditions for the application of optical measurement techniques. Consequently, the local flow structure was visualised with a high-speed video camera, delivering data with highresolution in space and time as needed for CFD code validation. Generic investigations were performed at atmospheric pressure and room temperature in two air/water channels made of acrylic glass. Divers preliminary experiments were conducted with various measuring systems in a test section mounted between two separators. The second test facility, the Horizontal Air/Water Channel (HAWAC), is dedicated to co-current flow investigations. The hydraulic jump as the quasi-stationary discontinuous transition between super- and subcritical flow was studied in this closed channel. Moreover, the instable wave growth leading to slug flow was investigated from the test section inlet. For quantitative analysis of the optical measurements, an algorithm was developed to recognise the stratified interface in the camera frames, allowing statistical treatments for comparison with CFD calculation results. The third test apparatus was installed in the pressure chamber of the TOPFLOW test facility in order to be operated at reactor typical conditions under pressure equilibrium with the vessel atmosphere. The test section representing a flat model of the hot leg of the German Konvoi pressurised water reactor (PWR) scaled at 1:3 is equipped with large glass side walls in the region of the elbow and of the steam generator inlet chamber to allow visual observations. The experiments were conducted with air and water at room temperature and maximum pressures of 3 bar as well as with steam and water at

Dynamics of the free surface of stratified two-phase flows in channels with rectangular cross-sections

Energy Technology Data Exchange (ETDEWEB)

Vallee, Christophe

2012-08-22

Stratified two-phase flows were investigated at different test facilities with horizontal test sections in order to provide an experimental database for the development and validation of computational fluid dynamics (CFD) codes. These channels were designed with rectangular cross-sections to enable optimal observation conditions for the application of optical measurement techniques. Consequently, the local flow structure was visualised with a high-speed video camera, delivering data with highresolution in space and time as needed for CFD code validation. Generic investigations were performed at atmospheric pressure and room temperature in two air/water channels made of acrylic glass. Divers preliminary experiments were conducted with various measuring systems in a test section mounted between two separators. The second test facility, the Horizontal Air/Water Channel (HAWAC), is dedicated to co-current flow investigations. The hydraulic jump as the quasi-stationary discontinuous transition between super- and subcritical flow was studied in this closed channel. Moreover, the instable wave growth leading to slug flow was investigated from the test section inlet. For quantitative analysis of the optical measurements, an algorithm was developed to recognise the stratified interface in the camera frames, allowing statistical treatments for comparison with CFD calculation results. The third test apparatus was installed in the pressure chamber of the TOPFLOW test facility in order to be operated at reactor typical conditions under pressure equilibrium with the vessel atmosphere. The test section representing a flat model of the hot leg of the German Konvoi pressurised water reactor (PWR) scaled at 1:3 is equipped with large glass side walls in the region of the elbow and of the steam generator inlet chamber to allow visual observations. The experiments were conducted with air and water at room temperature and maximum pressures of 3 bar as well as with steam and water at

Statistical Mechanics of the Geometric Control of Flow Topology in Two-Dimensional Turbulence

Science.gov (United States)

Nadiga, Balasubramanya; Loxley, Peter

2013-04-01

We apply the principle of maximum entropy to two dimensional turbulence in a new fashion to predict the effect of geometry on flow topology. We consider two prototypical regimes of turbulence that lead to frequently observed self-organized coherent structures. Our theory predicts bistable behavior that exhibits hysteresis and large abrupt changes in flow topology in one regime; the other regime is predicted to exhibit monstable behavior with a continuous change of flow topology. The predictions are confirmed in fully nonlinear numerical simulations of the two-dimensional Navier-Stokes equation. These results suggest an explanation of the low frequency regime transitions that have been observed in the non-equilibrium setting of this problem. Following further development in the non-equilibrium context, we expect that insights developed in this problem should be useful in developing a better understanding of the phenomenon of low frequency regime transitions that is a pervasive feature of the weather and climate systems. Familiar occurrences of this phenomenon---wherein extreme and abrupt qualitative changes occur, seemingly randomly, after very long periods of apparent stability---include blocking in the extra-tropical winter atmosphere, the bimodality of the Kuroshio extension system, the Dansgaard-Oeschger events, and the glacial-interglacial transitions.

Alignment dynamics of diffusive scalar gradient in a two-dimensional model flow

Science.gov (United States)

Gonzalez, M.

2018-04-01

The Lagrangian two-dimensional approach of scalar gradient kinematics is revisited accounting for molecular diffusion. Numerical simulations are performed in an analytic, parameterized model flow, which enables considering different regimes of scalar gradient dynamics. Attention is especially focused on the influence of molecular diffusion on Lagrangian statistical orientations and on the dynamics of scalar gradient alignment.

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

International Nuclear Information System (INIS)

Minato, Akihiko; Naitoh, Masanori

1987-01-01

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

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

International Nuclear Information System (INIS)

Raymond, P.; Toumi, I.

1992-01-01

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

Two-dimensional dynamics of elasto-inertial turbulence and its role in polymer drag reduction

Science.gov (United States)

Sid, S.; Terrapon, V. E.; Dubief, Y.

2018-02-01

The goal of the present study is threefold: (i) to demonstrate the two-dimensional nature of the elasto-inertial instability in elasto-inertial turbulence (EIT), (ii) to identify the role of the bidimensional instability in three-dimensional EIT flows, and (iii) to establish the role of the small elastic scales in the mechanism of self-sustained EIT. Direct numerical simulations of viscoelastic fluid flows are performed in both two- and three-dimensional straight periodic channels using the Peterlin finitely extensible nonlinear elastic model (FENE-P). The Reynolds number is set to Reτ=85 , which is subcritical for two-dimensional flows but beyond the transition for three-dimensional ones. The polymer properties selected correspond to those of typical dilute polymer solutions, and two moderate Weissenberg numbers, Wiτ=40 ,100 , are considered. The simulation results show that sustained turbulence can be observed in two-dimensional subcritical flows, confirming the existence of a bidimensional elasto-inertial instability. The same type of instability is also observed in three-dimensional simulations where both Newtonian and elasto-inertial turbulent structures coexist. Depending on the Wi number, one type of structure can dominate and drive the flow. For large Wi values, the elasto-inertial instability tends to prevail over the Newtonian turbulence. This statement is supported by (i) the absence of typical Newtonian near-wall vortices and (ii) strong similarities between two- and three-dimensional flows when considering larger Wi numbers. The role of small elastic scales is investigated by introducing global artificial diffusion (GAD) in the hyperbolic transport equation for polymers. The aim is to measure how the flow reacts when the smallest elastic scales are progressively filtered out. The study results show that the introduction of large polymer diffusion in the system strongly damps a significant part of the elastic scales that are necessary to feed

A Numerical Scheme Based on an Immersed Boundary Method for Compressible Turbulent Flows with Shocks: Application to Two-Dimensional Flows around Cylinders

Directory of Open Access Journals (Sweden)

Shun Takahashi

2014-01-01

Full Text Available A computational code adopting immersed boundary methods for compressible gas-particle multiphase turbulent flows is developed and validated through two-dimensional numerical experiments. The turbulent flow region is modeled by a second-order pseudo skew-symmetric form with minimum dissipation, while the monotone upstream-centered scheme for conservation laws (MUSCL scheme is employed in the shock region. The present scheme is applied to the flow around a two-dimensional cylinder under various freestream Mach numbers. Compared with the original MUSCL scheme, the minimum dissipation enabled by the pseudo skew-symmetric form significantly improves the resolution of the vortex generated in the wake while retaining the shock capturing ability. In addition, the resulting aerodynamic force is significantly improved. Also, the present scheme is successfully applied to moving two-cylinder problems.

Influence of a flow obstacle on the occurrence of burnout in boiling two-phase upward flow within a vertical annular channel

Energy Technology Data Exchange (ETDEWEB)

Mori, S.; Fukano, T. [Kyushu Univ., Fukuoka (Japan)

2003-07-01

When a flow obstruction such as a cylindrical spacer is set in a boiling two-phase flow with-in an annular channel, the inner tube of which is used as a heater, the temperature on the surface of the heating tube is severely affected by its existence. In some cases the cylindrical spacer has a cooling effect, and in the other cases it causes the dryout of the cooling water film on the heating surface resulting in the burnout of the heating tube. In the present paper we have focused our attention on the influence of a flow obstacle on the occurrence of burnout of the heating tube in boiling two-phase flow.

A model for the two-point velocity correlation function in turbulent channel flow

International Nuclear Information System (INIS)

Sahay, A.; Sreenivasan, K.R.

1996-01-01

A relatively simple analytical expression is presented to approximate the equal-time, two-point, double-velocity correlation function in turbulent channel flow. To assess the accuracy of the model, we perform the spectral decomposition of the integral operator having the model correlation function as its kernel. Comparisons of the empirical eigenvalues and eigenfunctions with those constructed from direct numerical simulations data show good agreement. copyright 1996 American Institute of Physics

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

Directory of Open Access Journals (Sweden)

Shao-Yiu Hsu

2017-01-01

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

Computation of gradually varied flow in compound open channel ...

Indian Academy of Sciences (India)

The flow of water in an open channel can be treated as steady, gradually varied flow for ... channel between two nodes is treated as a single reach to calculate the loss ... dition at control points and (iii) critical depth is also required to verify the ...

Numerical investigation of thermal-hydraulic performance of channel with protrusions by turbulent cross flow jet

Science.gov (United States)

Sahu, M. K.; Pandey, K. M.; Chatterjee, S.

2018-05-01

In this two dimensional numerical investigation, small rectangular channel with right angled triangular protrusions in the bottom wall of test section is considered. A slot nozzle is placed at the middle of top wall of channel which impinges air normal to the protruded surface. A duct flow and nozzle flow combined to form cross flow which is investigated for heat transfer enhancement of protruded channel. The governing equations for continuity, momentum, energy along with SST k-ω turbulence model are solved with finite volume based Computational fluid dynamics code ANSYS FLUENT 14.0. The range of duct Reynolds number considered for this analysis is 8357 to 51760. The ratios of pitch of protrusion to height of duct considered are 0.5, 0.64 and 0.82. The ratios of height of protrusion to height of duct considered are 0.14, 0.23 and 0.29. The effect of duct Reynolds number, pitch and height of protrusion on thermal-hydraulic performance is studied under cross flow condition. It is found that heat transfer rate is more at relatively larger pitch and small pressure drop is found in case of low height of protrusion.

Flow field optimization for proton exchange membrane fuel cells with varying channel heights and widths

International Nuclear Information System (INIS)

Wang Xiaodong; Huang Yuxian; Cheng, C.-H.; Jang, J.-Y.; Lee, D.-J.; Yan, W.-M.; Su Ay

2009-01-01

The optimal cathode flow field design of a single serpentine proton exchange membrane fuel cell is obtained by adopting a combined optimization procedure including a simplified conjugate-gradient method (SCGM) and a completely three-dimensional, two-phase, non-isothermal fuel cell model. The cell output power density P cell is the objective function to be maximized with channel heights, H 1 -H 5 , and channel widths, W 2 -W 5 as search variables. The optimal design has tapered channels 1, 3 and 4, and diverging channels 2 and 5, producing 22.51% increment compared with the basic design with all heights and widths setting as 1 mm. Reduced channel heights of channels 2-4 significantly enhance sub-rib convection to effectively transport oxygen to and liquid water out of diffusion layer. The final diverging channel prevents significant leakage of fuel to outlet via sub-rib convection from channel 4. Near-optimal design without huge loss in cell performance but is easily manufactured is discussed.

Fluid Mechanics of Taylor Bubbles and Slug Flows in Vertical Channels

International Nuclear Information System (INIS)

Anglart, Henryk; Podowski, Michael Z.

2002-01-01

Fluid mechanics of Taylor bubbles and slug flows is investigated in vertical, circular channels using detailed, three-dimensional computational fluid dynamics simulations. The Volume of Fluid model with the interface-sharpening algorithm, implemented in the commercial CFX4 code, is used to predict the shape and velocity of Taylor bubbles moving along a vertical channel. Several cases are investigated, including both a single Taylor bubble and a train of bubbles rising in water. It is shown that the potential flow solution underpredicts the water film thickness around Taylor bubbles. Furthermore, the computer simulations that are performed reveal the importance of properly modeling the three-dimensional nature of phenomena governing the motion of Taylor bubbles. Based on the present results, a new formula for the evaluation of bubble shape is derived. Both the shape of Taylor bubbles and the bubble rise velocity predicted by the proposed model agree well with experimental observations. Furthermore, the present model shows good promise in predicting the coalescence of Taylor bubbles

NUMERICAL SIMULATION OF FLOW OVER TWO-DIMENSIONAL MOUNTAIN RIDGE USING SIMPLE ISENTROPIC MODEL

Directory of Open Access Journals (Sweden)

Siswanto Siswanto

2009-07-01

Full Text Available Model sederhana isentropis telah diaplikasikan untuk mengidentifikasi perilaku aliran masa udara melewati topografi sebuah gunung. Dalam model isentropis, temperature potensial θ digunakan sebagai koordinat vertikal dalam rezim aliran adiabatis. Medan angin dalam arah vertikal dihilangkan dalam koordinat isentropis sehingga mereduksi sistim tiga dimensi menjadi sistim dua dimensi lapisan θ. Skema komputasi beda hingga tengah telah digunakan untuk memformulasikan model adveksi. Paper ini membahas aplikasi sederhana dari model isentropis untuk mempelajari gelombang gravitasi dan fenomena angin gunung  dengan desain komputasi periodik dan kondisi batas lateral serta simulasi dengan topografi yang berbeda.   The aim of this work is to study turbulent flow over two-dimensional hill using a simple isentropic model. The isentropic model is represented by applying the potential temperature θ, as the vertical coordinate and is conversed in adiabatic flow regimes. This implies a vanishing vertical wind in isentropic coordinates which reduces the three dimensional system to a stack of two dimensional θ–layers. The equations for each isentropic layer are formally identical with the shallow water equation. A computational scheme of centered finite differences is used to formulate an advective model. This work reviews a simple isentropic model application to investigate gravity wave and mountain wave phenomena regard to different experimental design of computation and topographic height.

Investigation on flow patterns and transition characteristics in a tube-bundle channel

International Nuclear Information System (INIS)

Xiang Wenyuan; Lu Yonghong; Zhao Guisheng

2012-01-01

Tube-bundle channels have been widely used in condenser-evaporator and other industrial heat-exchange equipment. The characteristics of two-phase flow patterns and their transitions for refrigerant R-113 through a vertical tube-bundle channel are experimentally investigated using high-speed camera. Experiments show that there are four main flow patterns in the tube-bundle channel, which are bubbly flow, bubbly-churn flow, churn flow and annular flow. And in the same cross-section of tube- bundle channels, it is shown that there might be different flow patterns in different sub-channels. The flow pattern transitions exhibit unsynchronized in different sub-channels. On the basis of experimental research, the flow pattern map is drawn and analyses are made on the comparison of differences between boiling flow patterns in a circular tube and those in a tube-bundle channel. (authors)

Mixed convective thermally radiative micro nanofluid flow in a stretchable channel with porous medium and magnetic field

Energy Technology Data Exchange (ETDEWEB)

Rauf, A., E-mail: raufamar@ciitsahiwal.edu.pk; Shahzad, S. A.; Meraj, M. A. [Department of Mathematics, Comsats Institute of Information Technology, Sahiwal 57000 (Pakistan); Siddiq, M. K. [Department of CASPAM, Bahauddin Zakariya University, Multan 63000 (Pakistan); Raza, J. [School of Quantitative Sciences, Universiti Utara Malaysia, 06010, Sintok, Kedah (Malaysia)

2016-03-15

A numerical study is carried out for two dimensional steady incompressible mixed convective flow of electrically conductive micro nanofluid in a stretchable channel. The flow is generated due to the stretching walls of the channel immersed in a porous medium. The magnetic field is applied perpendicular to the walls. The impact of radiation, viscous dissipation, thermophoretic and Brownian motion of nanoparticles appear in the energy equation. A numerical technique based on Runge-Kutta-Fehlberg fourth-fifth order (RFK45) method is used to express the solutions of velocity, microrotation, temperature and concentration fields. The dimensionless physical parameters are discussed both in tabular and graphical forms. The results are also found in a good agreement with previously published literature work.

Nonlinear mechanisms of two-dimensional wave-wave transformations in the initially coupled acoustic structure

Science.gov (United States)

Vorotnikov, K.; Starosvetsky, Y.

2018-01-01

The present study concerns two-dimensional nonlinear mechanisms of bidirectional and unidirectional channeling of longitudinal and shear waves emerging in the locally resonant acoustic structure. The system under consideration comprises an oscillatory chain of the axially coupled masses. Each mass of the chain is subject to the local linear potential along the lateral direction and incorporates the lightweight internal rotator. In the present work, we demonstrate the emergence of special resonant regimes of complete bi- and unidirectional transitions between the longitudinal and the shear waves of the locally resonant chain. These regimes are manifested by the two-dimensional energy channeling between the longitudinal and the shear traveling waves in the recurrent as well as the irreversible fashion. We show that the spatial control of the two dimensional energy flow between the longitudinal and the shear waves is solely governed by the motion of the internal rotators. Nonlinear analysis of the regimes of a bidirectional wave channeling unveils their global bifurcation structure and predicts the zones of their spontaneous transitions from a complete bi-directional wave channeling to the one-directional entrapment. An additional regime of a complete irreversible resonant transformation of the longitudinal wave into a shear wave is analyzed in the study. The intrinsic mechanism governing the unidirectional wave reorientation is described analytically. The results of the analysis of both mechanisms are substantiated by the numerical simulations of the full model and are found to be in a good agreement.

Hydrodynamics of vapor-liquid annular dispersed flows in channels with heated rod clusters under unsteady conditions

International Nuclear Information System (INIS)

Kroshilin, A.E.; Kroshilin, V.E.; Nigmatulin, B.I.

1984-01-01

A one-dimensional unsteady hydrodynamic model of vapour-liquid disperse-annular flows in channels with heated fuel rod clusters has been constructed. Regularities in the appearance of critical heat transfer due to the dryout of a near-wall liquid film on rod surfaces in such channels are investigated. The model developed takes into account the main flow regularities in the channels with heated rod clusters. The calculations made have shown that the time before crisis appearance agrees satisfactorily with the experimental data

Experimental analysis of upward vertical two-phase flow in four-cusp channels simulating the conditions of a typical nuclear reactor channel, degraded by a loss of coolant accident

International Nuclear Information System (INIS)

Assad, A.C.A.

1984-01-01

The present work deals with an experimental analysis of upward vertical two-phase flow in channels with circular and four-cusp cross-sections. The latter simulates the conditions of a typical nuclear reactor channel, degraded by a loss of coolant accident. Simultaneous flow of air and water has been employed to simulate adiabatic steam-water flow. The installation of air-water separators helped eliminate instabilities during pressure-drop measurements. The gamma ray attenuation was utilized for the void fraction determination. For the four-cusp geommetry, new criteria for two-phase flow regime transitions have been determined, as well as new correlatins for pressure drop and void fraction, as function of the Lockhart-Martinelli factor and vapour mass-fraction, respectively. (Author) [pt

Direct numerical simulation of turbulent velocity-, pressure- and temperature-fields in channel flows

International Nuclear Information System (INIS)

Goetzbach, G.

1977-10-01

For the simulation of non stationary, three-dimensional, turbulent flow- and temperature-fields in channel flows with constant properties a method is presented which is based on a finite difference scheme of the complete conservation equations for mass, momentum and enthalpie. The fluxes of momentum and heat within the grid cells are described by sub-grid scale models. The sub-grid scale model for momentum introduced here is for the first time applicable to small Reynolds-numbers, rather coarse grids, and channels with space dependent roughness distributions. (orig.) [de

Multi-dimensional two-phase flow measurements in a large-diameter pipe using wire-mesh sensor

International Nuclear Information System (INIS)

Kanai, Taizo; Furuya, Masahiro; Arai, Takahiro; Shirakawa, Kenetsu; Nishi, Yoshihisa; Ueda, Nobuyuki

2011-01-01

The authors developed a method of measurement to determine the multi-dimensionality of two phase flow. A wire-mesh sensor (WMS) can acquire a void fraction distribution at a high temporal and spatial resolution and also estimate the velocity of a vertical rising flow by investigating the signal time-delay of the upstream WMS relative to downstream. Previously, one-dimensional velocity was estimated by using the same point of each WMS at a temporal resolution of 1.0 - 5.0 s. The authors propose to extend this time series analysis to estimate the multi-dimensional velocity profile via cross-correlation analysis between a point of upstream WMS and multiple points downstream. Bubbles behave in various ways according to size, which is used to classify them into certain groups via wavelet analysis before cross-correlation analysis. This method was verified by air-water straight and swirl flows within a large-diameter vertical pipe. A high-speed camera is used to set the parameter of cross-correlation analysis. The results revealed that for the rising straight and swirl flows, large scale bubbles tend to move to the center, while the small bubble is pushed to the outside or sucked into the space where the large bubbles existed. Moreover, it is found that this method can estimate the rotational component of velocity of the swirl flow as well as measuring the multi-dimensional velocity vector at high temporal resolutions of 0.2 s. (author)

Micro-channel convective boiling heat transfer with flow instabilities

International Nuclear Information System (INIS)

Consolini, L.; Thome, J.R.

2009-01-01

Flow boiling heat transfer in micro-channels has attracted much interest in the past decade, and is currently a strong candidate for high performance compact heat sinks, such as those required in electronics systems, automobile air conditioning units, micro-reactors, fuel cells, etc. Currently the literature presents numerous experimental studies on two-phase heat transfer in micro-channels, providing an extensive database that covers many different fluids and operating conditions. Among the noteworthy elements that have been reported in previous studies, is the sensitivity of micro-channel evaporators to oscillatory two-phase instabilities. These periodic fluctuations in flow and pressure drop either result from the presence of upstream compressibility, or are simply due to the interaction among parallel channels in multi-port systems. An oscillating flow presents singular characteristics that are expected to produce an effect on the local heat transfer mechanisms, and thus on the estimation of the two-phase heat transfer coefficients. The present investigation illustrates results for flow boiling of refrigerants R-134a, R-236fa, and R-245fa in a 510 μm circular micro-channel, exposed to various degrees of oscillatory compressible volume instabilities. The data describe the main features of the fluctuations in the temperatures of the heated wall and fluid, and draw attention to the differences in the measured unstable time-averaged heat transfer coefficients with respect to those for stable flow boiling. (author)

Micro-channel convective boiling heat transfer with flow instabilities

Energy Technology Data Exchange (ETDEWEB)

Consolini, L.; Thome, J.R. [Ecole Polytechnique Federale de Lausanne (Switzerland). Lab. de Transfert de Chaleur et de Masse], e-mail: lorenzo.consolini@epfl.ch, e-mail: john.thome@epfl.ch

2009-07-01

Flow boiling heat transfer in micro-channels has attracted much interest in the past decade, and is currently a strong candidate for high performance compact heat sinks, such as those required in electronics systems, automobile air conditioning units, micro-reactors, fuel cells, etc. Currently the literature presents numerous experimental studies on two-phase heat transfer in micro-channels, providing an extensive database that covers many different fluids and operating conditions. Among the noteworthy elements that have been reported in previous studies, is the sensitivity of micro-channel evaporators to oscillatory two-phase instabilities. These periodic fluctuations in flow and pressure drop either result from the presence of upstream compressibility, or are simply due to the interaction among parallel channels in multi-port systems. An oscillating flow presents singular characteristics that are expected to produce an effect on the local heat transfer mechanisms, and thus on the estimation of the two-phase heat transfer coefficients. The present investigation illustrates results for flow boiling of refrigerants R-134a, R-236fa, and R-245fa in a 510 {mu}m circular micro-channel, exposed to various degrees of oscillatory compressible volume instabilities. The data describe the main features of the fluctuations in the temperatures of the heated wall and fluid, and draw attention to the differences in the measured unstable time-averaged heat transfer coefficients with respect to those for stable flow boiling. (author)

Study of the distributions of flow rate and enthalpy in the sub-channels of a bundle geometry of nuclear reactors in one and two-phase flow

International Nuclear Information System (INIS)

Bayoumi, M.A.A.

1976-10-01

A bibliographic study shows that the experimental studies examined, have been developed to understand the phenomenon acting on the mixing between the sub-channels of which geometries are such these of rod bundles used in some nuclear reactors. Experimental devices and tests have been developed to study the influence of the following parameters, operating conditions, pressure, flow rate, power brought to the bundle and inlet temperature on the distribution of flow rates and vapor content among the different sub-channels. By means of non isokinetic sampling, one has determined the enthalpy of the fluid participating to the mixing between the communicating sub-channels and it has been shown that the value of this enthalpy depends strongly on the type of fluid flow and that this enthalpy cannot be either the enthalpy of one of the two sub-channels, nor (always) an average of these two enthalpies. The experimental results have been compared with calculations developed with the code FLICA, concerning the mass velocity distribution, the exchange term of linear momentum, and the variation of the transversal enthalpy with regard to the type of fluid flow. A study of local void ratio measurement, by means of optical probes, has been proposed. The present study has been carried out with a smooth geometry [fr

Numerical studies of unsteady coherent structures and transport in two-dimensional flows

Energy Technology Data Exchange (ETDEWEB)

Hesthaven, J.S.

1995-08-01

The dynamics of unsteady two-dimensional coherent structures in various physical systems is studied through direct numerical solution of the dynamical equations using spectral methods. The relation between the Eulerian and the Lagrangian auto-correlation functions in two-dimensional homogeneous, isotropic turbulence is studied. A simple analytic expression for the Eulerian and Lagrangian auto-correlation function for the fluctuating velocity field is derived solely on the basis of the one-dimensional power spectrum. The long-time evolution of monopolar and dipolar vortices in anisotropic systems relevant for geophysics and plasma physics is studied by direct numerical solution. Transport properties and spatial reorganization of vortical structures are found to depend strongly on the initial conditions. Special attention is given to the dynamics of strong monopoles and the development of unsteady tripolar structures. The development of coherent structures in fluid flows, incompressible as well as compressible, is studied by novel numerical schemes. The emphasis is on the development of spectral methods sufficiently advanced as to allow for detailed and accurate studies of the self-organizing processes. (au) 1 ill., 94 refs.

Drag reduction statistics in a channel flow

Energy Technology Data Exchange (ETDEWEB)

Jimenez-Bernal, Jose A. [Instituto Politecnico Nacional, LABINTHAP-SEPI-ESIME, Edif. 5, 3er piso Col. Lindavista, Mexico DF 07738 (Mexico); Hassan, Yassin A.; Gutierrez-Torres, Claudia del C. [Nuclear Engineering Department, Texas A and M University, College Station, TX 77843-3123 (United States)

2005-07-01

Full text of publication follows: Methods to reduce the drag have been studied for many years because of the promising payoffs that can be attained. In this investigation, the evaluation of statistics such as skewness, flatness, spectra of the stream-wise velocity fluctuations is performed for single phase flow and for two phase flow. Micro-bubbles with an average diameter of 30 {mu}m and a local void fraction of 4.8 % were produced by electrolysis and injected inside the boundary layer. This value of void fraction produced a 38.45 % decrease of the drag. The experiments were conducted in a channel flow at a Reynolds number Re 5128 (considering half height of the channel, the bulk velocity and the kinematics viscosity of the water). The channel was made of acrylic due to the optical properties of this material; its dimensions are 3.85 m long, 0.206 m wide and 0.056 m high. A pressure transducer that ranges from 0 to 35 Pa is located in the test station to measure the pressure drop in single phase flow; this pressure value is used to calculate the shear wall stress. The shear wall stress of two phase flow was measured from the velocity fields obtained from Particle Image Velocimetry (PIV) technique. PIV was utilized to measure instantaneous velocity fields in the stream-wise-normal (x-y) plane. The use of low-local values of void fraction caused a reduction of undesirable speckles effects and an absence of extreme brightness provoked by high bubble saturation. The measurements were carried out in the upper wall of the channel at 3.15 m downstream the inlet's channel. The PIV system is formed by a CCD camera with a resolution of 1008 x 1018 pixels and a double pulse laser with a maximum power 400 mJ and a wavelength of 532 nm (green light). The laser beam was transformed into a sheet of light by an array of cylindrical lenses. Two hundred frames with an area of 1.28 cm{sup 2} were recorded to obtain one hundred velocity fields. The time separation between

Drag reduction statistics in a channel flow

International Nuclear Information System (INIS)

Jimenez-Bernal, Jose A.; Hassan, Yassin A.; Gutierrez-Torres, Claudia del C.

2005-01-01

Full text of publication follows: Methods to reduce the drag have been studied for many years because of the promising payoffs that can be attained. In this investigation, the evaluation of statistics such as skewness, flatness, spectra of the stream-wise velocity fluctuations is performed for single phase flow and for two phase flow. Micro-bubbles with an average diameter of 30 μm and a local void fraction of 4.8 % were produced by electrolysis and injected inside the boundary layer. This value of void fraction produced a 38.45 % decrease of the drag. The experiments were conducted in a channel flow at a Reynolds number Re 5128 (considering half height of the channel, the bulk velocity and the kinematics viscosity of the water). The channel was made of acrylic due to the optical properties of this material; its dimensions are 3.85 m long, 0.206 m wide and 0.056 m high. A pressure transducer that ranges from 0 to 35 Pa is located in the test station to measure the pressure drop in single phase flow; this pressure value is used to calculate the shear wall stress. The shear wall stress of two phase flow was measured from the velocity fields obtained from Particle Image Velocimetry (PIV) technique. PIV was utilized to measure instantaneous velocity fields in the stream-wise-normal (x-y) plane. The use of low-local values of void fraction caused a reduction of undesirable speckles effects and an absence of extreme brightness provoked by high bubble saturation. The measurements were carried out in the upper wall of the channel at 3.15 m downstream the inlet's channel. The PIV system is formed by a CCD camera with a resolution of 1008 x 1018 pixels and a double pulse laser with a maximum power 400 mJ and a wavelength of 532 nm (green light). The laser beam was transformed into a sheet of light by an array of cylindrical lenses. Two hundred frames with an area of 1.28 cm 2 were recorded to obtain one hundred velocity fields. The time separation between consecutive pulses

Streamline topologies near simple degenerate critical points in two-dimensional flow away from boundaries

DEFF Research Database (Denmark)

Brøns, Morten; Hartnack, Johan Nicolai

1998-01-01

Streamline patterns and their bifurcations in two-dimensional incompressible flow are investigated from a topological point of view. The velocity field is expanded at a point in the fluid, and the expansion coefficients are considered as bifurcation parameters. A series of non-linear coordinate c...

Streamline topologies near simple degenerate critical points in two-dimensional flow away from boundaries

DEFF Research Database (Denmark)

Brøns, Morten; Hartnack, Johan Nicolai

1999-01-01

Streamline patterns and their bifurcations in two-dimensional incompressible flow are investigated from a topological point of view. The velocity field is expanded at a point in the fluid, and the expansion coefficients are considered as bifurcation parameters. A series of nonlinear coordinate ch...

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

International Nuclear Information System (INIS)

Rasouli, F.

1981-01-01

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

Flow and Jamming of Granular Materials in a Two-dimensional Hopper

Science.gov (United States)

Tang, Junyao

Flow in a hopper is both a fertile testing ground for understanding fundamental granular flow rheology and industrially highly relevant. Despite increasing research efforts in this area, a comprehensive physical theory is still lacking for both jamming and flow of granular materials in a hopper. In this work, I have designed a two dimensional (2D) hopper experiment using photoelastic particles (particles' shape: disk or ellipse), with the goal to build a bridge between macroscopic phenomenon of hopper flow and microscopic particle-scale dynamics. Through synchronized data of particle tracking and stress distributions in particles, I have shown differences between my data of the time-averaged velocity/stress profile of 2D hopper flow with previous theoretical predictions. I have also demonstrated the importance of a mechanical stable arch near the opening on controlling hopper flow rheology and suggested a heuristic phase diagram for the hopper flow/jamming transition. Another part of this thesis work is focused on studying the impact of particle shape of particles on hopper flow. By comparing particle-tracking and photoelastic data for ellipses and disks at the appropriate length scale, I have demonstrated an important role for the rotational freedom of elliptical particles in controlling flow rheology through particle tracking and stress analysis. This work has been supported by International Fine Particle Research Institute (IFPRI) .

Three-dimensional visualization of material flow during friction stir welding by two pairs of X-ray transmission systems

International Nuclear Information System (INIS)

Morisada, Y.; Fujii, H.; Kawahito, Y.; Nakata, K.; Tanaka, M.

2011-01-01

Material flow during friction stir welding is crucial to obtaining sound joints. However, this phenomenon is still not fully understood despite many investigations and numerous models. In this study, the material flow is three-dimensionally visualized by X-ray radiography using a tiny spherical tungsten tracer. The movement of the tracer during the friction stir welding is observed by two pairs of X-ray transmission real-time imaging systems. The three-dimensional material flow is obtained by following the locus of the tracer.

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

International Nuclear Information System (INIS)

Srinivasan, J.

1978-05-01

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

Numerical solution of viscous and viscoelastic fluids flow through the branching channel by finite volume scheme

Science.gov (United States)

Keslerová, Radka; Trdlička, David

2015-09-01

This work deals with the numerical modelling of steady flows of incompressible viscous and viscoelastic fluids through the three dimensional channel with T-junction. The fundamental system of equations is the system of generalized Navier-Stokes equations for incompressible fluids. This system is based on the system of balance laws of mass and momentum for incompressible fluids. Two different mathematical models for the stress tensor are used for simulation of Newtonian and Oldroyd-B fluids flow. Numerical solution of the described models is based on cetral finite volume method using explicit Runge-Kutta time integration.

Experimental study on two-phase flow natural circulation in a core catcher cooling channel for EU-APR1400 using air-water system

Energy Technology Data Exchange (ETDEWEB)

Song, Ki Won [Division of Advanced Nuclear Engineering, POSTECH, Pohang 790-784 (Korea, Republic of); Korea Atomic Energy Research Institute, Daejeon 34057 (Korea, Republic of); Nguyen, Thanh Hung [School of Nuclear Engineering, Purdue University, West Lafayette, IN 47906 (United States); Ha, Kwang Soon; Kim, Hwan Yeol; Song, Jinho [Korea Atomic Energy Research Institute, Daejeon 34057 (Korea, Republic of); Park, Hyun Sun [Division of Advanced Nuclear Engineering, POSTECH, Pohang 790-784 (Korea, Republic of); Revankar, Shripad T., E-mail: shripad@postech.ac.kr [Division of Advanced Nuclear Engineering, POSTECH, Pohang 790-784 (Korea, Republic of); School of Nuclear Engineering, Purdue University, West Lafayette, IN 47906 (United States); Kim, Moo Hwan [Division of Advanced Nuclear Engineering, POSTECH, Pohang 790-784 (Korea, Republic of); Korea Institute of Nuclear Safety, Daejeon 305-338 (Korea, Republic of)

2017-05-15

Highlights: • Two-phase flow regimes and transition behavior were observed in the coolant channel. • Test were conducted for natural circulation with air-water. • Data were obtained on flow regime, void fraction, flow rates and re-wetting time. • The data were related to a cooling capability of core catcher system. - Abstract: Ex-vessel core catcher cooling system driven by natural circulation is designed using a full scaled air-water system. A transparent half symmetric section of a core catcher coolant channel of a pressurized water reactor was designed with instrumentations for local void fraction measurement and flow visualization. Two designs of air-water top separator water tanks are studied including one with modified ‘super-step’ design which prevents gas entrainment into down-comer. In the experiment air flow rates are set corresponding to steam generation rate for given corium decay power. Measurements of natural circulation flow rate, spatial local void fraction distribution and re-wetting time near the top wall are carried out for various air flow rates which simulate boiling-induced vapor generation. Since heat transfer and critical heat flux are strongly dependent on the water mass flow rate and development of two-phase flow on the heated wall, knowledge of two-phase flow characteristics in the coolant channel is essential. Results on flow visualization showing two phase flow structure specifically near the high void accumulation regions, local void profiles, rewetting time, and natural circulation flow rate are presented for various air flow rates that simulate corium power levels. The data are useful in assessing the cooling capability of and safety of the core catcher system.

Numerical study of the bubbly flow regime in micro-channel flow boiling

Science.gov (United States)

Bhuvankar, Pramod; Dabiri, Sadegh

2017-11-01

Two-phase flow accompanied by boiling in micro-channel heat sinks is an effective means for heat removal from computer chips. We present a numerical study of flow boiling in micro-channels with conjugate heat transfer with a focus on the bubbly flow regime. The bubbles are assumed to nucleate at a pre-determined location and frequency. The Navier Stokes equations are solved using a single fluid formulation with the Front tracking method. Phase change is implemented using the deficit in heat flux across the bubble interface. The analytical solution for bubble growth in a superheated liquid is used as a benchmark to validate the mentioned numerical method. Water and FC-72 are studied as the operating fluids in a micro-channel made of Copper with a focus on hotspot mitigation. The micro-channel of cross-section 231 μm × 1000 μm , is used to study the effects of vertical up-flow, vertical down-flow and horizontal flow of the mentioned fluids on the heat transfer coefficients. A simple film model accounting for mass and energy conservation is applied wherever the bubble approaches closer than a cell width to the wall. The results of the simulation are compared with existing experimental data for bubble growth rates and heat transfer coefficients.

Experimental visualization of temperature fields and study of heat transfer enhancement in oscillatory flow in a grooved channel

Energy Technology Data Exchange (ETDEWEB)

Herman, C.; Kang, E. [Dept. of Mechanical Engineering, Johns Hopkins Univ., Baltimore, MD (United States)

2001-01-01

An experimental study was conducted of incompressible, moderate Reynolds number flow of air over heated rectangular blocks in a two-dimensional, horizontal channel. Holographic interferometry combined with high-speed cinematography was used to visualize the unsteady temperature fields in self- sustained oscillatory flow. Experiments were conducted in the laminar, transitional and turbulent flow regimes for Reynolds numbers in the range from Re = 520 to Re = 6600. Interferometric measurements were obtained in the thermally and fluiddynamically periodically fully developed flow region on the ninth heated block. Flow oscillations were first observed between Re = 1054 and Re = 1318. The period of oscillations, wavelength and propagation speed of the Tollmien-Schlichting waves in the main channel were measured at two characteristic flow velocities, Re = 1580 and Re = 2370. For these Reynolds numbers it was observed that two to three waves span one geometric periodicity length. At Re = 1580 the dominant oscillation frequency was found to be around 26 Hz and at Re = 2370 the frequency distribution formed a band around 125 Hz. Results regarding heat transfer and pressure drop are presented as a function of the Reynolds number, in terms of the block-average Nusselt number and the local Nusselt number as well as the friction factor. Measurements of the local Nusselt number together with visual observations indicate that the lateral mixing caused by flow instabilities is most pronounced along the upstream vertical wall of the heated block in the groove region, and it is accompanied by high heat transfer coefficients. At Reynolds numbers beyond the onset of oscillations the heat transfer in the grooved channel exceeds the performance of the reference geometry, the asymmetrically heated parallel plate channel. (orig.)

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

International Nuclear Information System (INIS)

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

2002-01-01

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

Unsteady two-dimensional potential-flow model for thin variable geometry airfoils

DEFF Research Database (Denmark)

Gaunaa, Mac

2010-01-01

In the present work, analytical expressions for distributed and integral unsteady two-dimensional forces on a variable geometry airfoil undergoing arbitrary motion are derived under the assumption of incompressible, irrotational, inviscid flow. The airfoil is represented by its camber line...... in their equivalent state-space form, allowing for use of the present theory in problems employing the eigenvalue approach, such as stability analysis. The analytical expressions for the integral forces can be reduced to Munk's steady and Theodorsen's unsteady results for thin airfoils, and numerical evaluation shows...

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

International Nuclear Information System (INIS)

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

1997-01-01

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

Using artificial intelligence to improve identification of nanofluid gas–liquid two-phase flow pattern in mini-channel

Directory of Open Access Journals (Sweden)

Jian Xiao

2018-01-01

Full Text Available This work combines fuzzy logic and a support vector machine (SVM with a principal component analysis (PCA to create an artificial-intelligence system that identifies nanofluid gas-liquid two-phase flow states in a vertical mini-channel. Flow-pattern recognition requires finding the operational details of the process and doing computer simulations and image processing can be used to automate the description of flow patterns in nanofluid gas-liquid two-phase flow. This work uses fuzzy logic and a SVM with PCA to improve the accuracy with which the flow pattern of a nanofluid gas-liquid two-phase flow is identified. To acquire images of nanofluid gas-liquid two-phase flow patterns of flow boiling, a high-speed digital camera was used to record four different types of flow-pattern images, namely annular flow, bubbly flow, churn flow, and slug flow. The textural features extracted by processing the images of nanofluid gas–liquid two-phase flow patterns are used as inputs to various identification schemes such as fuzzy logic, SVM, and SVM with PCA to identify the type of flow pattern. The results indicate that the SVM with reduced characteristics of PCA provides the best identification accuracy and requires less calculation time than the other two schemes. The data reported herein should be very useful for the design and operation of industrial applications.

Using artificial intelligence to improve identification of nanofluid gas-liquid two-phase flow pattern in mini-channel

Science.gov (United States)

Xiao, Jian; Luo, Xiaoping; Feng, Zhenfei; Zhang, Jinxin

2018-01-01

This work combines fuzzy logic and a support vector machine (SVM) with a principal component analysis (PCA) to create an artificial-intelligence system that identifies nanofluid gas-liquid two-phase flow states in a vertical mini-channel. Flow-pattern recognition requires finding the operational details of the process and doing computer simulations and image processing can be used to automate the description of flow patterns in nanofluid gas-liquid two-phase flow. This work uses fuzzy logic and a SVM with PCA to improve the accuracy with which the flow pattern of a nanofluid gas-liquid two-phase flow is identified. To acquire images of nanofluid gas-liquid two-phase flow patterns of flow boiling, a high-speed digital camera was used to record four different types of flow-pattern images, namely annular flow, bubbly flow, churn flow, and slug flow. The textural features extracted by processing the images of nanofluid gas-liquid two-phase flow patterns are used as inputs to various identification schemes such as fuzzy logic, SVM, and SVM with PCA to identify the type of flow pattern. The results indicate that the SVM with reduced characteristics of PCA provides the best identification accuracy and requires less calculation time than the other two schemes. The data reported herein should be very useful for the design and operation of industrial applications.

Self-oscillations of a two-dimensional shear flow with forcing and dissipation

Science.gov (United States)

López Zazueta, A.; Zavala Sansón, L.

2018-04-01

Two-dimensional shear flows continuously forced in the presence of dissipative effects are studied by means of numerical simulations. In contrast with most previous studies, the forcing is confined in a finite region, so the behavior of the system is characterized by the long-term evolution of the global kinetic energy. We consider regimes with 1 limited to develop only one vortical instability by choosing an appropriate width of the forcing band. The most relevant regime is found for Reλ > 36, in which the energy maintains a regular oscillation around a reference value. The flow configuration is an elliptical vortex tilted with respect to the forcing axis, which oscillates steadily also. Second, the flow is allowed to develop two Kelvin-Helmholtz billows and eventually more complicated structures. The regimes of the one-vortex case are observed again, except for Reλ > 135. At these values, the energy oscillates chaotically as the two vortices merge, form dipolar structures, and split again, with irregular periodicity. The self-oscillations are explained as a result of the alternate competition between forcing and dissipation, which is verified by calculating the budget terms in the energy equation. The relevance of the forcing-vs.-dissipation competition is discussed for more general flow systems.

STRUYA a code for two-dimensional fluid flow analysis with and without structure coupling

International Nuclear Information System (INIS)

Katz, F.W.; Schlechtendahl, E.G.; Stoelting, K.

1979-11-01

STRUYA is a code for two-dimensional subsonic and supersonic flow analysis. Both Eulerian and Lagrangian grids are allowed. In the third dimension the flow domain may be bounded by a moving wall. The wall movement may be prescribed in a time-and space varying way or computed by a structural model. STRUYA offers a general scheme for adapting various structural models. As a standard feature it includes a cylindrical shell model (CYLDY2). (orig.) [de

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

International Nuclear Information System (INIS)

Moesinger, H.

1979-08-01

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

LLUVIA-II: A program for two-dimensional, transient flow through partially saturated porous media

International Nuclear Information System (INIS)

Eaton, R.R.; Hopkins, P.L.

1992-08-01

LLUVIA-II is a program designed for the efficient solution of two- dimensional transient flow of liquid water through partially saturated, porous media. The code solves Richards equation using the method-of-lines procedure. This document describes the solution procedure employed, input data structure, output, and code verification

High-frame rate imaging of two-phase flow in a thin rectangular channel using fast neutrons.

Science.gov (United States)

Zboray, R; Mor, I; Dangendorf, V; Stark, M; Tittelmeier, K; Cortesi, M; Adams, R

2014-08-01

We have demonstrated the feasibility of performing high-frame-rate, fast neutron radiography of air-water two-phase flows in a thin channel with rectangular cross section. The experiments have been carried out at the accelerator facility of the Physikalisch-Technische Bundesanstalt. A polychromatic, high-intensity fast neutron beam with average energy of 6 MeV was produced by 11.5 MeV deuterons hitting a thick Be target. Image sequences down to 10 ms exposure times were obtained using a fast-neutron imaging detector developed in the context of fast-neutron resonance imaging. Different two-phase flow regimes such as bubbly slug and churn flows have been examined. Two phase flow parameters like the volumetric gas fraction, bubble size and mean bubble velocities have been measured. The first results are promising, improvements for future experiments are also discussed. Copyright © 2014 Elsevier Ltd. All rights reserved.

Time evolution of the eddy viscosity in two-dimensional navier-stokes flow

Science.gov (United States)

Chaves; Gama

2000-02-01

The time evolution of the eddy viscosity associated with an unforced two-dimensional incompressible Navier-Stokes flow is analyzed by direct numerical simulation. The initial condition is such that the eddy viscosity is isotropic and negative. It is shown by concrete examples that the Navier-Stokes dynamics stabilizes negative eddy viscosity effects. In other words, this dynamics moves monotonically the initial negative eddy viscosity to positive values before relaxation due to viscous term occurs.

The Influence of Non-Equilibrium Excitation on the Electron Density in One-Dimensional MFD Channel Flow

Energy Technology Data Exchange (ETDEWEB)

Eichert, K.; Kaeppeler, H. J. [Institut fuer Plasmaforschung der Technischen Hochschule Stuttgart, Federal Republic of Germany (Germany)

1966-10-15

In previous publications, a system of equations was derived from the gas-kinetic description of a multi-component reacting plasma and employed for the calculation of one-dimensional subsonic flows. This system is now extended to include non-equilibrium excitation. No thermal or chemical equilibrium between the various components of the plasma is assumed. The components of the plasma considered are a non-reacting working fluid, an alkali metal vapour as a seeding material, ions of this seeding substance, and electrons. Three levels for the excited states are introduced. The reactions considered are excitation and ionization by electron collisions, and photo-ionization, as well as the corresponding reverse processes. For the reaction velocities, analytical equations are introduced permitting insertion of any excitation or ionization cross-sections of either experimental or theoretical origin. The method employed had been previously suggested by one of the authors. As examples, the degrees of excitation and ionization in the flow of a helium working fluid with 1% caesium seeding through a channel against transverse magnetic fields of 15 and 40 kg at Mach numbers of 0.7 and 0.8, respectively, were calculated. The results of the calculations show that for relatively small magnetic fields there is no rapid rise of the ionization to Saha-equilibrium as a function of electron temperature. A comparison with the results of a calculation neglecting excitation shows that especially for relatively large magnetic fields non-equilibrium excitation has an essential influence on the electron density and its approach to equilibrium. Neglecting excitation, there results a nearly frozen behaviour of the degree of ionization within channel lengths of technical interest for small magnetic fields. (author)

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

DEFF Research Database (Denmark)

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

2014-01-01

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

Analysis of flow distribution instability in parallel thin rectangular multi-channel system

Energy Technology Data Exchange (ETDEWEB)

Xia, G.L. [School of Nuclear Science and Technology, Xi’an Jiaotong University, Xi’an City 710049 (China); Fundamental Science on Nuclear Safety and Simulation Technology Laboratory, Harbin Engineering University, Harbin City 150001 (China); Su, G.H., E-mail: ghsu@mail.xjtu.edu.cn [School of Nuclear Science and Technology, Xi’an Jiaotong University, Xi’an City 710049 (China); Peng, M.J. [Fundamental Science on Nuclear Safety and Simulation Technology Laboratory, Harbin Engineering University, Harbin City 150001 (China)

2016-08-15

Highlights: • Flow distribution instability in parallel thin rectangular multi-channel system is studied using RELAP5 codes. • Flow excursion may bring parallel heating channel into the density wave oscillations region. • Flow distribution instability is more likely to happen at low power/flow ratio conditions. • The increase of channel number will not affect the flow distribution instability boundary. • Asymmetry inlet throttling and heating will make system more unstable. - Abstract: The flow distribution instability in parallel thin rectangular multi-channel system has been researched in the present study. The research model of parallel channel system is established by using RELAP5/MOD3.4 codes. The transient process of flow distribution instability is studied at imposed inlet mass flow rate and imposed pressure drop conditions. The influence of heating power, mass flow rate, system pressure and channel number on flow distribution instability are analyzed. Furthermore, the flow distribution instability of parallel two-channel system under asymmetric inlet throttling and heating power is studied. The results show that, if multi-channel system operates at the negative slope region of channel ΔP–G curve, small disturbance in pressure drop will lead to flow redistribution between parallel channels. Flow excursion may bring the operating point of heating channel into the density-wave oscillations region, this will result in out-phase or in-phase flow oscillations. Flow distribution instability is more likely to happen at low power/flow ratio conditions, the stability of parallel channel system increases with system pressure, the channel number has a little effect on system stability, but the asymmetry inlet throttling or heating power will make the system more unstable.

Flow study in the formatted channel for two disks in rotation; Estudo do escoamento no canal formado por dois discos em rotacao

Energy Technology Data Exchange (ETDEWEB)

Barbosa, Marcos Pinotti

1992-07-01

Flow study in the formatted channel for two disks in rotation is discussed including the following main issues: flow description between disks in rotation; computational model; and numerical results. The parametric studies accomplished of the spacing between disks, of the diameter, and of the rotor angular speed allowed the influence analysis of these variables in the flow inside the channel and in the generated pressure difference. The disks rotation, evaluated through Reynolds' rotational number, is the main parameter that influences the pressure difference between channel entrance and exit. It verified although how much larger the rotation larger the pressure difference.

Flow study in the formatted channel for two disks in rotation; Estudo do escoamento no canal formado por dois discos em rotacao

Energy Technology Data Exchange (ETDEWEB)

Barbosa, Marcos Pinotti

1992-07-01

Flow study in the formatted channel for two disks in rotation is discussed including the following main issues: flow description between disks in rotation; computational model; and numerical results. The parametric studies accomplished of the spacing between disks, of the diameter, and of the rotor angular speed allowed the influence analysis of these variables in the flow inside the channel and in the generated pressure difference. The disks rotation, evaluated through Reynolds' rotational number, is the main parameter that influences the pressure difference between channel entrance and exit. It verified although how much larger the rotation larger the pressure difference.

Generalisation of two-layer turbulent model for passive cooling in a channel

International Nuclear Information System (INIS)

Bennacer, R.; Hammami, T.; Mohamad, A.A.; Beji, H.

2003-01-01

Turbulent natural convection still under improvement and no perfect compromise exist. The near wall region modelisation poses numerical difficulties and current modeling are either expensive or lack universality. Uncertainness in evaluating the good heat transfer rate can be catastrophically in causing local overheat and materials destruction which can be of heavy consequence as cooling nuclear component (rodes). Using the recent DNS done on natural convection flow in an infinite channel differentially heated for (10 4 6 ) a scaling analysis is developed and a one-equation near-wall turbulence model is deduced (inner layer). The inner model is coupled with a Low Reynolds Model (LRM) in the outer region (second layer) and applied to calculate natural flow for different Ra numbers. It yields good performance, computation time reduction and much better heat transfer prediction compared to the diffusive Jones Launder LRM. The efficiency is tested in one-dimensional and two-dimensional case. (author)

NUMERICAL INVESTIGATION OF CURVATURE AND TORSION EFFECTS ON WATER FLOW FIELD IN HELICAL RECTANGULAR CHANNELS

Directory of Open Access Journals (Sweden)

A. H. ELBATRAN

2015-07-01

Full Text Available Helical channels have a wide range of applications in petroleum engineering, nuclear, heat exchanger, chemical, mineral and polymer industries. They are used in the separation processes for fluids of different densities. The centrifugal force, free surface and geometrical effects of the helical channel make the flow pattern more complicated; hence it is very difficult to perform physical experiment to predict channel performance. Computational Fluid Dynamics (CFD can be suitable alternative for studying the flow pattern characteristics in helical channels. The different ranges of dimensional parameters, such as curvature and torsion, often cause various flow regimes in the helical channels. In this study, the effects of physical parameters such as curvature, torsion, Reynolds number, Froude number and Dean Number on the characteristics of the turbulent flow in helical rectangular channels have been investigated numerically, using a finite volume RANSE code Fluent of Ansys workbench 10.1 UTM licensed. The physical parameters were reported for range of curvature (δ of 0.16 to 0.51 and torsion (λ of 0.032 to 0.1 .The numerical results of this study showed that the decrease in the channel curvature and the increase in the channel torsion numbers led to the increase of the flow velocity inside the channel and the change in the shape of water free surface at given Dean, Reynolds and Froude numbers.

Influential factors on debris flow events and hillslope-channel connectivity in Alpine regions: case studies from two Alpine regions in Styria, Austria

Science.gov (United States)

Traper, Sandra; Pöppl, Ronald; Rascher, Eric; Sass, Oliver

2016-04-01

In recent times different types of natural disasters like debris flow events have attracted increasing attention worldwide, since they can cause great damage and loss of infrastructure or even lives is not unusual when it comes to such an event. The engagement with debris flows is especially important in mountainous areas like Austria, since Alpine regions have proved to be particularly prone to the often harmful consequences of such events because of increasing settlement of previously uninhabited regions. Due to those frequently damaging effects of debris flows, research on this kind of natural disaster often focuses on mitigation and recovery measures after an event and on how to restore the initial situation. However, a view on the situation of an area, where severe debris flows recently occurred and are well documented, before the actual event can aid in discovering important preparatory factors that contribute to initiating debris flows and hillslope-channel connectivity in the first place. Valuable insights into the functioning and preconditions of debris flows and their potential connectivity to the main channel can be gained. The study focuses on two geologically different areas in the Austrian Alps, which are both prone to debris flows and have experienced rather severe events recently. Based on data from debris flow events in two regions in Styria (Austria), the Kleinsölk and the Johnsbach valleys, the aim of the study is to identify factors which influence the development of debris flows and the potential of such debris flows to reach the main channel potentially clogging up the river (hillslope-channel connectivity). The degree of hillslope-channel coupling was verified in extensive TLS and ALS surveys, resulting in DEMs of different resolution and spatial extension. Those factors are obtained, analyzed and evaluated with DEM-based GIS- and statistical analyses. These include factors that are attributed to catchment topography, such as slope angle

Numerical simulation in a two dimensional turbulent flow over a backward-facing step

International Nuclear Information System (INIS)

Silveira Neto, A. da; Grand, D.

1991-01-01

Numerical simulations of turbulent flows in complex geometries are generally restricted to the prediction of the mean flow and use semi-empirical turbulence models. The present study is devoted to the simulation of the coherence structures which develop in a flow submitted to a velocity change, downstream of a backward facing step. Two aspect ratios (height of the step over height of the channel) have been explored and the values of the Reynolds number vary from (6000 to 90000). In the isothermal case coherent structures have been obtained by the numerical simulation in the mixing layer downstream of the step. The numerical simulations provides results in fairly good agreement with available experimental results. In a second step a thermal stratification is imposed on this flow for one value of Richardson number (0.5) the coherent structures disappear downstream for increasing values of Richardson number. (author)

Numerical simulations of carbon monoxide poisoning in high temperature proton exchange membrane fuel cells with various flow channel designs

International Nuclear Information System (INIS)

Jiao, Kui; Zhou, Yibo; Du, Qing; Yin, Yan; Yu, Shuhai; Li, Xianguo

2013-01-01

Highlights: ► Simulations of CO poisoning in HT-PEMFC with different flow channels are conducted. ► Parallel and serpentine designs result in least and most CO effects, respectively. ► General CO distributions in CLs are similar with different flow channel designs. - Abstract: The performance of high temperature proton exchange membrane fuel cell (HT-PEMFC) is significantly affected by the carbon monoxide (CO) in hydrogen fuel, and the flow channel design may influence the CO poisoning characteristics by changing the reactant flow. In this study, three-dimensional non-isothermal simulations are carried out to investigate the comprehensive flow channel design and CO poisoning effects on the performance of HT-PEMFCs. The numerical results show that when pure hydrogen is supplied, the interdigitated design produces the highest power output, the power output with serpentine design is higher than the two parallel designs, and the parallel-Z and parallel-U designs have similar power outputs. The performance degradation caused by CO poisoning is the least significant with parallel flow channel design, but the most significant with serpentine and interdigitated designs because the cross flow through the electrode is stronger. At low cell voltages (high current densities), the highest power outputs are with interdigitated and parallel flow channel designs at low and high CO fractions in the supplied hydrogen, respectively. The general distributions of absorbed hydrogen and CO coverage fractions in anode catalyst layer (CL) are similar for the different flow channel designs. The hydrogen coverage fraction is higher under the channel than under the land, and is also higher on the gas diffusion layer (GDL) side than on the membrane side; and the CO coverage distribution is opposite to the hydrogen coverage distribution

Chemical control of channel interference in two-photon absorption processes.

Science.gov (United States)

Alam, Md Mehboob; Chattopadhyaya, Mausumi; Chakrabarti, Swapan; Ruud, Kenneth

2014-05-20

The two-photon absorption (TPA) process is the simplest and hence the most studied nonlinear optical phenomenon, and various aspects of this process have been explored in the past few decades, experimentally as well as theoretically. Previous investigations have shown that the two-photon (TP) activity of a molecular system can be tuned, and at present, performance-tailored TP active materials are easy to develop by monitoring factors such as length of conjugation, dimensionality of charge-transfer network, strength of donor-acceptor groups, polarity of solvents, self-aggregation, H-bonding, and micellar encapsulation to mention but a few. One of the most intriguing phenomena affecting the TP activity of a molecule is channel interference. The phrase "channel interference" implies that if the TP transition from one electronic state to another involves more than one optical pathway or channel, characterized by the corresponding transition dipole moment (TDM) vectors, the channels may interfere with each other depending upon the angles between the TDM vectors and hence can either increase (constructive interference) or decrease (destructive interference) the overall TP activity of a system to a significant extent. This phenomenon was first pointed out by Cronstrand, Luo, and Ågren [Chem. Phys. Lett. 2002, 352, 262-269] in two-dimensional systems (i.e., only involving two components of the transition moment vectors). For three-dimensional molecules, an extended version of this idea was required. In order to fill this gap, we developed a generalized model for describing and exploring channel interference, valid for systems of any dimensionality. We have in particular applied it to through-bond (TB) and through-space (TS) charge-transfer systems both in gas phase and in solvents with different polarities. In this Account, we will, in addition to briefly describing the concept of channel interference, discuss two key findings of our recent work: (1) how to control the

Sub-Nanometer Channels Embedded in Two-Dimensional Materials

KAUST Repository

Han, Yimo; Li, Ming-yang; Jung, Gang-Seob; Marsalis, Mark A.; Qin, Zhao; Buehler, Markus J.; Li, Lain-Jong; Muller, David A.

2017-01-01

Two-dimensional (2D) materials are among the most promising candidates for next-generation electronics due to their atomic thinness, allowing for flexible transparent electronics and ultimate length scaling1. Thus far, atomically-thin p-n junctions2

Strongly coupled fluid-particle flows in vertical channels. I. Reynolds-averaged two-phase turbulence statistics

International Nuclear Information System (INIS)

Capecelatro, Jesse; Desjardins, Olivier; Fox, Rodney O.

2016-01-01

Simulations of strongly coupled (i.e., high-mass-loading) fluid-particle flows in vertical channels are performed with the purpose of understanding the fundamental physics of wall-bounded multiphase turbulence. The exact Reynolds-averaged (RA) equations for high-mass-loading suspensions are presented, and the unclosed terms that are retained in the context of fully developed channel flow are evaluated in an Eulerian–Lagrangian (EL) framework for the first time. A key distinction between the RA formulation presented in the current work and previous derivations of multiphase turbulence models is the partitioning of the particle velocity fluctuations into spatially correlated and uncorrelated components, used to define the components of the particle-phase turbulent kinetic energy (TKE) and granular temperature, respectively. The adaptive spatial filtering technique developed in our previous work for homogeneous flows [J. Capecelatro, O. Desjardins, and R. O. Fox, “Numerical study of collisional particle dynamics in cluster-induced turbulence,” J. Fluid Mech. 747, R2 (2014)] is shown to accurately partition the particle velocity fluctuations at all distances from the wall. Strong segregation in the components of granular energy is observed, with the largest values of particle-phase TKE associated with clusters falling near the channel wall, while maximum granular temperature is observed at the center of the channel. The anisotropy of the Reynolds stresses both near the wall and far away is found to be a crucial component for understanding the distribution of the particle-phase volume fraction. In Part II of this paper, results from the EL simulations are used to validate a multiphase Reynolds-stress turbulence model that correctly predicts the wall-normal distribution of the two-phase turbulence statistics.

Strongly coupled fluid-particle flows in vertical channels. I. Reynolds-averaged two-phase turbulence statistics

Science.gov (United States)

Capecelatro, Jesse; Desjardins, Olivier; Fox, Rodney O.

2016-03-01

Simulations of strongly coupled (i.e., high-mass-loading) fluid-particle flows in vertical channels are performed with the purpose of understanding the fundamental physics of wall-bounded multiphase turbulence. The exact Reynolds-averaged (RA) equations for high-mass-loading suspensions are presented, and the unclosed terms that are retained in the context of fully developed channel flow are evaluated in an Eulerian-Lagrangian (EL) framework for the first time. A key distinction between the RA formulation presented in the current work and previous derivations of multiphase turbulence models is the partitioning of the particle velocity fluctuations into spatially correlated and uncorrelated components, used to define the components of the particle-phase turbulent kinetic energy (TKE) and granular temperature, respectively. The adaptive spatial filtering technique developed in our previous work for homogeneous flows [J. Capecelatro, O. Desjardins, and R. O. Fox, "Numerical study of collisional particle dynamics in cluster-induced turbulence," J. Fluid Mech. 747, R2 (2014)] is shown to accurately partition the particle velocity fluctuations at all distances from the wall. Strong segregation in the components of granular energy is observed, with the largest values of particle-phase TKE associated with clusters falling near the channel wall, while maximum granular temperature is observed at the center of the channel. The anisotropy of the Reynolds stresses both near the wall and far away is found to be a crucial component for understanding the distribution of the particle-phase volume fraction. In Part II of this paper, results from the EL simulations are used to validate a multiphase Reynolds-stress turbulence model that correctly predicts the wall-normal distribution of the two-phase turbulence statistics.

Integral model of linear momentum for one-dimensional two-phase flows

International Nuclear Information System (INIS)

Kuznetsov, Yu.A.; Sabaev, E.F.

1976-01-01

''An integrated momentum model'' obtained by Meyer-Rose and widely applicable in calculations of dynamics of the thermal power systems is generalized for a case of flow of a vapour-liquid mixture with phase creep and pressure variation in the heated channel. Pressure distribution along the channel length is shown for a number of cases to be negligible. The obtained equations are found as well applicable in case pressure greatly though slowly varies in the system

Heat transfer from two-side heated helical channels

International Nuclear Information System (INIS)

Shimonis, V.; Ragaishis, V.; Poshkas, P.

1995-01-01

Experimental results are presented on the heat transfer from two-side heated helical channels to gas (air) flows. The study covered six configurations and wide ranges of geometrical (D/h=5.5 to 84.2) and performance (Re=10 3 to 2*10 5 ) parameters. Under the influence of Re and of the channel curvature, the heat transfer from both the convex and the concave surfaces for two-side heating (q w1 ≅ q w2 ) is augmented by 20-30% over one-side heating. Improved relations to predict the critical values of Reynolds Re cr1 and Re cr2 are suggested. They enable more exact predictions of the heat transfer from convex surface in transient flows for one-side heating. The relation for annular channels is suggested for the turbulent heat transfer from the convex and concave surfaces of two-side heated helical channels. It can be adapted by introducing earlier expresions for one-side heated helical channels. (author). 6 refs., 2 tabs., 3 figs

Two-dimensional model for subthreshold current and subthreshold swing of graded-channel dual-material double-gate (GCDMDG) MOSFETs

Science.gov (United States)

Goel, Ekta; Kumar, Sanjay; Singh, Balraj; Singh, Kunal; Jit, Satyabrata

2017-06-01

The subthreshold performance of graded-channel dual-material double-gate (GCDMDG) MOSFETs is examined through two-dimensional (2D) analytical modeling of subthreshold-current (SC) and subthreshold-swing (SS). The potential function obtained by using the parabolic approach to solve the 2D Poisson's equation, has been used to formulate SC and SS characteristics of the device. The variations of SS against different device parameters have been obtained with the help of effective conduction path parameter. The SC and SS characteristics of the GCDMDG MOS transistor have been compared with those of the dual-material double-gate (DMDG) and simple graded-channel double-gate (GCDG) MOS structures to show its better subthreshold characteristics over the latter two devices. The results of the developed model are well-agreed with the commercially available SILVACO ATLAS™ simulator data.

Triple-porosity/permeability flow in faulted geothermal reservoirs: Two-dimensional effects

Energy Technology Data Exchange (ETDEWEB)

Cesar Suarez Arriaga, M. [Michoacan Univ. & CFE, Mich. (Mexico); Samaniego Verduzco, F. [National Autonomous Univ. of Mexico, Coyoacan (Mexico)

1995-03-01

An essential characteristic of some fractured geothermal reservoirs is noticeable when the drilled wells intersect an open fault or macrofracture. Several evidences observed, suggest that the fluid transport into this type of systems, occurs at least in three stages: flow between rock matrix and microfractures, flow between fractures and faults and flow between faults and wells. This pattern flow could define, by analogy to the classical double-porosity model, a triple-porosity, triple-permeability concept. From a mathematical modeling point of view, the non-linearity of the heterogeneous transport processes, occurring with abrupt changes on the petrophysical properties of the rock, makes impossible their exact or analytic solution. To simulate this phenomenon, a detailed two-dimensional geometric model was developed representing the matrix-fracture-fault system. The model was solved numerically using MULKOM with a H{sub 2}O=CO{sub 2} equation of state module. This approach helps to understand some real processes involved. Results obtained from this study, exhibit the importance of considering the triple porosity/permeability concept as a dominant mechanism producing, for example, strong pressure gradients between the reservoir and the bottom hole of some wells.

Study on characteristics of vapor-liquid two-phase flow in mini-channels

International Nuclear Information System (INIS)

Guo Lei; Zhang Shusheng; Cheng Lin

2011-01-01

Highlights: → I- and Z-shaped mini-channels are studied with water and ethanol respectively. → The smaller the wall contact angles are, the less time the bubbles take to depart. → Bubble bottom micro-layer can enhance heat transfer. → Z-shaped channels show a higher heat transfer coefficient but a larger pressure drop. → Water reflects a higher heat transfer coefficient than ethanol in the same channels. - Abstract: To explore the mechanism of boiling bubble dynamics in narrow channels, two types of channels are investigated which have I- and Z-shaped with width of 2 mm. Using VOF model and self-programming, the whole flow field is simulated with two different kinds of media, namely, water and ethanol. The influence of wall contact angle on the process of bubble generating and growth is studied, and the relationship between different channel shapes and the pressure drop is also investigated taking into account the effects of gravity, viscosity, surface tension and wall adhesion. The bubble generation, growth and departure processes are analyzed through numerical simulation and self-programming, and the influence of interface movements and changes on internal pressure difference and average surface heat transfer coefficient is investigated by using geometry reconstruction and interface tracking. It is found that wall contact angle has a great influence on the morphology of bubble. The smaller the wall contact angles are, the more round the bubbles are, and the less time the bubbles take to depart from the wall. The variation of contact angle also has effect upon the heat transfer coefficient. The greater the wall contact angle is, the larger the bubble-covered area is, thus the wall thermal resistance gets higher, and bubble nucleation is suppressed, and the heat transfer coefficient becomes lower. The role of surface tension in the process of boiling heat transfer is much more important than the gravity in narrow channels. The generation of bubbles

Spatio-temporal organization of dynamics in a two-dimensional periodically driven vortex flow: A Lagrangian flow network perspective.

Science.gov (United States)

Lindner, Michael; Donner, Reik V

2017-03-01

We study the Lagrangian dynamics of passive tracers in a simple model of a driven two-dimensional vortex resembling real-world geophysical flow patterns. Using a discrete approximation of the system's transfer operator, we construct a directed network that describes the exchange of mass between distinct regions of the flow domain. By studying different measures characterizing flow network connectivity at different time-scales, we are able to identify the location of dynamically invariant structures and regions of maximum dispersion. Specifically, our approach allows us to delimit co-existing flow regimes with different dynamics. To validate our findings, we compare several network characteristics to the well-established finite-time Lyapunov exponents and apply a receiver operating characteristic analysis to identify network measures that are particularly useful for unveiling the skeleton of Lagrangian chaos.

Completely two-dimensional model for analysis of characteristics of linear induction cylindrical pump

International Nuclear Information System (INIS)

Kirillov, I.R.; Obukhov, D.M.

2005-01-01

One introduces a completely two-dimensional mathematical model to calculate characteristics of induction magnetohydrodynamic (MHD) machines with a cylindrical channel. On the basis of the numerical analysis one obtained a pattern of liquid metal flow in a electromagnetic pump at presence of the MHD-instability characterized by initiation of large-scale vortices propagating longitudinally and azimuthally. Comparison of the basic calculated characteristics of pump with the experiment shows their adequate qualitative and satisfactory quantitative coincidence [ru

Dual curved photonic crystal ring resonator based channel drop filter using two-dimensional photonic crystal structure

Energy Technology Data Exchange (ETDEWEB)

Chhipa, Mayur Kumar, E-mail: mayurchhipa1@gmail.com [Deptt. of Electronics and Communication Engineering, Government Engineering College Ajmer Rajasthan INDIA (India); Dusad, Lalit Kumar [Rajasthan Technical University Kota, Rajasthan (India)

2016-05-06

In this paper channel drop filter (CDF) is designed using dual curved photonic crystal ring resonator (PCRR). The photonic band gap (PBG) is calculated by plane wave expansion (PWE) method and the photonic crystal (PhC) based on two dimensional (2D) square lattice periodic arrays of silicon (Si) rods in air structure have been investigated using finite difference time domain (FDTD) method. The number of rods in Z and X directions is 21 and 20 respectively with lattice constant 0.540 nm and rod radius r = 0.1 µm. The channel drop filter has been optimized for telecommunication wavelengths λ = 1.591 µm with refractive indices 3.533. In the designed structure further analysis is also done by changing whole rods refractive index and it has been observed that this filter may be used for filtering several other channels also. The designed structure is useful for CWDM systems. This device may serve as a key component in photonic integrated circuits. The device is ultra compact with the overall size around 123 µm{sup 2}.

Topology of Flow Separation on Three-Dimensional Bodies

Science.gov (United States)

Chapman, Gary T.; Yates, Leslie A.

1991-01-01

In recent years there has been extensive research on three-dimensional flow separation. There are two different approaches: the phenomenological approach and a mathematical approach using topology. These two approaches are reviewed briefly and the shortcomings of some of the past works are discussed. A comprehensive approach applicable to incompressible and compressible steady-state flows as well as incompressible unsteady flow is then presented. The approach is similar to earlier topological approaches to separation but is more complete and in some cases adds more emphasis to certain points than in the past. To assist in the classification of various types of flow, nomenclature is introduced to describe the skin-friction portraits on the surface. This method of classification is then demonstrated on several categories of flow to illustrate particular points as well as the diversity of flow separation. The categories include attached, two-dimensional separation and three different types of simple, three-dimensional primary separation, secondary separation, and compound separation. Hypothetical experiments are utilized to illustrate the topological terminology and its role in characterizing these flows. These hypothetical experiments use colored oil injected onto the surface at singular points in the skin-friction portrait. Actual flow-visualization information, if available, is used to corroborate the hypothetical examples.

Effect of void fraction correlations on two-phase pressure drop during flow boiling in narrow rectangular channel

International Nuclear Information System (INIS)

Huang, Dong; Gao, Puzhen; Chen, Chong; Lan, Shu

2013-01-01

Highlights: • Most of the slip ratio models and the Lockhart–Martinelli parameter based models give similar results. • The drift flux void fraction models give relatively small values. • The effect of void fraction correlations on two-phase friction pressure drop is inconspicuous. • The effect of void fraction correlations on two-phase acceleration pressure drop is significant. - Abstract: The void fraction of water during flow boiling in vertical narrow rectangular channel is experimentally investigated. The void fraction is indirectly determined using the present experimental data with various void fraction correlations or models published in the open literature. The effects of mass flux, mass quality, system pressure and inlet subcooling on the void fraction and pressure drop are discussed in detail. In addition, comparison and discussion among the numerous void fraction correlations are carried out. The effect of void fraction correlations on two-phase pressure drop is presented as well. The results reveal that most of the slip ratio correlations and the Lockhart–Martinelli parameter based void fraction correlations have results close to each other at mass quality higher than 0.2. The drift flux void fraction correlations give small values which are incompatible with other models making it inapplicable for narrow rectangular channel. The alteration of void fraction correlations has an inconspicuous effect on two-phase frictional pressure drop, while an obvious effect on two-phase accelerational pressure drop during flow boiling in narrow rectangular channel

Wake structure and thrust generation of a flapping foil in two-dimensional flow

DEFF Research Database (Denmark)

Andersen, Anders Peter; Bohr, Tomas; Schnipper, Teis

2017-01-01

We present a combined numerical (particle vortex method) and experimental (soap film tunnel) study of a symmetric foil undergoing prescribed oscillations in a two-dimensional free stream. We explore pure pitching and pure heaving, and contrast these two generic types of kinematics. We compare...... measurements and simulations when the foil is forced with pitching oscillations, and we find a close correspondence between flow visualisations using thickness variations in the soap film and the numerically determined vortex structures. Numerically, we determine wake maps spanned by oscillation frequency...

Finite volume model for two-dimensional shallow environmental flow

Science.gov (United States)

Simoes, F.J.M.

2011-01-01

This paper presents the development of a two-dimensional, depth integrated, unsteady, free-surface model based on the shallow water equations. The development was motivated by the desire of balancing computational efficiency and accuracy by selective and conjunctive use of different numerical techniques. The base framework of the discrete model uses Godunov methods on unstructured triangular grids, but the solution technique emphasizes the use of a high-resolution Riemann solver where needed, switching to a simpler and computationally more efficient upwind finite volume technique in the smooth regions of the flow. Explicit time marching is accomplished with strong stability preserving Runge-Kutta methods, with additional acceleration techniques for steady-state computations. A simplified mass-preserving algorithm is used to deal with wet/dry fronts. Application of the model is made to several benchmark cases that show the interplay of the diverse solution techniques.

Fluid flow and fuel-air mixing in a motored two-dimensional Wankel rotary engine

Science.gov (United States)

Shih, T. I.-P.; Nguyen, H. L.; Stegeman, J.

1986-01-01

The implicit-factored method of Beam and Warming was employed to obtain numerical solutions to the conservation equations of mass, species, momentum, and energy to study the unsteady, multidimensional flow and mixing of fuel and air inside the combustion chambers of a two-dimensional Wankel rotary engine under motored conditions. The effects of the following engine design and operating parameters on fluid flow and fuel-air mixing during the intake and compression cycles were studied: engine speed, angle of gaseous fuel injection during compression cycle, and speed of the fuel leaving fuel injector.

Large Eddy Simulation of turbulent flows in compound channels with a finite element code

International Nuclear Information System (INIS)

Xavier, C.M.; Petry, A.P.; Moeller, S.V.

2011-01-01

This paper presents the numerical investigation of the developing flow in a compound channel formed by a rectangular main channel and a gap in one of the sidewalls. A three dimensional Large Eddy Simulation computational code with the classic Smagorinsky model is introduced, where the transient flow is modeled through the conservation equations of mass and momentum of a quasi-incompressible, isothermal continuous medium. Finite Element Method, Taylor-Galerkin scheme and linear hexahedrical elements are applied. Numerical results of velocity profile show the development of a shear layer in agreement with experimental results obtained with Pitot tube and hot wires. (author)

A three-dimensional stratigraphic model for aggrading submarine channels based on laboratory experiments, numerical modeling, and sediment cores

Science.gov (United States)

Limaye, A. B.; Komatsu, Y.; Suzuki, K.; Paola, C.

2017-12-01

Turbidity currents deliver clastic sediment from continental margins to the deep ocean, and are the main driver of landscape and stratigraphic evolution in many low-relief, submarine environments. The sedimentary architecture of turbidites—including the spatial organization of coarse and fine sediments—is closely related to the aggradation, scour, and lateral shifting of channels. Seismic stratigraphy indicates that submarine, meandering channels often aggrade rapidly relative to lateral shifting, and develop channel sand bodies with high vertical connectivity. In comparison, the stratigraphic architecture developed by submarine, braided is relatively uncertain. We present a new stratigraphic model for submarine braided channels that integrates predictions from laboratory experiments and flow modeling with constraints from sediment cores. In the laboratory experiments, a saline density current developed subaqueous channels in plastic sediment. The channels aggraded to form a deposit with a vertical scale of approximately five channel depths. We collected topography data during aggradation to (1) establish relative stratigraphic age, and (2) estimate the sorting patterns of a hypothetical grain size distribution. We applied a numerical flow model to each topographic surface and used modeled flow depth as a proxy for relative grain size. We then conditioned the resulting stratigraphic model to observed grain size distributions using sediment core data from the Nankai Trough, offshore Japan. Using this stratigraphic model, we establish new, quantitative predictions for the two- and three-dimensional connectivity of coarse sediment as a function of fine-sediment fraction. Using this case study as an example, we will highlight outstanding challenges in relating the evolution of low-relief landscapes to the stratigraphic record.

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

Directory of Open Access Journals (Sweden)

Vincent Casseau

2016-12-01

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

Unsteady flow around a two-dimensional section of a vertical axis turbine for tidal stream energy conversion

Directory of Open Access Journals (Sweden)

Hyun Ju Jung

2009-12-01

Full Text Available The two-dimensional unsteady flow around a vertical axis turbine for tidal stream energy conversion was investigated using a computational fluid dynamics tool solving the Reynolds-Averaged Navier-Stokes equations. The geometry of the turbine blade section was NACA653-018 airfoil. The computational analysis was done at several different angles of attack and the results were compared with the corresponding experimental data for validation and calibration. Simulations were then carried out for the two-dimensional cross section of a vertical axis turbine. The simulation results demonstrated the usefulness of the method for the typical unsteady flows around vertical axis turbines. The optimum turbine efficiency was achieved for carefully selected combinations of the number of blades and tip speed ratios.

Generalization of Spatial Channel Theory to Three-Dimensional x-y-z Transport Computations

International Nuclear Information System (INIS)

Abu-Shumays, I. K.; Hunter, M. A.; Martz, R. L.; Risner, J. M.

2002-01-01

Spatial channel theory, initially introduced in 1977 by M. L. Williams and colleagues at ORNL, is a powerful tool for shield design optimization. It focuses on so called ''contributon'' flux and current of particles (a fraction of the total of neutrons, photons, etc.) which contribute directly or through their progeny to a pre-specified response, such as a detector reading, dose rate, reaction rate, etc., at certain locations of interest. Particles that do not contribute directly or indirectly to the pre-specified response, such as particles that are absorbed or leak out, are ignored. Contributon fluxes and currents are computed based on combined forward and adjoint transport solutions. The initial concepts were considerably improved by Abu-Shumays, Selva, and Shure by introducing steam functions and response flow functions. Plots of such functions provide both qualitative and quantitative information on dominant particle flow paths and identify locations within a shield configuration that are important in contributing to the response of interest. Previous work was restricted to two dimensional (2-D) x-y rectangular and r-z cylindrical geometries. This paper generalizes previous work to three-dimensional x-y-z geometry, since it is now practical to solve realistic 3-D problems with multidimensional transport programs. As in previous work, new analytic expressions are provided for folding spherical harmonics representations of forward and adjoint transport flux solutions. As a result, the main integrals involve in spatial channel theory are computed exactly and more efficiently than by numerical quadrature. The analogy with incompressible fluid flow is also applied to obtain visual qualitative and quantitative measures of important streaming paths that could prove vital for shield design optimization. Illustrative examples are provided. The connection between the current paper and the excellent work completed by M. L. Williams in 1991 is also discussed

Flow Interactions of Two- and Three-Dimensional Networked Bio-Inspired Control Elements in an In-Line Arrangement.

Science.gov (United States)

Kurt, Melike; Moored, Keith

2018-04-19

We present experiments that examine the modes of interaction, the collective performance and the role of three-dimensionality in two pitching propulsors in an in-line arrangement. Both two-dimensional foils and three-dimensional rectangular wings of $AR = 2$ are examined. \\kwm{In contrast to previous work, two interaction modes distinguished as the coherent and branched wake modes are not observed to be directly linked to the propulsive efficiency, although they are linked to peak thrust performance and minimum power consumption as previously described \\cite[]{boschitsch2014propulsive}.} \\kwm{In fact, in closely-spaced propulsors peak propulsive efficiency of the follower occurs near its minimum power and this condition \\kwm{ reveals a} branched wake mode. Alternatively, for propulsors spaced far apart peak propulsive efficiency of the follower occurs near its peak thrust and this condition \\kwm{reveals a} coherent wake mode.} By examining the collective performance, it is discovered that there is an optimal spacing between the propulsors to maximize the collective efficiency. For two-dimensional foils the optimal spacing of $X^* = 0.75$ and the synchrony of $\\phi = 2\\pi /3$ leads to a collective efficiency and thrust enhancement of 50\\% and 32\\%, respectively, as compared to two isolated foils. In comparison, for $AR = 2$ wings the optimal spacing of $X^* = 0.25$ and the synchrony of $\\phi = 7\\pi /6$ leads to a collective efficiency and thrust enhancement of 30\\% and 22\\%, respectively. In addition, at the optimal conditions the collective lateral force coefficients in both the two- and three-dimensional cases are negligible, while operating off these conditions can lead to non-negligible lateral forces. Finally, the peak efficiency of the collective and the follower are shown to have opposite trends with increasing spacing in two- and three-dimensional flows. This is correlated to the breakdown of the impinging vortex on the follower wing in three

Quantum ballistic transport by interacting two-electron states in quasi-one-dimensional channels

Energy Technology Data Exchange (ETDEWEB)

Huang, Danhong [Air Force Research Laboratory, Space Vehicles Directorate, Kirtland Air Force Base, New Mexico 87117 (United States); Center for High Technology Materials, University of New Mexico, 1313 Goddard St SE, Albuquerque, New Mexico 87106 (United States); Gumbs, Godfrey [Center for High Technology Materials, University of New Mexico, 1313 Goddard St SE, Albuquerque, New Mexico 87106 (United States); Abranyos, Yonatan [Department of Physics and Astronomy, Hunter College of the City University of New York, 695 Park Avenue, New York, New York 10065 (United States); Pepper, Michael; Kumar, Sanjeev [Department of Electronic and Electrical Engineering, University College London, London, WC1E 7JE (United Kingdom); London Centre for Nanotechnology, 17-19 Gordon Street, London, WC1H 0AH (United Kingdom)

2015-11-15

For quantum ballistic transport of electrons through a short conduction channel, the role of Coulomb interaction may significantly modify the energy levels of two-electron states at low temperatures as the channel becomes wide. In this regime, the Coulomb effect on the two-electron states is calculated and found to lead to four split energy levels, including two anticrossing-level and two crossing-level states. Moreover, due to the interplay of anticrossing and crossing effects, our calculations reveal that the ground two-electron state will switch from one anticrossing state (strong confinement) to a crossing state (intermediate confinement) as the channel width gradually increases and then back to the original anticrossing state (weak confinement) as the channel width becomes larger than a threshold value. This switching behavior leaves a footprint in the ballistic conductance as well as in the diffusion thermoelectric power of electrons. Such a switching is related to the triple spin degeneracy as well as to the Coulomb repulsion in the central region of the channel, which separates two electrons away and pushes them to different channel edges. The conductance reoccurrence region expands from the weak to the intermediate confinement regime with increasing electron density.

Visualization of an air-water interface on superhydrophobic surfaces in turbulent channel flows

Science.gov (United States)

Kim, Hyunseok; Park, Hyungmin

2017-11-01

In the present study, three-dimensional deformation of air-water interface on superhydrophobic surfaces in turbulent channel flows at the Reynolds numbers of Re = 3000 and 10000 is measured with RICM (Reflection Interference Contrast Microscopy) technique. Two different types of roughness feature of circular hole and rectangular grate are considered, whose depth is 20 μm and diameter (or width) is varied between 20-200 μm. Since the air-water interface is always at de-pinned state at the considered condition, air-water interface shape and its sagging velocity is maintained to be almost constant as time goes one. In comparison with the previous results under the laminar flow, due to turbulent characteristics of the flow, sagging velocity is much faster. Based on the measured sagging profiles, a modified model to describe the air-water interface dynamics under turbulent flows is suggested. Supported by City of Seoul through Seoul Urban Data Science Laboratory Project (Grant No 0660-20170004) administered by SNU Big Data Institute.

Two-phase flow patterns in horizontal rectangular minichannel

Directory of Open Access Journals (Sweden)

Ron’shin Fedor

2016-01-01

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

Measurement system of bubbly flow using Ultrasonic Velocity Profile Monitor and Video Data Processing Unit. 3. Comparison of flow characteristics between bubbly cocurrent and countercurrent flows

International Nuclear Information System (INIS)

Zhou, Shirong; Suzuki, Yumiko; Aritomi, Masanori; Matsuzaki, Mitsuo; Takeda, Yasushi; Mori, Michitsugu

1998-01-01

The authors have developed a new measurement system which consisted of an Ultrasonic Velocity Profile Monitor (UVP) and a Video Data Processing Unit (VDP) in order to clarify the two-dimensional flow characteristics in bubbly flows and to offer a data base to validate numerical codes for two-dimensional two-phase flow. In the present paper, the proposed measurement system is applied to fully developed bubbly cocurrent flows in a vertical rectangular channel. At first, both bubble and water velocity profiles and void fraction profiles in the channel were investigated statistically. In addition, the two-phase multiplier profile of turbulence intensity, which was defined as a ratio of the standard deviation of velocity fluctuation in a bubbly flow to that in a water single phase flow, were examined. Next, these flow characteristics were compared with those in bubbly countercurrent flows reported in our previous paper. Finally, concerning the drift flux model, the distribution parameter and drift velocity were obtained directly from both bubble and water velocity profiles and void fraction profiles, and their results were compared with those in bubbly countercurrent flows. (author)

Numerical investigation of fluid mud motion using a three-dimensional hydrodynamic and two-dimensional fluid mud coupling model

Science.gov (United States)

Yang, Xiaochen; Zhang, Qinghe; Hao, Linnan

2015-03-01

A water-fluid mud coupling model is developed based on the unstructured grid finite volume coastal ocean model (FVCOM) to investigate the fluid mud motion. The hydrodynamics and sediment transport of the overlying water column are solved using the original three-dimensional ocean model. A horizontal two-dimensional fluid mud model is integrated into the FVCOM model to simulate the underlying fluid mud flow. The fluid mud interacts with the water column through the sediment flux, current, and shear stress. The friction factor between the fluid mud and the bed, which is traditionally determined empirically, is derived with the assumption that the vertical distribution of shear stress below the yield surface of fluid mud is identical to that of uniform laminar flow of Newtonian fluid in the open channel. The model is validated by experimental data and reasonable agreement is found. Compared with numerical cases with fixed friction factors, the results simulated with the derived friction factor exhibit the best agreement with the experiment, which demonstrates the necessity of the derivation of the friction factor.

Lagrangian statistics and flow topology in forced two-dimensional turbulence.

Science.gov (United States)

Kadoch, B; Del-Castillo-Negrete, D; Bos, W J T; Schneider, K

2011-03-01

A study of the relationship between Lagrangian statistics and flow topology in fluid turbulence is presented. The topology is characterized using the Weiss criterion, which provides a conceptually simple tool to partition the flow into topologically different regions: elliptic (vortex dominated), hyperbolic (deformation dominated), and intermediate (turbulent background). The flow corresponds to forced two-dimensional Navier-Stokes turbulence in doubly periodic and circular bounded domains, the latter with no-slip boundary conditions. In the double periodic domain, the probability density function (pdf) of the Weiss field exhibits a negative skewness consistent with the fact that in periodic domains the flow is dominated by coherent vortex structures. On the other hand, in the circular domain, the elliptic and hyperbolic regions seem to be statistically similar. We follow a Lagrangian approach and obtain the statistics by tracking large ensembles of passively advected tracers. The pdfs of residence time in the topologically different regions are computed introducing the Lagrangian Weiss field, i.e., the Weiss field computed along the particles' trajectories. In elliptic and hyperbolic regions, the pdfs of the residence time have self-similar algebraic decaying tails. In contrast, in the intermediate regions the pdf has exponential decaying tails. The conditional pdfs (with respect to the flow topology) of the Lagrangian velocity exhibit Gaussian-like behavior in the periodic and in the bounded domains. In contrast to the freely decaying turbulence case, the conditional pdfs of the Lagrangian acceleration in forced turbulence show a comparable level of intermittency in both the periodic and the bounded domains. The conditional pdfs of the Lagrangian curvature are characterized, in all cases, by self-similar power-law behavior with a decay exponent of order -2.

A heat transfer model for evaporating micro-channel coalescing bubble flow

International Nuclear Information System (INIS)

Consolini, L.; Thome, J.R.

2009-01-01

The current study presents a one-dimensional model of confined coalescing bubble flow for the prediction of micro-channel convective boiling heat transfer. Coalescing bubble flow has recently been identified as one of the characteristic flow patterns to be found in micro-scale systems, occurring at intermediate vapor qualities between the isolated bubble and the fully annular regimes. As two or more bubbles bond under the action of inertia and surface tension, the passage frequency of the bubble liquid slug pair declines, with a redistribution of liquid among the remaining flow structures. Assuming heat transfer to occur only by conduction through the thin evaporating liquid film surrounding individual bubbles, the present model includes a simplified description of the dynamics of the thin film evaporation process that takes into account the added mass transfer by breakup of the bridging liquid slugs. The new model has been confronted against experimental data taken within the coalescing bubble flow mode that have been identified by a diabatic micro-scale flow pattern map. The comparisons for three different fluids (R-134a, R-236fa and R-245fa) gave encouraging results with 83% of the database predicted within a ± 30% error band. (author)

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

International Nuclear Information System (INIS)

Park, Chan Wook; Lee, Sung Su

2008-01-01

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

Evolution of the vorticity-area density during the formation of coherent structures in two-dimensional flows

NARCIS (Netherlands)

Capel, H.W.; Pasmanter, R.A.

2000-01-01

It is shown: (1) that in two-dimensional, incompressible, viscous flows the vorticity-area distribution evolves according to an advection-diffusion equation with a negative, time dependent diffusion coefficient and (2) how to use the vorticity-stream function relations, i.e., the so-called

Experimental Research into the Two-Level Cylindrical Cyclone with a Different Number of Channels

Directory of Open Access Journals (Sweden)

Egidijus Baliukas

2014-10-01

Full Text Available The multichannel two-level cyclone has been designed for separating solid particles from airflow and built at the Laboratory of Environmental Protection Technologies of Vilnius Gediminas Technical University. The conducted research is aimed at determining air flow distribution at two levels and channels of the multichannel cyclone. The multifunctional meter Testo-400 and the dynamic Pitot tube have been used form easuring air flow rates in the channels. The obtained results show that the equal volume of air gets into two levels installed inside the cyclone, and rates are distributed equally in the channels of these levels. The maximum air flow rate is recorded in the first channel and occurs when half-rings are set in such positions so that 75% of air flow returns to the previous channel. The biggest aerodynamic resistance is 1660 Pa and has been recorded in the cyclone having eight channels under air flow distribution ratio 75/25. The highest air purification efficiency has been observed in the two-level six-channel cyclone under air flow distribution ratio 75/25. The effectiveness of separating granite particles is 92.1% and that of wood particles – 91.1 when the particles are up to 20 μm in diameter.

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)

Three-dimensional attached viscous flow basic principles and theoretical foundations

CERN Document Server

Hirschel, Ernst Heinrich; Kordulla, Wilhelm

2014-01-01

Viscous flow is usually treated in the frame of boundary-layer theory and as a two-dimensional flow. At best, books on boundary layers provide the describing equations for three-dimensional boundary layers, and solutions only for certain special cases.   This book presents the basic principles and theoretical foundations of three-dimensional attached viscous flows as they apply to aircraft of all kinds. Though the primary flight speed range is that of civil air transport vehicles, flows past other flying vehicles up to hypersonic speeds are also considered. Emphasis is put on general three-dimensional attached viscous flows and not on three-dimensional boundary layers, as this wider scope is necessary in view of the theoretical and practical problems that have to be overcome in practice.   The specific topics covered include weak, strong, and global interaction; the locality principle; properties of three-dimensional viscous flows; thermal surface effects; characteristic properties; wall compatibility con...

On the existence of steady flow in a channel with one porous wall or two accelerating walls

Directory of Open Access Journals (Sweden)

Chunqing Lu

1998-11-01

Full Text Available channel either with no-slip at one wall and constant uniform suction or injection through another wall, or with two accelerating walls. The flows are governed by the fourth order nonlinear differential equation $F^{iv}+R(FF'''-F'F''=0$. In the former case, the flow is subject to the boundary conditions $F(-1=F'(-1=F'(1=0$, $F(1=-1$. In the latter case, the boundary conditions are $F(-1=F(1=0$, $F'(-1=-1$, $F'(1 = 1$.

Effects of couple stresses in MHD channel flow

International Nuclear Information System (INIS)

Soundalgekar, V.M.; Aranake, R.N.

1977-01-01

An analysis of fully developed MHD channel flow of an electrically conducting incompressible fluid, taking into account the couple stresses, is carried out. Exact solutions are derived for velocity profiles, current density, skin-friction and coefficient of mass flux. They are influenced by the magnetic field, the loading parameter k, and the non-dimensional parameter (a=b 1 /lambda). Their variations with respect to M, k and a are represented graphically, this is followed by a physical discussion. It is observed that the couple stresses are more effective in the presence of a very weak magnetic field. (Auth.)

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

Evaluation of correlations of flow boiling heat transfer of R22 in horizontal channels.

Science.gov (United States)

Zhou, Zhanru; Fang, Xiande; Li, Dingkun

2013-01-01

The calculation of two-phase flow boiling heat transfer of R22 in channels is required in a variety of applications, such as chemical process cooling systems, refrigeration, and air conditioning. A number of correlations for flow boiling heat transfer in channels have been proposed. This work evaluates the existing correlations for flow boiling heat transfer coefficient with 1669 experimental data points of flow boiling heat transfer of R22 collected from 18 published papers. The top two correlations for R22 are those of Liu and Winterton (1991) and Fang (2013), with the mean absolute deviation of 32.7% and 32.8%, respectively. More studies should be carried out to develop better ones. Effects of channel dimension and vapor quality on heat transfer are analyzed, and the results provide valuable information for further research in the correlation of two-phase flow boiling heat transfer of R22 in channels.

Differential geometric structures of stream functions: incompressible two-dimensional flow and curvatures

International Nuclear Information System (INIS)

Yamasaki, K; Iwayama, T; Yajima, T

2011-01-01

The Okubo-Weiss field, frequently used for partitioning incompressible two-dimensional (2D) fluids into coherent and incoherent regions, corresponds to the Gaussian curvature of the stream function. Therefore, we consider the differential geometric structures of stream functions and calculate the Gaussian curvatures of some basic flows. We find the following. (I) The vorticity corresponds to the mean curvature of the stream function. Thus, the stream-function surface for an irrotational flow and that for a parallel shear flow correspond to the minimal surface and a developable surface, respectively. (II) The relationship between the coherency and the magnitude of the vorticity is interpreted by the curvatures. (III) Using the Gaussian curvature, stability of single and double point vortex streets is analyzed. The results of this analysis are compared with the well-known linear stability analysis. (IV) Conformal mapping in fluid mechanics is the physical expression of the geometric fact that the sign of the Gaussian curvature does not change in conformal mapping. These findings suggest that the curvatures of stream functions are useful for understanding the geometric structure of an incompressible 2D flow.

Flow Oriented Channel Assignment for Multi-radio Wireless Mesh Networks

Directory of Open Access Journals (Sweden)

Niu Zhisheng

2010-01-01

Full Text Available We investigate channel assignment for a multichannel wireless mesh network backbone, where each router is equipped with multiple interfaces. Of particular interest is the development of channel assignment heuristics for multiple flows. We present an optimization formulation and then propose two iterative flow oriented heuristics for the conflict-free and interference-aware cases, respectively. To maximize the aggregate useful end-to-end flow rates, both algorithms identify and resolve congestion at instantaneous bottleneck link in each iteration. Then the link rate is optimally allocated among contending flows that share this link by solving a linear programming (LP problem. A thorough performance evaluation is undertaken as a function of the number of channels and interfaces/node and the number of contending flows. The performance of our algorithm is shown to be significantly superior to best known algorithm in its class in multichannel limited radio scenarios.

Experimental Investigation of Pressure Drop and Pressure Distribution Along a Heated Channel in Subcooled Flow Boiling

International Nuclear Information System (INIS)

Aharon, Y.; Hochbaum, I.; Shai, I.

2002-01-01

The state of knowledge relating to pressure drop in subcooled boiling region is very unsatisfactory. That pressure drop is an important factor in considering the design of nuclear reactors because of the possibility of flow excursion during a two phase flow in the channels. In operational systems with multiple flow channels, an increase in pressure drop in one flow channel, can cause the flow to be diverted to other channels. A burnout can occur in the unstable channel

VNAP2: a computer program for computation of two-dimensional, time-dependent, compressible, turbulent flow

Energy Technology Data Exchange (ETDEWEB)

Cline, M.C.

1981-08-01

VNAP2 is a computer program for calculating turbulent (as well as laminar and inviscid), steady, and unsteady flow. VNAP2 solves the two-dimensional, time-dependent, compressible Navier-Stokes equations. The turbulence is modeled with either an algebraic mixing-length model, a one-equation model, or the Jones-Launder two-equation model. The geometry may be a single- or a dual-flowing stream. The interior grid points are computed using the unsplit MacCormack scheme. Two options to speed up the calculations for high Reynolds number flows are included. The boundary grid points are computed using a reference-plane-characteristic scheme with the viscous terms treated as source functions. An explicit artificial viscosity is included for shock computations. The fluid is assumed to be a perfect gas. The flow boundaries may be arbitrary curved solid walls, inflow/outflow boundaries, or free-jet envelopes. Typical problems that can be solved concern nozzles, inlets, jet-powered afterbodies, airfoils, and free-jet expansions. The accuracy and efficiency of the program are shown by calculations of several inviscid and turbulent flows. The program and its use are described completely, and six sample cases and a code listing are included.

A nonlinear model of flow in meandering submarine and subaerial channels

Science.gov (United States)

Imran, Jasim; Parker, Gary; Pirmez, Carlos

1999-12-01

A generalized model of flow in meandering subaqueous and subaerial channels is developed. The conservation equations of mass and momentum are depth/layer integrated, normalized, and represented as deviations from a straight base state. This allows the determination of integrable forms which can be solved at both linear and nonlinear levels. The effects of various flow and geometric parameters on the flow dynamics are studied. Although the model is not limited to any specific planform, this study focuses on sine-generated curves. In analysing the flow patterns, the turbidity current of the subaqueous case is simplified to a conservative density flow with water entrainment from above neglected. The subaqueous model thus formally corresponds to a subcritical or only mildly supercritical mud-rich turbidity current. By extension, however the analysis can be applied to a depositional or erosional current carrying sand that is changing only slowly in the streamwise direction. By bringing the subaqueous and subaerial cases into a common form, flow behaviour in the two environments can be compared under similar geometric and boundary conditions. A major difference between the two cases is the degree of superelevation of channel flow around bends, which is modest in the subaerial case but substantial in the subaqueous case. Another difference concerns Coriolis effects: some of the largest subaqueous meandering systems are so large that Coriolis effects can become important. The model is applied to meander bends on the youngest channel in the mid-fan region of the Amazon Fan and a mildly sinuous bend of the North-West Atlantic Mid-Ocean Channel. In the absence of specific data on the turbid flows that created the channel, the model can be used to make inferences about the flow, and in particular the range of values of flow velocity and sediment concentration that would allow the growth and downfan migration of meander bends.

Bubbly flows around a two-dimensional circular cylinder

Science.gov (United States)

Lee, Jubeom; Park, Hyungmin

2016-11-01

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

An investigation of turbulent catalytically stabilized channel flow combustion of lean hydrogen - air mixtures

Energy Technology Data Exchange (ETDEWEB)

Mantzaras, I; Benz, P; Schaeren, R; Bombach, R [Paul Scherrer Inst. (PSI), Villigen (Switzerland)

1999-08-01

The catalytically stabilised thermal combustion (CST) of lean hydrogen-air mixtures was investigated numerically in a turbulent channel flow configuration using a two-dimensional elliptic model with detailed heterogeneous and homogeneous chemical reactions. Comparison between turbulent and laminar cases having the same incoming mean properties shows that turbulence inhibits homogeneous ignition due to increased heat transport away from the near-wall layer. The peak root-mean-square temperature and species fluctuations are always located outside the extent of the homogeneous reaction zone indicating that thermochemical fluctuations have no significant influence on gaseous combustion. (author) 4 figs., 6 refs.

Two dimensional analytical model for a reconfigurable field effect transistor

Science.gov (United States)

Ranjith, R.; Jayachandran, Remya; Suja, K. J.; Komaragiri, Rama S.

2018-02-01

This paper presents two-dimensional potential and current models for a reconfigurable field effect transistor (RFET). Two potential models which describe subthreshold and above-threshold channel potentials are developed by solving two-dimensional (2D) Poisson's equation. In the first potential model, 2D Poisson's equation is solved by considering constant/zero charge density in the channel region of the device to get the subthreshold potential characteristics. In the second model, accumulation charge density is considered to get above-threshold potential characteristics of the device. The proposed models are applicable for the device having lightly doped or intrinsic channel. While obtaining the mathematical model, whole body area is divided into two regions: gated region and un-gated region. The analytical models are compared with technology computer-aided design (TCAD) simulation results and are in complete agreement for different lengths of the gated regions as well as at various supply voltage levels.

Investigation of magneto-hemodynamic flow in a semi-porous channel using orthonormal Bernstein polynomials

Science.gov (United States)

Hosseini, E.; Loghmani, G. B.; Heydari, M.; Rashidi, M. M.

2017-07-01

In this paper, the problem of the magneto-hemodynamic laminar viscous flow of a conducting physiological fluid in a semi-porous channel under a transverse magnetic field is investigated numerically. Using a Berman's similarity transformation, the two-dimensional momentum conservation partial differential equations can be written as a system of nonlinear ordinary differential equations incorporating Lorentizian magneto-hydrodynamic body force terms. A new computational method based on the operational matrix of derivative of orthonormal Bernstein polynomials for solving the resulting differential systems is introduced. Moreover, by using the residual correction process, two types of error estimates are provided and reported to show the strength of the proposed method. Graphical and tabular results are presented to investigate the influence of the Hartmann number ( Ha) and the transpiration Reynolds number ( Re on velocity profiles in the channel. The results are compared with those obtained by previous works to confirm the accuracy and efficiency of the proposed scheme.