Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “continuous flow”

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 217 records · Page 12
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Best practices in software development for robust and reproducible geoscientific models based on insights from the Global Carbon Budget's dynamic vegetation models

Computational models play an increasingly vital role in scientific research by enabling the numerical simulation of complex processes. Such models are also fundamental in geosciences. For instance, they offer critical insights into the impacts of global change on the Earth system today and in the future. Beyond their value as research tools, models are also software products and should therefore adhere to certain established software engineering standards. However, scientists are rarely trained as software developers, which can lead to potential deficiencies in software quality like unreadable, inefficient, or erroneous code. The complexity of models, coupled with their integration into broader workflows, also often makes it challenging to reproduce results, evaluate processes, and build upon them. In this paper, we review the state and current practices of the development processes of the state-of-the-art land surface models used by the Global Carbon Budget. We combine the experience of modelers from the respective research groups with the expertise of software engineers from tech companies to outline key principles and tools for improving software quality in research. We explore four main areas: (1) model testing and validation, (2) scientific, technical, and user documentation, (3) version control, continuous integration, and code review, and (4) the portability and reproducibility of workflows. Our review reveals that while modeling communities are incorporating many best practices, significant room for improvement remains in areas such as automated testing, automated documentation, and reproducibility. Therefore, we here identify and promote essential software engineering practices, including numerous examples of practices from within the community that can serve as guidelines for other models and could help streamline processes across the entire community. We conclude with an open-source example implementation of these principles, demonstrating portable and reproducible data flows, a continuous integration setup, and web-based visualizations. This example may serve as a practical resource for model developers, users, and all scientists engaged in scientific programming.

Gregor, Konstantin [Technical Univ. of Munich (Ger↗
DOE OSTI2020

Retention of Long-period Gas Giant Planets: Type II Migration Revisited

During their formation, emerging protoplanets tidally interact with their natal disks. Proto–gas giant planets, with Hill radii larger than the disk thickness, open gaps and quench gas flow in the vicinity of their orbits. It is usually assumed that their type II migration is coupled to the viscous evolution of the disk. Although this hypothesis provides an explanation for the origin of close-in planets, it also encounters a predicament on the retention of long-period orbits for most gas giant planets. Moreover, numerical simulations indicate that the planets' migrations are not solely determined by the viscous diffusion of their natal disk. In this work, we carry out a series of hydrodynamic simulations combined with analytic studies to examine the transition between different paradigms of type II migration. We find a range of planetary mass for which gas continues to flow through a severely depleted gap so that the surface density distribution in the disk region beyond the gap is maintained in a quasi-steady state. The associated gap profile modifies the location of corotation and Lindblad resonances. In the proximity of the planet's orbit, high-order Lindblad and corotation torque are weakened by the gas depletion in the gap, while low-order Lindblad torques near the gap walls preserve their magnitude. Consequently, the intrinsic surface density distribution of the disk delicately determines both the pace and direction of the planets' type II migration. We show that this effect might stall the inward migration of giant planets and preserve them in disk regions where the surface density is steep.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Efficient many-jet event generation with flow matching

We apply for the first time, to the best of our knowledge, the flow matching method to the problem of phase-space sampling for event generation in high-energy collider physics. By training the model to remap the random numbers used to generate the momenta and helicities of the scattering matrix elements as implemented in the portable partonic event generator pepper, we find substantial efficiency improvements in the studied processes. We focus our study on the highest final-state multiplicities in Drell-Yan and top-antitop pair production used in simulated samples for the Large Hadron Collider, which computationally are the most relevant ones. We find that the unweighting efficiencies improve by factors of 184 and 25, respectively, when compared to the standard approach of using a vegas-based optimization. We also compare continuous normalizing flows trained with flow matching against the previously studied normalizing flows based on coupling layers and find that the former leads to better results, faster training and a better scaling behavior across the studied multiplicity range, while the latter evaluate faster. When combining the advantages of both methods using the regflow approach, we find parton-level unweighted event generation walltime gains of about a factor of 10 at the highest final-state multiplicities.

Bothmann, E. [CERN; Gottingen U.] (ORCID:000000016↗
DOE OSTI2026

Refining a Novel Process Monitoring Method to Safeguard Continuously Cycling Designs Using Isotopic Ratios

In advanced reactor (AR) designs, a common feature is continuous chemical processing and circulation of the nuclear material. This work bridges a significant measurement gap in safeguarding reactors with circulating fuel or continuous refueling by leveraging and building on the isotope ratio method first developed by our team under an FY21 Advanced Reactors International Safeguards Engagement (ARISE) project (Uribe et al. 2021). In circulating fuel designs, the radioisotope inventory changes from traditional effects (e.g., radioactive decay, fission) but also includes material transport due to pressure and temperature gradients. Such designs may also require regular or continuous additions or removals during operation, which significantly increases the rate of inventory change compared to traditional pressurized water reactor (PWR)s. Thus, directly tracking the nuclear inventory is ineffective since the isotopes are continuously added and removed. The isotope ratio method instead focuses on detecting changes to the input and output flows of radioisotopes. Previous work showed that for well-chosen pairs of isotopes, the isotopic ratio provides a sensitive and lasting indicator of deviation from normal conditions (e.g., startup, shutdown, diversion). The isotope ratio method is a process monitoring method with potential for application in for forward-looking approaches to International Atomic Energy Agency (IAEA) safeguards. The original process monitoring method was developed for a specific case—the decay tank of a thorium-fueled molten salt breeder reactor. In this expanded work, we explored other types of reactors and processes with nonstationary (e.g., flowing) nuclear material, which are difficult to safeguard with traditional methods because of the transient nature of the systems. The goal of the present work is to generalize the isotope ratio method for use in processes with continuously flowing nuclear material. All continuous processes have an average time for isotopes to be replaced in the system. The isotope ratio method works by choosing radioisotopes with half-lives both above and below the processing time. The present work seeks to explore which isotopes are suitable for the method by simulating the nuclear inventory, radioisotope emissions, and detector responses for several classes of advanced reactors. While the method can in principle be applied to other processes (e.g., enrichment or reprocessing facilities), the present work limits scope to ARs with continuously flowing fuel. Section 2 details the mathematics supporting the isotope ratio method, and Section 3 introduces the representative ARs selected for this work. Section 4 discusses how each reactor was analyzed, and Section 5 showcases the results for each representative reactor. Finally, Section 6 provides concluding remarks and suggests pathways for further analysis.

98 NUCLEAR DISARMAMENT, SAFEGUARDS, AND PHYSICAL P↗
DOE OSTI2020

Results of Multiple Tracer Injections into Fractures in the EGS Collab Testbed-1

The EGS Collab project constructed an intermediate scale (~10-20 m) testbed at the 4850 level of the Stanford Underground Research Facility (SURF) in South Dakota for testing and validating fracture stimulation and flow/transport models. This testbed consists of eight ~200 ft (~60 m) HQ-diameter (9.6 cm) boreholes that are drilled into the crystalline rocks of the Poorman Formation from the West Access Drift tunnel. Of the eight boreholes, one borehole is used as an injection/stimulation well, while another sub-parallel borehole located about 10 m away from the injection well is used as a production well, and rest of the other boreholes are used as geophysical/fluid sampling monitoring wells. Hydraulic stimulation activities were conducted at three locations along the injection hole in an attempt to create direct fracture connections to the production hole. A flow system has been established between injection and production boreholes through a set of hydraulically stimulated fractures propagated from a notch located at 164 ft in the injection hole. Although we planned for a single production well, a set of natural fractures in the testbed is believed to have intersected the stimulated hydraulic fractures and provided additional flow paths for water to be transported to the drift through multiple monitoring boreholes and weep zones. As the flow tests continued after stimulation activities, single or a combination of two or three producers become dominant producers at different times, mostly as a response to intersection of hydraulic fracture and testbed wells, activation of natural fractures and making new leak points to the monitoring wells, resealing of leaky wells, and so on. Since late October 2018, multiple tracers were injected into the fracture system at the 164 ft location that involves both stimulated and natural fractures, and tracers were recovered from multiple locations in nearby wells and weep. The cumulative water recoveries over the time have ranged from 50 to 90%; however, the injected tracer recoveries from various production sources are much less (ranging from a few percentages to 38%). Changes in water/tracer recoveries, shifting of major producing wells from one well to the others, and other observations (e.g., microseismic, electric resistivity tomography, etc.) indicate a testbed that has undergone several changes since October 2018. In this paper, we present tracer recovery data accumulated during several tracer campaigns, provide conceptual testbed flow pathways, and provide observations that suggest a evolutionary nature of fracture volume and fracture geometry in the testbed.

15 GEOTHERMAL ENERGY↗
DOE OSTI2022

An index to characterize gas‐solid and solid‐solid mixing from average volume fraction fields

Abstract A mixing index based on solid volume fraction fields is developed for gas‐solid flows. Conventional mixing indices are based on particle realizations of granular mixing and are applicable to experimental data or discrete element method simulations. However, these indices cannot be used as‐is for multifluid models, and an index for characterizing mixing in gas‐solid flows from continuous fields is needed. The performance of the new mixing index is tested in two applications. The first is a 3D simulation of the mixing of biomass and sand in a fluidized bed reactor, and the second is a 2D simulation of binary particle segregation in a fluidized bed. The simulations are performed using OpenFOAM®. The mixing index is used to quantify gas‐solid mixing using solid volume fractions and solid‐solid mixing using solid fractions. The formulation of conventional mixing indices is extended to be used with solid volume fractions fields, and methods for performance improvement are presented.

09 BIOMASS FUELS↗
DOE OSTI2021

Effect of operating parameters on H 2 /CO 2 conversion to methanol over Cu-Zn oxide supported on ZrO 2 polymorph catalysts: Characterization and kinetics

Kinetic aspects of the operating parameters for the catalytic conversion of H 2 /CO 2 to methanol over two novel catalysts were evaluated to understand the effect of the polymorphic ZrO 2 phase composed of Cu 0/+ -ZnO sites at the atomic level and its impact on the reaction mechanism. The catalysts were characterized by in situ and ex-situ XRD, N 2 adsorption/desorption isotherms, FRX, TPR, TPD-N 2 O, in situ XANES, TPD-CO 2 , and in situ DRIFTS techniques. The influence of different reaction variables such as the GHSV, temperature, pressure, and H 2 /CO 2 ratio were studied using a fixed bed continuous plug flow reactor. The Cu-ZnO catalyst supported on the tetragonal zirconia polymorph exhibited the highest methanol yield due to the lower activation energy when compared to the catalyst with a greater amount of the monoclinic phase. In addition, the catalysts were reused for 8 cycles of 6 hours to evaluate their stability, which can translate into lower costs for large-scale methanol production. As a result, the estimation of the kinetic parameters over Cu-Zn oxide supported on ZrO 2 polymorphs was significant for understanding the reaction mechanism, as well as to provide a pathway for scaling-up of the process.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Understanding the Dissolution and Passivation of an Aluminum Electrode during Electrocoagulation of Groundwater Using Neutron and X-ray Reflectometry

An aluminum (Al)-based electrocoagulation (EC) system can effectively remove dissolved silica and hardness in groundwater. The effectiveness of Al-EC in terms of pollutant removal, Faradaic efficiency, and energy consumption depends on the interfacial electrolysis or passivation of the electrode in water. Thus, understanding the electrolysis reaction at the liquid/electrode interface during operation is important for sustainable EC deployment. Here, a continuous flow-through Al-EC system was tested with various groundwater simulants, i.e., chloride (Cl – )-based, sulfate (SO 4 2– )-based, and mixed solutions. High pollutant removal with low energy consumption was observed in Cl – -based groundwater treatment, while low pollutant removal with high energy consumption was observed in SO 4 2– -based groundwater. For example, the required energy per unit mass of Al dosing in SO 4 2– -based groundwater is three times higher than that in Cl – -based groundwater at 10 mA/cm 2 . However, increasing the Cl – concentration significantly reduces this energy demand. In SO 4 2– -based groundwater, the silicate removal efficiency drops from 85.1% to 24.0% compared to that for Cl – -based groundwater, while Mg 2+ and Ca 2+ removal efficiencies decrease to 0.6% from 15.8% and 5.7% from 44.8%, respectively. To better understand this EC performance, we used in situ neutron reflectometry (NR) to examine the interfacial dynamics of Al dissolution and passivation at a 100 nm scale occurring on the surface of the sacrificial Al electrodes during EC. Ex situ X-ray reflectometry (XRR) was also used to support the in situ NR results. Both NR and XRR results revealed that Al dissolution is influenced by the presence of Cl – in the simulants, while a passivating layer forms on the electrode in a SO 4 2– -based solution. In the Cl – -based solution, anodic Al dissolution occurred locally and inhomogeneously across the surface of the Al anode film, resulting in a localized thickness reduction over time. In the SO 4 2– -based solution, no apparent dissolution of the Al anode was identified. Instead, Al underwent oxidation, forming an amorphous Al 2 O 3 surface layer within the Al electrode film that increased in thickness over time. In the mixed solution, both anodic Al dissolution and surface Al 2 O 3 layer formation occurred, indicating that Al dissolution and surface Al 2 O 3 layer formation are attributable to the Cl – and SO 4 2– ions, respectively.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Mechanocatalytic Ammonia Synthesis over TiN in Transient Microenvironments

Environmentally friendly and energy-efficient ways to produce ammonia are essential to meet global food demands. In this study, a new approach for ammonia production at nominally ambient conditions is introduced. As proof of concept, ammonia is synthesized mechanocatalytically by ball milling titanium in a continuous gas flow. The ammonia synthesis reaction is proposed to follow a transient Mars–van Krevelen mechanism under mechanically activated conditions, where molecular nitrogen incorporation into the titanium lattice and titanium nitride hydrogenation occur in thermodynamically distinct environments. X-ray powder diffraction and X-ray absorption spectroscopy confirm the formation of titanium nitride from titanium and N2. The reactivity of nitrided titanium supports that lattice nitrogen plays a role in ammonia formation. The in situ formed titanium nitride is catalytically active, and the nitride regeneration reaction is determined to be the rate-limiting step. A preliminary technoeconomic analysis shows that this approach could be feasible for distributed ammonia production.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Porous Iron Electrodes Reduce Energy Consumption During Electrocoagulation of a Virus Surrogate: Insights into Performance Enhancements Using Three-Dimensional Neutron Computed Tomography

Electrocoagulation has attracted significant attention as an alternative to conventional chemical coagulation because it is capable of removing a wide range of contaminants and has several potential advantages. In contrast to most electrocoagulation research that has been performed with nonporous electrodes, in this study, we demonstrate energy-efficient iron electrocoagulation using porous electrodes. In batch operation, investigation of the external pore structures through optical microscopy suggested that a low porosity electrode with sparse connection between pores may lead to mechanical failure of the pore network during electrolysis, whereas a high porosity electrode is vulnerable to pore clogging. Electrodes with intermediate porosity, instead, only suffered a moderate surface deposition, leading to electrical energy savings of 21% and 36% in terms of electrocoagulant delivery and unit log virus reduction, respectively. Neutron computed tomography revealed the critical role of electrode porosity in utilizing the electrode’s internal surface for electrodissolution and effective delivery of electrocoagulant to the bulk. Energy savings of up to 88% in short-term operation were obtained with porous electrodes in a continuous flow-through system. Further investigation on the impact of current density and porosity in long-term operation is desired as well as the capital cost of porous electrodes.

42 ENGINEERING↗
DOE OSTI2025

Ultracold cryogenic TEM with liquid helium and high stability

Cryogenic transmission electron microscopy has revolutionized structural biology and materials science. To image below liquid nitrogen temperatures, various liquid helium stages have been constructed but have proven to be complex and unstable, making high-resolution imaging challenging. This problem is even more pronounced in side-entry specimen holders common on modern transmission electron microscopes. Here, we introduce an ultracold liquid helium transmission electron microscope side-entry specimen holder, featuring continuous cryogen flow and vibration decoupling. This instrument is compatible with modern aberration-corrected microscopes and achieves sub-25 K base temperature, ±2 mK thermal stability over many hours, and atomic resolution—setting the stage for a new era of cryogenic electron microscopy.

42 ENGINEERING↗
DOE OSTI2024

Calibrating Bayesian generative machine learning for Bayesiamplification

Recently, combinations of generative and Bayesian deep learning have been introduced in particle physics for both fast detector simulation and inference tasks. These neural networks aim to quantify the uncertainty on the generated distribution originating from limited training statistics. The interpretation of a distribution-wide uncertainty however remains ill-defined. We show a clear scheme for quantifying the calibration of Bayesian generative machine learning models. For a Continuous Normalizing Flow applied to a low-dimensional toy example, we evaluate the calibration of Bayesian uncertainties from either a mean-field Gaussian weight posterior, or Monte Carlo sampling network weights, to gauge their behaviour on unsteady distribution edges. Well calibrated uncertainties can then be used to roughly estimate the number of uncorrelated truth samples that are equivalent to the generated sample and clearly indicate data amplification for smooth features of the distribution.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Coordinated time variability of multi-phase ultra-fast outflows in J132216.25 + 052446.3

ABSTRACT We present a time variability analysis of broad absorption lines (BAL; spread over the velocity range of 5800–29 000 km s−1) seen in the spectrum of J132216.25 + 052446.3 (zem = 2.04806) at ten different epochs spanning over 19 yr. The strongest absorption component (BAL-A; spread over 5800–9900 km s−1) is made up of several narrow components having velocity separations close to C iv doublet splitting. The C iv, N v, and Si iv absorption from BAL-A show correlated optical depth variability without major changes in the velocity structure. A very broad and shallow absorption (BAL-C; spread over the velocity range 15 000–29 000 km s−1) emerged during our monitoring period coinciding with a dimming episode of J1322 + 0524. All the identified absorption lines show correlated variability with the equivalent widths increasing with decreasing flux. This together with the C iv emission line variability is consistent with ionization being the main driver of the correlated variability. The observed UV-continuum variations are weaker than what is required by the photoionization models. This together with a scatter in the C iv equivalent width at a given continuum flux can be understood if variations of the C iv ionizing photons are much larger than that of the UV continuum, the variations in the ionizing photon and UV fluxes are not correlated and/or the covering factor of the flow varies continuously. We suggest BAL-A is produced by a stable clumpy outflow located beyond the broad emission line region and BAL-C is a newly formed wind component located near the accretion disc and both respond to changes in the ionizing continuum.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Relativistic gas accretion onto supermassive black hole binaries from inspiral through merger

Accreting supermassive black hole binaries are powerful multimessenger sources emitting both gravitational and electromagnetic (EM) radiation. Understanding the accretion dynamics of these systems and predicting their distinctive EM signals is crucial to informing and guiding upcoming efforts aimed at detecting gravitational waves produced by these binaries. To this end, accurate numerical modeling is required to describe both the spacetime and the magnetized gas around the black holes. In this paper, we present two key advances in this field of research. First, we have developed a novel 3D general relativistic magnetohydrodynamics (GRMHD) framework that combines multiple numerical codes to simulate the inspiral and merger of supermassive black hole binaries starting from realistic initial data and running all the way through merger. Throughout the evolution, we adopt a simple but functional prescription to account for gas cooling through photon emission. Next, we have applied our new computational method to follow the time evolution of a circular, equal-mass, nonspinning black hole binary for ∼200 orbits, starting from a separation of 20⁢𝑟 𝑔 and reaching the postmerger evolutionary stage of the system. We have shown how mass continues to flow toward the binary even after the binary “decouples” from its surrounding disk, but the accretion rate onto the black holes diminishes. We have identified how the minidisks orbiting each black hole are slowly drained and eventually dissolve as the binary compresses. We confirm previous findings that the system’s luminosity decreases by a factor of a few during inspiral; however, we observe an abrupt increase by ∼50% in this quantity at the time of merger, likely accompanied by an equally abrupt change in spectrum. Lastly, we have demonstrated that during the inspiral, fluid ram pressure regulates the fraction of the magnetic flux transported to the binary that attaches to the black holes’ horizons.

Accretion disk & black-hole plasma↗
DOE OSTI2022

Load Margin Constrained Moving Target Defense against False Data Injection Attacks

Cyber physical security of power systems with high penetration of renewable generation has attracted attention from researchers. One critical issue is that cyber-physical attacks, disguised as uncertain renewable generation, can target conventional power system state estimation (SE). Moving target defense (MTD) is a promising defense strategy to detect stealthy false data injection (FDI) attacks against SE. However, all existing studies myopically perturb the reactance of transmission lines equipped with distributed flexible AC transmission system (D-FACTS) devices without adequately considering the system voltage stability. Exacerbated by the renewable generation uncertainty, existing MTD may cause voltage instability when the power grid is under stress. To address this issue, we propose a novel MTD framework that explicitly considers system voltage stability by using continuation power flow. We utilize the sensitivity matrix of power injection to line impedance, on which an optimization problem for maximizing load margin is formulated. This framework is validated on the IEEE 14-bus system and the IEEE 118-bus system, in which net load redistribution attacks are launched by sophisticated attackers. Steady-state simulations and dynamic simulations on PSS/E show the effectiveness of the proposed framework in circumventing the voltage instability while maintaining the detection effectiveness of MTD. The impact of the proposed method on attack detection effectiveness is also revealed.

Zhang, Hang↗
DOE OSTI2023

Wassersplines for Neural Vector Field-Controlled Animation

Much of computer-generated animation is created by manipulating meshes with rigs. While this approach works well for animating articulated objects like animals, it has limited flexibility for animating less structured free-form objects. Here we introduce Wassersplines, a novel trajectory inference method for animating unstructured densities based on recent advances in continuous normalizing flows and optimal transport. The key idea is to train a neurally-parameterized velocity field that represents the motion between keyframes. Trajectories are then computed by advecting keyframes through the velocity field. We solve an additional Wasserstein barycenter interpolation problem to guarantee strict adherence to keyframes. Our tool can stylize trajectories through a variety of PDE-based regularizers to create different visual effects. We demonstrate our tool on various keyframe interpolation problems to produce temporally-coherent animations without meshing or rigging.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Numerical Modeling of Atmospheric Rime Ice Accretion on an Airfoil Using an Eulerian Approach

This paper presents an approach to numerically simulate the inherently unsteady rime ice accretion problem on a two-dimensional airfoil and elucidate the associated variations under different icing conditions. The airflow field and the water impingement on the airfoil are obtained based on an Eulerian two-phase model. A dynamic mesh strategy is employed to unsteadily account for the changes in the ice profile and its impact on the air and droplet flow by continuously reconstructing the computational grid at each time-step through smoothing and layering mechanisms. All main icing modules including the airflow field, droplet trajectory, icing thickness profile, and mesh management are fully coupled within the same computational framework without resorting to any external tools. Classical icing theory is employed to model the rime ice roughness, and it is assumed that the ice accretes in a direction normal to the airfoil surface. The governing Reynolds-averaged Navier–Stokes (RANS) conservation equations along with the energy and continuity equations are solved to produce the velocity and temperature fields. A convective film heat transfer coefficient is computed based on the surface heat flux and a recovery temperature which takes into account the dissipative heat release in the boundary layer in the vicinity of the airfoil surface. With the implemented strategy and calculating the convective heat transfer coefficient, the water film thickness is also calculated along with the ice shape. Furthermore, the model is validated by comparing the local collection efficiency distribution and ice shape with experimental data, and the results show that the implemented approach provides acceptable predictions of ice accretion profiles and rates.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗