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 775 records · Page 43

Realistic fuel spray modeling for gasoline direct injection engine applications

Fuel spray modeling plays a critical role during modern gasoline direct injection (GDI) engine development due to fuel injection’s dominant impact on engine performance and emissions as well as the complex physical processes involved. In engineering three-dimensional (3D) computational fluid dynamics (CFD) simulations, the liquid-phase fuel atomization, evaporation, and mixing are usually modeled with the discrete droplet model (DDM) adopting a Lagrangian approach for multiphase CFD simulations. To this end, general practices heavily depend on the reduced order characterization of the injector nozzle flow. However, such simplified injector modeling may lead to insufficient representations of the complex spray dynamics. To tackle this problem, this study proposes a novel workflow to numerically evaluate GDI sub-cooled and flash-boiling sprays under engine-relevant conditions using a side-mounted GDI injector together with real gasoline fuel properties. The workflow introduces a one-way coupling (OWC) method leveraging high-fidelity nozzle flow simulations to provide realistic boundary conditions to the Lagrangian injector model. The proposed workflow was first verified in a constant volume chamber (CVC) environment and then implemented in a practical GDI engine setup to study spray morphology, fuel-air mixing, and wall-wetting propensity. In addition, detailed comparison was performed between the OWC method and the conventional rate of injection (ROI) routine. In conclusion, quantitative analysis of spray characteristics was conducted to highlight possible source of discrepancies of the conventional ROI method.

33 ADVANCED PROPULSION SYSTEMS↗

Solving multiphysics-based inverse problems with learned surrogates and constraints

Abstract Solving multiphysics-based inverse problems for geological carbon storage monitoring can be challenging when multimodal time-lapse data are expensive to collect and costly to simulate numerically. We overcome these challenges by combining computationally cheap learned surrogates with learned constraints. Not only does this combination lead to vastly improved inversions for the important fluid-flow property, permeability, it also provides a natural platform for inverting multimodal data including well measurements and active-source time-lapse seismic data. By adding a learned constraint, we arrive at a computationally feasible inversion approach that remains accurate. This is accomplished by including a trained deep neural network, known as a normalizing flow, which forces the model iterates to remain in-distribution, thereby safeguarding the accuracy of trained Fourier neural operators that act as surrogates for the computationally expensive multiphase flow simulations involving partial differential equation solves. By means of carefully selected experiments, centered around the problem of geological carbon storage, we demonstrate the efficacy of the proposed constrained optimization method on two different data modalities, namely time-lapse well and time-lapse seismic data. While permeability inversions from both these two modalities have their pluses and minuses, their joint inversion benefits from either, yielding valuable superior permeability inversions and CO 2 plume predictions near, and far away, from the monitoring wells.

Yin, Ziyi (ORCID:0000000250248771)↗

Probing interfacial momentum closures in two-phase bubbly flow with machine learning-aided methods

Computational fluid dynamics (CFD) approach has already reached a high level of maturity for single-phase flows, however the development of closure models for two-phase flow requires additional attention. Multiphase CFD (M-CFD) methods resolve the conservation equations for mass, momentum and energy while differing in the approaches and strategies adopted in the physical closure models. The most widely adopted framework for M-CFD is the Eulerian-Eulerian two-fluid approach which assumes that all phases are co-existing inside each computational cell. For each fluid, the full set of conservation equations is solved; therefore, each fluid has a different velocity field. For adiabatic two-phase flow, the mechanisms of the interfacial momentum transfer are modeled by the interfacial forces representing different physical mechanisms. One of the crucial issues in the development and application of two-fluid model is the understanding of the interfacial momentum closures which determines the bubble distribution and migration behaviors. Dedicated experiments are performed to support the physical understanding and drive the closures’ development. However, limitations exist due to the uncertainties in the experimental measurement and the simplified analytical assumptions which have difficulties on representing the complex non-linear flow fields. In this paper, a data-driven approach, Feature Similarity Measurement (FSM), is developed and proposed to resolve the challenges of modeling the interfacial forces closures. Case study is performed with two-phase flow scenarios where the high-fidelity experimental data is available. Within the Eulerian-Eulerian two-fluid framework, only momentum equations for gas and liquid phases are solved and reduced-order interfacial momentum closures are aided with FSM. Predictions of void fraction and velocity fields are analyzed and demonstrate the potential of machine learning-driven interfacial forces closures.

97 MATHEMATICS AND COMPUTING↗

Evaluation of liquid-phase thermometry in impinging jet sprays using synchrotron x-ray scattering

Liquid thermometry during primary and secondary breakup of liquid sprays is challenging due to the presence of highly dynamic, optically complex flow features. This work evaluates the use of x-ray scattering from a focused, monochromatic beam of the Advanced Photon Source at Argonne National Laboratory for the measurement of liquid temperatures within the mixing zone of an impinging jet spray. The measured scattering profiles are converted to temperature through a previously developed two-component partial least squares (PLS) regression model. Transmitive mixing during jet merging is inferred through spatial mapping of temperatures within the impingement region. In this work, the technique exhibits uncertainties of ±2K in temperature and 2% in capturing the correct scattering profile, showing its potential utility for probing liquid temperature distributions in multiphase flows.

42 ENGINEERING↗

Challenges in correlating oxygen stable isotope ratios of hydrates on uranium ore concentrates to process waters

Exchange of oxygen stable isotopes (δ 18 O values) between precipitation waters and uranium oxides is governed by thermodynamics or kinetics. It has been assumed that meteoric waters can be related to precipitation waters in uranium ore concentrates and their calcination and reduced uranium oxide products. With this assumption, the δ 18 O values of uranium materials could provide forensic signatures that identify the production history and geolocation of nuclear materials. To further exploit the potential of δ 18 O values in nuclear material analysis, this study examines the oxygen stable isotope exchange in two UOCs, magnesium diuranate (MDU) and sodium diuranate (SDU). MDU and SDU were synthesized from solutions of uranyl nitrate hexahydrate using precipitation waters with unique oxygen isotope compositions. The structures of the MDU and SDU were analyzed using powder X-ray diffraction (p-XRD) and thermal mass loss curves, while the δ 18 O values of waters generated during thermal decomposition were analyzed using a thermogravimetric analyzer coupled to an isotope ratio infrared spectrometer (TGA-IRIS). By p-XRD, the MDU was uniform and amorphous across all syntheses with residual crystalline material incorporated as a minor component. Combined with the TGA results, all of the MDU is likely amorphous MgU 2 O 7 ·3H 2 O with MgO impurities present throughout. In contrast, the SDU synthesis resulted in multiple phases with many samples exhibiting crystalline phases including a combination of Na(UO 2 ) 4 O 2 (OH) 5 ·5H 2 O and Na 2 (UO2) 6 O 4 (OH) 6 ·8H 2 O with a Na 2 U 2 O 7 minor phase. A small fraction of the SDU samples were amorphous with no crystalline XRD peaks observed. Mass loss curves of the SDU samples revealed that the amorphous samples contained inclusions of similar crystalline phases compared to the crystalline materials. The uniformity of the MDU samples enabled highly reproducible measurements of δ 18 O values of the water vapor yielded for two dehydration events at 170 °C and 500 °C. In contrast, the multiphase composition of the SDU samples resulted in poor reproducibility in δ 18 O values. In conclusion, neither system revealed any correlation between the δ 18 O values of precipitation water, and the waters released during dehydration of the UOCs.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗

Microstructural characterization of phase-separated co-deposited Cu–Ta immiscible alloy thin films

Elevated temperature co-sputtering of immiscible elements ends in a variety of self-organized morphologies due to phase separation. Cu–Ta is used as a model system to understand the evolution of phase-separated microstructural morphologies by co-sputtering thin films with nominal 50–50 at.% composition at four temperatures: 25, 400, 600, and 800 °C. Scanning/transmission electron microscopy of the film cross sections showed the microstructure morphology varied from nanocrystalline Cu–Ta at 25 °C to a wavy ribbon-like structure at 400 °C, to Cu-rich agglomerates surrounded by Ta-rich veins at 600 and 800 °C. In the agglomerate-vein morphology, microstructural features were present on two length scales, from a few nanometers to a few tens of nanometers, thus making the structures hierarchical. On the nanoscale, the Cu-rich agglomerates contained Ta precipitates, whereas the Ta-rich veins had embedded Cu nanocrystals. The various microstructures can be attributed to the highly disparate constituent element interdiffusion at the deposition temperatures with the Cu having orders of magnitude higher mobility than Ta at the deposition temperatures. This research of processing–microstructure relationship will be useful in guiding the design of hierarchical multiphase microstructures in binary or multicomponent thin films with tailored mechanical properties.

36 MATERIALS SCIENCE↗

NUMERICAL INVESTIGATION OF AIR DEHUMIDIFICATION THROUGH WATER DROPLETS DIELECTROPHORESIS

This paper numerically investigated the dehumidification of moist airflow in a converging rectangular duct with electrically enhanced condensation. The aim was to develop a condensation model that predicted water vapor condensation by employing the dielectrophoresis principle. The charged droplets were injected into the computational domain to produce a gradient electric field. The dielectrophoretic interaction between the droplets and vapor molecules of the humid air resulted in local water vapor condensation at the vapor and liquid droplet interface. This phenomenon is described in detail in this paper. A hybrid Eulerian-Lagrangian solver SprayFoam was developed based on OpenFOAM® to simulate condensation in two-phase gas-liquid mixtures. The following developments were made based on the standard compressible multiphase flow solver SprayFoam in OpenFOAM®: (1) Eulerian solver for gas phase (2) multi-component species transport, and (3) Lagrangian solver for gas-droplet two-phase flows and sub-models for liquid droplets. The newly developed numerical model was experimentally validated with data from a series of tests conducted in the authors' laboratory. The results showed that the simulations followed the same trends as the data, and the model predicted the condensation due to the electro-spray injection process in the air stream. The injected water droplets, which are electrically charged to their Rayleigh limits, increased in size while wiping out the humidity from the air. The simulations indicated that the dehumidification was about 1% for 0.5 cubic feet per minute airflow rate. Scaling up to larger flows is a future follow-up work.

Yel Mahi, *Maliha↗

Supo Thermal Model Development IV

This report describes the continuation of the Computational Fluid Dynamics (CFD) model of the Supo cooling system described in the reports, Supo Thermal Model Development, II, and III. In the previous reports, the solution natural convection of Supo was simulated using a 2D, axisymmetric model that altered the actual geometry. A 3D geometry of Supo will more accurately estimate the natural convection heat transfer coefficient (HTC) of the system since Supo contains several non-symmetric internal components that can only be modeled in 3D. This, in turn, will improve the natural convection HTC results for comparison purposes with past operational data and previous 2D CFD results. Unfortunately, the complexities of performing multiphase 3D CFD simulations has proven too large of a task for our current capabilities. Instead, this report will focus on the 3D model development for the 2.39 kW/L (31.3 kW) case and limitations and challenges associated with the 3D simulations.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Stress from long-range interactions in particulate system

In particulate systems, the ensemble average accounting for effects from all particles is expressed as the expected value related to the nearest particle probability. Using this expression, a stress is defined even for long-range particle interactions without the divergence difficulty because of the rapid far field decay of the probability. Similar to the potential part of the virial stress in molecular systems, this stress represents particle-particle interactions through the surrounding field. As examples, the stresses are calculated for systems of charged particles and disperse multiphase flows.

97 MATHEMATICS AND COMPUTING↗

Verification and Validation of MFiX-PIC

This work is aimed at verification and validation (V&V) of particle-in-cell (PIC) approach in the National Energy Technology Laboratory’s (NETL) Multiphase Flow with Interphase eXchanges (MFiX-PIC) which was officially released in April 2019. Under verification, the results from MFiX-PIC are compared with exact or analytical solutions, while experimental results are used for benchmarking under validation. The test suite consists of the following cases: Terminal velocity. Particle advection in time-varying flow field. Particle settling in a dense medium. Evaporation. Rayleigh-Taylor instability. Minimum fluidization. Mini-circulating fluidized bed. The cases reported in a format consistent with the existing Verification and Validation Manual (Musser et al., 2018) explore a wide range of flow velocities and concentration regimes. Only one case considered in this work verifies the implementation of heat and mass transfer routines in MFiX-PIC, while a few other studies are currently underway to further expand the validation space in reacting flows. This study highlights the need for a rigorous analysis to objectively quantify the effect of model parameters. In addition, there is a scarcity of high-quality experimental data in the literature suitable for validating PIC methodology. Future efforts would include performing experiments which provide results for a valuable database. While designing such experiments, particle count should be significant to have sufficient particles per parcel while maintaining an optimal parcel to cell volume ratio. Also, datasets should have objectively assessed experimental uncertainty whereby confidence interval in those measurements become more explicit.

42 ENGINEERING↗

Aftershocks of a Chemical Explosion in Granite from the Source Physics Experiment Phase I

The Source Physics Experiment (SPE) is a multiphase experiment to better understand explosion source physics and thereby improve explosion monitoring. The rst phase (Phase I) examined the explosion source in a hard-rock medium and took place in a granite outcrop called the Climax Stock section of the Nevada National Security Site (NNSS). This series of chemical explosions are called the Wet Granite Geology (WGG) events, but are usually referred to simply as the SPE events (i.e., SPE-1, SPE-2, etc.). Event information for the chemical explosions is given in Table 1. The study of explosion aftershocks can aid in understanding of in-situ stresses and be used as a discriminant for earthquakes [Ford and Walter, 2010]. Sweeney and Harben reported no observed aftershocks for SPE-1 and -2, respectively. They calculated a minimum magnitude of observation based on the noise oor of the seismic array and predicted at least one or two events in the week after the events using the hard-rock aftershock model from Ford and Labak. Possible explanations for the null observation ranged from instrumental - poor seismometer coupling, to physical - the shallow depth of burial was unable to access the deep tectonic stress that causes some explosion aftershocks. We extend the analysis of Sweeney and Harben to the largest of the SPE Phase I chemical explosions, SPE-5, which was approximately a 5 tonne TNT-equivalent shot at 76.5 m depth. This event registered as an M L 2 on the University of Nevada, Reno seismic network and should produce a measurable aftershock sequence.

58 GEOSCIENCES↗

UNESE Data Analysis - Disko Elm Gas Transport Characteristics (NA-22 Quarterly Report)

Two research papers describing gas transport studies in subsurface containment environments appropriate for underground nuclear explosions were revised and accepted for publication in the peer-reviewed journals, Nature – Scientific Reports and Journal of Geophysical Research. The Nature – Scientific Reports paper is a primary deliverable and describes the gas-transport experiment performed in P Tunnel (NNSS) and 1) subsequent analysis of field data as part of the UNESE program. The Journal of Geophysical Research paper describes an experiment, supported mainly by DARPA and led by researchers at Weston Geophysical and New Mexico Tech., with 2) analyses supported by this project to evaluate the ability of SF6 gas tracer to track xenon gas migration under conditions related to a previous UNESE experiment. Two new tasks consider 3) evaluation of detonation heating and multiphase flow on arrival of xenon at the surface and 4) the impact of surface or Langmuir adsorption of gases being transported through pores and along fractures.

58 GEOSCIENCES↗

UNESE Data Analysis-Disko Elm Gas Transport Characteristics (NA-22 Quarterly Report)

Recently published papers have demonstrated the impact of sorption mechanisms on gas transport in the subsurface following an underground nuclear explosion. To include sorption processes in the multi-physics models of noble gas production and transport, researchers have conducted experiments of gas sorption in geologic materials and calibrated a linear (Henry’s law) sorption model for all gases in various geologic materials. As a part of Task #4, we requested the experimental data from Sandia Laboratories and derived a dual-mode (Henry’s absorption and Langmuir adsorption) sorption model. The derived model better describes the nonlinear sorption mechanism with much higher fitness with the experimental data. When gas concentrations are relatively low (in the conditions of field experiments or underground nuclear explosions) compared to those used sorption experiments, the nonlinear sorption model can be approximated as a linear sorption model. The equivalent Henry’s constant, which equals the product of Langmuir capacity and affinity, is systematically higher than the literature values. Further analyses and comparison of sorption models are planned in FY2020 Q4 and the sorption model will be coupled in models of nonisothermal multiphase transport for studying noble gas detectability.

58 GEOSCIENCES↗

World Congress on High Entropy Alloys (HEA 2019) (Final Technical Report)

The World Congress on High Entropy Alloys (HEA 2019) was a new, cross-disciplinary technical forum designed to share the latest research advances in materials with high configurational entropy, including both single-phase and multiphase (compositionally complex) alloys. This inaugural World Congress on High Entropy Alloys (HEA 2019) was held in Seattle, Washington from November 17-20, 2019 with the goal of HEA 2019 to assemble the international community of HEA stakeholders – including researchers, educators, and engineers – to capture the current state-of-the-art, identify current gaps in HEA efforts, and to set the stage for the future growth and application of HEAs worldwide. HEA 2019 was the environment for leaders in this field to showcase their researc, and network with their peers to informally discuss the latest developments in high entropy alloys, while also making valuable professional connections. The goal of the organizers was to ensure a high-quality conference that held a strong international presence and aimed to highlight the important advances in HEA research on a global scale. Student participation was greatly encouraged in order to give students exposure to recent advances in this field. This meeting was open to all interested individuals and attempted to provide the most advanced technical meeting possible, with a mixture of attendees from research and development, industry, and education.

36 MATERIALS SCIENCE↗

Computational Fluid Dynamics Analysis of a 12 MW Circulating Fluidized Bed Rise

A three-dimensional (3D) two fluid model (TFM) model of the 12 MW circulating fluidized bed (CFB) riser at Chalmers University was developed using the open-source Multiphase Flow with Interphase eXchanges (MFiX) Software Suite. The hydrodynamic behavior of the riser was validated by comparing axial pressure profiles from experimental measurements with the results from the simulations. Under some of the operating conditions where there was no external circulation of solids, the pressure profiles from the simulation matched well with the experimental measurements. For the cases where there was external circulation of solids, the pressure profiles with the monodisperse solid phase modeling underpredicted the solids concentration in the freeboard region. With the inclusion of polydispersity in the model, the predicted pressure profiles agreed well with the experimental measurements. A two-dimensional (2D) model of the riser was used to perform simulations of coal combustion. The predicted flue gas outlet compositions of oxygen and carbon dioxide compared reasonably to experiments. The experiments also reported some concentration of carbon monoxide which was overpredicted in the simulations. Additional experimental data such as axial oxygen concentration profile and pressure profiles were also compared with the experiment. It was observed that the oxygen concentration in the bottom bed from the simulations was higher than in the experiments, which suggested that most of the oxidation reactions occurred in the bottom bed. The concentration of volatiles and char in the simulations was lower in the bottom bed than in the experiments, which suggest that there is not sufficient penetration of the solid phase species into the bed.

20 FOSSIL-FUELED POWER PLANTS↗

SFAI and Related Technologies

Researchers at Los Alamos have developed Swept-Frequency Acoustic Interferometry (SFAI) a technology that enables businesses to have greater control over their industrial processes, allowing for significantly improved efficiency by giving operators vital composition and flow rate information in real-time. This technology fills the information gap by providing this data non-invasively and at a fraction of the cost of traditional composition meters. SFAI uses engineered ultrasound and advanced signal processing to characterize multiphase flow (e.g., liquids, gasses, multi-phase mixtures) and other flows inside sealed containers and flow pipes. By being both low-cost and non-invasive, this technology is set to disrupt the current oil field monitoring state-of-the-art by providing individual well-level composition data, compared to the multi-well macroscopic information that is currently available. We are seeking a commercialization partner to license the technology or fund us through a Cooperative Research and Development Agreement (CRADA) to a further develop a usable product for commercial purposes.

03 NATURAL GAS↗

Fast Computational Algorithms for Partial Differential Equations and Uncertainty Quantifications

This project concerned the construction, testing and analysis of computational algorithms for solving parameterized and stochastic partial differential equations. The study and understanding of equations of this type is of fundamental importance in numerous engineering and scientific applications. Examples include simulation of plasma dynamics in models of electric propulsion and nuclear fusion, simulation of multiphase flows, such as the flow of water, gas and oil in reservoirs, and structural analysis of the dependence of structures on materials. Parametrization is used in such settings when properties of the models such as viscosity of fluids or electric resistivity of materials are not precisely understood and instead are treated as random variables. The resulting solutions are themselves random, and having such solutions will enable engineers to use probabilistic methods to assess the likelihood of events, for example, whether a pollutant in a liquid will exceed a limit, and to use such analyses to develop ways to ensure positive outcomes. Construction of accurate (high resolution) computational solutions is expensive, requiring significant computer time and computational resources, and there is need to reduce computational cost to make simulation useful and effective. The aim of the project was to construct algorithms to efficiently compute surrogate solutions to parameterized problems to allow for efficient and accurate simulation. The technical approach used focused on two related strategies, based on rank-reduction methods and reduced-order models. These methods construct surrogate solutions of parameter-dependent models by projection or interpolation into low-dimensional approximation spaces. Cost savings are achieved if the low-dimensional spaces can be identified and constructed efficiently and if the resulting low-dimensional algebraic systems can be solved cheaply. Accomplishments include: Theoretical and empirical demonstration of the effectiveness of fast multigrid solution strategies for computing low-rank representations of parameter-dependent solutions to discrete partial differential equations, including the first proof establishing so-called textbook convergence properties for low-rank methods. Development of efficient solution algorithms for solving nonlinear parameter-dependent partial differential equations used in models of fluid dynamics. Developent of efficient algorithms for low-rank representation of solutions of time-dependent simulations of fluid dynamics using multi-dimensional tensor representations of solutions.

97 MATHEMATICS AND COMPUTING↗

Deterministic Calibration of MFiX-PIC, Part 1: Settling Bed

The Particle-in-cell (PIC) numerical approach for modeling granular solids in fluid flow has gained significant interest in recent years. Valued for its often shorter time-to-solution, the PIC formulation relies on modeling statistical groupings of particles called parcels in cooperation with a solids stress model to affect local solids velocity. This is in contrast to the discrete element model (DEM) where every particle in a system is modelled individually and directly coupled to local solids velocity through Newtonian mechanics. The U.S. Department of Energy (DOE), National Energy Technology Laboratory (NETL) develops and maintains Multiphase Flow with Interphase eXchanges (MFiX), a collection of open-source computational fluid dynamics (CFD) solvers. Included in the MFiX suite are traditional two-fluid model (TFM) and DEM solvers, and a recently added PIC solver (NETL, 2021). In general, PIC methodologies offer an accuracy trade-off in lieu of computational speed; and therefore, it is important to assess the credibility of MFiX-PIC simulations. For this purpose, a systematic verification, validation and uncertainty quantification (VVUQ) effort was initiated at NETL to assess the new PIC solver

42 ENGINEERING↗