Efficient inverse design optimization through multi-fidelity simulations, machine learning, and boundary refinement strategies
Not Available
97 MATHEMATICS AND COMPUTING↗
SEARCH · Search NASA
Results for “differentiable simulation”
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.
Improvements of Nuclear Data Evaluations for Lead Isotopes in Support of Next Generation Lead-Cooled Fast Systems
The neutron evaluation of the isotopes that comprise natural lead were undertaken as a part of DOE-NEUP Project #19-16739. The goal of the project was to update the neutron cross sections to account for new differential measurements and incorporate the most up-to-date physics. Shortcomings in the lead cross sections was made known by several independent reports. The work performed here repeated the simulation of all the “benchmark” validation systems and concluded that the major issue is the scattering cross sections in the major lead isotopes above 100 keV. Re-evaluation of 206,207,208 Pb included both the resolved resonance region and fast region evaluations of the cross sections. Combined these regions cover energies from thermal to 20 MeV. The most drastic improvement is the resolved resonance region evaluation of 208 Pb which is now extended from 1.0 to 1.5 MeV. Parameterization of these resonances in the R-matrix provides a superior reconstruction of not only the experimental cross section but also the scattering distributions via the Blatt-Biedenharn formalism. Fast region evaluations of the three major isotopes were done to include new inelastic experimental data from the neutron Time-of-Flight facility at CERN. The culmination of all the changes to the cross section is a drastic improvement in the scattering kernel as shown in Rensselaer Polytechnic Institute (RPI) Quasi-Differential scattering measurements and improved prediction of keff for fast integral experiments. Alongside the new cross sections, new nuclear data covariance (uncertainties) have been computed and are included in the evaluation. Little is changed in the magnitude of the uncertainties but the correlations within the covariance display non-trivial changes. The new evaluations of 206,207,208 Pb have been submitted to the National Nuclear Data Center to be included in the ENDF/B-VIII.1 library. While the cross sections have been updated extensively, knowledge and modeling of the cross sections between 1.0 and 3.0 MeV for the isotopes remain a challenge. Most notably is the double differential elastic cross section for 208 Pb and 206 Pb. New quasi-differential measurements for pure samples of these two nuclei would go a long way in reducing compensating errors in the evaluations. Qualitatively, the project has produced a prosperous collaboration between the nuclear data group at RPI with staff scientists at Brookhaven National Laboratory, Naval Nuclear Laboratory, Oak Ridge National Laboratory, Los Alamos National Laboratory, Lawrence Livermore National Laboratory, and Sandia National Laboratories. Quantitatively this is reflected in seven conference presentations, at least one journal submission, and one doctoral thesis.
73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
Measurement of the muon neutrino charged-current mesonless cross section in the NOvA near detector
NOvA is a long-baseline accelerator neutrino experiment at Fermilab. Its physics goals include precision neutrino oscillation measurements, neutrino interaction cross-section measurements and beyond Standard Model explorations. We present a measurement of muon neutrino charged-current cross section with zero mesons in the final state at the NOvA near detector. This measurement is performed as a function of the kinematics of the final state muon. Our chosen interaction channel is especially sensitive to quasielastic and meson exchange current interactions and it provides handles for constraining the cross section systematic uncertainties in oscillation analyses in present and future experiments. For particle identification, we use a convolutional neural network (CNN) trained on individual particles simulated in the NOvA Near detector. This allows us to select the desired signal while reducing the potential bias from neutrino interaction modeling. We study strategies for constraining the remaining charged-pion background via Michel electron information in a template fitting approach. The main experimental result is a two-dimensional differential cross section as a function of final-state muon kinetic energy and polar angle. The parameters of this measurement, including binning and unfolding, were optimized to reduce the expected systematic uncertainty in the total cross section. The final result shows good agreement with the main GENIE-based simulation framework that was independently fine-tuned in NOvA. We finally propose improvements and subsequent steps that build on this analysis and further dissect the final states of neutrino interactions. This work has been supported by US DOE grant DE-SC0015684.
Sánchez Falero, Sebastián Jesús [Iowa State U.]↗
Project Report: High Frequency Axion LDRD - A Multiport Approach to Thermal Noise and Scattering Parameter Simulation of Cryogenic Experiments
In this report, a simple algorithm for detailed system-level thermal noise analysis is developed, demonstrated, and verified. This method uses noise-wave theory and noise covariance matrices to cascade noise and scattering parameters of multiport devices at different temperatures. This method addresses the effects of return loss, multiport isolation/coupling, and static temperature differentials between components, and will work in cases where the noise temperature is at or near the quantum noise limit. An ideal multidevice network will first be demonstrated to show that this method’s results are consistent with the Friis cascade when component parameters such as return loss and isolation are ideal. Following the ideal multidevice example, a cryogenic experiment is conducted to demonstrate that the proposed simulation method is successful when real data are used.
42 ENGINEERING↗
Designing observables for measurements with deep learning
Many analyses in particle and nuclear physics use simulations to infer fundamental, effective, or phenomenological parameters of the underlying physics models. When the inference is performed with unfolded cross sections, the observables are designed using physics intuition and heuristics. We propose to design targeted observables with machine learning. Unfolded, differential cross sections in a neural network output contain the most information about parameters of interest and can be well-measured by construction. The networks are trained using a custom loss function that rewards outputs that are sensitive to the parameter(s) of interest while simultaneously penalizing outputs that are different between particle-level and detector-level (to minimize detector distortions). We demonstrate this idea in simulation using two physics models for inclusive measurements in deep inelastic scattering. We find that the new approach is more sensitive than classical observables at distinguishing the two models and also has a reduced unfolding uncertainty due to the reduced detector distortions.
72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
The Computational Analysis of a Generic Hypersonic Store Separation
Many current hypersonic vehicles involve a main vehicle from which components separate during flight and presents simulation challenges such as shock-shock interactions, highly nonlinear interactions, and shock wave boundary layer interactions. These aerodynamic phenomena influence the attitude and trajectory of each object, which needs to be understood before large-scale experiments can be run. To gain a generic understanding of the process, high-fidelity Reynolds-averaged Navier-Stokes solutions are performed on a conical vehicle geometry passing through an oblique shock wave that is representative of different configurations and scenarios. Trajectory and applied forces are tracked and show the vehicle’s dynamics are predominantly the result of differential flow incidence angles causing a strong shock; this leads to a large pressure increase over a fraction of the vehicle, which influences pitch. The vehicle appears to follow conventional stability theory with detached eddy simulation and small disturbances in initial attitude shown to have minimal influence on the scenario.
Savery, Ryan↗
Convective carbon dioxide dissolution in a closed porous medium at high-pressure real-gas conditions
Not provided.
Water Resources↗
Joint Resource Modeling and Assessment for Hybrid Distributed Solar and Wind Systems
The inherent variability and uncertainty in distributed energy resources can presents myriad challenges to the planning and operations of power systems. These risks are poised to become larger as the penetration of renewable energy sources rises in the power generation mix. Hybrid solar-wind energy systems are able to mitigate some of these risks by their complementary resource availability. Surface solar and wind fields are coupled and correlated in both space and time. Appropriately estimating the hybrid solar wind energy system requires simulating the spatio-temporal structure of these fields that can be produced for each time horizon. We introduce a novel joint spatio-temporal stochastic differential equation (SPDE) approach that captures the spatio-temporal dynamics of solar and wind fields and their joint dependency over a domain for each time step. In the case study on Colorado, we consider nonstationary three-level hierarchical spatio temporal models for both hourly solar irradiance data and wind speed data in Colorado. Dependence between the solar irradiance data and wind speed data is captured by a shared spatio-temporal random effect. Our approach performs well in terms of the prediction score criterion.
joint modeling↗
A Bayesian nonparametric analysis for zero-inflated multivariate count data with application to microbiome study
High-throughput sequencing technology has enabled researchers to profile microbial communities from a variety of environments, but analysis of multivariate taxon count data remains challenging. Here, we develop a Bayesian nonparametric (BNP) regression model with zero inflation to analyse multivariate count data from microbiome studies. A BNP approach flexibly models microbial associations with covariates, such as environmental factors and clinical characteristics. The model produces estimates for probability distributions which relate microbial diversity and differential abundance to covariates, and facilitates community comparisons beyond those provided by simple statistical tests. We compare the model to simpler models and popular alternatives in simulation studies, showing, in addition to these additional community-level insights, it yields superior parameter estimates and model fit in various settings. The model's utility is demonstrated by applying it to a chronic wound microbiome data set and a Human Microbiome Project data set, where it is used to compare microbial communities present in different environments.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Mathematical Foundations for Nonlocal Interface Problems: Multiscale Simulations of Heterogeneous Materials (Final LDRD Report)
Nonlocal models provide a much-needed predictive capability for important Sandia mission applications, ranging from fracture mechanics for nuclear components to subsurface flow for nuclear waste disposal, where traditional partial differential equations (PDEs) models fail to capture effects due to long-range forces at the microscale and mesoscale. However, utilization of this capability is seriously compromised by the lack of a rigorous nonlocal interface theory, required for both application and efficient solution of nonlocal models. To unlock the full potential of nonlocal modeling we developed a mathematically rigorous and physically consistent interface theory and demonstrate its scope in mission-relevant exemplar problems.
97 MATHEMATICS AND COMPUTING↗
Simulating single-particle dynamics in magnetized plasmas: The RMF code
The RMF (Rotating Magnetic Field) code is designed to calculate the motion of a charged particle in a given electromagnetic field. It integrates Hamilton’s equations in cylindrical coordinates using an adaptive predictor-corrector double-precision variable-coefficient ordinary differential equation solver for speed and accuracy. RMF has multiple capabilities for the field. Particle motion is initialized by specifying the position and velocity vectors. Here, the six-dimensional state vector and derived quantities are saved as functions of time. A post-processing graphics code, XDRAW, is used on the stored output to plot up to 12 windows of any two quantities using different colors to denote successive time intervals. Multiple cases of RMF may be run in parallel and perform data mining on the results. Recent features are a synthetic diagnostic for simulating the observations of charge-exchange-neutral energy distributions and RF grids to explore a Fermi acceleration parallel to static magnetic fields.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Dynamic Simulation Modeling and Control of a Desiccant Assisted Direct-expansion Air Handling Unit
Desirable built environments demand simultaneous regulation of thermal comfort and indoor air quality (IAQ) with energy-efficient operation of heating, ventilation and air conditioning (HV AC) systems, which involves controls of temperature, humidity and airborne contaminants simultaneously. This paper presents the efforts of dynamic modeling and initial development control strategy for a desiccant-assisted multi-functional air handling unit (AHU) coupled with a direct-expansion rooftop unit (RTU) system, which aims to achieve multiple functions for indoor environment conditioning with energy efficient control. The RTU-AHU system includes a desiccant wheel for dehumidification and a conceptual direct air capture (DAC) filtering device for CO2 regulation. A Modelica-based dynamic model is developed for this conceptual system, and a simple decentralized control strategy is designed, which combines a differential-enthalpy based AHU return-air ratio control, a demand-controlled ventilation, and supply-air temperature humidity control via the RTU and DW controls. The proposed control method is evaluated with the Modelica simulation model for a selected set of scenarios
Pan, Chao↗
High-Order Mesh r-Adaptivity with Tangential Relaxation and Guaranteed Mesh Validity
High-order meshes are crucial for achieving optimal convergence rates in curvilinear domains, preserving symmetry, and aligning with key flow features in moving mesh simulations [1], but their quality is challenging to control. In prior work, we have developed techniques based on Target-Matrix Optimization Paradigm (TMOP) to adapt a given high-order mesh to the geometry and solution of the partial differential equation (PDE) [2, 3]. Here, we extend this framework to address two key gaps in the literature for highorder mesh 𝑟-adaptivity. First, we introduce tangential relaxation on curved surfaces using solely the discrete mesh representation, eliminating the need for access to underlying geometry (e.g., CAD model). Second, we ensure a continuously positive Jacobian determinant throughout the domain. This determinant positivity is essential for using the high-order mesh resulting from 𝑟-adaptivity with arbitrary quadrature schemes in simulations. The proposed approach is demonstrated to be robust using a variety of numerical experiments.
Mathematics and Computing↗
Evaluation of current density measurement methods for high-current, low-beta electron beams
Recent current density measurements with a low-beta electron beam (β=0.5-0.75) are strongly affected by electron scatter and Cherenkov limits. These invasive measurements are heavily influenced by the electron energy and the intercepting material properties: density, index of refraction, and atomic number. We evaluate various measurement methods, including x-ray scintillation and Cherenkov emission. To optimize the x-ray scintillation technique for this energy regime, metal foils of differing atomic number (Z) and thickness were placed upstream of a scintillation screen. We have selected Cherenkov emitters, based on a material refractive index, in order to differentiate the contributions from electron scatter and produced fluorescence, Cherenkov emission, and total internal reflection. The observed distributions for both approaches are verified through MCNP6®. Additionally, the simulations provide insight as to what physical processes dominate the measured distribution. We conclude by determining the limits and optimal measurement range for each measurement technique.
72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
Measurement of the Ar ( e , e ' p ) and Ti ( e , e ' p ) cross sections in Jefferson Lab Hall A
The E12-14-012 experiment, performed in Jefferson Lab Hall A, has collected exclusive electron-scattering data ( e , e ' p ) in parallel kinematics using natural argon and natural titanium targets. In this study we report the first results of the analysis of the data set corresponding to beam energy 2222 GeV, electron scattering angle 21 . 5 ° , and proton emission angle – 50 ° . The differential cross sections, measured with ≈ 4 % uncertainty, have been studied as a function of missing energy and missing momentum, and compared to the results of Monte Carlo simulations, obtained from a model based on the distorted-wave impulse approximation.
73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
Binary Quantum Control Optimization with Uncertain Hamiltonians
Optimizing the controls of quantum systems plays a crucial role in advancing quantum technologies. The time-varying noises in quantum systems and the widespread use of inhomogeneous quantum ensembles raise the need for high-quality quantum controls under uncertainties. In this paper, we consider a stochastic discrete optimization formulation of a discretized binary optimal quantum control problem involving Hamiltonians with predictable uncertainties. We propose a sample-based reformulation that optimizes both risk-neutral and risk-averse measurements of control policies, and solve these with two gradient-based algorithms using sum-up-rounding approaches. Furthermore, we discuss the differentiability of the objective function and prove upper bounds of the gaps between the optimal solutions to binary control problems and their continuous relaxations. We conduct numerical simulations on various sized problem instances based on two applications of quantum pulse optimization; we evaluate different strategies to mitigate the impact of uncertainties in quantum systems. In conclusion, we demonstrate that the controls of our stochastic optimization model achieve significantly higher quality and robustness compared with the controls of a deterministic model.
conditional value-at-risk (CVaR)↗
A Robust Segmented Mixed Effect Regression Model for Baseline Electricity Consumption Forecasting
Renewable energy production has been surging around the world in recent years. To mitigate the increasing uncertainty and intermittency of the renewable generation, proactive demand response algorithms and programs are proposed and developed to further improve the utilization of load flexibility and increase the efficiency of power system operation. One of the biggest challenges to efficient control and operation of demand response resources is how to forecast the baseline electricity consumption and estimate the load impact from demand response resources accurately. In this paper, we propose a mixed effect segmented regression model and a new robust estimate for forecasting the baseline electricity consumption in Southern California, USA, by combining the ideas of random effect regression model, segmented regression model, and the least trimmed squares estimate. Since the log-likelihood of the considered model is not differentiable at breakpoints, we propose a new backfitting algorithm to estimate the unknown parameters. The estimation performance of the new estimation procedure has been demonstrated with both simulation studies and the real data application for the electric load baseline forecasting in Southern California.
42 ENGINEERING↗
Bulk Spray and Individual Plume Characterization of LPG and Iso-Octane Sprays at Engine-Like Conditions
This study presents experimental and numerical examination of directly injected (DI) propane and iso-octane, surrogates for liquified petroleum gas (LPG) and gasoline, respectively, at various engine like conditions with the overall objective to establish the baseline with regards to fuel delivery required for future high efficiency DI-LPG fueled heavy-duty engines. Sprays for both iso-octane and propane were characterized and the results from the optical diagnostic techniques including high-speed Schlieren and planar Mie scattering imaging were applied to differentiate the liquid-phase regions and the bulk spray phenomenon from single plume behaviors. The experimental results, coupled with high-fidelity internal nozzle-flow simulations were then used to define best practices in CFD Lagrangian spray models. Optical imaging revealed that unlike iso-octane, propane’s spray propagation was fed by its flash boiling, spray collapse, and high degree of vaporization, resulting in a direct proportionality of propane’s penetration length to temperature. These unique features of propane and its variation from iso-octane’s spray pattern, contributed to its classification as an unconventional spray. Appropriate corrections to the injection and breakup models were developed to reproduce the under-expanded jet dynamics and to mimic the flash boiling-driven spray development observed with propane sprays. The simulation results were found to be sensitive to cone and inclusion angles of the blob injector. The current work represents a first assessment of the capability of the commonly available models for engine-spray simulations and highlights the fact that, despite the reasonable agreement obtained in the fuel vapor morphology, the representation of the liquid phase lacks accuracy and requires further model development.
02 PETROLEUM↗