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.
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.
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.
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.
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.
Identification of Climatological Representative Days in the Mid-Atlantic for High-Fidelity Offshore Wind Energy Modeling
The goal of reaching 30 GW of offshore wind energy by 2030 becomes more realistic with the continued approval of offshore wind energy areas by the Biden Administration. In the Mid-Atlantic, where wind energy projects are in the most advanced stages of development, there is increased research focus on the eventual interaction of these wind farms. These interactions, in the form of wakes and cluster wakes, or wakes from multiple wind farms, could have detrimental effects on power production and forecastability for downwind wind farms (Pryor et al. 2022, Golbazi et al. 2022, Rosencrans et al. 2023). To help alleviate these issues, numerical simulations in the form of numerical weather prediction (NWP) and large eddy simulations (LES) can provide insight into when cluster wake situations may occur, but running such simulations can be expensive and difficult to run for multiple years. In this study, we leverage and build upon existing techniques in the literature (Fischereit et al. 2022) to identify climatologically representative days for wind energy areas in the Mid-Atlantic where conditions would promote cluster wake situations. We select meteorological variables (wind speed, wind direction, atmospheric stability, boundary-layer height, TKE) critical to understanding wind energy production and wake propagation. We then consider two different NWP datasets of varying spatial and temporal resolution: ERA5 provides data at hourly intervals from 1940 to present at 0.25 deg (31 km) spatial resolution (Hersbach et al. 2020), and the NOW-23 dataset provides data at 5-minute resolution for 21 years at 2-km spatial resolution (Bodini et al. 2020). Our first step is to compare these two datasets for an overlapping 21-year time period. Initial results show that the required number of days to represent the long-term climate increases with each additional variable considered. In their study of the German Bight, Fischereit et al. (2022) found that they could represent the long-term wind and wave climate in a "near-perfect" way with -180 days, by reaching a Perkins Skill Score (PSS) of 0.9; our investigation of the mid-Atlantic wind resource region with ERA5 and NOW-23 data suggests that we will need -100 days to reach a PSS of 0.9. As we expand our parameter space to include multiple variables, the number of required days will likely grow. These results will ultimately be used to select case studies to best represent cluster wake conditions that apply to this region for the lifetime of likely wind farms in this mid-Atlantic region.
Quality-Controlled Meteorological Data from the Flood Control District of Maricopa County (FCDMC) Network, Phoenix, Arizona (1987-2024)
This dataset contains 15- or 30-minute interval meteorological data from the Flood Control District of Maricopa County (FCDMC), Arizona, USA, covering eight key variables across multiple sensor stations between 1987 and 2024. Each variable is stored as a separate CSV file, containing time-series data that have undergone rigorous quality control (QC) procedures and, where appropriate, short-gap interpolation for consistency. The quality control (QC) pipeline consisted of four sequential tests: (1) a range test to ensure all values fall within physically realistic limits, (2) a step test to identify abrupt and implausible changes between consecutive records, (3) a proximity test that validates flagged values from step test using data from nearby stations and exceedance probability thresholds, and (4) a persistence test to detect and remove periods of unrealistically constant readings. These thresholds were calibrated to Arizona’s environmental conditions and sensor specifications. After QC, short gaps (≤2 hours) were linearly interpolated to ensure consistent temporal resolution, except for wind variables. Due to a major upgrade in FCDMC’s data transmission system, only ALERT-2 protocol data (2016–2024) for wind variables are included; earlier ALERT-1 data were excluded because of irregular sampling and high missing rates. This dataset supports regional climate and infrastructure resilience studies by providing standardized, high-resolution meteorological data for the greater Phoenix metropolitan area.
Low-Temperature Plasma-Based Metrology of Lithium-Ion Battery Electrode Materials (CRADA Final Report)
As part of the Cyclotron Road program, SirenOpt Inc. evaluated its low-temperature plasma-based metrology sensor prototype for measuring multiple critical properties of lithium-ion battery electrode materials in parallel and in real-time. Cost-effective, minimal-waste manufacturing of high-performance battery electrode materials will be vital for achieving society’s net-zero carbon emission goals. Because existing electrode metrology sensors cannot operate within most sections of manufacturing lines, manufacturers often complete hundreds of processing steps before they can test their products and detect problems. When manufacturers perform these offline tests, they typically only test a small portion of the manufactured products. Current electrode manufacturing thus often yields many low-quality products, or off-spec products that must be thrown away all together. For example, at least 6% of the total lithium-ion battery manufacturing cost (i.e., over $250 million/year for the average gigafactory) is devoted to processing defective electrodes that are not scrapped until performance tests are failed during late-stage quality control checks. Electrode variability also leads manufacturers to build extra cells into battery packs to reduce the risk of poor performance. For example, many electric vehicle (EV) manufacturers include up to 10% more cells than needed, which substantially increases the cost and weight of the final EV product. The SirenOpt sensor can potentially enable early detection of poorly manufactured electrodes and allow them to be removed earlier from manufacturing lines, which can save battery manufacturers (hundreds of) millions of dollars per year. The sensor can further be used to improve product quality by accelerating R&D and process optimization, improving quality control, and enabling real-time process control. Overall, a real-time, in-situ metrology strategy can create unprecedented opportunities for implementation of smart manufacturing practices and advanced quality and process control solutions to realize higher battery electrode throughput and performance.
Measurement of Neutron Multiplicity in Charged Current Neutrino Interactions on Oxygen
The Accelerator Neutrino Neutron Interaction Experiment (ANNIE) is a 26-ton gadolinium-doped water Cherenkov detector located 100~meters downstream in the Booster Neutrino Beam (BNB) at the Fermi National Accelerator Laboratory (Fermilab). Its primary goals are to (1) measure the neutron yield from $\nu_\mu$ interactions as a function of momentum transfer $Q^2$ so that neutrino-nucleus interaction models can be better constrained, and (2) demonstrate the power of novel, fast-timing detectors with the first deployment of Large Area Picosecond PhotoDetectors (LAPPDs). Current knowledge of neutrino-nucleus interactions fall short in modeling the topologies of such interactions, leading to inaccurate reconstruction of event kinematics such as particle energy, direction, and vertex. As a consequence, the accuracy of cross section measurements is impacted, which is necessary for precise physics measurements. Neutrons are an indication of inelasticity and affect the determination of the energy of the parent neutrino. Quantifying the neutron yield is a step towards reducing the associated uncertainties, and thus improving our understanding of these complex interactions and benefiting the next generation of long-baseline neutrino experiments. ANNIE will make use of LAPPDs to measure neutron multiplicity of CC-0$\pi$ $\nu_\mu$ interactions on oxygen, making it the first experiment to deploy an array of these photodetectors. Because the LAPPD is a novel photodetection technology, much customization is required to integrate it into existing electronics. The first half of this thesis covers the significant progress made towards the first deployment of the LAPPD system. From its test stand at Fermilab, the LAPPD system was systematically tested and put together until deployment readiness was achieved. I present my contributions to the design, fabrication, and testing of the waterproof housing and cables, and the commissioning of the LVHV board that powers the LAPPD and its readout electronics. These efforts brought the LAPPD system significantly closer to deployment, and eventually first data. The second half of this thesis presents the vertex and energy reconstruction algorithms developed to analysis the beam data with PMT-only information. While much progress has been made towards the deployment of LAPPDs, with several in the detector tank, efforts to integrate the LAPPD datastream are in progress. Thus, I developed a ring edge detection technique using PMT data to fit the muon vertex and determine its energy. The analysis in this thesis finds average neutron yields of $\Bar{n}_{data} = 0.452 \pm 0.039 (\text{stat}) \pm 0.27 (\text{sys})$ for a selection of muon neutrino candidates in the fiducial volume of ANNIE and corrected with an averaged neutron detection efficiency. An equivalent analysis for simulated beam data results in an average neutron yield of $\Bar{n}_{MC} = 0.582 \pm 0.018 (\text{stat}) \pm 0.25 (\text{sys})$. Future work includes application of efficiency corrections at a positional level, quantification of all systematic uncertainties, and neutron multiplicity measurements with other event topologies.
Machine Learning-Based Prediction of Distribution Network Voltage and Sensors Allocation: Preprint
Increasing penetration of fast-varying energy resources may negatively affect power systems' operation. At the same time, sensor deployment throughout distribution networks improves system awareness and enables the development of new and advanced voltage control solutions. Such control techniques rely on accurate prediction in anticipation to voltage violation scenarios. This paper analyzes various approaches for voltage prediction in a distribution system; it is shown that combining multiple techniques into a single regressor improves its predictive power. Moreover, a two-step regressor is proposed, where initial predictions based on a global regressor are refined by local regressors; in this case, prediction errors decrease significantly. Additionally, a clustering approach is employed for performing sensor allocation, so that only the most influential buses are selected for monitoring without diminishing prediction accuracy.
Machine Learning-Based Prediction of Distribution Network Voltage and Sensor Allocation
Increasing penetration of fast-varying energy resources may negatively affect power systems' operation. At the same time, sensor deployment throughout distribution networks improves system awareness and enables the development of new and advanced voltage control solutions. Such control techniques rely on accurate prediction in anticipation to voltage violation scenarios. This paper analyzes various approaches for voltage prediction in a distribution system; it is shown that combining multiple techniques into a single regressor improves its predictive power. Moreover, a two-step regressor is proposed, where initial predictions based on a global regressor are refined by local regressors; in this case, prediction errors decrease significantly. Additionally, a clustering approach is employed for performing sensor allocation, so that only the most influential buses are selected for monitoring without diminishing prediction accuracy.
Production of levulinic acid and biocarbon electrode material from corn stover through an integrated biorefinery process
To overcome the inefficient biomass conversion, waste generation, and lack of co-production in biorefineries, an integrated process was proposed for the conversion of corn stover into levulinic acid and biocarbon electrode material. Corn stover was pretreated through hydrothermal process using 0.45 wt% K 2 CO 3 which removed 76 wt% lignin and 85 wt% xylan while preserving 83 wt% glucan. This was followed by acid hydrolysis to produce levulinic acid at varying H 2 SO 4 concentrations and reaction time in a batch reactor at 190 °C. At a reaction time of 5 min in 2 wt% H 2 SO 4 , 35.8 wt% and 30 wt% glucan in raw and pretreated corn stover was converted to levulinic acid, respectively. The residue from acid hydrolysis was converted into biocarbon for supercapacitor electrodes via a two-step thermal activation process which showed a specific capacitance of 120 F g -1 . The proposed integrated biorefinery concept provides multiple value-added products for a greater financial and environmental sustainability.
Numerical schemes for 3-wave kinetic equations: A complete treatment of the collision operator
In our previous work Walton and Tran (2023), numerical schemes for a simplified version of 3-wave kinetic equations, in which only the simple forward-cascade terms of the collision operators are kept, have been successfully designed, especially to capture the long time dynamics of the equation given the multiple blow-up time phenomenon. In this second work in the series, we propose numerical treatments for the complete 3-wave kinetic equations, in which the complete, much more complicated collision operators are fully considered based on a novel conservative form of the equation. Here we then derive an implicit finite volume scheme to solve the equation. The new discretization uses an adaptive time-stepping method which allows for the simulations to be carried to very long times. Our computed solutions are compared with previously derived long-time asymptotic estimates for the decay rate of total energy of time-dependent solutions of 3-wave kinetic equations and found to be in excellent agreement.
Machine Learning-Based Prediction of Distribution Network Voltage and Sensor Allocation
Increasing penetration levels of fast-varying energy resources might negatively affect power system operation. At the same time, sensor deployment throughout distribution networks improves system awareness and enables the development of new and advanced voltage control solutions. Such control techniques rely on accurate prediction in anticipation of voltage violation scenarios. This paper analyzes various approaches to voltage prediction in a distribution system, and it is shown that combining multiple techniques into a single regressor improves its predictive power. Moreover, a two-step regressor is proposed in which initial predictions based on a global regressor are refined by local regressors; in this case, prediction errors decrease significantly. Additionally, a clustering approach is employed to perform sensor allocation so that only the most influential buses are selected for monitoring without diminishing prediction accuracy.
Gamma Driven Catalysis of Ammonia
Experiments were conducted to investigate a passive production mechanism for the world’s most energy intensive commodity, ammonia. A novel method, gamma catalyzed ammonia production at ambient conditions, was investigated. Ammonia is currently produced through the highly energy intensive Haber-Fritz process, which requires an operation pressure of 400 atmosphere and 600 degrees Celsius. Due to the high demand and need for ammonia, the Haber-Bosch process consumes 25% all energy produced globally. Reported herein was an attempt to produce ammonia at ambient temperature (20 C) and ambient pressure (1 atm), through a novel process developed at PNNL, gamma driven catalysis of ammonia. Although the measurements of the ammonia production suggest wild success, reports in the literature by Gao et.al. suggest an experimental positive bias in the results. To rule out the potential positive bias, multiple additional production campaigns would be needed to with an alternate analysis technique such as ion chromatography, as suggested by Gao et.al. Unfortunately, due to this late determination of potential positive bias, the results of this study remain inconclusive to the feasibility of gamma driven catalysis of ammonia and more work is needed to describe the chemical evolution with time. The results are a first step and demonstrate that gamma-catalyst mediated reactions are possible. This represents a key opportunity to explore the fundamental chemistry of high band gap catalysts that can change the paradigm of radiation, transforming it from a waste to a valuable energy source.
From lab to lamp: Understanding downconverter degradation in LED packages
Downconverters, primarily inorganic phosphors, are critical components in white solid-state LED-based lighting and liquid crystal display backlights. Research efforts have led to a fundamental understanding of a downconverter's absorption, photoluminescence, and efficiency as a function of composition, structure, and processing conditions. However, considerably less work has focused on the reliability of phosphors once they are incorporated into LED packages. Solving these issues is often the final step before the commercialization of new materials, but the significant resources and time required to evaluate and mitigate materials failure are rarely discussed in the literature. In this Perspective, we discuss the need for conducting downconverter reliability testing and the potential of accelerating, screening, and understanding downconverter failure modes. Our focus highlights the mechanisms of failure and discusses how this influences materials selection and the design of different LED packages. We also stress the potential for accelerated reliability testing protocols and note the potential role first-principles calculations and data-driven models could play in establishing the compositional-processing trends for different aspects of downconverter reliability. In conclusion, we close with possible research directions that could improve downconverter reliability and emphasize the importance of assessing a material's (chemical) stability where multiple manufacturing and processing steps can dictate system performance.
Modeling Time-Dependent Surrogates of Additive-Manufactured Nuclear Fuels Processes
Additive manufacturing (AM) technology is being increasingly adopted in a wide variety of application areas for its ability to rapidly produce, prototype, and customize designs. Recently, a hybrid AM technique was successfully developed at Idaho National Laboratory (INL) to manufacture nuclear fuels [1]. Despite the advantages, this AM technique needs optimization due to defects from a highly complex melting and sintering process. The complex metallurgical phenomena during AM processes are strongly related to parameters such as applied laser power, traveling speed, and scan style, which could lead to differences in density, residual stress, crystallographic texture, and mechanical properties. In addition, stochastic variations in laser energy interaction and associated multiscale/multiphysics phenomena cause variations in microstructure evolution and mechanical properties. Currently, researchers at INL are focusing on developing a comprehensive modeling framework, leveraging INL’s simulation tools MOOSE/MARMOT/BISON/RAVEN [2-4] to describe all steps of this AM process across multiple length scales. Although this advanced framework plays a critical role in enabling enhancements to traditional trial and error approaches for design and optimization of nuclear fuel materials, it remains computationally intense, limiting its use in sensitivity and optimization analysis. In this case, an accurate and inexpensive surrogate becomes an effective tool for providing a tractable approximation of the underlying underline physics. Surrogate models generally not based on the physics of a system are purely mathematical models used to capture the relationships between specific system inputs and outputs. Popular approaches, including neural networks [5], response surfaces [6], and subspace-based reduced order models [7], have been applied to a wide range of disciplines, such as nuclear reactor design, aerospace design and automotive design. In this summary, we employ advanced time-dependent surrogate models such as high-dimensional model representation (HDMR) [8] and physics-informed deep neural network (PINNs) [9] to accelerate the design and optimization of AM process.
Integrative Multi-PTM Proteomics Reveals Dynamic Global, Redox, Phosphorylation, and Acetylation Regulation in Cytokine-treated Pancreatic Beta Cells
Studying regulation of protein function at a systems level necessitates an understanding of the interplay among diverse post-translational modifications (PTMs). A variety of proteomics sample processing workflows are currently used to study specific PTMs but rarely characterize multiple types of PTMs from the same sample inputs. Method incompatibilities and laborious sample preparation steps complicate large-scale physiological investigations and can lead to variations in results. The single-pot, solid-phase-enhanced sample preparation (SP3) method for sample cleanup is compatible with different lysis buffers and amenable to automation, making it attractive for high-throughput multi-PTM profiling. Herein, we describe an integrative SP3 workflow for multiplexed quantification of protein abundance, cysteine thiol oxidation, phosphorylation, and acetylation. The broad applicability of this approach is demonstrated using cell and tissue samples, and its utility for studying interacting regulatory networks is highlighted in a time-course experiment of cytokine-treated ß-cells. We observed a swift response in global regulation of protein abundances consistent with rapid activation of JAK-STAT and NF-?B signaling pathways. Regulators of these pathways as well as proteins involved in their target processes displayed multi-PTM dynamics indicative of a complex cellular response stages: acute, adaptation, and chronic (prolonged stress). PARP14, a negative regulator of JAK-STAT, had multiple co-localized PTMs that may be involved in intraprotein regulatory crosstalk. Our workflow provides a high-throughput platform that can profile multi-PTMomes from the same sample set, which is valuable in unraveling the functional roles of PTMs and their co-regulation.