Search NASA⌕ Search

SEARCH · Search NASA

Results for “transportation simulations”

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 199 records · Page 11

Analytic Modeling of a Deep Shielding Problem

Previous generations of scientists would make tremendous efforts to simplify non-tractable problems and generate simpler models that preserved the fundamental physics. This process involved applying assumptions and simplifications to reduce the complexity of the problem until it reached a solvable form. Each assumption and simplification was chosen and applied with the intent to preserve the essential physics of the problem, since, if the core physics of the problem were eliminated, the simplified model served no purpose. Moreover, if done correctly, solutions to the reduced model would serve as useful approximations to the original problem. In a sense, solving the simple models laid the ground-work for and provided insight into the more complex problem. Today, however, the affordability of high performance computing has essentially replaced the process for analyzing complex problems. Rather than "building up" a problem by understanding smaller, simpler models, a user generally relies on powerful computational tools to directly arrive at solutions to complex problems. As computational resources grow, users continue trying to simulate new, more complex, or more detailed problems, resulting in continual stress on both the code and computational resources. When these resources are limited, the user will have to make concessions by simplifying the problem while trying to preserve important details. In the context of the Monte Carlo N-Particle radiation transport simulation tool, simplifications typically come as reductions in geometry, or by using variance reduction techniques. Both approaches can influence the physics of the problem, leading to potentially inaccurate or non-physical results. Errors can also be introduced as a result of faulty input into a computational tool: something as simple as transposing numbers in a tally input can result in incorrect answers. In this paradigm, reduced complexity computational and analytical models still have an important purpose. The explicit form of an analytic solution is arguably the best way to understand the qualitative properties of simple models. In contrast to "building up" a complex problem through understanding simpler problems, results from detailed computational scenarios can be better explained by "building down" the complex model through simple models rooted in the fundamental or essential phenomenology. Simplified analytic and computational models can be used to 1) increase a user's confidence in the computational solution of a complex model, 2) confirm there are no user input errors, and 3) ensure essential assumptions of the simulation tool are preserved. This process of using analytic models to develop a more valuable analysis of simulation results is named the results assessment methodology. The utility of the results assessment methodology and a complimentary sensitivity analysis is exemplified through the analysis of the neutron flux in a dry used fuel storage cask. This application was chosen due to current scientific interest in used nuclear fuel storage.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Analytic Modeling of a Deep Shielding Problem

Previous generations of scientists would make tremendous efforts to simplify non-tractable problems and generate simpler models that preserved the fundamental physics. This process involved applying assumptions and simplifications to reduce the complexity of the problem until it reached a solvable form. Each assumption and simplification was chosen and applied with the intent to preserve the essential physics of the problem, since, if the core physics of the problem were eliminated, the simplified model served no purpose. Moreover, if done correctly, solutions to the reduced model would serve as useful approximations to the original problem. In a sense, solving the simple models laid the ground-work for and provided insight into the more complex problem. Today, however, the affordability of high performance computing has essentially replaced the process for analyzing complex problems. Rather than "building up" a problem by understanding smaller, simpler models, a user generally relies on powerful computational tools to directly arrive at solutions to complex problems. As computational resources grow, users continue trying to simulate new, more complex, or more detailed problems, resulting in continual stress on both the code and computational resources. When these resources are limited, the user will have to make concessions by simplifying the problem while trying to preserve important details. In the context of the Monte Carlo N-Particle radiation transport simulation tool, simplifications typically come as reductions in geometry, or by using variance reduction techniques. Both approaches can influence the physics of the problem, leading to potentially inaccurate or non-physical results. Errors can also be introduced as a result of faulty input into a computational tool: something as simple as transposing numbers in a tally input can result in incorrect answers. In this paradigm, reduced complexity computational and analytical models still have an important purpose. The explicit form of an analytic solution is arguably the best way to understand the qualitative properties of simple models. In contrast to "building up" a complex problem through understanding simpler problems, results from detailed computational scenarios can be better explained by "building down" the complex model through simple models rooted in the fundamental or essential phenomenology. Simplified analytic and computational models can be used to 1) increase a user's confidence in the computational solution of a complex model, 2) con firm there are no user input errors, and 3) ensure essential assumptions of the simulation tool are preserved. This process of using analytic models to develop a more valuable analysis of simulation results is named the results assessment methodology. The utility of the results assessment methodology and a complimentary sensitivity analysis is exempli fied through the analysis of the neutron flux in a dry used fuel storage cask. This application was chosen due to current scientific interest in used nuclear fuel storage.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Performance Analysis of the Expected Track Length Estimator

This work investigates whether performance of the MCNP code is improved by the implementation of the expected track length estimator. The expected track length estimator is known to decrease variance when compared to the track length estimator in many classes of Monte Carlo particle transport simulations. However, the computational cost of computing an exponent for every particle collision within a tallied cell has previously been considered prohibitively expensive. Therefore, the goal of this work is to determine whether the variance improvements seen by the expected track length estimator in most cases makes up for the increase in transport time when implemented in the MCNP code.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Application of Scalar Filtered Density Function to Turbulent Flows Under Supercritical Condition

Abstract The scalar filtered density function (FDF) methodology is extended and employed for large eddy simulation (LES) of turbulent flows under supercritical condition. To describe real fluid behavior, the extended methodology incorporates the generalized heat and mass diffusion models along with real fluid thermodynamic relations which are derived using the cubic Peng–Robinson equation of state. These models are implemented within the stochastic differential equations comprising the scalar FDF transport. Simulations are conducted of a temporally developing mixing layer under supercritical condition and the results are assessed by comparing with data generated by direct numerical simulation (DNS) of the same layer. The consistency of the proposed FDF methodology is assessed. The LES-FDF predictions are shown to agree favorably with the DNS data and exhibit several key features pertaining to supercritical turbulent flows.

Energy & Fuels↗

Precision pulse shape simulation for proton detection at the Nab experiment

The Nab experiment at Oak Ridge National Laboratory, USA, aims to measure the beta-antineutrino angular correlation following neutron β decay to an anticipated precision of approximately 0.1%. The proton momentum is reconstructed through proton time-of-flight measurements, and potential systematic biases in the timing reconstruction due to detector effects must be controlled at the nanosecond level. In conclusion, we present a thorough and detailed semiconductor and quasiparticle transport simulation effort to provide precise pulse shapes, and report on relevant systematic effects and potential measurement schemes.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

A line-binned treatment of opacities for the spectra and light curves from neutron star mergers

ABSTRACT The electromagnetic observations of GW170817 were able to dramatically increase our understanding of neutron star mergers beyond what we learned from gravitational waves alone. These observations provided insight on all aspects of the merger from the nature of the gamma-ray burst to the characteristics of the ejected material. The ejecta of neutron star mergers are expected to produce such electromagnetic transients, called kilonovae or macronovae. Characteristics of the ejecta include large velocity gradients, relative to supernovae, and the presence of heavy r-process elements, which pose significant challenges to the accurate calculation of radiative opacities and radiation transport. For example, these opacities include a dense forest of bound–bound features arising from near-neutral lanthanide and actinide elements. Here we investigate the use of fine-structure, line-binned opacities that preserve the integral of the opacity over frequency. Advantages of this area-preserving approach over the traditional expansion–opacity formalism include the ability to pre-calculate opacity tables that are independent of the type of hydrodynamic expansion and thus eliminate the computational expense of calculating opacities within radiation-transport simulations. Tabular opacities are generated for all 14 lanthanides as well as a representative actinide element, uranium. We demonstrate that spectral simulations produced with the line-binned opacities agree well with results produced with the more accurate continuous Monte Carlo Sobolev approach, as well as with the commonly used expansion–opacity formalism. The agreement between the line-binned and expansion–opacity results is explained as arising from the similarity in their opacities in the limit of low optical depth, where radiation transport is important in the ejecta. Additional investigations illustrate the convergence of opacity with respect to the number of included lines, and elucidate sensitivities to different atomic physics approximations, such as fully and semirelativistic approaches.

79 ASTRONOMY AND ASTROPHYSICS↗

Spatial Correlations of the Poisson Model for Radiation Transport

Characterizing the relationship between bulk physical properties and mixing in randomly heterogeneous media is a central challenge across many areas of science and engineering. A benchmark model for such studies is the Poisson model, a random tessellation of space by a Poisson process of hyperplanes. In radiation transport studies, the lack of exact expressions for the Poisson model’s spatial multipoint functions has led to approximate methods being used, introducing unquantified sources of error. Here, we recently introduced an exact solution for the Poisson model’s multipoint functions and closely related conditional probability functions (CPFs), providing a new opportunity to understand and reduce these sources of error. In this paper, we enable a more rigorous investigation of radiation transport in stochastic media by applying the recently introduced exact solution for the Poisson model’s CPFs. This paper consists of three main contributions. First, we introduce a unified framework for CPFs of the Poisson model, encompassing the recently introduced exact CPFs as well as the previously introduced atomic mix, nearest-neighbor, and combination CPFs. This framework also includes existing pruning techniques for the approximate CPFs, such as angular exclusion, as well as a novel form of angular exclusion suitable for the exact CPFs. Second, we use the exact CPFs to characterize the spatial regions where each approximate three-point CPF is most accurate, thereby explaining the observed hierarchy of accuracy among the approximate models. Finally, we evaluate material transmittance, reflectance, and flux in a three-dimensional test problem using conditional point sampling, demonstrating the relationship between CPF accuracy and transport simulation accuracy.

Poisson model↗

Development of a broadband hard x-ray radiography platform for pulsed-power experiments

In this article, we develop and demonstrate a broadband hard x-ray radiography platform at the Zebra Pulsed Power Laboratory that integrates point-projection radiography, bremsstrahlung measurements, and hard x-ray pinhole imaging, designed to diagnose current-driven, cylindrically compressed matter. Initial laser-pulsed-power coupled experiments revealed that intense background radiation generated during 1 MA Zebra current shots overwhelmed laser-produced hard x-rays, obscuring radiographic images. Using combined spectral and spatial diagnostics, we identify energetic electrons accelerated by return currents as the dominant source of background hard x-rays, with electron energies inferred to be 3–4 MeV based on Monte Carlo simulations, and demonstrate mitigation through modifications to the radiation shielding and return-current configuration. The diagnostic platform was validated using a wire-pinch hard x-ray source, allowing radiographs of static 1-mm-diameter aluminum wires to be obtained while simultaneously measuring x-ray source spectra and spatial emission distributions within a single shot. Measured wire transmission profiles were quantitatively reconstructed using radiation transport simulations that incorporate an experimentally inferred two-temperature exponential x-ray spectrum from bremsstrahlung signal analysis and spatially distributed emission sources identified by pinhole imaging. Agreement between measured and simulated transmission profiles demonstrates the validity of the radiographic and x-ray source characterization approach, establishing this diagnostic platform as a promising tool for diagnosing magnetically driven, high-density plasmas relevant to warm dense matter and inertial fusion energy research.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Emerging investigator series: ion diffusivities in nanoconfined interfacial water films contribute to mineral carbonation thresholds

The dynamics and reactivity of nanoconfined fluids play critical roles across a wide range of environmental and technological systems, though reaction mechanisms and kinetics are not well understood. Here, the carbonation kinetics of forsterite (Mg 2 SiO 4 ) exposed to 90 atm supercritical carbon dioxide at 35–65 °C and 85–100% relative humidity (RH) was monitored with in situ X-ray diffraction, and partner molecular dynamics simulations were used to describe the free energy landscape of Mg 2+ adsorption and diffusion on forsterite surfaces covered in water films 3–10 monolayers thick. The collective findings reveal how decreasing the water film thickness by ~1.4 monolayers, from ~0.92 to ~0.64 nm, inhibited reaction rates by up to 97%, promoted anhydrous Mg-carbonate (magnesite, MgCO 3 ) precipitation, and more than doubled the apparent activation energy of carbonation. The transport simulations suggest that four monolayers are required to enable sufficiently facile Mg 2+ diffusion, helping explain previously observed water film thickness-dependent reactivity thresholds.

54 ENVIRONMENTAL SCIENCES↗

Performance Assessment for the E-Area Low-Level Radioactive Waste Disposal Facility at the Savannah River Site: Appendix H

This section provides supporting material for the development of DU-specific final inventory limits for the GW pathways for every generic waste form and SWF parent radionuclide requiring an inventory limit. The final inventory limits are based on nominal PA transport simulations using PORFLOW as reported in Chapter 5. The nominal PA settings represent the compliance case where some modeling parameter settings are defined based on conservative (pessimistically leaning) arguments. In the overall computational approach employed in this PA, a multitiered radionuclide screening process is adopted as discussed in Chapter 2, Section 2.3.7. In the multitiered process, the initial list of 1,252 parent radionuclides is shortened substantially using conservative, but simple, transport models, along with a reasonably low cutoff criterion of 0.1% SOF value. Multidimensional PORFLOW flow and transport modeling is employed for every parent radionuclide that failed the GW screening. The generic waste form limits represent Tier-3 analyses, while Tier-4 analyses are employed for SWF limits, where warranted.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Unstructured mesh-based multi-physics reactor transient simulation with iMC code

In this study, we present a time-dependent Monte Carlo simulation result about a complex fuel-loaded light water reactor system with thermal-hydraulic feedback using the KAIST iMC code. The intra-pin temperature distribution is calculated based on Monte Carlo-tallied detailed power distribution and the finite element heat transfer method. The temperature effect of fuel, absorber, and coolant are adequately reflected into the dynamic response by the on-the-fly Doppler broadened cross-section correction. For the time-dependent Monte Carlo transport simulation, the predictor-corrector quasi-static Monte Carlo method is used for this study. We tested the coupled analysis framework with a 17-by-17 fuel assembly loaded with the centrally-shielded burnable absorber (CSBA) fuel element. The temperature feedback effect on the material cross-section and the moderator density incurred negative reactivity in response to the external positive reactivity insertion, showing the inherent stability of the reactor system. The pin-wise power and temperature change during the system transient are presented and discussed. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Windowed multipole representation of R -matrix cross sections

Nuclear cross sections are basic inputs to any nuclear computation. Campaigns of experiments are fitted with the parametric R-matrix model of quantum nuclear interactions, and the resulting cross sections are documented—both pointwise and as resonance parameters (with uncertainties)—in standard evaluated nuclear data libraries (ENDF, JEFF, BROND, JENDL, CENDL, TENDL): these constitute our common knowledge of fundamental low-energy nuclear cross sections. In the past decade, a collaborative effort has been deployed to establish a new nuclear cross-section library format—the Windowed Multipole Library—with the goal of considerably reducing the computational cost of cross-section calculations in nuclear transport simulations. This work lays the theoretical foundations underpinning these efforts. From general R-matrix scattering theory, we derive the windowed multipole representation of nuclear cross sections. Though physically and mathematically equivalent to R-matrix cross sections, the windowed multipole representation is particularly well suited for subsequent temperature treatment of angle-integrated cross sections, in particular Doppler broadening, which is the averaging of cross sections over the thermal motion of the target atoms. Doppler broadening is of critical importance in neutron transport applications, as it ensures the stability of many nuclear reactors (negative thermal reactivity). Yet, Doppler broadening of nuclear cross sections has been a considerable bottleneck for nuclear transport computations, often requiring memory-costly pretabulations. We show that the windowed multipole representation can perform accurate Doppler broadening analytically (up to the first reaction threshold), from which we derive cross-section temperature derivatives to any order—all computable on the fly (without precalculations stored in memory). Furthermore, we here establish a way of converting the R-matrix resonance parameters uncertainty (covariance matrices) into windowed multipole parameters uncertainty. We show that generating stochastic nuclear cross sections by sampling from the resulting windowed multipole covariance matrix can reproduce the cross-section uncertainty in the original nuclear data file. The windowed multipole representation is therefore a novel nuclear physics formalism able to generate Doppler broadened stochastic nuclear cross sections on the fly, unlocking breakthrough computational gains for nuclear computations. Through this foundational paper, we hope to make the windowed multipole representation accessible, reproducible, and usable for the nuclear physics community, as well as provide the theoretical basis for future research on expanding its capabilities.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Metrics for Interzonal Dispersion Assessment of Airborne SARS-CoV-2 Within Office Buildings

High ventilation rate is increasingly considered as a mitigation strategy in the rapid spread of COVID-19 in existing buildings in addition to other options like UVGI and added filtration. However, this and other similar measures may also impact energy use and the efficiency of the buildings in which they are installed. Without concrete metrics, the impact of all these potential changes to the system cannot be properly assessed relative to the other impacts the measures might have, including additional energy use, added operation and maintenance costs, and reduced system life. Multizone airflow and contaminant transport simulations involving parameters such as infiltration rate, system configuration, and level of occupancy can provide critical information on relative risk distribution within a building. In this study, prototype airflow models are used to assess the degree to which ventilation-related measures can mitigate the spread of a virus like SARS-CoV-2 in a similar way that prototype energy models are used to study energy efficiency measures. A new detailed medium-sized office building prototype model is developed using the CONTAM software to represent contaminant transport in buildings like those modeled by the U.S. Department of Energy’s detailed medium office prototype model. The goal of this study is to generate a set of suitable metrics that will give a whole-building picture of airborne SARS-CoV-2 virus distributions across the different zones of the building under different parameters and scenarios.

DeGraw, Jason↗

Machine learning-enhanced model-based scenario optimization for DIII-D

Abstract Scenario development in tokamaks is an open area of investigation that can be approached in a variety of different ways. Experimental trial and error has been the traditional method, but this required a massive amount of experimental time and resources. As high fidelity predictive models have become available, offline development and testing of proposed scenarios has become an option to reduce the required experimental resources. The use of predictive models also offers the possibility of using a numerical optimization process to find the controllable inputs that most closely achieve the desired plasma state. However, this type of optimization can require as many as hundreds or thousands of predictive simulation cases to converge to a solution; many of the commonly used high fidelity models have high computational burdens, so it is only reasonable to run a handful of predictive simulations. In order to make use of numerical optimization approaches, a compromise needs to be found between model fidelity and computational burden. This compromise can be achieved using neural networks surrogates of high fidelity models that retain nearly the same level of accuracy as the models they are trained to replicate while reducing the computation time by orders of magnitude. In this work, a model-based numerical optimization tool for scenario development is described. The predictive model used by the optimizer includes neural network surrogate models integrated into the fast Control-Oriented Transport simulation framework. This optimization scheme is able to converge to the optimal values of the controllable inputs that produce the target plasma scenario by running thousands of predictive simulations in under an hour without sacrificing too much prediction accuracy.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Space–time reduced order model for large-scale linear dynamical systems with application to Boltzmann transport problems

A classical reduced order model for dynamical problems involves spatial reduction of the problem size. However, temporal reduction accompanied by the spatial reduction can further reduce the problem size without losing much accuracy, which results in a considerably more speed-up than the spatial reduction only. Recently, a novel space–time reduced order model for dynamical problems has been developed [17], where the space–time reduced order model shows an order of a hundred speed-up with a relative error of 10 –4 for small academic problems. However, in order for the method to be applicable to a large-scale problem, an efficient space–time reduced basis construction algorithm needs to be developed. Here we present the incremental space–time reduced basis construction algorithm. The incremental algorithm is fully parallel and scalable. Additionally, the block structure in the space–time reduced basis is exploited, which enables the avoidance of constructing the reduced space–time basis. These novel techniques are applied to a large-scale particle transport simulation with million and billion degrees of freedom. The numerical example shows that the algorithm is scalable and practical. Also, it achieves a tremendous speed-up, maintaining a good accuracy. Finally, error bounds for space-only and space–time reduced order models are derived.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Breakdown of helical edge state topologically protected conductance in time-reversal-breaking excitonic insulators

Gapless helical edge modes are a hallmark of the quantum spin Hall effect. Protected by time-reversal symmetry, each edge contributes a quantized zero-temperature conductance quantum G 0 ≡ e 2 /h. However, the experimentally observed conductance in WT e2 decreases below G 0 per edge already at edge lengths around 100 nm, even in the absence of explicit time-reversal breaking due to an external field or magnetic impurities. Here we show how a time-reversal breaking excitonic condensate with a spin-spiral order that can form in WT e2 leads to the breakdown of conductance quantization. We perform Hartree-Fock calculations to compare time-reversal breaking and preserving excitonic insulators. Using these mean-field models we demonstrate via quantum transport simulations that weak nonmagnetic disorder reproduces the edge length scaling of resistance observed in the experiments. We complement this by analysis in the Luttinger liquid picture, shedding additional light on the mechanism behind the quantization breakdown.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Modeling and simulation study for the design of the Fuel Interrogation and Examination using Submersible Tomography Analysis Mk II instrument

Here, the design of a submersible, gamma-ray tomography system for imaging irradiated nuclear fuel is described. The system—named Fuel Interrogation and Examination using Submersible Tomography Analysis (FIESTA) Mk. II—is a variation on a previous Mk. I I design, which was developed to non-destructively image fuel capsules irradiated in the Advanced Test Reactor at the Idaho National Laboratory. The FIESTA system uses a combination of transmission computed tomography and emission computed tomography to image the restructuring and fission product migration at different points of burnup. Changes made to FIESTA Mk. I reflect the revised design requirements and a need to reduce background noise, largely originating from downscattered photons from fuel and transmission source. The computational design and radiation transport simulations are described.

46 - INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AN↗

Dimensionality effects on multicomponent ionic transport and surface complexation in porous media

Coulombic interactions between charged species in pore water and at surface/solution interfaces are of pivotal importance for multicomponent reactive transport in porous media. In this study, we investigate the impact of domain dimensionality on electrostatically coupled dispersion and surface-solution reactions during transport of acidic plumes and major ions in porous media. Column and quasi two-dimensional flow-through experiments were performed, with identical silica porous media and under the same advection-dominated conditions. Equal mass fluxes of different electrolyte solutions (i.e., HCl - pH ~ 2.8, NaBr - 100 mM, HCl - pH ~ 2.8 plus NaBr - 100 mM) were continuously injected in the 1-D and 2-D experiments and breakthrough curves of pH and major ions were measured at the outlet of the domains. The presence of pronounced ionic strength gradients in the transverse direction in the 2-D setup caused distinct retardation and transport behaviors of protons and major ions which were not observed in the one-dimensional column experiments. Furthermore, in the cases of salt electrolytes injection, considerably enhanced release of H + (>61%) from the quartz surface was observed in the multidimensional system compared to the one-dimensional setup. Reactive transport modeling was performed to reproduce the experimental outcomes and to analyse the coupling between transport processes, based on the Nernst-Planck formulation of diffusive/dispersive fluxes and on surface complexation reactions at the solid-solution interface. Electrostatic interactions between Na + , B – , and H + , and deprotonation of the quartz surface upon the formation of sodium outer-sphere complexes, are the primary controllers of the spatial and temporal features displayed by the pH and major ions measurements. The reactive transport simulations allowed us to interpret the experimental observations, to visualize the distribution and spatio-temporal evolution of dissolved and solid species, to identify a spatially heterogeneous zonation of Coulombic interactions with distinct behavior at the fringe and core of the injected plumes in the multidimensional setup, and to quantify the different components of the Nerst-Planck fluxes of the charged solutes. This study demonstrates that the domain dimensionality directly affects electrostatic interactions between charged aqueous species in the pore water and surface complexation reactions at the solid-solution interface. The non-trivial effects of dimensionality on multicomponent ionic transport result in a significantly different behavior in 1-D and 2-D systems.

58 GEOSCIENCES↗