SEARCH · Search NASA
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.
pyperc is a Python package for invasion percolation which models multiphase fluid migration in porous media. The model accounts for capillary pressure, buoyancy, and gravity. Regular and unstructured grids can be used. A tuning parameter can be set to adjust the stochastic nature of the process. Graphics depict fluid migration according to the invasion pressure.
TOUGH is a widely used general purpose numerical simulator designed for fluid and heat flows of multiphase, multicomponent fluid mixtures in porous and fractured media, which has been applied to many real-world problems such as underground geological storage, geothermal reservoir engineering or nuclear waste disposal, to name a few. `toughio` is a lightweight, object-oriented and vectorized Python library that aims to provide user-friendly routines to facilitate pre- and post-processing of a TOUGH simulation. Although its primary target is the latest version TOUGH3, `toughio` can also be used to read/write TOUGH2 input/output files. Currently, to the best of our knowledge, only PyTOUGH offers an exhaustive number of features to carry out a complete TOUGH simulation using a scripting language. `toughio` and PyTOUGH share the same objectives, yet with different approaches and philosophies. On the one hand, a PyTOUGH mesh is represented as a MULGRAPH geometry where elements can be unstructured horizontally but only layered vertically. It also checks that the model parameters are valid to ensure that a TOUGH simulation runs smoothly. On the other hand, although it provides basic meshing features, `toughio` mostly relies on common third-party softwares (e.g. Abaqus, FLAC3D, Gmsh, LaGriT) to generate the mesh by importing and converting it to a TOUGH mesh, which also conveniently facilitates the coupling of TOUGH with any other simulator that also supports the same mesh formats. Besides, `toughio` is designed to be just a TOUGH input and output file reader/writer. It assumes that the user knows what he is doing and does not verify the validity of the model parameters (`toughio` does verify the functional arguments).
PFLOTRAN is a GNU LGPL licensed code that simulates single phase water and multiphase (air, co2, water) flow, energy transport, geochemical reaction and biological reaction in the earth's subsurface environment. The problem may be discretized on structured and unstructured grids. The code is designed to run in parallel on supercomputers and leverages the open source PETSc library for solvers and data structures and open source HDF5 for binary file I/O. See www.pflotran.org and/or doc-dev.pflotran.org for further information. SAND2020-12161 M Sandia National Laboratories is a multimission laboratory managed and operated by National Technology & Engineering Solutions of Sandia, LLC, a wholly owned subsidiary of Honeywell International Inc., for the U.S. Department of Energy’s National Nuclear Security Administration under contract DE-NA0003525.
LLEPE (liquid-liquid equilibrium parameter estimator) is a tool for estimating thermodynamic parameters like standard enthalpies, standard molar volumes and Pitzer coefficients, by combining Gibbs minimization with regression to multiphase equilibrium isotherm data. It also includes visualization tools to access data quality and prediction model performance.
Grit is a software library designed specifically for conducting particle-based Lagrangian simulations with CPU/GPUperformance portability. This library allows researchers and engineers to perform a wide range of simulations, including multiphase flow, and more. Grit employs Message Passing Interface (MPI) for distributed memory parallelism and Kokkos programming model for on-node shared memory parallelism with performance portability across different architectures of GPUs and multi-core/manycore CPUs.
PFLOTRAN leverages massively parallel, high performance computing to simulate large-scale non-isothermal multiphase flow, multicomponent reactive transport and electrical resistivity tomography (ERT) problems in the subsurface environment. Researchers have employed PFLOTRAN to simulate these Earth system processes on leadership class supercomputers for over two decades. The code is designed to predict the future estate of environmental systems and better inform stakeholders in the regulatory decision making process (e.g., fate of contaminants, long-term stewardship for nuclear waste, impact of climate change, etc.). A diverse team of scientists oversees PFLOTRAN development and maintenance under an open-source licensing agreement and manages contributions from an international community of researchers.
SAND2022-3463 O This code is intended to accompany the journal manuscript, “Sub-resolution modeling of the apparent mass loss in quantitative broadband X-ray radiography." The manuscript covers in detail how to improve the quantitative mass distribution measurements made for optical diagnostics of multiphase flows. The code contains four separate script files in MATLAB format to support the objective of improving mass measurements. Sandia National Laboratories is a multimission laboratory managed and operated by National Technology & Engineering Solutions of Sandia, LLC, a wholly owned subsidiary of Honeywell International Inc., for the U.S. Department of Energy’s National Nuclear Security Administration under contract DE-NA0003525.
We have studied the hydrocarbon production from oil shale reservoirs filled with diverse initial saturations of fluid phases by implementing numerical simulations of various thermal in-situ upgrading processes. We use our in-house fully functional, fully implicit, and non-isothermal simulator, which describes the in-situ upgrading processes and hydrocarbon recovery by multiphase-multicomponent systems. We have conducted two sets of simulation cases—five-spot well pattern problems and Shell In-situ Conversion Process (ICP) problems. In the five-spot well pattern problems, we have analyzed the effects of initial fluid phase that fills the single-phase reservoir and thermal processes by four cases—electrical heating in the single-phase-aqueous reservoir, electrical heating in the single-phase-gaseous reservoir, hot water injection in the single-phase-aqueous reservoir, and hot CO 2 -injection in the single-phase-gaseous reservoir. In the ICP problems, we have analyzed the effects of initial saturations of fluid phases that fill two-phase-aqueous-and-gaseous reservoir by three cases—initial aqueous phase saturations of 0.16, 0.44, and 0.72. Through the simulation cases, system response and production behavior including temperature profile, kerogen fraction profile, evolution of effective porosity and absolute permeability, phase production, and product selectivity are analyzed. In the five-spot well pattern problems, it is found that the hot water injection in the aqueous phase reservoir shows the highest total hydrocarbon production, but also shows the highest water-oil-mass-ratio. Productions of phases and components show very different behavior in the cases of electrical heating in the aqueous phase reservoir and the gaseous phase reservoir. In the ICP problems, it is found that the speed of kerogen decomposition is almost identical in the cases, but the production behavior of phases and components is very different. It is found that more liquid organic phase has been produced in the case with the higher initial saturation of aqueous phase by the less production of gaseous phase.
The interrogation of materials with X-rays or neutrons to determine structure, energetics, and dynamics is fundamental to advancing physical and chemical materials science and enabling innovative material technologies. A persistent challenge in materials development is that progress depends on understanding structure and dynamics across multiple length and time scales in increasingly complex, multicomponent systems featuring interfaces, heterogeneity, and hierarchical organization. Despite rapidly growing demands on materials characterization, current experimental approaches are almost exclusively based on isolated X-ray or neutron scattering and spectroscopy, reflecting a paradigm largely unchanged for decades. To assess the scientific need for a new experimental paradigm, a 3-day workshop sponsored by the U.S. National Science Foundation (NSF) was held at the SpringHill Suites, San Jose, California, from June 2 to 4, 2022. The workshop brought together 70 national and international experts who critically evaluated opportunities enabled by concurrent neutron and X-ray (NeX) scattering, spectroscopy, and imaging experiments. The participants reached a clear consensus that establishing NeX capabilities is crucial for advancing the science of complex materials in the United States. This report illustrates the scientific drivers for NeX experiments through representative examples spanning biomaterials, energy materials, soft matter, nanomaterials, quantum materials, geoscience, and applied materials research. The complementarity of neutrons and X-rays is essential for robust model development and refinement, particularly in multiphase and multicomponent systems. While joint refinement of data from separate experiments is valuable, concurrent measurements uniquely eliminate uncertainties arising from sample evolution, environmental drift, and irreproducibility associated with experiments performed at different locations and times. Realizing NeX capabilities will require the development of new instrumentation, data analysis frameworks, and robust sample environments compatible with both neutron and X-ray probes. Addressing these challenges will enable unambiguous interpretation of complex materials behavior and open new frontiers in materials research.
Block-structured adaptive mesh refinement (AMR) provides the basis for the temporal and spatial discretization strategy for a number of Exascale Computing Project applications in the areas of accelerator design, additive manufacturing, astrophysics, combustion, cosmology, multiphase flow, and wind plant modeling. AMReX is a software framework that provides a unified infrastructure with the functionality needed for these and other AMR applications to be able to effectively and efficiently utilize machines from laptops to exascale architectures. AMR reduces the computational cost and memory footprint compared to a uniform mesh while preserving accurate descriptions of different physical processes in complex multiphysics algorithms. AMReX supports algorithms that solve systems of partial differential equations in simple or complex geometries and those that use particles and/or particle–mesh operations to represent component physical processes. In this article, we will discuss the core elements of the AMReX framework such as data containers and iterators as well as several specialized operations to meet the needs of the application projects. In addition, we will highlight the strategy that the AMReX team is pursuing to achieve highly performant code across a range of accelerator-based architectures for a variety of different 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.
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.
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.
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.
Structural rearrangements of HIV-1 glycoprotein Env promote viral entry through membrane fusion. Env is a symmetric homotrimer with each protomer composed of surface subunit gp120 and transmembrane subunit gp41. Cellular CD4- and chemokine receptor-binding to gp120 coordinate conformational changes in gp41, first to an extended prehairpin intermediate (PHI) and, ultimately, into a fusogenic trimer-of-hairpins (TOH). HIV-1 fusion inhibitors target gp41 in the PHI and block TOH formation. To characterize structural transformations into and through the PHI, we employed asymmetric Env trimers containing both high and low affinity binding sites for individual fusion inhibitors. Asymmetry was achieved using engineered Env heterotrimers composed of protomers deficient in either CD4- or chemokine receptor-binding. Linking receptor engagement to inhibitor affinity allowed us to assess conformational changes of individual Env protomers in the context of a functioning trimer. We found that the transition into the PHI could occur symmetrically or asymmetrically depending on the stoichiometry of CD4 binding. Sequential engagement of multiple CD4s promoted progressive exposure of individual fusion inhibitor binding sites in a CD4-dependent fashion. By contrast, engagement of only a single CD4 molecule led to a delayed, but symmetric, exposure of the gp41 trimer. This complex coupling between Env-CD4 interaction and gp41 exposure explained the multiphasic fusion-inhibitor titration observed for a mutant Env homotrimer with a naturally asymmetric gp41. Our results suggest that the spatial and temporal exposure of gp41 can proceed in a nonconcerted, asymmetric manner depending on the number of CD4s that engage the Env trimer. The findings have important implications for the mechanism of viral membrane fusion and the development of vaccine candidates designed to elicit neutralizing antibodies targeting gp41 in the PHI.
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.
Although bedrock weathering strongly influences water quality and global carbon and nitrogen budgets, the weathering depths and rates within subsurface are not well understood nor predictable. Determination of both porewater chemistry and subsurface water flow are needed in order to develop more complete understanding and obtain weathering rates. In a long-term field study, we applied a multiphase approach along a mountainous watershed hillslope transect underlain by marine shale. Here we report three findings. First, the deepest extent of the water table determines the weathering front, and the range of annually water table oscillations determines the thickness of the weathering zone. Below the lowest water table, permanently water-saturated bedrock remains reducing, preventing deeper pyrite oxidation. Secondly, carbonate minerals and potentially rock organic matter share the same weathering front depth with pyrite, contrary to models where weathering fronts are stratified. Thirdly, the measurements-based weathering rates from subsurface shale are high, amounting to base cation exports of about 70 kmolc ha−1 y−1, yet consistent with weathering of marine shale. Finally, by integrating geochemical and hydrological data we present a new conceptual model that can be applied in other settings to predict weathering and water quality responses to climate change.