Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “multiple time stepping”

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 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

3D Printing of Functional Hydrogel Devices for Screenings of Membrane Permeability and Selectivity

Developing a fundamental understanding of the effects of varying ligand chemistries on mass transport rates is key to designing membranes with solute-specific selectivity. While permeation cells offer a robust method to characterize membrane performance, they are limited to assessing a single membrane chemistry or salt solution per test. As a result, investigating the effects of varying ligand chemistries on membrane performance can be a tedious process, involving both the preparation of multiple samples and numerous, time-consuming permeation tests. This study uses digital light processing (DLP) 3D printing to fabricate a millifluidic flow-based permeation device made from a hydrogel active ester network that can be easily functionalized with ion-selective ligands. Without the need for bonding or assembly steps, ligands can be introduced and tested in the permeation device by simply injecting a small volume of a ligand solution. Various salt concentrations and molecular species can be cycled through a single device by switching the solution feeding into the salt reservoir, thereby reducing the number of samples needed for permeability and selectivity screenings. This research sets the groundwork for formulation development and postprocessing methods to 3D-print functional millifluidic devices capable of assessing solute selectivity in membranes and polymer adsorbents for aqueous separations. In this work, comparable salt permeability trends were observed with both 3D-printed devices and traditional assays. Devices were functionalized with an imidazole ligand to investigate salt permeability and selectivity of monovalent and divalent salts. Measurements showed increasing permeability for monovalent salts (NaCl) relative to divalent salts (MgCl 2 , CuCl 2 ) in functionalized membranes, with higher monovalent/divalent selectivity at increasing imidazole grafting densities. Here, the methods and findings described here represent a step toward developing higher-throughput methods with 3D-printed devices for screening the effects of ligand chemistry on mass transport rates in membrane materials.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Rapid assessment of SARS-CoV-2–evolved variants using virus-like particles

Efforts to determine why new severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) variants demonstrate improved fitness have been limited to analyzing mutations in the spike (S) protein with the use of S-pseudotyped particles. In this study, we show that SARS-CoV-2 virus-like particles (SC2-VLPs) can package and deliver exogenous transcripts, enabling analysis of mutations within all structural proteins and at multiple steps in the viral life cycle. In SC2-VLPs, four nucleocapsid (N) mutations found universally in more-transmissible variants independently increased messenger RNA delivery and expression ~10-fold, and in a reverse genetics model, the serine-202→arginine (S202R) and arginine-203→methionine (R203M) mutations each produced >50 times as much virus. SC2-VLPs provide a platform for rapid testing of viral variants outside of a biosafety level 3 setting and demonstrate N mutations and particle assembly to be mechanisms that could explain the increased spread of variants, including B.1.617.2 (Delta, which contains the R203M mutation).

60 APPLIED LIFE SCIENCES↗
DOE OSTI2025

Conversion of Compositionally Diverse Plastic Waste over Earth-Abundant Sulfides

Chemical deconstruction of polyolefin plastic wastes via hydroconversion is promising for mitigating plastic accumulation in landfills and the environment. However, hydroconversion catalysts cannot handle complex feedstocks containing multiple polymers, additives, and heteroatom impurities. Here, we report a single-step strategy using earth-abundant metal sulfide catalysts to deconstruct these wastes. We show that NiMoS x /HY catalysts deconstruct polyolefin feedstocks, achieving ~81–94% selectivity to liquid products. Postsynthetic zeolite modification enhances the catalyst’s activity by >2.5 times, achieving over 95% selectivity to liquid fuels with controllable product distribution in the naphtha, jet fuel, and diesel range. The catalyst is resilient to increasingly complex feedstocks, such as additive-containing polymers and mixed plastics composed of polyolefins and heteroatom-containing polymers, including poly(vinyl chloride). As a result, we extend the strategy to single-use polyolefin wastes that can generate toxic byproducts, such as HCl and NH 3 , and eliminate their emissions by integrating reaction and sorption in a one-step process.

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

Revealing Hidden Orbital Pseudospin Texture with Time-Reversal Dichroism in Photoelectron Angular Distributions

We performed angle-resolved photoemission spectroscopy (ARPES) of bulk 2H-WSe 2 for different crystal orientations linked to each other by time-reversal symmetry. We introduce a new observable called time-reversal dichroism in photoelectron angular distributions (TRDAD), which quantifies the modulation of the photoemission intensity upon effective time-reversal operation. We demonstrate that the hidden orbital pseudospin texture leaves its imprint on TRDAD, due to multiple orbital interference effects in photoemission. Our experimental results are in quantitative agreement with both the tight-binding model and state-of-the-art fully relativistic calculations performed using the one-step model of photoemission. While spin-resolved ARPES probes the spin component of entangled spin-orbital texture in multiorbital systems, we unambiguously demonstrate that TRDAD reveals its orbital pseudospin texture counterpart.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Modeled estimates of permafrost hydrology and related fields for pan-Arctic region over the period 1980-2100

This dataset contains gridded estimates of total runoff (mm/day); subsurface runoff (mm/day); snow water equivalent (SWE, mm/day); vertical upward water flux of combined surface evaporation, evapotranspiration, and snow sublimation (mm/day); and soil temperature (C) across the pan-Arctic basin from simulations of the Permafrost Water Balance Model (PWBM v4). These data were generated to allow researchers to advance understanding of the changes unfolding in the Arctic's terrestrial water cycle under the influences of climate warming, hydrological cycle intensification, and permafrost thaw. The dataset will enable researchers to address important questions spanning multiple disciplines. For example, how will rising river flows exported to the coast impact critical zone ecosystems, near-shore ice dynamics, and ocean biogeochemical processes. Spatial resolution is 25x25 km with N=39926 grid cells spanning the domain. Three forcing datasets were used to drive simulations at an implicit daily time step over the period 1980-2100. Daily air temperature, precipitation, and wind speed data from the WATCH Forcing Data methodology applied to ERA5 reanalysis (W5E5) were used for meteorological forcings over the period 1980-2019. Simulations for the period 1980-2100 used meteorological data from two global climate models; IPSL-CM6A-LR and MPI-ESM1-2-HR. The PWBM hydrological outputs are organized by grid and month in gzipped ASCII text files within each archived tar file. Soil temperature is organized by grid cell and day of year from DOY 205 to DOY 260. No special software is required to work with these data. The package contains spatial domain and river basin information files. These data are most useful for analyses of the spatial and temporal variations in terrestrial hydrology and soil thermal dynamics. A manuscript describing the model simulations and associated analysis has been accepted for publication in The Cryosphere (Rawlins and Karmalkar, 2024).

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Machine learning based rate optimization under geologic uncertainty

We propose a novel approach for rate optimization during a waterflood under geologic uncertainty in reservoir properties such as permeability and porosity. The traditional approach typically involves several runs of the forward simulator. This may not scale well when the optimization is to be performed at the full field-level and over multiple geologic realizations. A machine-learning (ML) based approach which is quick and scalable for rate optimization over multiple geologic realizations is proposed instead. The training data for the model is generated by running the forward simulator with randomly assigned well rates using multiple geologic realizations. A reduced order representation of the permeability heterogeneity in each of the realizations is derived using a grid connectivity transformation (GCT). This step involves finding basis functions corresponding to the different modal frequencies of the grid connectivity represented by the grid Laplacian. The projection of the heterogeneous property field along these basis functions gives the basis coefficients that form the reduced order representation. Subsequently, for each training datapoint, streamlines are traced and the minimum time of flight (TOF) representing the tracer breakthrough time at each producer is recorded. The basis coefficients and well rates are fed to a machine learning model as input and the minimum TOF at the producers forms the output of the model. This trained model can then be used along with an optimizer for computing the optimal injection rates to maximize the injection sweep efficiency. This corresponds to minimizing the variance in the minimum TOF within each well group. Different architectures of neural network are tested using 5-fold cross validation to decide the best ML model to compute the streamline time of flight. The trained model is used to perform well rate optimization over multiple realizations of geology by using a risk tolerance penalty. The optimal well rates thus obtained are compared with two cases: a) equal well rates assigned to all injectors and producers and b) well rates obtained by optimizing over a single realization without considering the uncertainty in geology. The optimal well rates are seen to offer better oil recovery and sweep efficiency than both cases.

02 PETROLEUM↗
DOE OSTI2020

Air-transfer production method for large-area picosecond photodetectors

We have designed and prototyped the process steps for the batch production of large-area micro-channel-plate photomultipliers (MCPPMT) using the “air-transfer” assembly process developed with single LAPPD TM modules. Results are presented addressing the challenges of designing a robust package that can transmit large numbers of electrical signals for pad or strip readout from inside the vacuum tube and of hermetically sealing the large-perimeter window–body interface. We have also synthesized a photocathode in a large-area low-aspect-ratio volume and have shown that the micro-channel plates recover their functionality after cathode synthesis. These steps inform a design for a multi-module batch facility employing dual nested low-vacuum and ultra-high-vacuum systems in a small-footprint. The facility design provides full access to multiple MCP-PMT modules prior to hermetic pinch-off for leak-checking and real-time photocathode optimization.

47 OTHER INSTRUMENTATION↗
DOE OSTI2025

Surrogate models for development of unconventional shale reservoirs by an integrated numerical approach of hydraulic fracturing, flow and geomechanics, and machine learning

We develop well-completion surrogate models by taking an integrated workflow of hydraulic fracturing, flow, geomechanics, and machine learning simulation. There are three steps in the proposed workflow. First, history-matching processes are conducted with the field data including pumping and production data for characterization. Second, full-physics simulation is performed with various parameters of the field development (e.g., cluster spacing, clusters per stage, pumping rates and times, amount of proppant, and well spacing) to generate multiple simulation results by changing the parameters of the completion design with well-known hydraulic fracturing, reservoir, geomechanics simulators to calculate fracture geometry, reservoir depressurization, induced stress changes. The workflow is demonstrated over a field in the Southern Midland Basin. Here, we take two completion scenarios: a single well case followed by a multi-well case. Finally, a Long Short-Term Memory (LSTM) machine learning algorithm is employed to create surrogate models that can replicate the full-physics simulation results. Furthermore, results show that the trained models applied in the single well and multi-well cases for a particular geological system can provide good accuracy close to those provided by full-physics simulations. Specifically, the site-specific surrogate models can predict fracture parameters (length, height, and surface area) and cumulative production accurately with computational efficiency, suggesting our proposed workflow can be used as a pragmatic tool for expediting the well completion optimization process.

Geomechanics↗
DOE OSTI2022

Energetic particle optimization of quasi-axisymmetric stellarator equilibria

Abstract An important goal of stellarator optimization is to achieve good confinement of energetic particles such as, in the case of a reactor, alphas created by deuterium–tritium fusion. In this work, a fixed-boundary stellarator equilibrium was re-optimized for energetic particle confinement via a two-step process: first, by minimizing deviations from quasi-axisymmetry (QA) on a single flux surface near the mid-radius, and secondly by maintaining this improved QA while minimizing the analytical quantity Γ C , which represents the angle between magnetic flux surfaces and contours of J | | , the second adiabatic invariant. This was performed multiple times, resulting in a group of equilibria with significantly reduced energetic particle losses, as evaluated by Monte Carlo simulations of alpha particles in scaled-up versions of the equilibria. This is the first time that energetic particle losses in a QA stellarator have successfully been reduced by optimizing Γ C . The relationship between energetic particle losses and metrics such as QA error ( E q a ) and Γ C in this set of equilibria were examined via statistical methods and a nearly linear relationship between volume-averaged Γ C and prompt particle losses was found.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

The LSST DESC DC2 Simulated Sky Survey

Here, we describe the simulated sky survey underlying the second data challenge (DC2) carried out in preparation for analysis of the Vera C. Rubin Observatory Legacy Survey of Space and Time (LSST) by the LSST Dark Energy Science Collaboration (LSST DESC). Significant connections across multiple science domains will be a hallmark of LSST; the DC2 program represents a unique modeling effort that stresses this interconnectivity in a way that has not been attempted before. This effort encompasses a full end-to-end approach: starting from a large N-body simulation, through setting up LSST-like observations including realistic cadences, through image simulations, and finally processing with Rubin's LSST Science Pipelines. This last step ensures that we generate data products resembling those to be delivered by the Rubin Observatory as closely as is currently possible. The simulated DC2 sky survey covers six optical bands in a wide-fast-deep area of approximately 300 deg 2 , as well as a deep drilling field of approximately 1 deg 2 . We simulate 5 yr of the planned 10 yr survey. The DC2 sky survey has multiple purposes. First, the LSST DESC working groups can use the data set to develop a range of DESC analysis pipelines to prepare for the advent of actual data. Second, it serves as a realistic test bed for the image processing software under development for LSST by the Rubin Observatory. In particular, simulated data provide a controlled way to investigate certain image-level systematic effects. Finally, the DC2 sky survey enables the exploration of new scientific ideas in both static and time domain cosmology.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Hardware-based Advanced Electromagnetic Transient Simulation for A Large-Scale PV Plant in Real Time Digital Simulator

Power electronics-based resources, such as high-voltage direct current (HVdc) substations, photovoltaic (PV) plants, wind plants, electric vehicle charging stations, and energy storage systems, are increasingly being integrated within the power grid. Recently, multiple reports have emphasized the necessity for high-fidelity electromagnetic transient (EMT) simulations of these large-scale power electronics-based resources to accurately understand their behavior in power grids. However, performing hardware-based EMT simulations with high-fidelity models for such large power electronics systems is challenging due to the small time-step requirements and the involvement of a large number of states. This paper presents the implementation of hardware-based high-fidelity EMT dynamic model of a large-scale PV plant, accomplished through custom model development using the specific-C language in real-time digital simulator hardware (RTDS) and software (RSCAD).

Choi, Jongchan↗
DOE OSTI2024

A Comprehensive Analysis of Uncertainties in Warm-Rain Parameterizations in Climate Models Based on In Situ Measurements

Abstract Because of the coarse grid size of Earth system models (ESMs), representing warm-rain processes in ESMs is a challenging task involving multiple sources of uncertainty. Previous studies evaluated warm-rain parameterizations mainly according to their performance in emulating collision–coalescence rates for local droplet populations over a short period of a few seconds. The representativeness of these local process rates comes into question when applied in ESMs for grid sizes on the order of 100 km and time steps on the order of 20–30 min. We evaluate several widely used warm-rain parameterizations in ESM application scenarios. In the comparison of local and instantaneous autoconversion rates, the two parameterization schemes based on numerical fitting to stochastic collection equation (SCE) results perform best. However, because of Jessen’s inequality, their performance deteriorates when grid-mean, instead of locally resolved, cloud properties are used in their simulations. In contrast, the effect of Jessen’s inequality partly cancels the overestimation problem of two semianalytical schemes, leading to an improvement in the ESM-like comparison. In the assessment of uncertainty due to the large time step of ESMs, it is found that the rainwater tendency simulated by the SCE is roughly linear for time steps smaller than 10 min, but the nonlinearity effect becomes significant for larger time steps, leading to errors up to a factor of 4 for a time step of 20 min. After considering all uncertainties, the grid-mean and time-averaged rainwater tendency based on the parameterization schemes is mostly within a factor of 4 of the local benchmark results simulated by SCE.

Meteorology & Atmospheric Sciences↗
DOE OSTI2026

Comparison of time-resolved photoluminescence and deep-level transient spectroscopy defect evaluations in an InAs nBn detector subjected to in situ and ex situ 63 MeV proton irradiation

Deep-level transient spectroscopy and temperature-dependent time-resolved photoluminescence experiments are performed on identical InAs nBn photodetector structures as a function of in situ and ex situ 63 MeV proton irradiation to assess their generation and recombination dynamics. Pre-irradiation, the n-type InAs absorbing region, exhibits a steadily increasing minority carrier lifetime with increasing temperature, providing evidence that excited minority carriers may be recombining via shallow defect levels. From deep-level transient spectroscopy, two features are found between 10 and 275 K: a low temperature broad “shoulder,” which suggests emission from multiple shallow electron defect levels with energies <29 meV and a high temperature minimum occurring at ∼230 K with an activation energy of 539 meV, which suggests a defect in the barrier layer in the device. Two similar nBn detectors are then subjected to 63 MeV proton irradiation in step doses and measured between steps. One experiment is performed in situ with an nBn held at ∼10 K during dosing, and the other experiment is performed ex situ with a similar nBn held at room temperature for dosing. The ex situ dosing results in an evaluation of the defect introduction rate that is three to four times lower than in situ due to partial annealing of the proton-induced displacement damage at room temperature. The results of these two experiments are then compared with the dose-dependent recombination rate analysis, resulting in an estimated recombination defect cross section of 1.6 × 10 −13 cm 2 for the shallow shoulder defect.

Carrasco, Rigo A. [Air Force Research Laboratory (↗
DOE OSTI2025

Subject-specific modeling framework for particle deposition using computational fluid dynamics

Quantifying particle deposition and dose in the respiratory tract requires a physiologically realistic representation and reproducible computational workflows. However, existing modeling frameworks, such as the International Commission on Radiological Protection (ICRP) compartmental models and the Multiple Path Particle Dosimetry (MPPD) tool, lack detailed deposition profiles and subject-specific capabilities. The combination of advances in computer vision algorithms applied to the respiratory tract and Computational Fluid and Particle Dynamics (CFPD) allows high-fidelity simulations of particle behavior in anatomically accurate geometries derived from individual CT scans. The segmentation, preprocessing, and file preparation task for a CFPD simulation was often time-consuming, and no prior studies to-date have yet presented a fully automated framework. This work presents a fully automated workflow to obtain individualized particle deposition profiles in the human respiratory tract. The pipeline starts with segmenting upper and lower airway geometries using morphological and deep learning-based methods, generating three-dimensional (3D) models from CT imaging data. Next, a series of algorithms are presented to quality check and prepare the 3D geometry for a CFD or CFPD simulation. The preprocessing step includes correcting geometric artifacts, enforcing a physically consistent mesh, and automatically identifying and capping multiple outlets, which is required for CFD/CFPD simulations. These processed models are then input into open-source (OpenFOAM) or commercial (StarCCM+) CFD solvers, where flow and transient particle transport equations — including turbulence and particle–wall interactions are solved under realistic breathing conditions. Finally, the resulting particle deposition profiles can be integrated with Monte Carlo radiation transport codes and state-of-the-art computational phantoms to assess organ-specific absorbed doses in scenarios of radioactive aerosol inhalation. The presented work streamlines respiratory tract segmentation, preprocessing for CFD/CFPD simulations, and integration with dose assessment workflows, reducing manual intervention and improving access to high-fidelity, subject-specific modeling. The high precision in predicted particle deposition and dose distributions can improve personalized treatment strategies in respiratory medicine and refine dose estimates for radiation protection.

AI↗
DOE OSTI2022

A techno-economic assessment framework for hydrogen energy storage toward multiple energy delivery pathways and grid services

Hydrogen energy storage (HES) transforms and stores electric energy from the grid into hydrogen, and supplements other energy storage and demand response resources in addressing challenges in renewable-intensive power systems. Understanding how to optimally utilize an HES system to maximize its economic benefits from stacked value streams is highly important to its development and deployment. Here, in this paper, we present a techno-economic assessment framework for an HES system considering three common energy delivery pathways and multiple grid and end-user services. Models are developed to capture the operational capability, flexibility, and constraints associated with hydrogen production, compression, storage, and utilization as well as different grid services in an economic assessment. To define the technically achievable benefits, an optimal dispatch formulation is proposed to maximize the economic benefits over a representative year with an hourly time step considering the trade-offs among different value streams. Representative case studies are designed and carried out to show how system configuration, energy delivery pathways, and grid services may affect economic benefits. It was found that value streams from bundling grid services account for up to 76% of the total benefits and are critical for an HES project to be financially viable.

25 ENERGY STORAGE↗
DOE OSTI2025

Coupling flux balance analysis with reactive transport modeling through machine learning for rapid and stable simulation of microbial metabolic switching

Integrating genome-scale metabolic networks with reactive transport models (RTMs) provides a detailed description of the dynamic changes in microbial growth and metabolism. Despite promising demonstrations in the past, computational inefficiency has been pointed out as a critical issue to overcome because it requires repeated application of linear programming (LP) to obtain flux balance analysis (FBA) solutions in every time step and spatial grid. To address this challenge, we propose a new simulation method where we train and validate artificial neural networks (ANNs) using randomly sampled FBA solutions and incorporate the resulting surrogate FBA model (represented as algebraic equations) into RTMs as source/sink terms. We demonstrate the efficiency of our method via a case study of Shewanella oneidensis MR-1. During aerobic growth on lactate, S. oneidensis produces metabolic byproducts (such as pyruvate and acetate), which are subsequently consumed as alternative carbon sources when the preferred nutrients are depleted. To effectively simulate these complex dynamics, we used a cybernetic approach that models metabolic switches as the outcome of dynamic competition among multiple growth options. In both zero-dimensional batch and one-dimensional column configurations, the ANN-based surrogate models achieved substantial reduction of computational time by several orders of magnitude compared to the original LP-based FBA models. Moreover, the ANN models produced robust solutions without any special measures to prevent numerical instability. These developments significantly promote our ability to utilize genome-scale networks in complex, multi-physics, and multi-dimensional ecosystem modeling.

59 BASIC BIOLOGICAL SCIENCES↗