Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Multiphase”

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 919 records · Page 51
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

An extended trajectory-mechanics approach for calculating two-phase flow paths

A technique originating in quantum dynamics is used to derive a trajectory-based, semi-analytical solution for two-phase flow. The partial differential equation governing the evolution of the aqueous phase is equivalent to a family of ordinary differential equations defined along a path through the porous medium. The trajectories may be found by solving the differential equations directly or by post-processing the output of a numerical solution to the full set of governing equations. The trajectories, which differ from conventional streamlines, are found to bend downward in response to gravitational forces. The curvature is more pronounced as the dip of the porous layer containing the flow increases. Subtle changes in the relative permeability curve can lead to significant variations in the trajectories. The ordinary differential equation for the trajectory provides an expression for the travel time along the path. The expression produces a semi-analytical approximation to the model parameter sensitivities, the partial derivatives of the travel times with respect to changes in the permeability model. The semi-analytical trajectory-based sensitivities generally agree with those computed using a numerical reservoir simulator and a perturbation approach. The sensitivities are useful in tomographic imaging algorithms designed to estimate the spatial variation in permeability within a porous medium using multiphase observations.

58 GEOSCIENCES↗
DOE OSTI2021

Investigation of theoretical scaling laws using large eddy simulations for airborne spreading of viral contagion from sneezing and coughing

Using a set of large eddy point-particle simulations, we explore the fluid dynamics of an ejected puff resulting from a cough/sneeze. The ejection contains over 61,000 potentially virus-laden droplets at an injection Reynolds number of about 46,000, comparable to an actual cough/ sneeze. We observe that global puff properties, such as centroid, puff volume, momentum, and buoyancy vary little across realizations. Other properties, such as maximum extent, shape, and edge velocity of the puff, may exhibit substantial variation. In many realizations, a portion of the puff splits off and advances along a random direction, while keeping airborne droplet nuclei afloat. This peeled-off portion provides a mechanism for virus-laden droplets to travel over large distances in a short amount of time. We also observe that the vast majority of droplets remain suspended within the puff after all liquid has evaporated. The main objectives of the study are to (i) evaluate assumptions of Balachandar’s et al. theory [Int. J. Multiphase Flow 132, 103439 (2020)], which include buoyancy effects, shape of the puff, and droplet evaporation rate, (ii) obtain values of closure parameters, which include location and time of the virtual origin, and puff entrainment and drag coefficients, and (iii) evaluate the accuracy of the theory in predicting the shape, size, and location of the puff, as well as droplet number density long after ejection. The theory adequately predicts global puff properties including size, velocity, and distance traveled, the largest size of droplets that exit the puff due to settling, and the droplet size distribution within the puff long after ejection.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Modeling atomically mixed graded density impactors

Graded density impactors (GDIs) are multi-material composite impactors used in gas gun experiments to tailor the drive conditions imparted to a sample test material. Previous graded density impactors generally rely on thin, but discrete, layers of different materials. The thinner and the greater number of layers will result in smoother compression. Taken to the limit of very thin layers would be pure material 1 at one surface, such as the front surface of an impactor, smoothly transitioning at the atomic scale to pure material 2 on the back surface. Such an impactor can initially shock, then smoothly compresses a material during a dynamic experiment. This type of experiment can serve to explore a larger region of thermodynamic space than a single or even multi-shock experiments. An overview of how graded density impactors are made is reviewed and sample results are given. Additionally, a strategy for modeling these kinds of impactors is presented. The length scales of constituent mixing are given from the experimental build through electrochemical-deposition. Equation of state models for pure constituents and their subsequent mixtures are presented. It is demonstrated that the time scales for pressure and temperature equilibration, for atomically mixed GDIs, are short enough to be a justifiable closure for the resulting multiphase flow. Furthermore, we present simulation results of dynamic shock followed by a ramp compression, utilizing a silver/gold graded density impactor, onto a tantalum sample.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Constraining the release of Sn to the ambient melting point following shock loading using time-resolved x-ray diffraction

The formation of liquid following release from a shocked state governs the transition from spall to cavitation and the formation of ejecta in metals. In order to build physics-based models of these processes, it is necessary to critically evaluate the relative importance of kinetics and entropy generation during the release along with the accuracy of multiphase equations of state. Tin (Sn) has served as a testbed for a variety of experiments examining strength and ejecta due to its accessible melt boundary and solid–solid phase transitions. This work presents experiments examining the phase evolution of high purity Sn following the shock and release to ambient stress near the melting point. Sn is found to release to states between its ambient solidus and liquidus from approximately 19 to 33 GPa under uniaxial loading, with the two-phase region being characterized by a reduction in the intensity of the (220), (211) β -Sn doublet. Jetting experiments performed at 27–28 GPa exhibit comparable diffraction patterns with what is observed following the uniaxial release. The solid fractions of β-Sn in the ambient mixed phase region are found to decrease linearly with increasing shock stress as increasing liquid Sn is formed. Here, the results provide much needed information for interpreting measurements of dynamic strength at a high strain rate and experiments examining cavitation and shallow bubble collapse in Sn.

, spall↗
DOE OSTI2019

Simulation of Boiling Two-Phase Flow in a Helical Coil Steam Generator Using the Spectral Element Code Nek-2P

To enable the design of a light water small modular reactor, the boiling flow inside a helical coil steam generator has been simulated with the two-fluid model in Nek-2P. Nek-2P is the multiphase branch of the spectral element code Nek5000. Details of the implementation of the two-fluid model and the included closure models are discussed. The presented closure models include interactions for momentum, heat, and mass transfer between phases. Models for the drag, lift, and turbulent dispersion forces are included. The complete model is fully consistent in the limits of both phasic volume fractions approaching zero and is able to simulate flows of dispersed vapor, continuous liquid, dispersed liquid, continuous vapor, or any combination thereof. The closure models and their implementation in Nek-2P have been validated by comparing to experimental data for a boiling flow, demonstrating excellent agreement. Results from the simulation of the helical coil indicate strong phasic separation driven by the effects of buoyancy and inertia. Finally, significant differences were observed in the results compared to simulations performed using Star-CCM+, although these differences were somewhat expected.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2020

Progress Toward Simulating Departure from Nucleate Boiling at High-Pressure Applications with Selected Wall Boiling Closures

Recently, a Eulerian-based two-fluid computational fluid dynamics (CFD) framework with a wall heat flux partitioning approach has been intensively investigated for departure from nucleate boiling (DNB) simulation under the U.S. Department of Energy–funded Consortium for Advanced Simulation of Light Water Reactors (CASL) program. Understanding of the DNB characteristics over a range of pressurized water reactor–like operating conditions and accurate prediction of boiling crisis in the nuclear power system have been grand challenges because of the large impact of DNB on reactor safety and operational economics. The ultimate goal of this task in the CASL program is to introduce a robust multiphase CFD–based DNB modeling framework that is capable of characterizing an entire boiling history in which the wall boiling mode experiences the following through multiple stages of heat transfer mode: (1) single-phase convective heat transfer, (2) nucleate boiling heat transfer, and (3) identification of the departure of nucleate boiling. To validate the CASL boiling model, we have benchmarked simulated DNB over three different flow channel configurations (pipe flow, 5 × 5 fuel bundle with mixing vane tests, and 5 × 5 fuel bundle without mixing vane tests) against experimental measurements, and the validation result with open literature is reported. The DNB detection criteria in the simulation are checked by monitoring the peak wall temperature, wall dryout factor, and net energy balance. In addition to the DNB performance test, some preliminary sensitivity results on closure model selection are reported to address the prediction capability of local void profile against measurements. The boiling simulation tested in this study exhibits a maximum deviation of 24% from the measured DNB value in a high-pressure (i.e., 138 bars) subcooled pipe flow test. The ranges of operating conditions are as follows: 1650 to 2650 kg/m 2 ·s for mass flux and 8.5 to 96 K for subcooled inlet temperature. The deviation is even reduced to 7% when the subcooled temperature is less than 40 K. Besides accuracy, base practice guidelines for DNB detection criteria are tested by monitoring three simulation variables: (1) maximum wall temperature, (2) wall dryout factor (i.e., K-value), and (3) energy balance. Numerical robustness of DNB simulation is largely achieved in most of the validation test except for a few high subcooled test cases.

42 ENGINEERING↗
DOE OSTI2017

Modernization and optimization of a legacy open-source CFD code for high-performance computing architectures

Legacy codes remain a crucial element of today's simulation-based engineering ecosystem due to the extensive validation process and investment in such software. The rapid evolution of high-performance computing architectures necessitates the modernization of these codes. One approach to modernization is a complete overhaul of the code. However, this could require extensive investments, such as rewriting in modern languages, new data constructs, etc., which will necessitate systematic verification and validation to re-establish the credibility of the computational models. The current study advocates using a more incremental approach and is a culmination of several modernization efforts of the legacy code MFIX, which is an open-source computational fluid dynamics code that has evolved over several decades, widely used in multiphase flows and still being developed by the National Energy Technology Laboratory. Two different modernization approaches,‘bottom-up’ and ‘top-down’, are illustrated. Here, preliminary results show up to 8.5x improvement at the selected kernel level with the first approach, and up to 50% improvement in total simulated time with the latter were achieved for the demonstration cases and target HPC systems employed.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2016

Plasma–liquid interactions: a review and roadmap

We present that plasma-liquid interactions represent a growing interdisciplinary area of research involving plasma science, fluid dynamics, heat and mass transfer, photolysis, multiphase chemistry and aerosol science. This review provides an assessment of the state-of-the-art of this multidisciplinary area and identifies the key research challenges. The developments in diagnostics, modeling and further extensions of cross section and reaction rate databases that are necessary to address these challenges are discussed. Finally, the review focuses on non-equilibrium plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Effect of the volume fraction gradient on the phase interaction force model for disperse two-phase flows

In this work, the effects of the particle volume fraction gradient on fluid-particle interactions are studied. The phase interaction force is decomposed into three terms. For the first term, namely the symmetrized force density, we present theoretical reasoning and numerical evidence to assume that it is independent of the particle volume fraction gradient. The second term is the particle volume fraction gradient times a newly introduced diffusion stress. The third term is the divergence of the particle-fluid-particle (PFP) stress. If this assumption of independence of the particle volume fraction gradient for the first term can be verified, to the first order of the ratio of the mean distance between particles to the macroscopic lengthscale, all three terms can be studied and modeled in flows with uniform particle distributions. Models thus obtained are applicable to statistically inhomogeneous flows, with the second and third terms accounting for statistical inhomogeneity. To verify this assumption, numerical simulations of flows passing fixed arrays of particles are performed. Both uniform and nonuniform particle volume fractions are studied and compared for disperse multiphase flows with the particle Reynolds numbers ranging from 1 to 100, and particle volume fraction ranging from 1% to 26% in statistically steady states. It is found that the symmetrized force (first) term can be well approximated by the drag force obtained from studies of uniform flows. The diffusion stress is positive along the flow direction and negative in the directions perpendicular to the flow. In the case of moving particles, this stress could potentially cause particle aggregation in the flow direction and dispersion in the directions perpendicular to the flow. Finally, the diffusion stress is only important when there is a volume fraction gradient, while the PFP stress can be important in inhomogeneous flows with either nonuniform particle concentrations or nonuniform average relative velocities between the phases.

42 ENGINEERING↗
DOE OSTI2020

Electronic Scanning Strategies in Adaptive Electrical Capacitance Volume Tomography: Tradeoffs and Prospects

Electrical Capacitance Volume Tomography (ECVT) has been applied for imaging of multiphase flows found in industrial applications. The ill-posed nature of the image reconstruction problem in ECVT and the consequent low resolution can be alleviated by employing electronic scanning enabled by electrode segmentation and reconfiguration during data acquisition. Here we study electronic scanning strategies that mimic physical rotation and shifting of the sensor along its symmetry axis. First, we study the feasibility of electronic scanning by analyzing the capacitance transducer circuit in SPICE. Then, we simulate electronic scanning using the finite element method for different electrode shapes, for which we compare the image reconstruction results and acquisition time. We find a noticeable improvement in image resolution for the scanning cases over conventional ECVT. Lastly, among the scanning cases, we emphasize a particular electrode shape that provides the best image resolution along with the minimal amount of acquisition time.

47 OTHER INSTRUMENTATION↗
DOE OSTI2017

The Fascinating and Complex Dynamics of Geyser Eruptions

Geysers episodically erupt liquid and vapor. Despite two centuries of scientific study, basic questions persist-why do geysers exist? What determines eruption intervals, durations, and heights? What initiates eruptions? Through monitoring eruption intervals, analyzing geophysical data, taking measurements within geyser conduits, performing numerical simulations, and constructing laboratory models, some of these questions have been addressed. Geysers are uncommon because they require a combination of abundant water recharge, magmatism, and rhyolite flows to supply heat and silica, and large fractures and cavities overlain by low-permeability materials to trap rising multiphase and multicomponent fluids. Eruptions are driven by the conversion of thermal to kinetic energy during decompression. Larger and deeper cavities permit larger eruptions and promote regularity by isolating water from weather variations. Furthermore, the ejection velocity may be limited by the speed of sound of the liquid + vapor mixture.

58 GEOSCIENCES↗
DOE OSTI2021

MFIX-Exa: A path toward exascale CFD-DEM simulations

MFIX-Exa is a computational fluid dynamics–discrete element model (CFD-DEM) code designed to run efficiently on current and next-generation supercomputing architectures. MFIX-Exa combines the CFD-DEM expertise embodied in the MFIX code—which was developed at NETL and is used widely in academia and industry—with the modern software framework, AMReX, developed at LBNL. The fundamental physics models follow those of the original MFIX, but the combination of new algorithmic approaches and a new software infrastructure will enable MFIX-Exa to leverage future exascale machines to optimize the modeling and design of multiphase chemical reactors.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2018

SECARB Pressure Monitoring Data

Bottom-hole, above zone monitoring interval, and injection zone pressure data collected during the SECARB project in Cranfield, Mississippi to assess the relationship between pressure field and multiphase field. Submission includes 10-second interval data from Detailed Area of Study wells: CFU31-F1 (injector), CFU31-F2 (observation), CFU31-F3 (observation) and Ella G Lees no. 7 (observation) well located west of the DAS. Associated Publications: Joy, C. A., 2011, The effects of pressure variation and chemical reactions on the elasticity of the lower Tuscaloosa sandstone of the Cranfield Field Mississippi, The University of Texas at Austin, Master’s thesis, 97 p. Kim, S., and Hosseini, S. A., 2013, Above-zone pressure monitoring and geomechanical analysis of a field scale CO2 injection, Cranfield Mississippi, Greenhouse Gases: Science and Technology, doi:10.1002/ghg.1388. Kim, S., and Hosseini, S. A., 2017, Study on the ratio of pore-pressure/stress changes during fluid injection and its implications for CO2 geologic storage: Journal of Petroleum Science and Engineering, v. 149, p. 138-150, doi:10.1016/j.petrol.2016.10.037. Mathias, S. A., Gluyas, J. G., Gonzalez Martinez de Miguel, G. J., and Hosseini, S. A., 2011, Role of partial miscibility on pressure buildup due to constant rate injection of CO2 into closed and open brine aquifers: Water Resources Research, v. 47, W12525, 11 p., doi:10.1029/2011WR011051. Meckel, T. A., Zeidouni, M., Hovorka, S. D., and Hosseini, S.A., 2013, Assessing sensitivity to well leakage from three years of continuous reservoir pressure monitoring during CO2 injection at Cranfield, MS, USA: International Journal of Greenhouse Gas Control, [insert volume no., page numbers], doi:10.1016/j.ijggc.2013.01.019. Nicot, J.-P., Oldenburg, C. M., Bryant, S. L., and Hovorka, S. D., 2009, Pressure perturbations from geologic carbon sequestration: area-of-review boundaries and borehole leakage driving forces, in Energy Procedia (v. 1, no.1), Proceedings of 9th International Conference on Greenhouse Gas Control Technologies, GHGT9, 16–20 November, Washington DC, p. 47–54. Tao, Q., Bryant, S. L., and Meckel, T. A., 2013, Modeling above-zone measurements of pressure and temperature for monitoring CCS sites: International Journal of Greenhouse Gas Control, v. 18, p. 523–530, doi:10.1016/j.ijggc.2012.08.011.

above zone monitoring interval↗
DOE OSTI2023

Optimization of a cyclone using MFIX and Nodeworks

Video depicting the optimization process of a cyclone on NETL's chemical looping reactor (CLR) using MFIX and Nodeworks. MFIX is used to model the cyclone using PIC. Nodeworks is then used to generate proposed geometry changes using a Latin hypercube. Each design is simulated, with an objective value being computed based on the cyclone efficiency and pressure drop. A Gaussian Process surrogate model is then constructed from the objective values. This surrogate model is then used by a differential evolution optimization algorithm to identify the optimal cyclone design. Details published here: Weber, J., Fullmer, W., Gel, A., and Musser, J. (February 4, 2020). "Optimization of a Cyclone Using Multiphase Flow Computational Fluid Dynamics." ASME. J. Fluids Eng. March 2020; 142(3): 031111. https://doi.org/10.1115/1.4045952 OSTI: https://www.osti.gov/pages/servlets/purl/1763893

cyclone↗
DOE OSTI2012

RELAP-7 Level 2 Milestone Report: Demonstration of a Steady State Single Phase PWR Simulation with RELAP-7

The document contains the simulation results of a steady state model PWR problem with the RELAP-7 code. The RELAP-7 code is the next generation nuclear reactor system safety analysis code being developed at Idaho National Laboratory (INL). The code is based on INL's modern scientific software development framework - MOOSE (Multi-Physics Object-Oriented Simulation Environment). This report summarizes the initial results of simulating a model steady-state single phase PWR problem using the current version of the RELAP-7 code. The major purpose of this demonstration simulation is to show that RELAP-7 code can be rapidly developed to simulate single-phase reactor problems. RELAP-7 is a new project started on October 1 st , 2011. It will become the main reactor systems simulation toolkit for RISMC (Risk Informed Safety Margin Characterization) and the next generation tool in the RELAP reactor safety/systems analysis application series (the replacement for RELAP5). The key to the success of RELAP-7 is the simultaneous advancement of physical models, numerical methods, and software design while maintaining a solid user perspective. Physical models include both PDEs (Partial Differential Equations) and ODEs (Ordinary Differential Equations) and experimental based closure models. RELAP-7 will eventually utilize well posed governing equations for multiphase flow, which can be strictly verified. Closure models used in RELAP5 and newly developed models will be reviewed and selected to reflect the progress made during the past three decades. RELAP-7 uses modern numerical methods, which allow implicit time integration, higher order schemes in both time and space, and strongly coupled multi-physics simulations. RELAP-7 is written with object oriented programming language C++. Its development follows modern software design paradigms. The code is easy to read, develop, maintain, and couple with other codes. Most importantly, the modern software design allows the RELAP-7 code to evolve with time. RELAP-7 is a MOOSE-based application. MOOSE (Multiphysics Object-Oriented Simulation Environment) is a framework for solving computational engineering problems in a well-planned, managed, and coordinated way. By leveraging millions of lines of open source software packages, such as PETSC (a nonlinear solver developed at Argonne National Laboratory) and LibMesh (a Finite Element Analysis package developed at University of Texas), MOOSE significantly reduces the expense and time required to develop new applications. Numerical integration methods and mesh management for parallel computation are provided by MOOSE. Therefore RELAP-7 code developers only need to focus on physics and user experiences. By using the MOOSE development environment, RELAP-7 code is developed by following the same modern software design paradigms used for other MOOSE development efforts. There are currently over 20 different MOOSE based applications ranging from 3-D transient neutron transport, detailed 3-D transient fuel performance analysis, to long-term material aging. Multi-physics and multiple dimensional analyses capabilities can be obtained by coupling RELAP-7 and other MOOSE based applications and by leveraging with capabilities developed by other DOE programs. This allows restricting the focus of RELAP-7 to systems analysis-type simulations and gives priority to retain and significantly extend RELAP5's capabilities.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2013

Annual Report: Innovative Process Technologies (30 September 2013)

Innovative Process Technologies is concerned with the development of innovative cost-effective technologies that promote efficiency, environmental performance, availability of advanced energy systems, and the development of computational tools that shorten development timelines of advanced energy systems. The project is focused on five distinct research tasks: (1) Sensors and Controls, (2) Power Electronics and Energetic Materials, (3) Innovative Energy Concepts (IEC), (4) Computational Materials, and (5) Multiphase Flow.

20 FOSSIL-FUELED POWER PLANTS↗
DOE OSTI2015

Fault-Related CO 2 Degassing, Geothermics, and Fluid Flow in Southern California Basins--Physiochemical Evidence and Modeling

Our studies have had an important impact on societal issues. Experimental and field observations show that CO 2 degassing, such as might occur from stored CO 2 reservoir gas, can result in significant stable isotopic disequilibrium. In the offshore South Ellwood field of the Santa Barbara channel, we show how oil production has reduced natural seep rates in the area, thereby reducing greenhouse gases. Permeability is calculated to be ~20-30 millidarcys for km-scale fault-focused fluid flow, using changes in natural gas seepage rates from well production, and poroelastic changes in formation pore-water pressure. In the Los Angeles (LA) basin, our characterization of formation water chemistry, including stable isotopic studies, allows the distinction between deep and shallow formations waters. Our multiphase computational-based modeling of petroleum migration demonstrates the important role of major faults on geological-scale fluid migration in the LA basin, and show how petroleum was dammed up against the Newport-Inglewood fault zone in a “geologically fast” interval of time (less than 0.5 million years). Furthermore, these fluid studies also will allow evaluation of potential cross-formational mixing of formation fluids. Lastly, our new study of helium isotopes in the LA basin shows a significant leakage of mantle helium along the Newport Inglewood fault zone (NIFZ), at flow rates up to 2 cm/yr. Crustal-scale fault permeability (~60 microdarcys) and advective versus conductive heat transport rates have been estimated using the observed helium isotopic data. The NIFZ is an important deep-seated fault that may crosscut a proposed basin decollement fault in this heavily populated area, and appears to allow seepage of helium from the mantle sources about 30 km beneath Los Angeles. The helium study has been widely cited in recent weeks by the news media, both in radio and on numerous web sites.

58 GEOSCIENCES↗
DOE OSTI2015

Experimental Results for Direct Electron Irradiation of a Uranyl Sulfate Solution: Bubble Formation and Thermal Hydraulics Studies

In support of the development of accelerator-driven production of fission product Mo-99 as proposed by SHINE Medical Technologies, a 35 MeV electron linac was used to irradiate depleted-uranium (DU) uranyl sulfate dissolved in pH 1 sulfuric acid at average power densities of 6 kW, 12 kW, and 15 kW. During these irradiations, gas bubbles were generated in the solution due to the radiolytic decomposition of water molecules in the solution. Multiple video cameras were used to record the behavior of bubble generation and transport in the solution. Seven six-channel thermocouples were used to record temperature gradients in the solution from self-heating. Measurements of hydrogen and oxygen concentrations in a helium sweep gas were recorded by a gas chromatograph to estimate production rates during irradiation. These data are being used to validate a computational fluid dynamics (CFD) model of the experiment that includes multiphase flow and a custom bubble injection model for the solution region.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗