Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Interfacial flows”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 271 records · Page 15
All sourcesNASA NTRSDOE OSTI
DOE OSTI2017

Computation of Air Discharge into a Water Pool by Using Volume-of-Fluid Method

Suppression tanks, in-containment refueling water storage tanks (IRWSTs) and filtered containment venting systems (FCVSs) of nuclear power plants use water pools to accommodate high-pressure and high-temperature discharge of non-condensable gases and/or steam. Dynamic behavior of gas bubbles inside a water pool determines pressure and temperature loading to the system and thus we need validated method for applications. The underlying physics relevant to these problems is still a difficult issue due to the complicated nature of gas- liquid interfacial transfers of mass, momentum, and energy in a two-phase flow. Therefore, notwithstanding many studies performed over the past two or three decades, two-fluid theory is still a formidable task and direct numerical simulation (DNS) needs huge computing cost. Pellegrini et al. analyzed steam condensation and chugging in the water pool of the PPOOLEX STB-24 test by using STAR-CCM+ code implemented with two-phase interface model based on the Rayleigh-Taylor instability theory (RTI) which triggers steam depressurization. For the solution of the fine interface, they used Eulerian-Eulerian model for the fluid domain since the volume-of-fluid (VOF) model requires larger mesh points. The RTI resulted in complete steam condensation consistent with the experiment. The fluid tank was open to the atmosphere. They reported the volume fraction distributions in the pool, not the pool temperature or pressure history. For non-condensing gas discharges, however, the VOF is still a good alternative to the Eulerian two-fluid or the DNS approaches. The spatial distribution of each phase is defined in terms of volume fractions and the cell physical properties are obtained from linear combinations of those of constituent phases weighted by each volume fraction. Therefore, conditions where thermal and phase change effects become important are less likely to be simulated realistically. The VOF model is suited for the simulation of flows where each phase constitutes a large structure with a relatively small total contact area between phases. Present work is also interested in this VOF method for application to air discharge into a water pool which occurs in the early operating stages of the said facilities. Our purpose is to validate the VOF method against an existing air discharge experiment of PPOOLEX by using a commercial CFD code, STAR-CCM+ version 10.04, in terms of macroscopic flow properties. Effect of turbulence is also discussed.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
NASA NTRS2004

Fixed Packed Bed Reactors in Reduced Gravity

We present experimental data on flow pattern transitions, pressure drop and flow characteristics for cocurrent gas-liquid flow through packed columns in microgravity. The flow pattern transition data indicates that the pulse flow regime exists over a wider range of gas and liquid flow rates under microgravity conditions compared to 1-g and the widely used Talmor map in 1-g is not applicable for predicting the transition boundaries. A new transition criterion between bubble and pulse flow in microgravity is proposed and tested using the data. Since there is no static head in microgravity, the pressure drop measured is the true frictional pressure drop. The pressure drop data, which has much smaller scatter than most reported 1-g data clearly shows that capillary effects can enhance the pressure drop (especially in the bubble flow regime) as much as 200% compared to that predicted by the single phase Ergun equation. The pressure drop data are correlated in terms of a two-phase friction factor and its dependence on the gas and liquid Reynolds numbers and the Suratman number. The influence of gravity on the pulse amplitude and frequency is also discussed and compared to that under normal gravity conditions. Experimental work is planned to determine the gas-liquid and liquid-solid mass transfer coefficients. Because of enhanced interfacial effects, we expect the gas-liquid transfer coefficients kLa and kGa (where a is the gas-liquid interfacial area) to be higher in microgravity than in normal gravity at the same flow conditions. This will be verified by gas absorption experiments, with and without reaction in the liquid phase, using oxygen, carbon dioxide, water and dilute aqueous amine solutions. The liquid-solid mass transfer coefficient will also be determined in the bubble as well as the pulse flow regimes using solid benzoic acid particles in the packing and measuring their rate of dissolution. The mass transfer coefficients in microgravity will be compared to those in normal gravity cocurrent flow to determine the mass transfer enhancement and propose new mass transfer correlations for two-phase gas-liquid flows through packed beds in microgravity.

Motil, Brian J.↗
DOE OSTI2016

Uncertainty Quantification with MCMC method and BFBT Data

The reliability of predictions of the system codes is closely related to the validation of their physical models. For example, the accuracy of void fraction prediction in a Boiling Water Reactor (BWR) using TRACE code depends on the uncertainties of closure relations developed for the two-phase flow models. This work focus on quantifying the uncertainties of two physical models, interfacial drag model and subcooled boiling model, using Markov Chain Monte Carlo (MCMC) method and BFBT experimental benchmark data. We applied the MCMC and MLE algorithms to BFBT benchmark data and estimated the uncertainty distribution of two model parameters: subcooled boiling heat transfer coefficient and interfacial drag (bubbly/slug rod bundle) coefficient. These algorithms tend to give various estimation results. Despite the differences in these algorithms, we find that the estimation results are consistent. After obtaining the uncertainty estimates, we validated the uncertainty estimates by applying the estimated distributions of model parameters to TRACE prediction. The new predictions succeed in reducing the void fraction prediction error, but it is not possible to eliminate all prediction error. We are only adjusting two model parameters, and there exists uncertainties in all model parameters. Application of the algorithms to other physical models and experimental data will be helpful for future work. (authors)

42 ENGINEERING↗
DOE OSTI2022

Heat transfer in directly-irradiated high-temperature solid–gas flows laden with polydisperse particles

Heat transfer in directly-irradiated high-temperature solid–gas flows laden with polydisperse particles is investigated using a novel transient three-dimensional computational fluid dynamics model. The model couples particle–gas hydrodynamics of solid–gas flows laden with polydisperse particles, radiative heat transfer in non-grey absorbing, emitting and anisotropically-scattering multi-component participating media, conduction heat transfer in the gas phase, and interfacial convection heat transfer. The multiphase particle-in-cell method is used to predict high-fidelity solid–gas flow characteristics, such as the local discrete particle size distribution, with increased computational efficiency by combining the advantages of both Eulerian and Lagrangian methods. The multi-component radiative transfer model is implemented using an advanced collision-based Monte Carlo ray-tracing method. The number of the prescribed discrete particle components is found to be the key parameter affecting the computational accuracy and efficiency, which primarily depends on the size distribution of the particles. For the model particle–gas flow featuring free-falling Gamma-distributed ceramic particles exposed to concentrated solar irradiation, the particle volume fraction, radiative, fluid flow and thermal characteristics appear to converge with the increasing number of the discrete particle components. Five particle components are sufficient to obtain physically meaningful results. A further increase in the number of the particle components only slightly increases the accuracy of the numerical predictions at the expense of a rapidly increasing computational time. For five particle components, the particle vertical velocity at the receiver exit for particles with the diameter of 43.4 μm is 57% of that for the particles with the diameter of 202.8 μm. The temperatures of these two particle components increase from the initial ambient values by factors of 2 and 1.2, respectively, during the simulation time. In conclusion, the model developed allows for increased fidelity of particle–gas flow simulations with significant radiative effects.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2015

Flow of immiscible ferrofluids in a planar gap in a rotating magnetic field

Analytical solutions are obtained for the steady, fully developed flow of two layers of immiscible ferrofluids of different thicknesses between two parallel plates. Interfacial linear and internal angular momentum balance relations are derived for the case when there is a ferrofluid-ferrofluid interface to obtain the translational and spin velocity profiles in the gap. As expected for the limit of low applied field amplitude, the magnitude of the translational velocity is directly proportional to the frequency of the applied magnetic field and to the square of the magnetic field amplitude. Expressions for the velocity profiles are obtained for the zero spin viscosity and non-zero spin viscosity cases and the effect of applied pressure gradient on the flows is studied. The spin velocity in both ferrofluid phases is in the direction of the rotating magnetic field, except for cases of extreme applied pressure gradients for which the fluid vorticity opposes the spin. We find that for the case of non-zero spin viscosity, flow reversals are predicted using representative ferrofluid property values and field conditions. The unique predictions of the solution with non-zero spin viscosity could be used to experimentally test the existence of couple stresses in ferrofluids and the validity of previously reported values of the so-called spin viscosity.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1983

Motion driven by the interface

Due to the reduction in buoyant forces aboard orbiting spacecraft such as the Space Shuttle, fluid motion driven by gradients in interfacial tension will be important in the processing of materials in space. In this paper, preliminary results from a study of surface tension driven flow in a pendant drop are reported. The drop is heated from above, and the resulting temperature gradients on the drop surface give rise to interfacial tension gradients. These, in turn, drive a circulation in the drop which is made visible by suitable tracers. The velocities are measured using a video technique, and the data on core velocities are found to agree well with results from a predictive theoretical model.

Jayaraj, K.↗
DOE OSTI2017

Simplified Pore Network Model for Analysis of Flow Capacity and Residual Oil Saturation Distribution Based on Computerized Tomography

The pore network (connected tunnels) can be used not only to describe pore geometry inside a rock but also to calculate various fluid flow properties on a micro scale. Based on a “digital core” extracted from micro computerized tomography (μ-CT), we have developed a simplified pore network model. This study simplifies the pore structure as a lattice of cubes. Pore size, pore distribution, tunnel connectivity, fluid viscosity, interfacial tension, and the external driving pressure gradient are considered for each adjacent cube to evaluate the flow capacity of the pore network. We calculated the permeability, effective permeability, and effective sweep efficiency. The results show that fluid flow initially starts in the largest tunnel, which corresponds to the threshold pressure of the rock sample (fluid flow onset pressure). Permeability changes nonlinearly before all the tunnels are involved in the flow. Changes in the microscopic residual oil saturation and the sweep efficiency were visualized based on the three-dimensional pore distribution. Visualization shows that more oil remains in the formation, the heterogeneity is greater, and the oil viscosity is higher.

58 GEOSCIENCES↗
DOE OSTI2017

Evaluation of the shear rate effect on the interfacial forces acting on a single bubble

Owing to the outstanding characteristics of mass and heat transfer efficiency, two-phase bubbly flows have wide applications in the nuclear engineering fields. However, the physics behind the bubbly flows phenomenon can be very complex due to the various parameters and phenomena involved. Thanks to the increasing computing power and advanced algorithms, computational fluid dynamics (CFD) approach has emerged as an effective tool for exploring the bubbly flows. In typical laminar channel or pipe flow, the flow field develops a parabolic velocity profile due to the wall friction. Assuming the shear rate across the bubble diameter is constant, the bubble will encounter various shear rates depending on its transverse position. Different shear rates may lead to the acceleration of transverse migration speed and influence the bubble rising velocity. In this summary, we evaluate the effect of shear rate on the interfacial forces (lift and drag) acting on a single bubble with the shear rate ranging from 3.8 to 8.3 s{sup -1} and the Eotvos number ranging from 2.58 to 4.45.

42 ENGINEERING↗
DOE OSTI2016

New analytical solutions to the two-phase water faucet problem

Here, the one-dimensional water faucet problem is one of the classical benchmark problems originally proposed by Ransom to study the two-fluid two-phase flow model. With certain simplifications, such as massless gas phase and no wall and interfacial frictions, analytical solutions had been previously obtained for the transient liquid velocity and void fraction distribution. The water faucet problem and its analytical solutions have been widely used for the purposes of code assessment, benchmark and numerical verifications. In our previous study, the Ransom’s solutions were used for the mesh convergence study of a high-resolution spatial discretization scheme. It was found that, at the steady state, an anticipated second-order spatial accuracy could not be achieved, when compared to the existing Ransom’s analytical solutions. A further investigation showed that the existing analytical solutions do not actually satisfy the commonly used two-fluid single-pressure two-phase flow equations. In this work, we present a new set of analytical solutions of the water faucet problem at the steady state, considering the gas phase density’s effect on pressure distribution. This new set of analytical solutions are used for mesh convergence studies, from which anticipated second-order of accuracy is achieved for the 2nd order spatial discretization scheme. In addition, extended Ransom’s transient solutions for the gas phase velocity and pressure are derived, with the assumption of decoupled liquid and gas pressures. Numerical verifications on the extended Ransom’s solutions are also presented.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2017

Implementation of a pressure correction term in the CTF sub-channel code

CTF [1] is a modernized and improved version of the legacy subchannel code, COBRA-TF [2], that is being jointly developed and maintained by Oak Ridge National Laboratory (ORNL) and North Carolina State University. The code was adopted by ORNL for use in the Consortium for Advanced Simulation of Light Water Reactors (CASL) in 2012 for aiding in addressing CASL challenge problems. Since that time, activities related to CTF have included implementing software quality assurance measures, implementing new features and models, performing validation and verification testing, establishing and supporting a CTF User Group, and developing the code for use in coupled applications [3, 4, 5]. CTF is applicable to single- and two-phase flows in light water reactor geometries at normal and accident operating conditions. The two-phase flow model consists of the six-equation two-phase flow model [6] augmented by equations for the droplet field and also includes appropriate source terms. This work focuses on the numerical stability of the six-equation two-phase flow model and is part of a recent eort to improve CTF capabilities for simulating flows in boiling nuclear reactors (BWRs). As it is of common knowledge among the thermodynamic community, the six-equation two-phase flow model is not strongly hyperbolic, meaning some of its eigenvalues are either not real or not distinct [7]. In a Cauchy sense and for an initial-value problem, strong hyperbolicity was found to be equivalent to well-posedness which ensures uniqueness of the solution in the space-time used. This concept is also of importance when designing high-order numerical methods and also when implementing boundary conditions that rely on real eigenvalues such as Riemann solvers. A tremendous amount of research is available in the literature on how to hyperbolize the six-equation two-phase flow model by adding ad hoc terms to the system of equations often referred to as virtual mass or interfacial pressure correction term. From a mathematical prospective, the main idea driving this approach is to add terms to the phasic momentum equation containing first-order partial derivatives of the unknowns to modify the Jacobian matrix and thus canceling the imaginary part of the eigenvalues.For this work it was chosen to implement a interfacial pressure correction (IPC) term on the same model as what was done for the CATHARE system code [8], that allows to recover a two-pressure formulation. In the following, the six-equation two-phase flow model is recalled along with the IPC term and its influence on the eigenvalues. Then, numerical results and convergence studies are presented for a water-faucet problem [9] at a pressure of 10 bar. Finally, conclusions and future work are addressed.

Delchini, Marco↗
DOE OSTI2014

Conjugate heat and mass transfer in the lattice Boltzmann equation method

An interface treatment for conjugate heat and mass transfer in the lattice Boltzmann equation method is proposed based on our previously proposed second-order accurate Dirichlet and Neumann boundary schemes. The continuity of temperature (concentration) and its flux at the interface for heat (mass) transfer is intrinsically satisfied without iterative computations, and the interfacial temperature (concentration) and their fluxes are conveniently obtained from the microscopic distribution functions without finite-difference calculations. The present treatment takes into account the local geometry of the interface so that it can be directly applied to curved interface problems such as conjugate heat and mass transfer in porous media. For straight interfaces or curved interfaces with no tangential gradient, the coupling between the interfacial fluxes along the discrete lattice velocity directions is eliminated and thus the proposed interface schemes can be greatly simplified. Several numerical tests are conducted to verify the applicability and accuracy of the proposed conjugate interface treatment, including (i) steady convection-diffusion in a channel containing two different fluids, (ii) unsteady convection-diffusion in the channel, (iii) steady heat conduction inside a circular domain with two different solid materials, and (iv) unsteady mass transfer from a spherical droplet in an extensional creeping flow. The accuracy and order of convergence of the simulated interior temperature (concentration) field, the interfacial temperature (concentration), and heat (mass) flux are examined in detail and compared with those obtained from the "half-lattice division" treatment in the literature. The present analysis and numerical results show that the half-lattice division scheme is second-order accurate only when the interface is fixed at the center of the lattice links, while the present treatment preserves second-order accuracy for arbitrary link fractions. For curved interfaces, the present treatment yields second-order accurate interior and interfacial temperatures (concentrations) and first-order accurate interfacial heat (mass) flux. An increase of order of convergence by one degree is obtained for each of these three quantities compared with the half-lattice division scheme. The surface-averaged Sherwood numbers computed in test (iv) agree well with published results.

Li, LK↗
NASA NTRS1989

The lubrication force between two viscous drops

The present determination of the hydrodynamic force resisting the relative motion of two unequal drops moving in Stokes flow conditions along their line of centers assumes that the drops are in near-contact, and possess sufficiently high interfacial tension to remain spherical. Depending on the ratio of drop viscosity to that of the continuous phase, as well as on the ratio of the distance between drops to their reduced radius, three possible flow situations arise: these correspond to nearly-rigid drops, drops with partially mobile interfaces, and drops with fully mobile interfaces.

Davis, Robert H.↗
NASA NTRS2003

Investigation of Critical Heat Flux in Reduced Gravity Using Photomicrographic Techniques

Experiments were performed to examine the effects of body force on flow boiling critical heat flux (CHF). FC-72 was boiled along one wall of a transparent rectangular flow channel that permitted photographic study of the vapor-liquid interface just prior to CHF. High-speed video imaging techniques were used to identify dominant CHF mechanisms corresponding to different flow orientations and liquid velocities. Six different CHF regimes were identified: Wavy Vapor Layer, Pool Boiling, Stratification, Vapor Counterflow, Vapor Stagnation, and Separated Concurrent Vapor Flow. CHF showed significant sensitivity to orientation for flow velocities below 0.2 m/s, where extremely low CHF values where measured, especially with downward-facing heated wall and downflow orientations. High flow velocities dampened the effects of orientation considerably. The CHF data were used to assess the suitability of previous CHF models and correlations. It is shown the Interfacial Lift-off Model is very effective at predicting CHF for high velocities at all orientations. The flooding limit, on the other hand, is useful at estimating CHF at low velocities and for downflow orientations. A new method consisting of three dimensionless criteria is developed for determining the minimum flow velocity required to overcome body force effects on near-saturated flow boiling CHF. Vertical upflow boiling experiments were performed in pursuit of identifying the trigger mechanism for subcooled flow boiling CHF. While virtually all prior studies on flow boiling CHF concern the prediction or measurement of conditions that lead to CHF, this study was focused on events that take place during the CHF transient. High-speed video imaging and photomicrographic techniques were used to record the transient behavior of interfacial features from the last steady-state power level before CHF until the moment of power cut-off following CHF. The video records show the development of a wavy vapor layer which propagates along the heated wall, permitting cooling prior to CHF only in wetting fronts corresponding to the wave troughs. Image analysis software was developed to estimate void fraction from the individual video images. The void fraction records for subcooled flow boiling show the CHF transient is accompanied by gradual lift-off of wetting fronts culminating in some maximum vapor layer mean thickness, following which the vapor layer begins to thin down as the transition to film boiling ensues. This study proves the Interfacial Lift-off Model, which has been validated for near-saturated flow boiling CHF, is equally valid for subcooled conditions.

Mudawar, Issam↗
DOE OSTI2019

Silica Surface States and their Wetting Characteristics

The wetting characteristics of silica (SiO 2 ) surfaces can be described by molecular dynamics (MD) and ab initio simulations, including comparison of silica surfaces (talc (001), siloxanated quartz, tridymite (001) and quartz (001)), some of which have not been considered previously in the literature. Classic MD and ab initio simulation methods have been used to determine the contact angle, interfacial water structure, hydroxylation reaction and hydration energy, the results of which are compared with experimental results reported in the literature. Additionally, wetting of silica surfaces depends on surface polarity and extent of hydroxylation. The non-polar siloxane surfaces are characterized by a contact angle of about 80°, an MD ‘water exclusion zone’ of about 3 Å, a relaxed interfacial water orientation, inertness to hydroxylation and minimal hydration energy. The polar silica surfaces can be wetted by water and have a more ordered interfacial water structure. Silanol groups form at the polar silica surface during hydroxylation reactions, and the calculated hydration energy of -1·2 and -1·6 eV matches the experimental heat of immersion measurements reported in the literature, which correspond to hydrogen (H) bonding with interfacial water. In conclusion, fundamental understanding of silica surfaces is important for understanding flotation phenomena and fluid flow in silica nanopores.

36 MATERIALS SCIENCE↗
NASA NTRS2001

The Dynamics of Miscible Fluids: A Space Flight Experiment (MIDAS)

We propose a space flight experiment to study the dynamics of miscible interfaces. A less viscous fluid displaces one of higher viscosity within a tube. The two fluids are miscible in all proportions. An intruding "finger" forms that occupies a fraction of the tube. As time progresses diffusion at the interface combined with flow induced straining between the two fluids modifies the concentration and velocity distributions within the whole tube. Also, under such circumstances it has been proposed that the interfacial stresses could depend on the local concentration gradients (Korteweg stresses) and that the divergence of the velocity need not be zero, even though the flow is incompressible. We have obtained reasonable agreement for the tip velocity between numerical simulations (that ignored the Korteweg stress and divergence effects) and physical experiments only at high Peclet Numbers. However at moderate to low Pe agreement was poor. As one possibility we attributed this lack of agreement to the disregard of these effects. We propose a space experiment to measure the finger shape, tip velocity, and the velocity and concentration fields. From intercomparisons between the experiment and the calculations we can then extract values for the coefficients of the Korteweg stress terms and confirm or deny the importance of these stresses.

Maxworthy, T.↗
NASA NTRS2001

The Dynamics of Miscible Fluids: A Space Flight Experiment (MIDAS)

We propose a space flight experiment to study the dynamics of miscible interfaces. A less viscous fluid displaces one of higher viscosity within a tube. The two fluids are miscible in all proportions. An intruding "finger" forms that occupies a fraction of the tube. As time progresses diffusion at the interface combined with flow induced straining between the two fluids modifies the concentration and velocity distributions within the whole tube. Also, under such circumstances it has been proposed that the interfacial stresses could depend on the local concentration gradients (Korteweg stresses) and that the divergence of the velocity need not be zero, even though the flow is incompressible. We have obtained reasonable agreement for the tip velocity between numerical simulations (that ignored the Korteweg stress and divergence effects) and physical experiments only at high Peelet Numbers. However at moderate to low Pe agreement was poor. As one possibility we attributed this lack of agreement to the disregard of these effects. We propose a space experiment to measure the finger shape, tip velocity, and the velocity and concentration fields. From intercomparisons between the experiment and the calculations we can then extract values for the coefficients of the Korteweg stress terms and confirm or deny the importance of these stresses.

Maxworthy, T.↗
DOE OSTI2014

Theory and Implementation of Nuclear Safety System Codes - Part II: System Code Closure Relations, Validation, and Limitations

This is Part II of two articles describing the details of thermal-hydraulic sys- tem codes. In this second part of the article series, the system code closure relationships (used to model thermal and mechanical non-equilibrium and the coupling of the phases) for the governing equations are discussed and evaluated. These include several thermal and hydraulic models, such as heat transfer coefficients for various flow regimes, two phase pressure correlations, two phase friction correlations, drag coefficients and interfacial models be- tween the fields. These models are often developed from experimental data. The experiment conditions should be understood to evaluate the efficacy of the closure models. Code verification and validation, including Separate Effects Tests (SETs) and Integral effects tests (IETs) is also assessed. It can be shown from the assessments that the test cases cover a significant section of the system code capabilities, but some of the more advanced reactor designs will push the limits of validation for the codes. Lastly, the limitations of the codes are discussed by considering next generation power plants, such as Small Modular Reactors (SMRs), analyz- ing not only existing nuclear power plants, but also next generation nuclear power plants. The nuclear industry is developing new, innovative reactor designs, such as Small Modular Reactors (SMRs), High-Temperature Gas-cooled Reactors (HTGRs) and others. Sub-types of these reactor designs utilize pebbles, prismatic graphite moderators, helical steam generators, in- novative fuel types, and many other design features that may not be fully analyzed by current system codes. This second part completes the series on the comparison and evaluation of the selected reactor system codes by discussing the closure relations, val- idation and limitations. These two articles indicate areas where the models can be improved to adequately address issues with new reactor design and development.

99 GENERAL AND MISCELLANEOUS↗