Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Simulation and optimization”

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 847 records · Page 47
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Accurate Force Field for Carbon Dioxide–Silica Interactions Based on Density Functional Theory

Fluid–silica interfaces are ubiquitous in chemistry, occurring in both natural geochemical environments and practical applications ranging from separations to catalysis. Simulations of these interfaces have been, and continue to be, a significant avenue for understanding their behavior. A constraining factor, however, is the availability of accurate force fields. Most simulations use traditional “mixing rules” to determine nonbonded dispersion interactions, an approach that has not been critically examined. Here, in this study, we present Lennard-Jones parameters for the interaction of carbon dioxide with silica interfaces that are optimized to reproduce density functional theory (DFT)-based binding energies. The modeling is based on the recently developed silica-DDEC force field, whose atomic charges are consistent with DFT calculations. Standard mixing rules are found to predict weaker CO 2 binding to silica than that obtained from DFT, an effect corrected by the optimized parameters given here. This behavior extends to other silica force fields (Clayff and Gulmen-Thompson), and the present Lennard-Jones parameters improve their performance as well. The effects of improved Lennard-Jones parameters on the structural and dynamical properties of condensed CO 2 in silica slit pores are also examined.

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

High-Capacity Enzymatic Degradation of Postconsumer Poly(ethylene terephthalate): Modeling and Experimental Investigations

Enzymatic degradation is considered as one of the key steps for biorecycling of polyethylene terephthalate (PET), a widely used plastic. Recently, we have developed an optimized leaf-branch compost cutinase PelB-LCC ICCG and demonstrated ∼80% degradation of 200 g L −1 recycled PET (RPET) within two days in bioreactors. In this study, continued research efforts were made to achieve a complete degradation of up to 300 g L −1 RPET in bioreactors. First, low-capacity reaction conditions with ≤5 g L −1 PET were examined to identify the key variables for PET degradation with PelB-LCC ICCG , including temperature, enzyme loading, and product inhibition. To better understand and further optimize the enzymatic degradation process, a mechanism-based model was established to describe the kinetics of PET degradation and formation of the main product terephthalic acid (TPA) and byproducts MHET and BHET in high-capacity bioreactors. Model simulation suggested that a minimal enzyme loading of ∼1.4 mg PelB-LCC ICCG g −1 RPET is required to achieve a nearly complete degradation of RPET within 48 h, which was used to guide more high-capacity experiments with 100−300 g L −1 RPET in fully controlled 1 L bioreactors. A higher temperature (≥65 °C) was found not only to enable rapid degradation in the beginning but also to induce a gradual increase in RPET’s crystallinity and significantly slow down the degradation after 48 h. A high loading of RPET solids and the accumulation of the produced insoluble TPA pose a big challenge on mixing and mass transfer in the stirred bioreactor, which can be addressed by increasing the stirring speed. The results pave the way toward biorecycling of PET at a large scale.

enzymatic degradation↗
DOE OSTI2025

Using the ATLAS experiment software on heterogeneous resources

With the large dataset expected from 2030 onwards by the HL-LHC at CERN, the ATLAS experiment is reaching the limits of the current data processing model in terms of traditional CPU resources based on x86_64 architectures and an extensive program for software upgrades towards the HL-LHC has been set up. The ARM CPU architecture is becoming a competitive and energy efficient alternative. Accelerators like GPUs are available in any recent HPC. In the past years ATLAS has successfully ported its full data processing and simulation software framework Athena to ARM and has invested significant effort in porting parts of the reconstruction and simulation algorithms to GPUs. We report on the successful usage of the ATLAS experiment offline and online software framework Athena on ARM and GPUs through the PanDA workflow management system at various WLCG sites. Furthermore we report on performance optimizations of the builds for ARM CPUs and the GPU integration efforts. We will discuss performance comparisons of different ARM and x86_64 architectures on WLCG resources and Cloud compute providers like GCP and AWS using ATLAS productions workflows as used in the Hep-Score23 benchmark suite.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Semi-analytical covariance matrices for two-point correlation function for DESI 2024 data

We present an optimized way of producing the fast semi-analytical covariance matrices for the Legendre moments of the two-point correlation function, taking into account survey geometry and mimicking the non-Gaussian effects. We validate the approach on simulated (mock) catalogs for different galaxy types, representative of the Dark Energy Spectroscopic Instrument (DESI) Data Release 1, used in 2024 analyses. We find only a few percent differences between the mock sample covariance matrix and our results, which can be expected given the approximate nature of the mocks, although we do identify discrepancies between the shot-noise properties of the DESI fiber assignment algorithm and the faster approximation (emulator) used in the mocks. Importantly, we find a close agreement (≤ 8% relative differences) in the projected errorbars for distance scale parameters for the baryon acoustic oscillation measurements. This confirms our method as an attractive alternative to simulation-based covariance matrices, especially for non-standard models or galaxy sample selections, making it particularly relevant to the broad current and future analyses of DESI data.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

PUMA:POWDER UTILIZATION MODELING APPLICATION

SF-25-084 PUMA a high performance modeling framework to simulate powder processing. It provides a scalable tool for manufacturers to simulate powder pre- and post-processing. The tool can predict the distortion, residual stress, and (for reactive processes) reaction completion fraction of complex parts after curing/debinding, sintering, and infiltration processes. These predictions are key metrics industry uses to optimize these processes to produce dense, defect-free, stable components.

HU, TIANCHEN (GARY)↗
DOE OSTI2024

AEOLUS: Advances in Experimental Design, Optimal Control, and Learning for Uncertain Complex Systems

The AEOLUS Center is dedicated to developing a unified optimization-under-uncertainty framework for (1) learning predictive models from data and (2) optimizing experiments, processes, and designs governed by these models, all driven by complex, uncertain energy systems. AEOLUS addressed the critical need for principled, rigorous, scalable, and structure-exploiting capabilities for exploring parameter and decision spaces of complex forward simulation models---the so-called outer loop. This report summarizes the work done under DE-SC0021077 on (1) nonlocal models for solidification problems, (2) a multifidelity method for a nonlocal diffusion model, and (3) multifidelity Monte Carlo methods.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Comparative Analysis of HEATNETS for Geothermal Network Performance

Thermal energy networks (TENs), also known as 5th generation district energy systems, or more specifically geothermal networks when exchanging heat with geothermal boreholes, are an important technology for decarbonization. In these networks an ambient loop connects buildings and thermal sources, such as a borehole field, to exchange energy and maintain a desired loop temperature. Water-source heat pumps are used at the buildings to connect to the ambient or thermal loop to meet to the building heating and cooling loads and maintain comfort. A semi-transient, reduced-order technical model and techno-economic model, called HEATNETS, has been developed at NREL that captures the flow of energy around a TEN. In this work, a comparison of the HEATNETS technical model and a well-known coding platform used for modeling geothermal networks, TRNSYS, has been completed for a proposed geothermal network as a verification and validation process. Hourly data provided from the TRNSYS simulation included building loads, pumping power, heat pump power, temperature entering and leaving the borehole field, and mass flow rates. The hourly borehole temperatures were used to create a linear regression model utilized in HEATNETS to estimate the borehole field heat exchange. The building loads and mass flow rates were direct inputs to HEATNETS while the pumping power, heat pump power, borehole temperatures, and coefficients of performance were all simulated and calculated by HEATNETS, allowing for direct comparison of the thermal energy transfer HEATNETS considers the full process from design inputs to economic outputs and can provide modeling options for high-level initial system design and operational optimization. This study focuses on a validation of HEATNETS using results from TRNSYS. HEATNETS is not intended to replace other modeling tools, but this work demonstrates, via a comparison with an industry standard code, that HEATNETS can be a unique, high-level and rapid modeling tool for estimating the performance of a full geothermal network system.

15 GEOTHERMAL ENERGY↗
DOE OSTI2026

Hydrogen Dispersion Modeling for Development of Smart Distributed Monitoring

Studying hydrogen dispersion is crucial for ensuring the safe and effective deployment of hydrogen as an energy carrier. This study presents a comprehensive CFD modeling framework for simulating hydrogen dispersion at a real-world hydrogen production, storage, and utilization facility. Utilizing the Hydrogen Research Facility under the Advanced Research on Integrated Energy Systems (ARIES) at the National Renewable Energy Laboratory's (NREL) Flatirons campus, controlled hydrogen releases at 27 kg-H2/hr were simulated. The model incorporated site-specific atmospheric conditions, including hourly wind speeds and temperatures recorded between 8 AM and 8 PM from October to December 2023. To reduce computational demands, a statistical reduction technique was applied to condense the dataset to 100 representative scenarios, validated by statistical tests for wind speeds and power law coefficients. Simulations were conducted using the Reynolds-Averaged Navier-Stokes equations. Results demonstrated that wind speed substantially influences hydrogen dispersion, with low wind conditions forming concentrated clouds and higher wind speeds stretching the plume. Additionally, clustering analysis informed optimal sensor placement at various elevations with up to 10 sensor locations on each elevation. This framework offers a robust approach for understanding hydrogen behavior in ambient conditions and informing detection strategies.

08 HYDROGEN↗
DOE OSTI2025

Analysis of the Trusted Inertial Terrain-Aided Navigation Measurement Function

The trusted inertial terrain-aided navigation (TITAN) algorithm leverages an airborne vertical synthetic aperture radar to measure the range to the closest ground points along several prescribed iso-Doppler contours. These TITAN minimum-range, prescribed-Doppler measurements are the result of a constrained nonlinear optimization problem whose optimization function and constraints both depend on the radar position and velocity. Owing to the complexity of this measurement definition, analysis of the TITAN algorithm is lacking in prior work. This publication offers such an analysis, making the following three contributions: (1) an analytical solution to the TITAN constrained optimization measurement problem, (2) a derivation of the TITAN measurement function Jacobian, and (3) a derivation of the Cramér-Rao lower bound on the estimated position and velocity error covariance. These three contributions are verified via Monte Carlo simulations over synthetic terrain, which further reveal two remarkable properties of the TITAN algorithm: (1) the along-track positioning errors tend to be smaller than the cross-track positioning errors, and (2) the cross-track positioning errors are independent of the terrain roughness.

TITAN↗
DOE OSTI2025

Error field measurements with rotating RMP fields for DIII-D H-mode

3D magnetic sensors are employed to identify the amplitude and toroidal phase of error fields (EF) by analyzing the torque balance for magnetic islands entrained by rotating resonant magnetic perturbations (RMPs) in DIII-D H-mode plasmas. This technique of torque balance allows for efficient error field identification, offering a valuable tool for scenario-specific and optimized error field compensation (EFC) and requires only magnetic diagnostics. The torque balance used in this work includes the contribution from electromagnetic torque due to error fields, wall response, RMP fields, and viscous. Results show that viscous torque plays a crucial role, particularly during locked modes and H-mode plasmas, ensuring accurate data fits with lower residuals. The torque balance technique reveals that the L- and H-mode plasmas have distinct EF configurations, and consistent with a model-based EF analysis including MHD response in IPEC and the SURFMN EF simulation. This technique shows great robustness in measuring the intrinsic EF amplitude regardless of its amplitude or toroidal phase. Repeated discharges with EFC disparities exhibit consistent results of intrinsic error field within a reasonable range near the “standard” error field compensation. Additionally, the use of a rotating n = 1 resonant magnetic perturbation offers the advantage of reducing disruption risks by entraining saturated magnetic islands. These findings are instrumental for optimizing EF correction in fusion devices, thereby enhancing tearing mode suppression and overall plasma stability.

3D magnetic sensors↗
DOE OSTI2025

Revealing the Crystallization Pathways of Mixed‐Halide Low‐Dimensional Perovskites: A First Step Toward Solar Cell Applications

Ruddlesden–Popper perovskites (RPPs) are promising materials for optoelectronic devices. While iodide‐based RPPs are well‐studied, the crystallization of mixed‐halide RPPs remains less explored. Understanding the factors affecting their formation and crystallization are vital for optimizing morphology, phase purity, and orientation, which directly impact device performance. Here, we investigate the crystallization and properties of mixed‐halide RPPs (PEA) 2 FA n−1 Pb n (Br 1/3 I 2/3 ) 3n + 1 (PEA = C 6 H 5 (CH 2 ) 2 NH 3 + and FA = CH(NH 2 ) 2 + ) (n = 1, 5, 10) using DMSO ((CH 3 ) 2 SO) or NMP (OC 4 H 6 NCH 3 ) as cosolvents and MACl (MA = CH 3 NH 3 + ) as an additive. For the first time, the presence of planar defects in RPPs is directly observed by in situ grazing‐incidence wide‐angle X‐ray scattering (GIWAXS) and confirmed through the simulation of the patterns that matched the experimental. GIWAXS data also reveals that DMSO promotes higher crystallinity and vertical orientation, while MACl enhances crystal quality but increases halide segregation, shown here by nano X‐ray fluorescence (nano‐XRF) experiments. For low‐n RPPs, orientation is crucial for solar cell efficiency, but its impact decreases with increasing n. Our findings provide insights into optimizing mixed‐halide RPPs, guiding strategies to improve crystallization, phase control, and orientation for better performance not only in solar cells but also in other potential optoelectronic devices.

Guaita, Maria G D↗
DOE OSTI2025

Efficient BiVO 4 /CoFeO x H y photoanodes using controlled annealing and conformal linear-sweep electrocatalyst photodeposition

Although monoclinic bismuth vanadate (BiVO 4 ) is a promising photoanode for solar water splitting, its practical use is hindered by imperfect photocurrent generation/collection, low photovoltage compared to the bandgap, and corrosion side reactions that limit durability. Here, we introduce a controlled-annealing sol–gel process for BiVO 4 thin-film photoanodes along with an optimized linear-sweep-voltammetry photodeposition of CoFeO x H y cocatalysts. The resulting BiVO 4 films annealed at 550 °C exhibited a photocurrent density of 4.1 mA/cm 2 at 1.23 V RHE under 1 sun AM 1.5G solar simulation and a low onset potential of 0.26 V RHE due to high majority carrier conductivity, a crystalline bulk with reduced defects as evidenced by x-ray photoelectron spectroscopy and photoluminescence lifetime analysis, and thus enhanced photocarrier collection. However, significant degradation in performance was found due to interfacial photocorrosion. To protect the surface and speed the oxygen-evolution reaction CoFeO x H y cocatalyst layers were deposited. By varying the number of consecutive sweeps and adjusting the applied bias range, an ultra-thin (~15 nm) CoFeO x H y cocatalyst layer was uniformly grown deposited over 30 cycles on the BiVO 4 surface. The resulting BiVO 4 /CoFeO x H y yielded 4.03 mA/cm 2 at 1.23 V RHE and onset potential of 0.24 V RHE , with stable operation (~15 % loss in photocurrent at 1.23 V RHE relative to ~60 % loss in the uncatalyzed control sample). These conformal CoFeO x H y catalytic layers function simultaneously to selectively collect photoexcited holes from the BiVO 4 , catalyze the water-oxidation reaction, and protect the BiVO 4 from photodegradation.

42 ENGINEERING↗
DOE OSTI2025

gRASPA

GPU Monte Carlo Simulation Code with a taste of RASPA We present enhancements in Monte Carlo simulation speed and functionality within an open-source code, gRASPA, which uses graphical processing units (GPUs) to achieve significant performance improvements compared to serial, CPU implementations of Monte Carlo. The code supports a wide range of Monte Carlo simulations, including canonical ensemble (NVT), grand canonical, NVT Gibbs, Widom test particle insertions, and continuous-fractional component Monte Carlo. Implementation of grand canonical transition matrix Monte Carlo (GC-TMMC) and a novel feature to allow different moves for the different components of metal-organic framework (MOF) structures exemplify the capabilities of gRASPA for precise free energy calculations and enhanced adsorption studies, respectively. The introduction of a High-Throughput Computing (HTC) mode permits many Monte Carlo simulations on a single GPU device for accelerated materials discovery. The code can incorporate machine learning (ML) potentials. The open-source nature of gRASPA promotes reproducibility and openness in science, and users may add features to the code and optimize it for their own purposes. The code is written in CUDA/C++ and SYCL/C++ to support different GPU vendors. The gRASPA code is publicly available at https://github.com/snurr-group/gRASPA.

Li, Zhao [Purdue/Northwestern/Notre Dame Universit↗
DOE OSTI2025

Observed declines in leaf nitrogen explained by photosynthetic acclimation to CO 2

Widespread evidence of decreasing leaf nutrients has raised concerns about ecosystem productivity under global change. Interpreting trends in leaf nutrients has important implications for the fate of ecosystem services, particularly the role of forests in mitigating climate change and sustaining quality food sources. Here, we challenge the common interpretation that decreasing leaf nitrogen concentration (LNC) is evidence of increasing nutrient limitations on ecosystem primary productivity. Instead, we show that declines in LNC (4% decrease per 50 ppm CO 2 increase), observed across 409 European forest plots over 22 y, can be explained by reduced photosynthetic nitrogen demand. This regional trend is consistent with leaf acclimation to increasing atmospheric CO 2 according to optimality theory. This finding suggests that enhanced photosynthetic nitrogen use efficiency due to CO 2 fertilization may lead to less nitrogen uptake and/or reallocation of nitrogen for plant growth and other functions. Our results have large implications for understanding and simulating interactions between ecosystem nitrogen and carbon cycles and suggest nitrogen requirements for terrestrial carbon uptake under elevated CO 2 may be lower than previously thought.

CO2 fertilization↗
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 OSTI2025

A Computational Procedure for Assessing I$_{c}(\varepsilon$) in Nb$_{3}$Sn/Bi-2212 Hybrid Magnets

The critical current of superconductors is commonly measured by testing unloaded wires under an external magnetic field. While stressed by intense Lorentz forces, the existing HTS/LTS superconductors are prone to a reduction in critical current before reaching their structural mechanical limit. Here, in this work, the magnetic and mechanical analysis of the FNAL 4-layer Bi-2212/Nb$_{3}$Sn hybrid dipole magnet is reported, aimed at predicting the critical current degradation for both the superconductors during powering at 16T. All the Rutherford cables in the coils of the hybrid magnet were modeled at the strand level in Ansys APDL with the heterogeneous cable model. Utilizing this detailed geometry, it was possible to evaluate the effects of strain on the critical current degradation for both the Nb$_{3}$Sn and Bi-2212 superconductors under the intense Lorentz forces. The analysis presented in this paper integrates strain-dependent critical current laws, with parameters derived from experimental data, to simulate the hybrid magnet's performance for all possible current-powering configurations. The proposed methodology enables a detailed assessment of conductor integrity and I$_{C}(\varepsilon$) reduction in existing hybrid magnet designs, providing a versatile and rigorous framework for optimizing future high-field hybrid magnets.

D'Agliano, A. [Lawrence Berkeley National Laborato↗
DOE OSTI2026

Studying the Impact of Momentary Cessation of IBRs Compliant With IEEE Std. 2800-2022 on Transmission Line Protection Elements

This paper examines the impact of inverter-based resource (IBR) momentary cessation timing on transmission system protection relays. IEEE Std. 2800 mandates voltage ride-through, requiring IBRs to remain online during voltage disturbances. However, for severe voltage dips (|V| = 0.1 p.u.), the standard permits current blocking. In practice, vendors implement momentary cessation with varying delays, potentially affecting relay operation if IBRs cease too early. This article studies the impact of the time IBRs enter momentary cessation on transmission line protective relays. Using an electromagnetic transient (EMT) simulation, this study models IBR momentary cessation with tunable delays in a realistic transmission network. COMTRADE files generated from EMT simulations are analyzed in MATLAB using detailed protective relay models. Results indicate that IBRs must sustain operation for at least one cycle to ensure proper relay response; otherwise, relays may fail to oper- ate. These findings guide IBR vendors in optimizing momentary cessation settings to enhance transmission line protection.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2026

Searching for the Most Harmful Field Errors in the HSR IR Superconducting Magnets

In this project, we improve beam stability for the Electron-Ion Collider. Magnetic field errors can reduce beam stability, making it essential to identify the field errors that have the greatest impact on accelerator performance. However, this is particularly challenging because beam stability depends on the complex interactions of many magnetic field errors, resulting in a high-dimensional and nonlinear optimization problem. We determine which field errors are the most influential for the large physical aperture superconducting magnet B2PF, a critical magnet in the Interaction Region (IR) in the Hadron Storage Ring (HSR). We complete and analyze nearly 30,000 simulations on the Brookhaven National Laboratory Linux Cluster by varying 18 nonlinear magnetic field errors. We evaluate beam stability using the dynamic aperture and the tune diffusion. We identify the field errors that most strongly influence beam stability and establish quantitative field error tolerances that improve accelerator performance.

43 PARTICLE ACCELERATORS↗