Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Sensitivity Study”

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 577 records · Page 32
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Super‐Droplet Method to Simulate Lagrangian Microphysics of Nuclear Fallout in a Homogeneous Cloud

Abstract Nuclear detonations produce hazardous local and global particles or fallout. Predicting fallout size, chemical components, and location is necessary to inform officials and determine immediate guidance for the public. However, existing nuclear detonation fallout models prescribe the particle size distributions based on limited observations. In this work, we apply the super‐droplet method, which is a numerical modeling technique developed for cloud microphysics, to simulate size distributions of particles in a mushroom cloud formed post‐detonation of a nuclear device. We model fallout formation and evolution with homogeneous nucleation and condensation of a single species and a Monte Carlo coagulation algorithm. We verify the numerical methods representing coagulation and condensation processes against analytical test problems. Additionally, we explore several scenarios for the integral system mass and yield in equivalent kilotons (kt) of TNT (trinitrotoluene). The fallout size distribution median diameter d pg follows a scaling law based on the integral system mass m v 0 kg and yield Y kt: nm. We test the effect of cloud turbulence, enhanced nucleation and growth, and vapor volatility with a sensitivity study. The range in median diameter predictions for simulations of historical tests performed over the Pacific encompass the measurements of particles sampled from the cloud caps. Predicted median particle size ranges up to 217, 123, 86, and 35 nm for historical tests with yields of 0.2, 0.7, 2, and 10 Mt, respectively. This work can be expanded in many different directions to build a more predictive model for fallout formation.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗
DOE OSTI2025

Global Simulations Suggest Biomass Burning Aerosol Emissions From Grassland Fires Could Be Important Ice Nucleating Particles

Ice nucleating particles (INP) capable of nucleating ice crystals via immersion freezing at temperatures above approximately −35°C may strongly influence cloud glaciation, with implications for global precipitation and climate feedback. In addition to mineral dust, soil dust, and marine organics, laboratory and field measurements suggest biomass burning aerosols (BBA) can act as immersion-mode INP between around −30°C and −15°C. However, the contribution of BBA to the global INP budget remains poorly understood due to poor knowledge of which fuels yield INPs, uncertainties in global coverage of those fuels, and unknown size distributions of the INPs in the BBA. Nonetheless, with some understanding of these uncertainties from sensitivity studies, the relative importance of ice nucleation activity of BBA compared to other INP sources can be quantified. In this work, we investigate the potential global importance of BBA as INP using a global aerosol-climate model, specifically the UK Met Office Unified Model (UM). We evaluate the model using field campaign data sets. We examine potential uncertainties in fuel types and particle sizes on BBA-based INP concentrations. Averaged over June–September between 15°S and 50°S, BBA is a more important INP than dust and marine INP about 30% of the time at altitudes with temperatures between −30°C and −20°C. Our simulations therefore suggest BBA INPs may be at least as important as mineral dust and marine INP over the atmospheric regions and seasons where grassland fires are frequent.

Gohil, Kanishk [Carnegie Mellon University, Pittsb↗
DOE OSTI2024

Enhanced light absorption for solid-state brown carbon from wildfires due to organic and water coatings

Abstract Wildfires emit solid-state strongly absorptive brown carbon (solid S-BrC, commonly known as tar ball), critical to Earth’s radiation budget and climate, but their highly variable light absorption properties are typically not accounted for in climate models. Here, we show that from a Pacific Northwest wildfire, over 90% of particles are solid S-BrC with a mean refractive index of 1.49 + 0.056 i at 550 nm. Model sensitivity studies show refractive index variation can cause a ~200% difference in regional absorption aerosol optical depth. We show that ~50% of solid S-BrC particles from this sample uptake water above 97% relative humidity. We hypothesize these results from a hygroscopic organic coating, potentially facilitating solid S-BrC as nuclei for cloud droplets. This water uptake doubles absorption at 550 nm and the organic coating on solid S-BrC can lead to even higher absorption enhancements than water. Incorporating solid S-BrC and water interactions should improve Earth’s radiation budget predictions.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Characteristics of Vertical Ground Motions and Their Effect on the Seismic Response of Bridges in the Near-Field: A State-of-the-Art Review

Despite the evidence from past earthquakes and several numerical investigations demonstrating the detrimental impact of vertical ground motions (VGMs) on the integrity of bridge structures, incorporating their effects into seismic assessment and design procedures has traditionally been given limited consideration. Current codes utilize rather simplistic approaches to account for the concurrent effects of vertical and horizontal motions in structural performance evaluations, potentially leading to unconservative estimates of structural demands. This paper reviews the main features of VGMs and their effect on the seismic response of bridges. The methods and empirical models available to estimate vertical motions for design purposes are discussed, and research gaps and related research needs are identified. Finally, the emerging role of physics-based ground-motion simulations, as well as their limitations, in supporting future research and informing the development of simplified design procedures is examined. The main areas of interest for future research are identified in the need to carry out systematic sensitivity studies to gain insight into the main earthquake parameters that influence key VGMs features, understand the influence of soil nonlinearities on VGMs amplitude and frequency content, inform the development of empirical models that cover a range of site conditions and source-to-site distances where current models are poorly constrained, generate arrays of motions to update coherency models to properly inform the analysis of distributed infrastructure, investigate the impulsive character of VGMs, and assess the approximations made in estimating VGMs with 1D site response analyses. Furthermore, specific focus is laid on large-magnitude earthquakes in the near-field.

Asynchronous ground motion↗
DOE OSTI2021

Skyshine Calculations for a Large Spent Nuclear Fuel Storage Facility with SCALE 6.2.3

The SCALE code system developed at Oak Ridge National Laboratory includes state-of-the-art capabilities for radiation source term and radiation transport simulations that can be used in numerous applications, including dose rate analyses of complex consolidated interim storage facilities (CISFs). A licensed CISF could be used to store tens of thousands of tonnes of spent nuclear fuel discharged from commercial power reactors using various cask and storage pad designs. A CISF design must comply with the regulatory requirements provided in 10 CFR Part 72, including requirements related to annual dose limits applicable to real individuals located beyond the area controlled by the licensee. Therefore, calculating a dose to the public is a necessary part of the licensing process for the construction of a CISF. These calculations are very challenging because of the complexity of the CISF design and the low magnitude of dose rate at large distances from the facility. This paper describes detailed far-field dose rate calculations performed for a proposed CISF using MAVRIC, the Monte Carlo radiation shielding sequence in SCALE 6.2.3, with automated variance reduction based on discrete ordinates calculations. The method presented in this paper uses a detailed Monte Carlo radiation transport simulation in one step from source to dose rate. A series of independent simulations was made using the complete site geometry (all casks present), but with only one cask containing radiation sources to obtain the dose rate maps produced by each storage cask. The CISF dose rate map was obtained by adding the dose rate maps produced by the independent individual cask simulations. Ample volumes of air and soil extending beyond the location of interest for dose rate calculation were included in the calculation model to properly simulate important radiation attenuation and scattering events that affect far-field dose rates. Furthermore, a comprehensive sensitivity study is included in this paper to illustrate the importance of selecting appropriate air volume, mass density, and composition for CISF skyshine dose rate calculations. Dry soil and soil containing water were analyzed to determine their effects on groundshine radiation.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2025

Advancing the STS Neutron Moderator Design with an Automated Optimization Workflow and Unstructured Mesh Modeling

With the Second Target Station approaching its final design phase, a detailed neutronics evaluation of its critical components is necessary. Optimizing the dimensions of the two cold-source moderators that are at the heart of this facility presents a multi-objective optimization problem for which an accurate geometric description is crucial. We have applied a fully automated optimization workflow in which a detailed unstructured mesh geometry is automatically generated with Attila4MC, starting from a parametrized CREO geometry followed by preprocessing with SpaceClaim. With this geometry, a MCNP run is performed to calculate the brightness metrics, which are subsequently provided to the optimization algorithm in DAKOTA that provides new parameters and drives the optimization loop until convergence. In this paper, we show the results of the analysis that are used for the final design of the cylindrical and tube moderator. The optimization simulations provide a refinement to and confirmation of the conclusions of the previous design iteration. Additional to the optimization, a sensitivity study is performed to study the effect of minor geometry changes, which is important for the final engineering design. In conclusion, with these studies, we demonstrate that the automated workflow and high-fidelity unstructured mesh modeling are efficient tools for a thorough design evaluation.

DAKOTA↗
DOE OSTI2024

Accuracy of kinetic equilibrium reconstruction of NSTX and NSTX-U plasmas and its impact on the transport and stability analysis

An accurate magnetohydrodynamic (MHD) equilibrium reconstruction is an essential starting point for stability and transport plasma analysis. Herein this work describes an approach for obtaining kinetic equilibrium reconstructions using the OMFIT framework, which has been applied for the first time to spherical tokamak data from NSTX and NSTX-U. The EFIT equilibrium solver is integrated with experimental data analysis procedures and subsequent TRANSP transport simulations to enhance the accuracy of the reconstruction, in particular, at the edge region, by adding constraints on the total pressure and current density profiles, based on the transport code solution. The accuracy of the equilibrium reconstruction depends on the uncertainty and number of constraints, as well as the choice of basis functions to represent the pressure and current density profiles. Improved fidelity of the equilibrium reconstruction is demonstrated by reducing the variability of the magnetic axis and boundary locations from several centimeters, for reconstructions based on magnetic and experimental pressure constraints, to only several millimeters, for kinetic reconstructions based on transport code constraints, when different representations of basis functions were tested. The variability of the safety factor on axis was reduced ten times in the same sensitivity study. The accuracy of the equilibrium reconstruction and subsequent mapping of the experimental kinetic profile data have a significant impact on the trapped gyro Landau fluid and linear CGYRO turbulence simulations, which predict different spectra of unstable modes and turbulent fluxes for cases with different numbers of constraints in the equilibrium reconstruction. Conversely, the stability analysis performed using the GATO code shows plasmas that are stable to n = 1 MHD modes in both equilibria using magnetic and experimental pressure constraints as well as the transport code constrained equilibrium. However, a scan of parameters away from these conditions shows considerable deviation in the threshold of unstable modes between these reconstructions. Therefore, for reliable plasma analysis and use in turbulence and stability calculations, a high-fidelity equilibrium reconstruction with accurate kinetic constraints based on transport code solutions is necessary.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Fast ion transport by sawtooth instability in the presence of ICRF-NBI synergy in JET plasmas

JET experiments have shown that the three-ion scenarios using waves in the ion cyclotron range of frequencies (ICRF) is an efficient way to build fast ion population through beam ion acceleration by radio frequency (RF) waves. Here, such a heating scheme is applied to plasmas with at least two thermal ion species. Analysis of mixed discharges with complex heating schemes requires a workflow that allows to model thermal and fast ion transport consistently. This paper is dedicated to modelling of a mixed plasma discharge with significant fraction of fast ions and contributes to development of fast ion transport models. For interpretive analysis with the TRANSP code a JET hydrogen-deuterium (HD) plasma discharge with neutral beam injection (NBI) and ICRF heating has been chosen. The task is complicated by NBI-ICRF synergy and plasma magnetohydrodynamic activity, like sawtooth crashes. D beam ions accelerated by RF waves form a high energy tail in fast ion distribution. Significant difference between the neutron rate computed by TRANSP and measured one is observed if the same diffusivity for electrons and ions is assumed. Sensitivity studies show that uncertainties in input plasma parameters and thermal ion transport models are crucial for modelling mixed plasma discharges and increased D transport is required to reach the plasma composition consistent with diagnostic measurements at the plasma edge. Fast ion redistribution by a sawtooth instability is characterised by non-resonant transport due to reconnection of magnetic field lines and resonant transport caused by resonance interaction between the instability and fast ions. With ORBIT simulations it has been shown that resonant interaction strongly affects fast ions of high energies, like beam ions accelerated by RF waves and fusion products. For the considered case, fast ion profiles simulated by ORBIT remain peaked after the sawtooth crashes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Artificial intelligence-based predictive modeling for imaging neutral particle analyzers on the DIII-D tokamak

The Imaging Neutral Particle Analyzer (INPA) at DIII-D is a diagnostic system used to accurately resolve the energy and spatial distributions of fast ions in fusion plasmas. A novel artificial intelligence (AI) technique named INPA-net is based on Reservoir Computing Networks and developed here to predict active and passive signals produced by charge-exchange reactions from injected and edge-cold neutrals, respectively, in magnetically confined fusion plasmas. This model is trained using a set of 21 time domain signals between 0 s to 3.35 s that includes injected beam and thermal plasma information, and 6444 real 2D experimental images of the INPA in 12 plasma discharges at DIII-D. The trained neural network is able to forecast experimental images in real-time. The model achieves an R-squared value of 0.91, which is higher than the 0.83 value achieved by a simple linear regression model. This improvement highlights the model's enhanced predictive accuracy for measured images from the validation set. This AI approach is valuable due to its rapid response times and potential for integration into real-time plasma control systems. A version of this model capable of generating syntehic images would be useful for the real-time monitoring of fast-ion transport. A comprehensive sensitivity study reveals that INPA-net maintains high performance even with variations in the input parameters, indicating the model's robustness and reliability. While developed for the INPA, the underlying architecture is adaptable and may be applied to various 2D imaging diagnostics in fusion research.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Cross-section measurements to low-lying excited final states in the 24 M(α,p) 27 Al*(γ) reaction as an energy source for x-ray bursts

Energy production in Type Ia x-ray bursts depends on a number of nuclear reactions that become efficient in a hot environment (up to 2 GK). Model sensitivity studies have been performed in an effort to better identify the reactions that have the largest effect, but these efforts are hampered by the high level of complexity of the astrophysical model and large nuclear physics uncertainties. In a recent study, the 24 Mg(α,p) 27 Al reaction was found to significantly affect the energy generation in x-ray bursts. This manuscript reports the first study of the 24 Mg(α,p 1,2 γ) 27 Al reaction at energies relevant for x-ray bursts. In conclusion, the branches to the 27 Al excited states increase to a small degree the estimates of the total astrophysical 24 Mg(α,p) 27 Al reaction rate.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Proton capture on 34 S in the astrophysical energy regime of O-Ne novae

Nuclear reaction sensitivity studies have shown that the final isotopic abundance of O-Ne nova nucleosynthesis is dependent on the 34 S(p,γ) 35 Cl reaction at astrophysical energies corresponding to peak nova burning temperatures of 0.1–0.4 GK. Isotopic ratios of the S, Cl, and Ar products are all used in various methods of cosmochemical analysis of presolar meteoritic grains. Due to the lack of direct experimental data, the 34 S+p reaction rate has been estimated using statistical modeling or information from indirect nucleon transfer experiments. In order to provide direct reaction information, here the resonance strengths of several low energy resonances, E c.m. = 272–495 keV, in the 34 S(p,γ) 35 Cl reaction were measured for the first time in inverse kinematics using the DRAGON recoil separator located at TRIUMF, Canada’s Particle Accelerator Centre in Vancouver.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Measurement of the Mo 100 ( α , x n ) cross section at weak r -process energies

Here, the weak r process in neutrino-driven winds following a core-collapse supernova is thought to contribute to the cosmic abundances of the first r-process peak elements between Se and Ag. Sensitivity studies have found that the early nucleosynthesis in the weak r process is primarily driven by (α, xn) reactions due to the high temperatures, and that current nuclear physics uncertainties in the (α, xn) rates result in significant uncertainties of the calculated abundances. The weak r-process path proceeds several nuclei away from stability where (α, xn) reaction cross sections have not yet been measured. In this paper we report the 100 Mo (α, xn) cross section (between 8.9 and 13.2 MeV in the center of mass, corresponding to 3.5–6.8 GK) in inverse kinematics using the Multi-Sampling Ionization Chamber (MUSIC) detector at the Argonne Tandem Linac Accelerator System (ATLAS) facility. With this first measurement of the 100 Mo (α, xn) cross section, we have demonstrated the ability of MUSIC to measure (α, xn) cross sections for A up to 100, therefore paving the way for further measurements with radioactive beams at ATLAS or the Facility for Rare Isotope Beams.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Isomeric yield ratios of fission products: A missing piece in reactor antineutrino summation calculations

The calculation of the spectrum of antineutrinos ($\bar{v}_e$) from a reactor is a complicated problem requiring several nuclear data and physics inputs. Many of these have been investigated thoroughly to improve calculations and properly account for uncertainties. One input which has heretofore escaped consideration is the fission-yield distribution between ground and isomeric states. Here, in this work, we explore the effect of incorporating newly evaluated isomeric yield ratios (IYR) for 43 fission products into summation calculations and identify the disproportionate importance of certain isotopes, particularly at higher energies. Our analysis indicates that updated IYRs contribute to a significant increase in the $\bar{v}_e$ spectrum around and above 7 MeV, with increases of more than 50% at higher energies. Through a detailed sensitivity study, we highlight a number of isotopes for which the IYR has a substantial effect on the $\bar{v}_e$ spectrum. This work stresses the critical role of isomeric yields in calculations of reactor $\bar{v}_e$ spectra and points to the necessity for their accurate experimental determination, especially for key fission products, in order to refine our understanding and address the observed anomalies between measured and calculated $\bar{v}_e$ spectra.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

A Nonlinear Regression Method for Composite Protection Modeling of Induction Motor Loads

Protection equipment are used to prevent damages to induction motor loads by isolating those from the power network in the event of severe faults. Modeling the response of induction motor loads and their protection is vital for power system planning and operation, especially in understanding system's response moments after a fault has occurred. This article proposes an optimization based framework to generate composite protection models for commercial building motor loads. Introducing a mathematical abstraction, the task of finding a suitable (simplified) model of the composite protection scheme is formulated as a nonlinear regression problem. Numerical examples are provided to illustrate the application of the framework.

Kundu, Soumya↗
DOE OSTI2021

The Sensitivity Analysis of Coil Misalignment for a 200-kW Dynamic Wireless Power Transfer System with an LCC-S and LCC-P Compensation

The dynamic wireless power transfer systems may reduce the battery size of electric vehicles while maintaining the travel range. Similar to the stationary wireless power transfer systems, the compensation networks are desired in dynamic wireless charging systems at the primary and secondary sides to reduce the reactive power requirement from the source and improve the efficiency. Consequently, it is essential to understand the behavior of the compensation networks, especially the sensitivity to misalignments. This paper presents a sensitivity analysis of LCC-S and LCC-P compensation networks for a 200-kW power transfer with variations in the coupling coefficient due to the electric vehicle travel and misalignments between the pads. The sensitivity study combines the electromagnetic finite element analysis and circuit analysis to determine the impact of misalignments on the wireless power transfer system characteristics. The theoretical analysis is verified by conducting circuit simulations at discrete points to verify the accuracy of the mathematical model.

Kavimandan, Utkarsh↗
DOE OSTI2025

Quench Simulation and Validation for EIC Direct Wind Magnets

The Electron-Ion Collider (EIC), a new facility to be built in the United States at the U.S. Department of Energy’s Brookhaven National Laboratory in collaboration with Thomas Jefferson National Accelerator Facility-will be the only new particle collider in the U.S. for the next 10∼20 years. The EIC project will investigate the origin of mass and building blocks of the Universe. It will accelerate and collide beams of electrons with protons and nuclei. Unlike any collider that has ever existed before, both EIC’s particle beams will be polarized, meaning their spins will be aligned in a controlled way to allow for precision measurements. Different types of new superconducting magnets will be required near the interaction region (IR) of EIC. A 6-around-1 NbTi cable will be used for some of those new superconducting magnets. Given the relatively smaller conductors used for those direct wind magnets, the overall engineering current density will be 2–3 times higher in comparison to traditional collared superconducting magnets in other accelerator magnets. This is technically challenging and potentially represents high risks for quench protection. Here, in this paper we use different tools and methods for quench simulation for several EIC direct wind magnets. We also present results of a sensitivity study, exploring different material properties and Jc curves. This paper also discusses and validates the quench simulation results with measurement results from a small-scale serpentine magnet tested at BNL in October 2024.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Hierarchical Convolutional Neural Networks for Event Classification on PMU Measurements

Event classification is one of the central components of automated disturbance analysis based on PMU measurements. Obtaining high-quality event labels remains a challenge for supervised learning-based classification of local and system-wide events in power grids due to its labor-intensive requirement. We present a sensitivity study considering rapidly refined, partially and fully inspected event labels that leads to evidence that hierarchical convolutional neural networks (HCNNs) outperform traditional classification models regardless of the quality of the available event labels. Furthermore, it is demonstrated that performance similar to the one obtained using entirely domain-driven labeling can be achieved as long as the involved expert does not mislabel more than ~5% of the event data captured by PMU measurements.

47 OTHER INSTRUMENTATION↗
DOE OSTI2024

Optimization of a High Pressure Turbine Blade and Sector-Based Annular Rig Design for Supercritical CO2 Power Cycle Representative Testing

Abstract As part of the ongoing research into the design of hardware for zero emission cycles, a first-stage high-pressure turbine (HPT) blade is optimized for a 300 MWe supercritical CO2 (sCO2) power cycle using the surrogate-assisted genetic algorithm optimizer in Numeca FINE/Design three-dimensional with objectives of increasing efficiency and decreasing heat load to the blade. Supercritical CO2 property tables are constructed from NIST REFPROP data for the condensable gas simulation in FINE/Turbo. A detailed mesh sensitivity study is performed for a baseline design to identify the proper-grid refinement and efficiently allocate resources for the optimization. Seventy design variables are selected for the initial population generation. Self-organizing maps are then used to focus the design variables on the most important ones affecting the objective functions. The optimization results in approximately 3000 three-dimensional Reynolds Averaged Navier Stokes simulations of different blade shapes with increases in efficiency of up to 0.85% and decreases in heat load of 14%. Families of blade shapes are identified for experimental testing in an annular rig at the Purdue Experimental Turbine Aerothermal Laboratory. A design to adapt the annular cascade for testing optimized geometries is introduced, which features eccentric radius sectors allowing for scaled-up geometries of sCO2 optimized blade profiles to be tested at design cycle representative conditions at high Reynolds numbers in dry air. Analysis into the effects of Reynolds number, working fluid, and geometric relations are presented to prove the efficacy of the test method.

Engineering↗