Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Model”

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 253 records · Page 14
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Estimates of Southern Hemispheric Gravity Wave Momentum Fluxes across Observations, Reanalyses, and Kilometer-Scale Numerical Weather Prediction Model

Abstract Gravity waves (GWs) are among the key drivers of the meridional overturning circulation in the mesosphere and upper stratosphere. Their representation in climate models suffers from insufficient resolution and limited observational constraints on their parameterizations. This obscures assessments of middle atmospheric circulation changes in a changing climate. This study presents a comprehensive analysis of stratospheric GW activity above and downstream of the Andes from 1 to 15 August 2019, with special focus on GW representation ranging from an unprecedented kilometer-scale global forecast model (1.4 km ECMWF IFS), ground-based Rayleigh lidar (CORAL) observations, modern reanalysis (ERA5), to a coarse-resolution climate model (EMAC). Resolved vertical flux of zonal GW momentum (GWMF) is found to be stronger by a factor of at least 2–2.5 in IFS compared to ERA5. Compared to resolved GWMF in IFS, parameterizations in ERA5 and EMAC continue to inaccurately generate excessive GWMF poleward of 60°S, yielding prominent differences between resolved and parameterized GWMFs. A like-to-like validation of GW profiles in IFS and ERA5 reveals similar wave structures. Still, even at ∼1 km resolution, the resolved waves in IFS are weaker than those observed by lidar. Further, GWMF estimates across datasets reveal that temperature-based proxies, based on midfrequency approximations for linear GWs, overestimate GWMF due to simplifications and uncertainties in GW wavelength estimation from data. Overall, the analysis provides GWMF benchmarks for parameterization validation and calls for three-dimensional GW parameterizations, better upper-boundary treatment, and vertical resolution increases commensurate with increases in horizontal resolution in models, for a more realistic GW analysis. Significance Statement Gravity wave–induced momentum forcing forms a key component of the middle atmospheric circulation. However, complete knowledge of gravity waves, their atmospheric effects, and their long-term trends are obscured due to limited global observations, and the inability of current climate models to fully resolve them. This study combines a kilometer-scale forecast model, modern reanalysis, and a coarse-resolution climate model to first compare the resolved and parameterized momentum fluxes by gravity waves generated over the Andes, and then evaluate the fluxes using a state-of-the-art ground-based Rayleigh lidar. Our analysis reveals shortcomings in current model parameterizations of gravity waves in the middle atmosphere and highlights the sensitivity of the estimated flux to the formulation used.

Meteorology & Atmospheric Sciences↗
DOE OSTI2023

A computational fluid dynamics model to estimate local quantities in firebrand char oxidation

Firebrand burning is a complex phenomenon that is influenced by several parameters which are difficult to fully explore experimentally. Computational fluid dynamics models capable of predicting local quantities are essential for accurate prediction of char oxidation in firebrands. This article presents a computational fluid dynamics model to estimate firebrand mass loss, diameter change, and surface temperature during char oxidation. The model was validated using previously conducted wind tunnel experiments. These experiments were conducted for firebrands of two different aspect ratios, which were arranged in three different configurations (single, horizontal array, and vertical array), and for four different wind speeds (0.5, 1, 1.5, and 2 m/s). The computational fluid dynamics results were compared with a previous 1 D model. In all the test cases, the computational fluid dynamics model predicted the physical phenomena with significantly improved accuracy compared to a 1 D model. The char oxidation model presented in this article can be coupled with other models to study firebrand generation and trajectory, biomass pyrolysis, fluidized bed reactors, and coal combustion.

Engineering↗
DOE OSTI2022

Important powertrain dynamics for developing models for control of connected and automated electrified vehicles

Connected and Automated Vehicles (CAV) technology presents significant opportunities for energy saving in the transportation sector. CAV technology forecasts vehicle and powertrain power needs under various terrain, ambient, and traffic conditions. Integration of the CAV technology in Hybrid Electric Vehicles (HEVs) provides the opportunity for optimal vehicle operation. Indeed, Hybrid Electric Vehicle powertrains present high degrees of flexibility and possibility for choosing optimum powertrain modes based on the predicted traction power needs. In modeling complex CAV powertrain dynamics, the modeler needs to consider short-time scale powertrain dynamics, such as engine transients, and hysteresis of mode-switching for a multi-mode HEV. Therefore, the powertrain dynamics essential for developing powertrain controllers for a class of connected HEVs is presented. To this end, control-oriented powertrain dynamic models for a test vehicle consisting of full electric, hybrid, and conventional engine operating modes are developed. The resulting powertrain model can forecast vehicle traction torque and energy consumption for the specified prediction horizon of the test vehicle. The model considers different operating modes and associated energy penalty terms for mode switching. Thus, the vehicle controller can determine the optimum powertrain mode, torque, and speed for forecasted vehicle operation via utilizing connectivity data. The powertrain model is validated against the experimental data and shows prediction error of less than 5% for predicting vehicle energy consumption. The model is used to create energy penalty maps that can be used for CAV control, for example fuel penalty map for engine torque changes (10–40 Nm) at each engine speed. The results of model-based optimization show optimum switching delays ranging from 0.4 to 1.4 s to avoid hysteresis in mode switching.

Engineering↗
DOE OSTI2022

A new dynamic zOnal model with air-diffuser (DOMA) - Application to thermal comfort prediction

A new Dynamic zOnal Model with Air-diffuser (DOMA) was developed. Several case studies were investigated and tested to evaluate and validate this program using measurement data. This new model was integrated into a TRaNsient SYstems Simulation program library and coupled with the multi-zone thermal model. The DOMA/TRNSYS coupled model was then used to predict room temperature distribution over an entire day of a single-zone building. The results show that increasing the heating outputs of the electric floor system, for example, from 75 to 200 W/m 2 , would not effectively improve the indoor thermal comfort, since the thermostat will reach the set point first and then turn off the system before the room gets enough heat and reach a comfortable level. This indicates the importance of selecting an appropriate location and set point for the thermostat when using a floor heating system. This potential thermal comfort issue can only be identified through the two-node model with a dynamic zonal model rather than the conventional PMV model, which thus suggests that for optimizing indoor thermal comfort of a building equipped with a time-sensitive control strategy and/or HVAC system, the TSENS results obtained from the two-node model integrated with DOMA are more appropriate than PMVs.

Construction & Building Technology↗
DOE OSTI2026

Data for Impact of Vertical and Seasonal Variation in Leaf Traits on Simulating Soybean Canopy Photosynthesis via 1D and 3D Modeling

Accurate modeling of photosynthesis is crucial for predicting crop productivity and quantifying the carbon cycle in agroecosystems. Leaf traits are essential inputs for modeling canopy photosynthesis. Yet, many existing models still use fixed plant functional type (PTF)-based values to parameterize leaf traits under a big-leaf or two-big-leaf assumption, neglecting their vertical profiles and seasonal changes. This simplification may introduce significant uncertainties in estimating gross primary productivity (GPP). In this study, we simulated soybean GPP and tested the effects of vertical and seasonal variation in three key leaf photosynthetic traits: the maximum carboxylation rate at 25 °C (Vcmax25), leaf chlorophyll content (LCC), and leaf mass per area (LMA) in the 1D-SCOPE and 3D-Helios models. Weekly field measurements were conducted during the growing season of 2024 to support the simulation. We designed ten leaf trait parameterization schemes by incorporating different combinations of vertical profiles and seasonal changes, while assuming homogeneous canopy architecture in both models. Our results revealed that Vcmax25 vertical and seasonal variation had the strongest influence on simulated GPP in both 1D and 3D models, while LCC and LMA effects were minimal. Particularly, the scheme with an empirically parameterized Vcmax25 profile achieved comparable performance to the scheme with the measured Vcmax25 profile. Both 1D-SCOPE and 3D-Helios accurately modeled GPP (SCOPE: R2 = 0.87, Bias = 0.55 µmol m⁻² s⁻¹; Helios: R2 = 0.9, Bias = 0.22 µmol m⁻² s⁻¹) under the most complex scheme, and their responses to vertical and seasonal variation in leaf traits were consistent, demonstrating the robustness of our findings. Based on our findings, we propose a scalable framework for parameterizing leaf traits to improve GPP simulations. This study contributes to improving the representation of leaf trait dynamics in canopy-level photosynthesis models, potentially enhancing our ability to predict crop productivity and understand agroecosystem carbon dynamics.

Photosynthesis↗
DOE OSTI2015

Determining Reduced Order Models for Optimal Stochastic Reduced Order Models

The use of parameterized reduced order models(PROMs) within the stochastic reduced order model (SROM) framework is a logical progression for both methods. In this report, five different parameterized reduced order models are selected and critiqued against the other models along with truth model for the example of the Brake-Reuss beam. The models are: a Taylor series using finite difference, a proper orthogonal decomposition of the the output, a Craig-Bampton representation of the model, a method that uses Hyper-Dual numbers to determine the sensitivities, and a Meta-Model method that uses the Hyper-Dual results and constructs a polynomial curve to better represent the output data. The methods are compared against a parameter sweep and a distribution propagation where the first four statistical moments are used as a comparison. Each method produces very accurate results with the Craig-Bampton reduction having the least accurate results. The models are also compared based on time requirements for the evaluation of each model where the Meta- Model requires the least amount of time for computation by a significant amount. Each of the five models provided accurate results in a reasonable time frame. The determination of which model to use is dependent on the availability of the high-fidelity model and how many evaluations can be performed. Analysis of the output distribution is examined by using a large Monte-Carlo simulation along with a reduced simulation using Latin Hypercube and the stochastic reduced order model sampling technique. Both techniques produced accurate results. The stochastic reduced order modeling technique produced less error when compared to an exhaustive sampling for the majority of methods.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2014

Breaking Computational Barriers: Real-time Analysis and Optimization with Large-scale Nonlinear Models via Model Reduction

Model reduction for dynamical systems is a promising approach for reducing the computational cost of large-scale physics-based simulations to enable high-fidelity models to be used in many- query (e.g., Bayesian inference) and near-real-time (e.g., fast-turnaround simulation) contexts. While model reduction works well for specialized problems such as linear time-invariant systems, it is much more difficult to obtain accurate, stable, and efficient reduced-order models (ROMs) for systems with general nonlinearities. This report describes several advances that enable nonlinear reduced-order models (ROMs) to be deployed in a variety of time-critical settings. First, we present an error bound for the Gauss-Newton with Approximated Tensors (GNAT) nonlinear model reduction technique. This bound allows the state-space error for the GNAT method to be quantified when applied with the backward Euler time-integration scheme. Second, we present a methodology for preserving classical Lagrangian structure in nonlinear model reduction. This technique guarantees that important properties--such as energy conservation and symplectic time-evolution maps--are preserved when performing model reduction for models described by a Lagrangian formalism (e.g., molecular dynamics, structural dynamics). Third, we present a novel technique for decreasing the temporal complexity --defined as the number of Newton-like iterations performed over the course of the simulation--by exploiting time-domain data. Fourth, we describe a novel method for refining projection-based reduced-order models a posteriori using a goal-oriented framework similar to mesh-adaptive h -refinement in finite elements. The technique allows the ROM to generate arbitrarily accurate solutions, thereby providing the ROM with a 'failsafe' mechanism in the event of insufficient training data. Finally, we present the reduced-order model error surrogate (ROMES) method for statistically quantifying reduced- order-model errors. This enables ROMs to be rigorously incorporated in uncertainty-quantification settings, as the error model can be treated as a source of epistemic uncertainty. This work was completed as part of a Truman Fellowship appointment. We note that much additional work was performed as part of the Fellowship. One salient project is the development of the Trilinos-based model-reduction software module Razor , which is currently bundled with the Albany PDE code and currently allows nonlinear reduced-order models to be constructed for any application supported in Albany. Other important projects include the following: 1. ROMES-equipped ROMs for Bayesian inference: K. Carlberg, M. Drohmann, F. Lu (Lawrence Berkeley National Laboratory), M. Morzfeld (Lawrence Berkeley National Laboratory). 2. ROM-enabled Krylov-subspace recycling: K. Carlberg, V. Forstall (University of Maryland), P. Tsuji, R. Tuminaro. 3. A pseudo balanced POD method using only dual snapshots: K. Carlberg, M. Sarovar. 4. An analysis of discrete v. continuous optimality in nonlinear model reduction: K. Carlberg, M. Barone, H. Antil (George Mason University). Journal articles for these projects are in progress at the time of this writing.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Modeling District Heating and Cooling Systems with URBANopt, GeoJSON to Modelica Translator, and the Modelica Buildings Library

The URBANopt project has successfully leveraged OpenStudio/EnergyPlus to model buildings and electrical systems at an urban scale; however, URBANopt has lacked the ability to model district thermal energy systems until recently. This paper will present the modeling infrastructure that was developed specifically for the analysis of district heating and cooling systems, and how it is integrated into the existing URBANopt framework. The paper also discusses the development of new models added to the Modelica Buildings Library to model various district energy system components including loads, energy transfer stations (ETS), distribution networks, and central plants. The paper describes how different building loads can be modeled including time series, TEASER reduced-order models, or Spawn of EnergyPlus models. URBANopt District Energy Systems allows the user to switch between the various configurations.

district heating and cooling↗
DOE OSTI2025

Comparing computational times for simulations when using PBPK model template and stand-alone implementations of PBPK models

Introduction We previously developed a PBPK model template that consists of a single model “superstructure” with equations and logic found in many physiologically based pharmacokinetic (PBPK) models. Using the template, one can implement PBPK models with different combinations of structures and features. Methods To identify factors that influence computational time required for PBPK model simulations, we conducted timing experiments using various implementations of PBPK models for dichloromethane and chloroform, including template and stand-alone implementations, and simulating four different exposure scenarios. For each experiment, we measured the required computational time and evaluated the impacts of including various model features (e.g., number of output variables calculated) and incorporating various design choices (e.g., different methods for estimating blood concentrations). Results We observed that model implementations that treat body weight and dependent quantities as constant (fixed) parameters can result in a 30% time savings compared with options that treat body weight and dependent quantities as time-varying. We also observed that decreasing the number of state variables by 36% in our PBPK model template led to a decrease of 20–35% in computational time. Other factors, such as the number of output variables, the method for implementing conditional statements, and the method for estimating blood concentrations, did not have large impacts on simulation time. In general, simulations with PBPK model template implementations of models required more time than simulations with stand-alone implementations, but the flexibility and (human) time savings in preparing and reviewing a model implemented using the PBPK model template may justify the increases in computational time requirements. Conclusion Our findings concerning how PBPK model design and implementation decisions impact computational speed can benefit anyone seeking to develop, improve, or apply a PBPK model, with or without the PBPK model template.

Bernstein, Amanda S.↗
DOE OSTI2016

Evaluation of observed and modelled aerosol lifetimes using radioactive tracers of opportunity and an ensemble of 19 global models

Aerosols have important impacts on air quality and climate, but the processes affecting their removal from the atmosphere are not fully understood and are poorly constrained by observations. This makes modelled aerosol lifetimes uncertain. In this study, we make use of an observational constraint on aerosol lifetimes provided by radionuclide measurements and investigate the causes of differences within a set of global models. During the Fukushima Dai-Ichi nuclear power plant accident of March 2011, the radioactive isotopes cesium-137 ( 137 Cs) and xenon-133 ( 133 Xe) were released in large quantities. Cesium attached to particles in the ambient air, approximately according to their available aerosol surface area. 137 Cs size distribution measurements taken close to the power plant suggested that accumulation-mode (AM) sulfate aerosols were the main carriers of cesium. Hence, 137 Cs can be used as a proxy tracer for the AM sulfate aerosol's fate in the atmosphere. In contrast, the noble gas 133 Xe behaves almost like a passive transport tracer. Global surface measurements of the two radioactive isotopes taken over several months after the release allow the derivation of a lifetime of the carrier aerosol. We compare this to the lifetimes simulated by 19 different atmospheric transport models initialized with identical emissions of 137 Cs that were assigned to an aerosol tracer with each model's default properties of AM sulfate, and 133 Xe emissions that were assigned to a passive tracer. We investigate to what extent the modelled sulfate tracer can reproduce the measurements, especially with respect to the observed loss of aerosol mass with time. Modelled 137 Cs and 133 Xe concentrations sampled at the same location and times as station measurements allow a direct comparison between measured and modelled aerosol lifetime. The e-folding lifetime τ e , calculated from station measurement data taken between 2 and 9 weeks after the start of the emissions, is 14.3 days (95 % confidence interval 13.1–15.7 days). The equivalent modelled τ e lifetimes have a large spread, varying between 4.8 and 26.7 days with a model median of 9.4 ± 2.3 days, indicating too fast a removal in most models. Because sufficient measurement data were only available from about 2 weeks after the release, the estimated lifetimes apply to aerosols that have undergone long-range transport, i.e. not for freshly emitted aerosol. However, modelled instantaneous lifetimes show that the initial removal in the first 2 weeks was quicker (lifetimes between 1 and 5 days) due to the emissions occurring at low altitudes and co-location of the fresh plume with strong precipitation. Deviations between measured and modelled aerosol lifetimes are largest for the northernmost stations and at later time periods, suggesting that models do not transport enough of the aerosol towards the Arctic. The models underestimate passive tracer ( 133 Xe) concentrations in the Arctic as well but to a smaller extent than for the aerosol ( 137 Cs) tracer. As a result, this indicates that in addition to too fast an aerosol removal in the models, errors in simulated atmospheric transport towards the Arctic in most models also contribute to the underestimation of the Arctic aerosol concentrations.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Advancing Building Energy Modeling with Large Language Models: Exploration and Case Studies

The rapid progression in artificial intelligence has facilitated the emergence of large language models like ChatGPT, offering potential applications extending into specialized engineering modeling, especially physics-based building energy modeling. This paper investigates the innovative integration of large language models with building energy modeling software, focusing specifically on the fusion of ChatGPT with EnergyPlus. A literature review is first conducted to reveal a growing trend of incorporating large language models in engineering modeling, albeit limited research on their application in building energy modeling. We underscore the potential of large language models in addressing building energy modeling challenges and outline potential applications including simulation input generation, simulation output analysis and visualization, conducting error analysis, co-simulation, simulation knowledge extraction and training, and simulation optimization. Three case studies reveal the transformative potential of large language models in automating and optimizing building energy modeling tasks, underscoring the pivotal role of artificial intelligence in advancing sustainable building practices and energy efficiency. The case studies demonstrate that selecting the right large language model techniques is essential to enhance performance and reduce engineering efforts. The findings advocate a multidisciplinary approach in future artificial intelligence research, with implications extending beyond building energy modeling to other specialized engineering modeling.

building energy modeling↗
DOE OSTI2021

High-resolution hourly surrogate modeling framework for physics-based large-scale building stock modeling

Surrogate modeling can play a key role in reducing high computational burdens for large-scale physics-based modeling and uncertainty quantification. With the rapid development of large-scale building stock energy modeling, surrogate modeling has also begun to be widely applied in this field; however, most existing surrogate models lack hourly time resolution for regional-scale modeling, which is essential for understanding building demand profiles and grid impacts. Further, there is generally a lack of necessary data and feature engineering frameworks specific to building modeling for efficiently managing large datasets and complex computations. Furthermore, this paper proposes a modeling framework for large-scale (city-/region-scale), high-resolution, high-fidelity surrogate building stock energy models. Our developed framework consists of six modules: (1) building stock energy modeling (ComStock TM and ResStock TM ), (2) data engineering for large simulation data, (3) high performance computing workflow, (4) feature engineering, (5) machine learning model development, and (6) model performance evaluation. Two case studies apply the developed framework in both residential and commercial building stock analysis to demonstrate its computational efficiency and surrogate modeling accuracies. Results show that surrogate models, when efficiently trained using the HPC workflow module, reach a high level of modeling accuracy for two case studies.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2017

Synthesis and Electrocatalytic Activity of [FeFe]‐Hydrogenase Model Complexes with Non‐Innocent Chelating Nitrogen‐Donor Ligands

To probe the influence of redox non‐innocent ligands on a well‐known class of [FeFe]‐hydrogenase models, three new asymmetrically disubstituted diiron complexes of the general formula (µ‐SRS)[Fe(CO) 3 ][Fe(CO)(N‐N)] {SRS = propane‐1,3‐dithiolate (pdt) or benzene‐1,2‐dithiolate (bdt), and N‐N = 2,2′‐bipyridine (bipy) or 2,2′‐bipyrimidine (bpym)} have been synthesized from their parent hexacarbonyls and characterized. The new complexes, (µ‐pdt)Fe 2 (CO) 4 (κ 2 ‐bpym) ( 2 ), (µ‐bdt)Fe 2 (CO) 4 (κ 2 ‐bipy) ( 3 ), and (µ‐bdt)Fe 2 (CO) 4 (κ 2 ‐bpym) ( 4 ), were fully characterized by spectroscopic and electrochemical techniques, and the results are compared to those of a related complex, (µ‐pdt)Fe 2 (CO) 4 (κ 2 ‐bipy) ( 1 ). The crystal structures of 2 – 4 show that, in each complex, the two iron units are in an eclipsed orientation, and the N‐N ligand lies in the basal plane. IR spectra and electrochemical analyses indicate that electron density at the iron centers decreases in the order 1 > 2 > 3 > 4 . Furthermore, compound 2 undergoes a ligand‐centered reduction at the same potential that the hexacarbonyl precursor undergoes its first reduction. However, unlike the 2,2′‐bipy derivatives 1 and 3 , the 2,2′‐bpym complexes 2 and 4 are not effective catalysts for electrochemical proton reduction from acetic acid.

Roy, Souvik↗
DOE OSTI2014

Electrochemical solvent reorganization energies in the framework of the polarizable continuum model

Electron transfer reactions at electrochemical interfaces play a critical role in a wide range of catalytic processes. A key parameter in the rate constant expressions for such processes is the reorganization energy, which reflects the energetic cost of the solute and solvent rearrangements upon electron transfer. In this paper, we present dielectric continuum methods for calculating the solvent reorganization energy for electrochemical processes. We extend the simple approach in which the solute is represented as a point charge located a specified distance from the electrode surface to the representation of the solute as a collection of point charges corresponding to the partial atomic charges of the molecule. We also develop a method for calculating the electrochemical solvent reorganization energies with molecular-shaped cavities within the framework of the polarizable continuum model (PCM). The electronic and inertial responses of the solvent are separated according to their respective time scales, and two limiting cases of the relation between the solute and solvent electrons are examined. The effects of the electrode are included with the integral equations formalism PCM (IEF-PCM), in which the molecule-solvent boundary is treated explicitly, but the effects of the electrode-solvent boundary are included through an external Green’s function. This approach accounts for the effects of detailed molecular charge redistribution in a molecular-shaped cavity, as well as the electronic and inertial solvent responses and the effects of the electrode. The calculated total reorganization energies are in reasonable agreement with experimental measurements for a series of electrochemical systems. Inclusion of the effects of the electrode is found to be essential for obtaining even qualitatively accurate solvent reorganization energies. These approaches are applicable to a wide range of systems and can be extended to include other types of boundaries, such as a self-assembled monolayer or double layer separating the electrode and the molecule. This research was supported as part of the Center for Molecular Electrocatalysis, an Energy Frontier Research Center funded by the U.S. Department of Energy, Office of Science, Office of Basic Energy Sciences.

Ghosh, Soumya↗
DOE OSTI2014

Evolution of helical perturbations in a thin-shell model of an imploding liner

A thin-shell model of the liner stability has been revisited and applied to the stability of the helical perturbations. Several stages of the implosion have been identified, starting from a long initial “latent” phase of an almost resting liner, continuing to the second stage of a rapid contraction and significant perturbation growth, and then transitioning to the third stage where perturbations become ballistic and highly non-linear. The stage of stagnation and rebound is beyond the scope of this paper. An importance of vorticity conservation during the late stages is emphasized. Nonlinear evolution of perturbations is followed up to the point of the formation of cusp structures. Effects of in-surface flows and of their enhancement due to the vorticity conservation are discussed. It is shown that the pre-machined perturbations created only on the outer surface of the liner grow much slower than one could anticipate. The limitations on the thin-shell description are discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2018

Modeling the one-dimensional imager of neutrons (ODIN) for neutron response functions at the Sandia Z facility

The one-dimensional imager of neutrons (ODIN) at the Sandia Z facility consists of a 10-cm block of tungsten with rolled edges, creating a slit imager with slit widths of either 250, 500, or 750 μm. Designed with a 1-m neutron imaging line of sight, we achieve about 4:1 magnification and 500-μm axial spatial resolution. The baseline inertial confinement fusion concept at Sandia is magnetized liner inertial fusion, which nominally creates a 1-cm line source of neutrons. ODIN was designed to determine the size, shape, and location of the neutron producing region, furthering the understanding of compression quality along the cylindrical axis of magnetized liner implosions. Challenges include discriminating neutrons from hard x-rays and gammas with adequate signal-to-noise in the 2 × 10 12 deuterium-deuterium (DD) neutron yield range, as well as understanding the point spread function of the imager to various imaging detectors (namely, CR-39). Modeling efforts were conducted with MCNP6.1 to determine neutron response functions for varying configurations in a clean DD neutron environment (without x-rays or gammas). Configuration alterations that will be shown include rolled-edge slit orientation and slit width, affecting the resolution and response function. The experiment to determine CR-39 neutron sensitivity, with and without a high density polyethylene (n, p) converter, an edge spread function, and resolution will be discussed.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗