Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Simulation Quality”

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 595 records · Page 33

Performance Understanding and Analysis for Exascale Data Management Workflows (Collaboration)

The general approach of the MONA project is depicted in Figure 1. The figure shows performance monitoring applied to an I/O workflow associated with a scientific simulation, with online measurements captured from this workflow informing methods for workflow reconfiguration and dynamic adjustment. The goal is to maintain suitable levels of Quality of Service for workflow execution, by understanding the underlying causes of workflow performance and behavior. One outcome will be workflow performance models able to characterize realistic workflows. Another outcome will be performance ‘mini-apps’ implementing such workflow behavior. Out of scope for this project are advanced methods for online workflow control.

97 MATHEMATICS AND COMPUTING↗

Identification and Photometric Classification of Extragalactic Transients in the Vera C. Rubin Observatory’s Data Preview 1

The Vera C. Rubin Observatory will soon survey the southern sky, delivering a depth and sky coverage that is unprecedented in time-domain astronomy. As part of commissioning, Data Preview 1 (DP1) has been released. It comprises a Legacy Survey of Space and Time (LSST) Commissioning Camera observing campaign between 2024 November and December with multiband imaging of seven fields, covering roughly 0.4 deg 2 each, providing a first glimpse into the data products that will become available once the LSST begins. In this work, we search three fields for extragalactic transients. We identify eight new likely supernovae (SNe), and three known ones from a sample of 369,644 difference image analysis objects. Photometric classification using Superphot+ assigns subclasses with >95% confidence to only one SN Ia and one SN II in this sample. Our findings are in agreement with SN detection rate predictions of 15 ± 4 SNe from simulations using simsurvey. The SN detection rate in the data is possibly affected by the lack of suitable templates. Nevertheless, this work demonstrates the quality of the data products delivered in DP1 and indicates that the Rubin Observatory’s LSST is well placed to fulfill its discovery potential in time-domain astronomy.

Freeburn, James [University of North Carolina, Cha↗

Systems, methods and computer program products for collaborative agent control

Systems, methods and unmanned agents for collaboratively controlling agents in a collaborative network by one or more agents continuously simulating numeric models of one or more other agents in the network to dramatically reduce the computational bandwidth required between agents, and improve the quality of shared estimates of the agent locations as well as the locations and characteristics of other objects of interest, e.g. targets. Bandwidth is reduced by using the models to intelligently filter data before communicating.

Buerger, Stephen↗

Active Learning Framework

Machine learning (ML) of interatomic potentials show great promise to accelerate scientific simulation, e.g., by emulating expensive computations at a high accuracy but much reduced computational cost. Training datasets are calculated from computationally expensive ab initio quantum mechanics methods, density functional theory (DFT). Trained on this data, an ML model can be very successful in predicting energy and forces for new atomic configurations. A critical factor is the quality and diversity of the training dataset. Thus, a highly automated approach to dataset construction based on active learning framework is designed suitable for material physics. The active learning scheme begins with fully randomized atomic configurations. Then, many Molecular Dynamics (MD) trajectories are simulated using current ML potentials, where each MD trajectory is initialized to a random disordered configuration. The temperature is varied in order to diversify the sampled configuration during these simulations. The variance of predictions for eight neural networks within an ensemble is analyzed to determine whether the model is operating as expected. This helps in determining whether collecting more data would be helpful to the model by checking the ensemble variance is greater than the threshold. In this case, the MD trajectory is terminated and the final atomic configuration is placed on a queue (SQL database) for DFT calculations and added to training dataset. Periodically, ML model is retrained to the updated training model. This Active Learning loop is iterated until the cost of MD simulations becomes prohibitively expensive. The MD simulations will hopefully be sufficiently robust to support nucleation after many active learning iterations. In this sense, active learning scheme must automatically discover the important low energy and nonequilibrium physics.

Nebgen, Benjamin↗

Software Quality Assurance for the MOOSE-Based Open-Source Multiphysics Code Cardinal - An Expanded CI Testing Suite

Cardinal is a wrapping of the GPU-oriented spectral element Computational Fluid Dynamics (CFD) code NekRS and the Monte Carlo particle transport code OpenMC within the Multiphysics Object-Oriented Simulation Environment (MOOSE). Cardinal provides high-resolution thermal-hydraulics and/or radiation transport feedback to MOOSE multiphysics simulations. Multiphysics feedback is implemented in a geometry-agnostic manner which eliminates the need for rigid one-to-one mappings. A generic data transfer implementation also allows NekRS and OpenMC to couple to any MOOSE application, enabling a broad set of multiphysics capabilities. Cardinal simulations can also leverage combinations of MPI, OpenMP, and GPU resources. Cardinal continuous development and improvement efforts have led to the software being considered as a high-fidelity design and licensing tool for key areas of nuclear reactor relevant physics, including neutron transport, fluid flow, heat transfer, and mechanical processes. The fast development and expansion of the software from a pure R&D framework towards its application in the nuclear industry and regulation require a focus on developing, enhancing and, maintaining Cardinal’s software quality through strict adherence to a Software Quality Assurance (SQA) framework and SQA program. To facilitate compliance with SQA standards, the Cardinal SQA Program has been initiated during Fiscal Year 2023 (FY23). During the development of the Cardinal SQA Program, multiple gaps have been identified. These gaps are primarily related to model verification and code pedigree as they relate to the use of Cardinal as a safety analysis tool. These gaps have been captured in a report published in 2023. A second report highlighted the progress made during Fiscal Year 2024 (FY24) and described Argonne’s effort to document and integrate software verification within Cardinal’s software development process. This report documents a snapshot of the verification test cases currently available for Cardinal and NekRS in their assimilation into a Continuous Integration (CI) platform. Following the CI practice permits the integrating of source code changes frequently and ensuring that the integrated codebase clears the verification testing for the software. It should be noted that the SQA program itself, including the program plans, procedures, configuration management, and testing strategies, need to be developed in a future step of this task.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Effects of residue removal and tillage on greenhouse gas emissions in continuous corn systems as simulated with RZWQM2

Agricultural production is a major source of carbon dioxide (CO 2 ) and nitrous oxide (N 2 O) globally. The effects of conservation practices on soil CO 2 and N 2 O emissions remain a high degree of uncertainty. In this study, soil CO 2 and N 2 O emissions under different residue and tillage practices in an irrigated, continuous corn system, were investigated using the Root Zone Water Quality Model (RZWQM2). Combinations of no/high stover removal (NR and HR, respectively) and no-till/conventional tillage (NT and CT, respectively) field experiments were tested over the four crop-years (Apr. 2011–Apr. 2015). The model was calibrated using the NRCT, and validated with other treatments. The simulation results showed that soil volumetric water content (VWC) in the NR treatments (i.e., NRCT and NRNT) was 1.3%–1.9% higher than that in the HR treatments (i.e., HRCT and HRNT) averaged across the four years. A higher amount of CO 2 and N 2 O emissions were simulated in the NRCT across the four years (annual average: 7034 kg C/ha/yr for CO 2 and 3.8 kg N/ha/yr for N 2 O), and lower emissions were in the HRNT (annual average: 6329 kg C/ha/yr and 3.7 kg N/ha/yr for N 2 O). A long-term simulation (2001–2015) suggested that the CO 2 and N 2 O emissions were closely correlated with the stover removal degree (SRD), tillage, VWC, soil temperature (ST), years in management (Y), and fertilizer application. Stover and tillage practices had cumulative effects on CO 2 emissions. The simulated annual CO 2 emissions in 1st year from NRCT, NRNT, and HRCT were 7.8%, 0.0%, and 7.7% higher than that from HRNT, respectively; then the emissions in 15th year were 63.6%, 47.7%, and 29.1% higher, respectively. Meanwhile, there were no cumulative effects on N 2 O emissions. Furthermore, the results also demonstrated that the RZWQM2 is a promising tool for evaluating the long-term effects of CO 2 and N 2 O emissions on different conservation practices.

54 ENVIRONMENTAL SCIENCES↗

Propagation of the 2014–2016 Northeast Pacific Marine Heatwave Through the Salish Sea

Effects and impacts of the Northeast Pacific marine heatwave of 2014–2016 on the inner coastal estuarine waters of the Salish Sea were examined using a combination of monitoring data and an established three-dimensional hydrodynamic and biogeochemical model of the region. The anomalous high temperatures reached the U.S. Pacific Northwest continental shelf toward the end of 2014 and primarily entered the Salish Sea waters through an existing strong estuarine exchange. Elevated temperatures up to + 2.3°C were observed at the monitoring stations throughout 2015 and 2016 relative to 2013 before dissipating in 2017. The hydrodynamic and biogeochemical responses to this circulating high-temperature event were examined using the Salish Sea Model over a 5-year window from 2013 to 2017. Responses of conventional water-quality indicator variables, such as temperature and salinity, nutrients and phytoplankton, zooplankton, dissolved oxygen, and pH, were evaluated relative to a baseline without the marine heatwave forcing. The simulation results relative to 2014 show an increase in biological activity (+14%, and 6% Δ phytoplankton biomass, respectively) during the peak heatwave year 2015 and 2016 propagating toward higher zooplankton biomass (+14%, +18% Δ mesozooplankton biomass). However, sensitivity tests show that this increase was a direct result of higher freshwater and associated nutrient loads accompanied by stronger estuarine exchange with the Pacific Ocean rather than warming due to the heatwave. Strong vertical circulation and mixing provided mitigation with only ≈+0.6°C domain-wide annual average temperature increase within Salish Sea, and served as a physical buffer to keep waters cooler relative to the continental shelf during the marine heatwave.

54 ENVIRONMENTAL SCIENCES↗

Structural determination of a full-length plant cellulose synthase informed by experimental and in silico methods

Three-dimensional structure determination and prediction of proteins with intrinsically disordered regions, unstructured regions, conformational flexibility, and lacking homologous structures are challenging. We previously predicted and refined an in silico structure of a plant cellulose synthase from cotton (GhCESA1), and more recently, cryo-electron microscopy (cryo-EM) has resolved a majority of the lengths of two CESA structures from poplar (PttCESA8) and cotton (GhCESA7). However, 26–30% of these cryo-EM structures remain unresolved, including the N-terminal domain, half of the class-specific region, the gating loop region, and the C-terminal domain. Here, we describe the generation and evaluation of a full-length hybrid GhCESA1 model based on this cryo-EM PttCESA8 structure, with unresolved regions completed using this in silico refined GhCESA1 model. All-atom molecular dynamics simulations and subsequent energy minimizations were performed for the in silico and hybrid GhCESA1 models in a lipid bilayer-water-ion environment, and structural stability, dynamics, energetics, contacts, and quality were evaluated. The unresolved regions were found to be the most dynamic, in agreement with their poor electron density with cryo-EM. The hybrid model exhibited a higher total secondary structure content, more favorable intra-protein and protein-lipid interaction energies, and improved quality metrics. Moreover, hydrogen bonding was revealed to be a primary mechanism for intra-protein and protein-lipid contacts. These results demonstrate that in silico structure prediction and refinement may be useful to augment experimental structure determination, especially for disordered and unstructured regions. Furthermore, this hybrid model can serve as a steppingstone to derive full-length homology models of other CESAs found in more experimentally tractable organisms.

59 BASIC BIOLOGICAL SCIENCES↗

Fast and accurate simulations of calorimeter showers with normalizing flows

In this study, we introduce caloflow, a fast detector simulation framework based on normalizing flows. For the first time, we demonstrate that normalizing flows can reproduce many-channel calorimeter showers with extremely high fidelity, providing a fresh alternative to computationally expensive geant4 simulations, as well as other state-of-the-art fast simulation frameworks based on generative adversarial networks (GANs) or variational autoencoders (VAEs). In addition to the usual histograms of physical features and images of calorimeter showers, we introduce a new metric for judging the quality of generative modeling: the performance of a classifier trained to differentiate real from generated images. We show that GAN-generated images can be identified by the classifier with nearly 100% accuracy, while images generated from caloflow are better able to fool the classifier. More broadly, normalizing flows offer several advantages compared to other state-of-the-art approaches (GANs and VAEs), including tractable likelihoods, stable and convergent training, and principled model selection. Normalizing flows also provide a bijective mapping between data and the latent space, which could have other applications beyond simulation, for example, to detector unfolding.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Deep Optical Switching on Subpicosecond Timescales in an Amorphous Ge Metamaterial

Active nanostructured optical components show promise as potential building blocks for novel light-based computing and data processing architectures. However, nanoscale all-optical switches that have low activation powers and high-contrast ultrafast switching have been elusive so far. Here, pump–probe measurements performed on amorphous-Ge-based micro-resonator metasurfaces that exhibit strong resonant modes in the mid-infrared are reported. Relative change is observed in transmittance of ΔT/T ≈ 1 with picosecond (down to τ ≈ 0.5 ps) free carrier relaxation rates, obtained with very low pump fluences of 50 μJ cm –2 . In this work, these observations are attributed to efficient free carrier promotion, affecting light transmittance via high quality-factor optical resonances, followed by an increased electron–phonon scattering of free carriers due to the amorphous crystal structure of Ge. Full-wave simulations based on a permittivity model that describes free-carrier damping through crystal structure disorder find excellent agreement with the experimental data. These findings offer an efficient and robust platform for all-optical switching at the nanoscale.

36 MATERIALS SCIENCE↗

Harmonic cancellation within AC low voltage distribution for a realistic office environment

An increase of non-linear loads, primarily from power electronics, has substantially increased current harmonics in commercial buildings, which contributes to decreased transformer efficiency / lifespan and poor power quality. This study uses recorded power consumption data from common miscellaneous electric loads (MELs) seen in offices, combined with detailed characterizations of example MELs, to simulate harmonic cancellation within building circuits. Typically, harmonic cancellation studies assume that AC converters operate across their rated power range. However, this study finds that common MELs operate below 40% of rated power the majority of the time when not quiescent; 89% of sampled devices never operated above 60% of rated power. Simulations using these more realistic power levels indicate current-harmonic cancellation (3rd to 13th harmonic) is significantly lower than that predicted when using full-range power assumptions, resulting in minor errors for low-order harmonics and larger errors for higher order harmonics. Furthermore, increased MELs load diversity increases harmonic cancellation, but insufficiently to eliminate errors. In contrast, blending lighting loads with MELs on the secondaries of distribution transformers improves harmonic cancellation to near those predicted by traditional methods. These results indicate that realistic power levels, as well as better characterization of harmonics from typical MELs, should be used to estimate harmonic cancellation.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Ab initio study of the structure and properties of amorphous silicon hydride from accelerated molecular dynamics simulations

This paper presents a large-scale ab initio simulation study of amorphous silicon hydride (a-Si 1-x H x ) with an emphasis on the structure and properties of the material across a range of hydrogen concentration by combining accelerated molecular dynamics (MD) simulations with first-principles density-functional calculations. The accelerated MD scheme relied on classical metadynamics, which enabled the development of 2500+ high-quality structural models of a-Si 1-x H x , with system sizes ranging from 150 to 6000 atoms and hydrogen concentrations vary from 6 to 20 at. %. The resulting amorphous networks were found to be completely free from any coordination defects and that they all exhibited a pristine band-gap in their electronic spectrum. The microstructural properties of hydrogen distributions were examined with an emphasis on the presence of isolated and clustered environments of hydrogen atoms. The results were compared with experimental data obtained from X-ray diffraction, infrared spectroscopy and nuclear magnetic resonance studies.

36 MATERIALS SCIENCE↗

Vertical van der Waals heterojunction diodes comprising 2D semiconductors on 3D β-Ga 2 O 3

Wide bandgap semiconductors such as gallium oxide (Ga 2 O 3 ) have attracted much attention for their use in next-generation high-power electronics. Although single-crystal Ga 2 O 3 substrates can be routinely grown from melt along various orientations, the influence of such orientations has been seldom reported. Further, making rectifying p–n diodes from Ga 2 O 3 has been difficult due to lack of p-type doping. In this study, we fabricated and optimized 2D/3D vertical diodes on β-Ga 2 O 3 by varying the following three factors: substrate planar orientation, choice of 2D material and metal contacts. Furthermore, the quality of our devices was validated using high-temperature dependent measurements, atomic-force microscopy (AFM) techniques and technology computer-aided design (TCAD) simulations. Our findings suggest that 2D/3D β-Ga 2 O 3 vertical heterojunctions are optimized by substrate planar orientation (–201), combined with 2D WS 2 exfoliated layers and Ti contacts, and show record rectification ratios (>10 6 ) concurrently with ON-Current density (>10 3 A cm –2 ) for application in power rectifiers.

36 MATERIALS SCIENCE↗

Spatiotemporally Adaptive Compression for Scientific Dataset with Feature Preservation – A Case Study on Simulation Data with Extreme Climate Events Analysis

Scientific discoveries are increasingly constrained by limited storage space and I/O capacities. For time-series simulations and experiments, their data often need to be decimated over timesteps to accommodate storage and I/O limitations. In this paper, we propose a technique that addresses storage costs while improving post-analysis accuracy through spatiotemporal adaptive, error-controlled lossy compression. We investigate the trade-off between data precision and temporal output rates, revealing that reducing data precision and increasing timestep frequency lead to more accurate analysis outcomes. Additionally, we integrate spatiotemporal feature detection with data compression and demonstrate that performing adaptive error-bounded compression in higher dimensional space enables greater compression ratios, leveraging the error propagation theory of a transformation-based compressor. To evaluate our approach, we conduct experiments using the well-known E3SM climate simulation code and apply our method to compress variables used for cyclone tracking. Our results show a significant reduction in storage size while enhancing the quality of cyclone tracking analysis, both quantitatively and qualitatively, in comparison to the prevalent timestep decimation approach. Compared to three state-of-the-art lossy compressors lacking feature preservation capabilities, our adaptive compression framework improves perfectly matched cases in TC tracking by 26.4-51.3% at medium compression ratios and by 77.3-571.1% at large compression ratios, with a merely 5–11% computational overhead.

Gong, Qian↗

Ab initio study of the structure and properties of amorphous silicon hydride from accelerated molecular dynamics simulations

This paper presents a large-scale ab initio simulation study of amorphous silicon hydride (a-Si 1-x H x ) with an emphasis on the structure and properties of the material across a range of hydrogen concentration by combining accelerated molecular dynamics (MD) simulations with first-principles density-functional calculations. The accelerated MD scheme relied on classical metadynamics, which enabled the development of 2600+ high-quality structural models of a-Si 1-x H x , with system sizes ranging from 150 to 6,000 atoms and hydrogen concentrations vary from 6 to 20 at. %. The resulting amorphous networks were found to be completely free from any coordination defects and that they all exhibited a pristine band-gap in their electronic spectrum. The microstructural properties of hydrogen distributions were examined with great emphasis on the presence of isolated and clustered environments of hydrogen atoms. The results were compared with a suite of experimental data obtained from x-ray diffraction, infrared spectroscopy, spectroscopic ellipsometry and nuclear magnetic resonance studies.

36 MATERIALS SCIENCE↗

Exploring the Intersection of AI and Visualization in the Nuclear Industry

This presentation explores the impact of AI and visualization in advancing the nuclear industry by improving safety, operational efficiency, and decision-making processes. It highlights key applications such as real-time monitoring, predictive maintenance, and immersive training, while addressing challenges like data quality, regulatory hurdles, and the need for explainable AI. Additionally, the presentation outlines future directions, emphasizing the potential of AI-driven reactor design, advanced simulation tools, and ethical considerations to drive innovation and sustainability in nuclear operations.

99 GENERAL AND MISCELLANEOUS↗

Simulating organic aerosol in Delhi with WRF-Chem using the volatility-basis-set approach: exploring model uncertainty with a Gaussian process emulator

The nature and origin of organic aerosol in the atmosphere remain unclear. The gas–particle partitioning of semi-volatile organic compounds (SVOCs) that constitute primary organic aerosols (POAs) and the multigenerational chemical aging of SVOCs are particularly poorly understood. The volatility basis set (VBS) approach, implemented in air quality models such as WRF-Chem (Weather Research and Forecasting model with Chemistry), can be a useful tool to describe emissions of POA and its chemical evolution. However, the evaluation of model uncertainty and the optimal model parameterization may be expensive to probe using only WRF-Chem simulations. Gaussian process emulators, trained on simulations from relatively few WRF-Chem simulations, are capable of reproducing model results and estimating the sources of model uncertainty within a defined range of model parameters. In this study, a WRF-Chem VBS parameterization is proposed; we then generate a perturbed parameter ensemble of 111 model runs, perturbing 10 parameters of the WRF-Chem model relating to organic aerosol emissions and the VBS oxidation reactions. This allowed us to cover the model's uncertainty space and to compare outputs from each run to aerosol mass spectrometer observations of organic aerosol concentrations and O:C ratios measured in New Delhi, India. The simulations spanned the organic aerosol concentrations measured with the aerosol mass spectrometer (AMS). However, they also highlighted potential structural errors in the model that may be related to unsuitable diurnal cycles in the emissions and/or failure to adequately represent the dynamics of the planetary boundary layer. While the structural errors prevented us from clearly identifying an optimized VBS approach in WRF-Chem, we were able to apply the emulator in the following two periods: the full period (1–29 May) and a subperiod period of 14:00–16:00 h LT (local time) on 1–29 May. The combination of emulator analysis and model evaluation metrics allowed us to identify plausible parameter combinations for the analyzed periods. We demonstrate that the methodology presented in this study can be used to determine the model uncertainty and to identify the appropriate parameter combination for the VBS approach and hence to provide valuable information to improve our understanding of OA production.

54 ENVIRONMENTAL SCIENCES↗

Effects of temperature inversion on densification in chemical vapor infiltration

In a classical chemical vapor infiltration (CVI) process, the competing effects of chemical kinetics and reagent gas transport lead to non–uniform depositions such that outer layers of a preform densify faster leaving the core highly porous. Currently, CVI must be performed at a sufficiently low temperature to achieve good densification quality which leads to high processing time and cost. Volumetric heating of the preform, especially through microwaves, can create temperature inversion such that the core is hotter than the outer surface and potentially, overcome the challenges associated with isothermal CVI. Direct numerical simulations (DNS) of densification under various such temperature distributions indicate that microwave heating in CVI processing can lead to better (uniform) densification of porous preforms. Here the role of key parameters describing the temperature distributions on the densification behavior is investigated. Strategic temporal control of the temperature distribution shows that processing times can be reduced by almost half while maintaining a good densification quality similar to that of low–temperature isothermal processing. Inside–out densification due to the inverted temperature profile is a key distinguishing characteristic of microwave assisted CVI.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗