Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Multiphase 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 361 records · Page 20
All sourcesNASA NTRSDOE OSTI
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 OSTI2017

Simulating Bubble Deformations Using the Intersection Marker (ISM) Front-Tracking method

Surface and volume tracking methods have been under continuous development over the last 50+ years since the first MAC (Marker and Cell) and SLIC (Simple Line Interface Calculation) methods appeared in the 1960's. The sustained interest in simulating multiphase flow stem from their many practical applications in the physical world; from biological systems, to chemical systems, to boiling dynamics and atmospheric heat and mass transfer. The ability to capture the details of interface evolution as well as the physics of interfacial heat and mass transfer has thus been a tempting goal for many numerical researchers. Numerical analysts, apart from the well-established and popular VOF PLIC (Volume of Fluid, Piecewise Linear Interface Construction) method, can now call upon even more computationally intensive options such as the LS (Level Set) method, the FT (Front Tracking) method, the MoF (Moment of Fluid) method and hybrid methods.

42 ENGINEERING↗
DOE OSTI2022

A particle resolved simulation approach for studying shock interactions with moving, colliding solid particles

This work applies a new combination of techniques for the fully resolved simulation of compressible, gas–particle multiphase flows. Here, the adaptive wavelet collocation method is used to dynamically, and efficiently, adapt the computational grid to localized flow features and the particles. A characteristic-based volume penalization method that imposes arbitrary Dirichlet, Neumann, or Robin-type immersed boundary conditions, is used to enforce the no-slip condition at particle surfaces. A hard-sphere collision model is applied to capture the particle–particle collisions. Proof of concept test cases are presented, showcasing the dynamic grid adaptation and fully resolved two-way coupling between the phases that is possible with this approach. Results for a shock-driven single cylinder under viscous and inviscid conditions are presented along with a demonstration of a shock interacting with a cloud of randomly distributed cylinders and spheres.

42 ENGINEERING↗
DOE OSTI2015

Development of a multiphysics analysis system for sodium-water reaction phenomena in steam generators of sodium-cooled fast reactors

A multiphysics analysis system for sodium-water reaction phenomena in a steam generator of sodium-cooled fast reactors was newly developed. The analysis system consists of the mechanistic numerical analysis codes, SERAPHIM, TACT, and RELAP5. The SERAPHIM code calculates the multicomponent multiphase flow and sodium-water chemical reaction caused by discharging of pressurized water vapor. Applicability of the SERAPHIM code was confirmed through the analyses of the experiment on water vapor discharging in liquid sodium. The TACT code was developed to calculate heat transfer from the reacting jet to the adjacent tube and to predict the tube failure occurrence. The numerical models integrated into the TACT code were verified through some related experiments. The RELAP5 code evaluates thermal hydraulic behavior of water inside the tube. The original heat transfer correlations were corrected for the tube rapidly heated by the reacting jet. The developed system enables evaluation of the wastage environment and the possibility of the failure propagation.

42 ENGINEERING↗
DOE OSTI2016

Resolved-particle simulation by the Physalis method: Enhancements and new capabilities

We present enhancements and new capabilities of the Physalis method for simulating disperse multiphase flows using particle-resolved simulation. The current work enhances the previous method by incorporating a new type of pressure-Poisson solver that couples with a new Physalis particle pressure boundary condition scheme and a new particle interior treatment to significantly improve overall numerical efficiency. Further, we implement a more efficient method of calculating the Physalis scalar products and incorporate short-range particle interaction models. We provide validation and benchmarking for the Physalis method against experiments of a sedimenting particle and of normal wall collisions. We conclude with an illustrative simulation of 2048 particles sedimenting in a duct. In the appendix, we present a complete and self-consistent description of the analytical development and numerical methods.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Radionuclide Transport Simulations Supporting Proposed Borehole Waste Disposal in Israel

A scientific collaboration between the U.S. and Israel is underway to assess the suitability of a potential site for subsurface radioactive waste disposal in the Negev Desert, Israel. The Negev Desert has several favorable attributes for geologic disposal, including an arid climate, a deep vadose zone, interlayered low-permeability lithologies, and carbonate rocks with high uranium-sorption potential. These features may provide a robust natural barrier to radionuclide migration. Geologic and laboratory characterization data from the Negev Desert are incorporated into multiphase flow and transport models, solved using PFLOTRAN, to aid in site characterization and risk analysis that will support decision-making for waste disposal in an intermediate-depth borehole design. The lithology with the greatest uranium sorption potential at the site is phosphorite. We use modeling to evaluate the ability of this layer to impact uranium transport around a proposed disposal borehole. The current objective of the simulations is focused on characterizing hypothetical leakage from waste canisters and subsequent uranium migration under three infiltration scenarios. Here, we describe a hydrogeologic model based on data from a local exploratory borehole and present results for uranium flow and transport simulations under varying infiltration scenarios. We find that under the current climate conditions, it is likely that uranium will remain in the near-field of the borehole for thousands of years. However, under a hypothesized extreme climate scenario representing an increase in infiltration by a factor of 300x above present-day values, uranium may break through the phosphorite layer and exit the base of the model domain (~200 m above the water table) within 1000 years. Simulation results have direct implications for the planning of nuclear waste disposal in the Negev Desert, and specifically in intermediate-depth boreholes.

58 GEOSCIENCES↗
DOE OSTI2025

Accurate models of the added mass force of a uniform random distribution of spherical particles or bubbles

The added mass force resulting from the acceleration of a body in a fluid is of fundamental and practical interest in dispersed multiphase flows. Euler–Lagrange (EL) and Euler–Euler (EE) simulations require closure terms for the added mass force in order to accurately couple the conserved variables between phases. Presently, a more thorough understanding of the added mass force in a multi-particle system is developed based on potential flow resulting in a resistance matrix formulation analogous to Stokesian dynamics. This formulation is then used to generate a dataset of added mass resistance matrices for large systems of randomly generated particles. This methodology is used to create a volume fraction corrected binary model for predicting the added mass force in large systems as well as generate statistics of the added mass force in such systems. This work provides clarification to the theory of the added mass force for particle clouds, and modelling options that may be implemented in existing EL and EE codes.

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

A multiphase three-dimensional multi-relaxation time (MRT) lattice Boltzmann model with surface tension adjustment

The interdependence of surface tension and density ratio is a weakness of pseudo-potential based lattice Boltzmann models (LB). In this paper, we propose a 3D multi-relaxation time (MRT) model for multiphase flows at large density ratios. The proposed model is capable of adjusting the surface tension independently of the density ratio. We also present the 3D macroscopic equations recovered by the proposed forcing scheme. A high order of isotropy for the interaction force is used to reduce the amplitude of spurious currents. The proposed 3D-MRT model is validated by verifying Laplace's law and by analyzing its thermodynamic consistency and the oscillation period of a deformed droplet. The model is then applied to the simulation of the impact of a droplet on a dry surface. Impact dynamics are determined and the maximum spread factor calculated for different Reynolds and Weber numbers. The numerical results are in agreement with data published in the literature. The influence of surface wettability on the spread factor is also investigated. Finally, our 3D-MRT model is applied to the simulation of the impact of a droplet on a wet surface. The propagation of transverse waves is observed on the liquid surface.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Comparison of Surface Tension Generation Methods in Smoothed Particle Hydrodynamics for Dynamic Systems

Developing robust numerical models of dynamic surface tension dominated multiphase systems is an ongoing challenge, especially in scenarios with large density and viscosity ratios. This is critical to the design and understanding of various physical and engineering systems, such as fluidized beds, fuel injectors, and drug delivery schemes. Much of the computational work in surface tension dominated multiphase flows has employed the continuum surface force method (CSF) of Brackbill et al. [1], which recasts surface tension from a surface force to a volumetric force that can be imposed in the vicinity of an interface. The CSF method produces accurate results across a variety of systems, however it relies on the identication of surface normals, which can be unreliable under certain conditions. Alternative methods of simulating surface tension have been proposed. Here the advantages and disadvantages of the CSF method in comparison to a pairwise forces (PF) method proposed by Tartakovsky and Panchenko [2] are explored. The CSF and PF methods are used in a smoothed particle hydrodynamics (SPH) framework to model dynamic systems. Results are compared to existing test cases from the literature and to analytic solutions derived from fundamental normal mode behavior of bubbles and droplets. A physical system for which the PF method is more stable and physically appropriate than the CSF method is identied.

smoothed particle hydrodynamics, multiphase flow↗
DOE OSTI2016

Rise of an argon bubble in liquid steel in the presence of a transverse magnetic field

The rise of gaseous bubbles in viscous liquids is a fundamental problem in fluid physics, and it is also a common phenomenon in many industrial applications such as materials processing, food processing, and fusion reactor cooling. In this work, the motion of a single argon gas bubble rising in quiescent liquid steel under an external magnetic field is studied numerically using a Volume-of-Fluid method. To mitigate spurious velocities normally generated during numerical simulation of multiphase flows with large density differences, an improved algorithm for surface tension modeling, originally proposed by Wang and Tong [“Deformation and oscillations of a single gas bubble rising in a narrow vertical tube,” Int. J. Therm. Sci. 47, 221–228 (2008)] is implemented, validated and used in the present computations. The governing equations are integrated by a second-order space and time accurate numerical scheme, and implemented on multiple Graphics Processing Units with high parallel efficiency. The motion and terminal velocities of the rising bubble under different magnetic fields are compared and a reduction in rise velocity is seen in cases with the magnetic field applied. The shape deformation and the path of the bubble are discussed. An elongation of the bubble along the field direction is seen, and the physics behind these phenomena is discussed. The wake structures behind the bubble are visualized and effects of the magnetic field on the wake structures are presented. A modified drag coefficient is obtained to include the additional resistance force caused by adding a transverse magnetic field.

42 ENGINEERING↗
DOE OSTI2019

Characterizing the Impact of Fractured Caprock Heterogeneity on Supercritical CO 2 Injection

We present a set of multiphase flow simulations where supercritical CO 2 (scCO 2 ) displaces water at hydrostatic conditions within three-dimensional discrete fracture networks that represent paths for potential leakage through caprock above CO 2 storage reservoirs. The simulations are performed to characterize and compare the relative impact of hydraulic and structural heterogeneity in fractured media on the initial movement of scCO 2 through these caprock formations. In one scenario, intrinsic fracture permeabilities are varied stochastically within a fixed network structure. In another scenario, we generate multiple independent, identically distributed network realizations with varying fracture network densities to explore a wide range of geometric and topological configurations. Analysis of the simulations indicates that network structure, specifically connectivity and the presence of hanging fractures, plays a larger role in controlling the displacement of water by scCO 2 than variations in local hydraulic properties. Finally, we identify active surface area of the network as a single-phase feature that could provide a lower bound on the percentage of the network surface area reached by scCO 2 .

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2018

Three-Dimensional Mathematical Model of Oxygen Transport Behavior in Electroslag Remelting Process

A transient three-dimensional model has been proposed to investigate the oxygen transport behavior in electroslag remelting process. The electromagnetism, heat transfer, multiphase flow, and species transport were calculated simultaneously by finite volume method. The volume of fluid approach was adopted to trace the metal–slag–air three-phase flow. Based on the necessary thermodynamics of oxygen transport behavior, a kinetic model was established to predict the mass source terms in species transport equation. The kinetic correction factor was proposed to account for the effect of the oxide scale formed on the electrode on the FeO content in slag. Finally, the effect of applied current on the oxygen transfer was studied. The predicted result agrees well with the measured data when the kinetic correction factor is set to be 0.5. The temperature distribution that affects the thermodynamics differs at the interfaces. The oxygen in air is absorbed into slag due to the oxidation at the slag/air interface. The Fe{sub 2}O{sub 3} in slag and the oxide scale contribute to the increase of FeO content in slag, and the latter one plays the leading role. The oxygen transfer from slag to metal mainly occurs during the formation of the droplet at the slag/metal droplet interface. With the current increasing from 1200 to 1800 A, the oxygen content increases from 76.4 to 89.8 ppm, and then slightly declines to 89.2 ppm when the current increases to 2100 A.

36 MATERIALS SCIENCE↗
DOE OSTI2015

Influence of dynamic factors on nonwetting fluid snap-off in pores

Snap-off is an important dynamic multiphase flow phenomenon which occurs in porous media. It plays a dominant role in the residual trapping and mobilization/immobilization of nonwetting fluids such as hydrocarbons or CO 2 . Current studies, applications, and threshold criteria of snap-off are mostly based on static or equilibrium conditions. Thus, the dynamics of snap-off which is relevant for many real world applications has rarely been systematically studied. While a static criterion indicates the snap-off potential for nonwetting fluids, the competition between the time required for snap-off and the local pore throat capillary number determines whether snap-off actually occurs. Using a theoretical model to couple the wetting film thickness to the local capillary number at the pore throat, we analyzed the dynamics of the wetting/nonwetting interface instability in sinusoidally constricted capillary tubes. The influence of dynamic factors as encapsulated by the effect of local capillary number on nonwetting fluid snap-off time were investigated for varying pore throat to pore body aspect ratio and pore body distances. Finally, the analysis showed that snap-off can be inhibited by a sufficiently large local capillary number even in cases where the static snap-off criterion has been met.

58 GEOSCIENCES↗
DOE OSTI2019

Role of Well Operations and Multiphase Geomechanics in Controlling Fault Stability During CO 2 Storage and Enhanced Oil Recovery

CO 2 injection in active oil fields is a technology proposed for industrializing carbon storage operations. However, monitoring CO 2 migration, oil recovery, and mechanical deformation within caprock and on faults is a challenging problem in large-scale storage-enhanced oil recovery operations due to intersecting faults, multiphase flow effects, and wells with complex production-injection schedule. We develop and demonstrate a methodology based on our coupled modeling framework to monitor the movement of CO 2 , hydrocarbons and water, and the associated evolution in mechanical stability of faults during CO 2 storage-enhanced oil recovery in the Farnsworth Unit oil field in Texas, United States. The methodology honors geological, geophysical, and production-injection data acquired in the field over six decades. Differential depletion from hydraulically isolated fault compartments followed by water and CO 2 injection-induced overpressure causes volumetric contraction and expansion of the reservoir, and changes in the total and effective stresses in the overburden-reservoir-underburden complex. CO 2 migrates upward to accumulate near top of the geologic structure, and water migrates downward to pressurize the faults. Three-dimensional changes in the pressure and stress fields in the system lead to changes in the shear and effective normal tractions on three major faults compartmentalizing the field. Evolution in fault tractions is used to compute the evolution in the Coulomb failure function of the faults to quantify induced-slip tendency under production and injection. We explain the spatial heterogeneity and time variability of Coulomb failure function in terms of well location heterogeneity and well rate variability.

58 GEOSCIENCES↗
DOE OSTI2017

Effect of Slotted Anode on Gas Bubble Behaviors in Aluminum Reduction Cell

In the aluminum reduction cells, gas bubbles are generated at the bottom of the anode which eventually reduces the effective current contact area and the system efficiency. To encourage the removal of gas bubbles, slotted anode has been proposed and increasingly adopted by some industrial aluminum reduction cells. Nonetheless, the exact gas bubble removal mechanisms are yet to be fully understood. A three-dimensional (3D) transient, multiphase flow mathematical model coupled with magnetohydrodynamics has been developed to investigate the effect of slotted anode on the gas bubble movement. The Eulerian volume of fluid approach is applied to track the electrolyte (bath)–molten aluminum (metal) interface. Meanwhile, the Lagrangian discrete particle model is employed to handle the dynamics of gas bubbles with considerations of the buoyancy force, drag force, virtual mass force, and pressure gradient force. The gas bubble coalescence process is also taken into account based on the O’Rourke’s algorithm. The two-way coupling between discrete bubbles and fluids is achieved by the inter-phase momentum exchange. Numerical predictions are validated against the anode current variation in an industrial test. Comparing the results using slotted anode with the traditional one, the time-averaged gas bubble removal rate increases from 36 to 63 pct; confirming that the slotted anode provides more escaping ways and shortens the trajectories for gas bubbles. Furthermore, the slotted anode also reduces gas bubble’s residence time and the probability of coalescence. Moreover, the bubble layer thickness in aluminum cell with slotted anode is reduced about 3.5 mm (17.4 pct), so the resistance can be cut down for the sake of energy saving and the metal surface fluctuation amplitude is significantly reduced for the stable operation due to the slighter perturbation with smaller bubbles.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2018

Burial-driven methane recycling in marine gas hydrate systems

Natural gas hydrate may be buried with sediments until it is no longer stable at a given pressure and temperature, resulting in conversion of hydrate into free gas. This gas may migrate upward and recycle back into the hydrate stability zone to form hydrate. As of yet, however, no quantitative description of the methane recycling process has been developed using multiphase flow simulations to model burial-driven gas hydrate recycling. In this study, we present a series of 1D multiphase transport simulations to investigate the methane recycling process in detail. By invoking the effects of capillary phenomena on hydrate and gas formation in pores of varying size, we find that a free gas phase can migrate a significant distance above the bulk base of hydrate stability. Since the top of the free gas occurrence is often identified as the base of the hydrate stability zone from seismic data, our results demonstrate that not only could this assumption mischaracterize a hydrate system, but that under recycling conditions the highest hydrate saturations can occur beneath the top of the free gas occurrence. We show that the presence of pore size distributions requires a replacement zone through which hydrate saturations progressively decrease with depth and are replaced with free gas. This replacement zone works to buffer against significant gas buildup that could lead to fracturing of overlying sediments. In conclusion, this work provides a framework for simulating flow and transport of methane within the 3-phase stability zone from a mass conservation perspective.

58 GEOSCIENCES↗