Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Deterministic transport”

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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
NASA NTRS1988

Benchmark solutions for the galactic ion transport equations with spatial and energy coupling

In order to anticipate future space shielding requirements, NASA has initiated an effort to formulate computational methods to simulate radiation effects in space. As part of the program, numerical transport algorithms have been developed for the deterministic Boltzman equation describing galactic cosmic ray (GCR) interactions with matter. It thus becomes necessary to assess the accuracy of proposed deterministic algorithms. For this reason, analytical benchmark solutions to mathematically tractable galactic cosmic ray equations have recently been obtained. Even though these problems involve simplifying assumptions of the associated physics, they still contain the essential features of the basic transport processes. The solutions obtained are features of the basic transport processes. The solutions obtained are compared to results from numerical algorithms in order to ensure proper coding and to provide a measure of the accuracy of the numerical methods used in the algorithm. For the first time, mathematical methods have been applied to the galactic ion transport (GIT) equations in the straight ahead approximation with constant nuclear properties. The approach utilizes a Laplace transforms inversion yielding a closed form benchmark solution which is also computationally efficient.

Ganapol, Barry D.↗
DOE OSTI2026

Quantum statistical plasmonic metacrystals

Engineering materials that control quantum many-body dynamics remains challenging, as multiparticle interactions typically produce complex emergent behaviour that is difficult to predict. Here we introduce quantum statistical plasmonic metacrystals, structures in which the multiparticle dynamics mediated by optical near fields produce forbidden quantum statistical bands that enable selective transmission of different types of light. This functionality arises from a plasmonic structure composed of nanoantennas acting as meta-atoms. Multiphoton fields with statistics within the allowed bands propagate without distortion, whereas fields in forbidden bands are suppressed or driven towards the nearest accessible statistical state. We show that these bands are determined by the geometry and collective arrangement of the meta-atoms, providing a deterministic route to engineering quantum statistical transport. This platform establishes a room-temperature quantum material intrinsically sensitive to the quantum coherence of many-body photonic systems, enabling their robust manipulation and transport. Our results have implications for coherence-sensitive photonic materials for energy harvesting and scalable many-body quantum technologies.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Historical review and proof-of-concept future method demonstration of adaptive mesh refinement in nuclear engineering for increased fidelity and computational efficiency

As the nuclear industry's use of computational tool increases, the need for increased fidelity and computational efficiency is well known. While most approaches to increased fidelity rely on applying a fine mesh over the problem domain, a more efficient method is to apply an adaptive mesh refinement (AMR) algorithm to the mesh definition. In the field of nuclear engineering, AMR has previously been used in conjunction with deterministic methods, including: S{sub N} transport methods, Lattice Boltzmann Methods, and COMSOL. The future of AMR in nuclear engineering is to couple it to a Monte Carlo code with the goal of reducing calculation time. A proof-of-concept example yielded positive results for using the gradient of the flux as a refinement criteria. The refinement criteria was varied from 0.01 to 0.10, which yielded a recommended range of 0.01 to 0.04, and the number of refinement iterations was varied from 0 to 7, with diminishing returns seen after 5 iterations. After the success of the proof-of-concept exercise, work began on creating a full program coupling MCNP6.2 and the AMR algorithm in the deal.II library. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2020

SCALE Code System

The SCALE Code System is a widely used modeling and simulation suite for nuclear safety analysis and design that is developed, maintained, tested, and managed by the Reactor and Nuclear Systems Division (RNSD) of Oak Ridge National Laboratory (ORNL). SCALE provides a comprehensive, verified and validated, user-friendly tool set for criticality safety, reactor and lattice physics, radiation shielding, spent fuel and radioactive source term characterization, and sensitivity and uncertainty analysis. Since 1980, regulators, licensees, and research institutions around the world have used SCALE for safety analysis and design. SCALE provides an integrated framework with dozens of computational modules, including three deterministic and three Monte Carlo radiation transport solvers that are selected based on the desired solution strategy. SCALE includes current nuclear data libraries and problem-dependent processing tools for continuous-energy (CE) and multigroup (MG) neutronics and coupled neutron-gamma calculations, as well as activation, depletion, and decay calculations. SCALE includes unique capabilities for automated variance reduction for shielding calculations, as well as sensitivity and uncertainty analysis. SCALE’s graphical user interfaces assist with accurate system modeling, visualization of nuclear data, and convenient access to desired results.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
NASA NTRS1973

Further examination of pilot instrument scanning data and development of a new link value estimator

Pilot instrument scanning data collected during simulated transport instrument landing approaches is examined. Deterministic features of random pilot instrument scanning behavior are confirmed to be: transitions in point of eye fixation which originate and terminate on the same instrument are rare; and transitions in point of eye fixation which originate on one secondary instrument and terminate on another secondary instrument are rare. Link value estimators are developed using statistics and these experimental facts. The result has a special significance when there is but a single primary instrument, a flight director. The result is used to simplify the iterative computational procedure of STI's display theory to a non-iterative procedure for the flight director case.

Hofmann, L. G.↗
NASA NTRS1975

Mathematical formulation of scatter-free propagation of solar cosmic rays

The observations of strong persistent velocity anisotropies in solar flare events demand a mathematical theory closer to the extreme of scatter-free (deterministic) propagation rather than diffusive (stochastic) transport, since the latter breaks down as inferred mean-free-paths exceed 0.1 AU. Equations are derived for the time-dependent phase-space density, and Laplace transform techniques are used to obtain solutions under rather general conditions. The case of an Archimedean spiral field has been solved numerically, and the results compared with observations from Mariner and Explorer spacecraft of nearly 0.4 MeV proton intensity and anisotropy histories. These can both be replicated if the inner boundary of the modulation region is placed beyond 2 AU.

Nolte, J. T.↗
NASA NTRS2000

On Geomagnetism and Paleomagnetism I

A partial description of Earth's broad scale, core-source magnetic field has been developed and tested three ways. The description features an expected, or mean, spatial magnetic power spectrum that is approximately inversely proportional to horizontal wavenumber atop Earth's core. This multipole spectrum describes a magnetic energy range; it is not steep enough for Gubbins' magnetic dissipation range. Temporal variations of core multipole powers about mean values are to be expected and are described statistically, via trial probability distribution functions, instead of deterministically, via trial solution of closed transport equations. The distributions considered here are closed and neither require nor prohibit magnetic isotropy. The description is therefore applicable to, and tested against, both dipole and low degree non-dipole fields. In Part 1, a physical basis for an expectation spectrum is developed and checked. The description is then combined with main field models of twentieth century satellite and surface geomagnetic field measurements to make testable predictions of the radius of Earth's core. The predicted core radius is 0.7% above the 3480 km seismological value. Partial descriptions of other planetary dipole fields are noted.

Voorhies, Coerte V.↗
NASA NTRS2010

Proton Lateral Broadening Distribution Comparisons Between GRNTRN, MCNPX, and Laboratory Beam Measurements

Recent developments in NASA s deterministic High charge (Z) and Energy TRaNsport (HZETRN) code have included lateral broadening of primary ion beams due to small-angle multiple Coulomb scattering, and coupling of the ion-nuclear scattering interactions with energy loss and straggling. This new version of HZETRN is based on Green function methods, called GRNTRN, and is suitable for modeling transport with both space environment and laboratory boundary conditions. Multiple scattering processes are a necessary extension to GRNTRN in order to accurately model ion beam experiments, to simulate the physical and biological-effective radiation dose, and to develop new methods and strategies for light ion radiation therapy. In this paper we compare GRNTRN simulations of proton lateral broadening distributions with beam measurements taken at Loma Linda University Proton Therapy Facility. The simulated and measured lateral broadening distributions are compared for a 250 MeV proton beam on aluminum, polyethylene, polystyrene, bone substitute, iron, and lead target materials. The GRNTRN results are also compared to simulations from the Monte Carlo MCNPX code for the same projectile-target combinations described above.

Mertens, Christopher J.↗
DOE OSTI2026

Score-based deterministic density sampling

We propose a deterministic sampling framework using Score-Based Transport Modeling for sampling an unnormalized target density π given only its score ∇ log π. Our method approximates the Wasserstein gradient flow on KL($f_t$∥π) by learning the time-varying score ∇ log $f_t$ on the fly using score matching. While having the same marginal distribution as Langevin dynamics, our method produces smooth deterministic trajectories, resulting in monotone noise-free convergence. We prove that our method dissipates relative entropy at the same rate as the exact gradient flow, provided sufficient training. Numerical experiments validate our theoretical findings: our method converges at the optimal rate, has smooth trajectories, and is often more sample efficient than its stochastic counterpart. Experiments on high-dimensional image data show that our method produces high-quality generations in as few as 15 steps and exhibits natural exploratory behavior. The memory and runtime scale linearly in the sample size.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

A nodal transport method for coupled fast-thermal reactor analysis

We report that for efficient core analyses of coupled fast-thermal reactors, a new deterministic method has been developed based on the variational nodal transport method of VARIANT. A new multigroup cross section generation procedure was devised by combining Monte Carlo lattice calculations for thermal assemblies and the two-step procedure of the MC 2 -3 code for fast assemblies. To reduce assembly homogenization errors in nodal transport calculations, a new nodal equivalence method was developed based on the partial current discontinuity factor (PCDF) and incorporated in VARIANT. A practical procedure to calculate PCDFs with fixed source supercell calculations was also devised. The performance of the proposed method was investigated using a test problem derived from the versatile coupled test reactor (VCTR) design. The new procedure produced multigroup cross sections accurately for all assemblies. VARIANT transport calculations with PCDFs produced accurate multiplication factor and power distribution compared to reference Serpent-2 Monte Carlo solutions.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS1992

Probabilistic design applications for the Space Transportation Main Engines

All of the current rocket engine designs are based on deterministic design analysis methods. The most recently designed engine, Space Shuttle Main Engines (SSME), was designed and developed using these techniques in the 1970's and 80's. The Space Transportation Main Engine (STME) design objective is to use a robust design methodology and to pro-actively assess and improve the design during the design process. The design approach is a combination of the traditional deterministic methodologies and recently developed probabilistic methodologies. The designer must use a more rigorous structural criteria, conduct pro-active reliability assessments (reliability allocations, failure modes and effects analysis, reliability growth studies, and design reliability assessment), use a damage tolerant design philosophy and perform probabilistic analyses. Rigorous application of this structured methodology will develop a better engine product.

Ohara, Kevin↗
DOE OSTI2021

TNSL Overview

Thermal neutron scattering law (TNSL) data describe low-energy neutrons scattering off of bound materials, and can have a significant impact on modeling any system with slow neutrons, including nuclear reactors. Previous work to introduce TNSL data to neutron transport codes at LLNL focused on COG and TART [1], with the limitation that these codes require highly specialized data processing and formatting. We have recently increased efforts to process TNSL data with the central LLNL nuclear data processing code FUDGE, to be stored in the generalized nuclear database structure (GNDS) for use in any general transport code with the ability to read GNDS data. The first step in this effort is to verify that the TNSL processing with FUDGE yields results comparable to results obtained using the LANL nuclear data processing code NJOY. The next step is to verify the transport of thermal neutrons in Mercury (a Monte Carlo code) and Ardra (a deterministic code) against one another, as well as against the LANL Monte Carlo neutron transport code MCNP. This verification step has not been completed, due to a number of discrepancies between results obtained using differently processed data. There is ongoing effort to understand differences between FUDGE and NJOY. Finally, we map out our current capability to validate TNSL data against benchmark systems.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Importance of 3-D S{sub N} depletion in non-proliferation using BSOLVE

We present a single Pressurized Water Reactor (PWR) 3-D fuel rod design for depletion analysis using BSOLVE, our newly developed Runge-Kutta-Fehlberg based depletion code. BSOLVE is coupled with the deterministic 3-D S{sub N} particle transport code, PENTRAN, applied here with a 4-neutron energy group comparison to Continuous Energy (C/E) SERPENT2 Monte Carlo results. Differences are expected, as PENTRAN+BSOLVE retains full (multi-group) energy information for reactions, nuclide specific fission contributions, and energy dependent fission yields, using the latest available ENDF-BVIII data, important to retain accurate burned fuel inventories; SERPENT2 collapses burnup reactions to a single energy value. For depletion times up to ∼ 700 days and typical PWR power densities, relative differences between multigroup 3-D S{sub N} with full energy data and Monte Carlo one group burnup for trans-uranium nuclide concentrations and fission products are up to ∼20%. System eigenvalues are consistent, but with differences early and late in the cycle attributed to multigroup vs. C/E Monte Carlo cross sections. This work highlights the importance of low variance transport driven burnup for non-proliferation concerns, since plutonium quality varies significantly along axial lengths, and is more challenging to converge using Monte Carlo; details of depletion steps with spatial/zone dependent plutonium quality are provided. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2025

Incorporating geological structure into sensitivity analysis of subsurface contaminant transport

Simulating subsurface contaminant transport at the kilometer-scale often entails modeling reactive flow and transport within and through complex geologic structures. These structures are typically meshed by hand and as a result geologic structure is usually represented by one or a few deterministically generated geological models for uncertainty studies of flow and transport in the subsurface. Uncertainty in geologic structure can have a significant impact on contaminant transport. In this study, the impact of geologic structure on contaminant tracer transport in a shale formation is investigated for a simplified generic deep geologic repository for permanent disposal of spent nuclear fuel. An open-source modeling framework is used to perform a sensitivity analysis study on transport of two tracers from a generic spent nuclear fuel repository with uncertain location of the interfaces between the stratum of the geologic structure. The automated workflow uses sampled realizations of the geological structural model in addition to uncertain flow parameters in a nested sensitivity analysis. Concentration of the tracers at observation points within, in line with, and downstream of the repository are used as the quantities of interest for determining model sensitivity to input parameters and geological realization. Finally, the results of the study indicate that the location of strata interfaces in the geological structure has a first-order impact on tracer transport in the example shale formation, and that this impact may be greater than that of the uncertain flow parameters.

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

Software Control Program For Transportable Microgrid State-of-charge Balancing And Frequency Stability Controls

A deterministic state-of-charge (SOC) balancing approach software control code is introduced as an integral secondary management to primary control layer of an islanded small microgrid or nanogrid system made up of multiple grid-forming inverter/battery/solar combination systems, where each set of batteries with each inverter are on independent DC buses (i.e. non-paralleled on the DC sides). A DERMS-level control approach, algorithm and automation controller program was developed to improve coordination and enable microgrid asset compliance and SOC balancing, enabling provision of a system-level power stability support architecture, load support, and asset scalability. The architecture is configured to treat each unit or micro/nano-grid as a node in a microgrid network, allowing for autonomous DERMS control regarding load and SOC balancing and power stability. As the network grows with the addition of units, greater coordination efforts may be required. The ideal small network microgrid ranges from 2-10 inverter/battery units before additional control parameters must be considered in the existing architecture. The control approach focuses on a deterministic state-of-charge analysis as the primary level control process followed by a secondary control loop using a forced frequency-watt droop strategy to conform off-the-shelf components into behaving under a leader-follower configuration. Adopting this control scheme has been shown to allow for a balanced, unit-coordinated microgrid network, enabling stable power flow. The deterministic state-of-charge approach is introduced as an integral primary control layer of an islanded small network microgrid. A standard strategy for SOC balancing is implementing a battery management system (BMS) to control SOC on the DC side. An alternative approach is to determine how to coordinate sending and receiving power on the AC side with multiple units. The latter approach assesses all the integrated units in the microgrid network. Once the individual units are identified, further system data is required to calculate each unit's total kWh, provided information about its capability to supply or consume kWh and availability. The secondary control layer in the multi-layered small network microgrid methodology uses the primary layer’s decision to initiate frequency setpoint changes, initializing the SOC balancing. The secondary control layer considers numerous system-dependent variables to enable a charging and discharging profile based on adjustable frequency setpoints. The combined architecture will result in stable, coordinated power flow enhancing an AC microgrid's functionalities.

Myers, KurtS [Idaho National Laboratory (INL), Ida↗
DOE OSTI2020

Learning and meta-learning of stochastic advection–diffusion–reaction systems from sparse measurements

Physics-informed neural networks (PINNs) were recently proposed in [18] as an alternative way to solve partial differential equations (PDEs). A neural network (NN) represents the solution, while a PDE-induced NN is coupled to the solution NN, and all differential operators are treated using automatic differentiation. Here, we first employ the standard PINN and a stochastic version, sPINN, to solve forward and inverse problems governed by a non-linear advection–diffusion–reaction (ADR) equation, assuming we have some sparse measurements of the concentration field at random or pre-selected locations. Subsequently, we attempt to optimise the hyper-parameters of sPINN by using the Bayesian optimisation method (meta-learning) and compare the results with the empirically selected hyper-parameters of sPINN. In particular, for the first part in solving the inverse deterministic ADR, we assume that we only have a few high-fidelity measurements, whereas the rest of the data is of lower fidelity. Hence, the PINN is trained using a composite multi-fidelity network, first introduced in [12], that learns the correlations between the multi-fidelity data and predicts the unknown values of diffusivity, transport velocity and two reaction constants as well as the concentration field. For the stochastic ADR, we employ a Karhunen–Loève (KL) expansion to represent the stochastic diffusivity, and arbitrary polynomial chaos (aPC) to represent the stochastic solution. Correspondingly, we design multiple NNs to represent the mean of the solution and learn each aPC mode separately, whereas we employ a separate NN to represent the mean of diffusivity and another NN to learn all modes of the KL expansion. For the inverse problem, in addition to stochastic diffusivity and concentration fields, we also aim to obtain the (unknown) deterministic values of transport velocity and reaction constants. The available data correspond to 7spatial points for the diffusivity and 20 space–time points for the solution, both sampled 2000 times. We obtain good accuracy for the deterministic parameters of the order of 1–2% and excellent accuracy for the mean and variance of the stochastic fields, better than three digits of accuracy. In the second part, we consider the previous stochastic inverse problem, and we use Bayesian optimisation to find five hyper-parameters of sPINN, namely the width, depth and learning rate of two NNs for learning the modes. Here, we obtain much deeper and wider optimal NNs compared to the manual tuning, leading to even better accuracy, i.e., errors less than 1% for the deterministic values, and about an order of magnitude less for the stochastic fields.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Periodic Porous 3D Electrodes Mitigate Gas Bubble Traffic during Alkaline Water Electrolysis at High Current Densities

We report alkaline water electrolysis at high current densities is plagued by gas bubble generation and trapping in stochastic porous electrodes (e.g., Ni foams), which causes a significant reduction in the number of electrolyte accessible catalyst active sites. Here, 3D printed Ni (3DPNi) electrodes with highly controlled, periodic structures are reported that suppress gas bubble coalescence, jamming, and trapping and, hence, result in rapid bubble release. The 3DPNi electrodes decorated with carbon-doped NiO achieve a high current density of 1000 mA cm -2 in 1.0 m KOH electrolyte at hydrogen evolution reaction and oxygen evolution reaction overpotentials of 245 and 425 mV, respectively. This work demonstrates a new approach to the deterministic design of 3D electrodes to facilitate rapid bubble transport and release to enhance the total electrode catalytic activity at commercially relevant current densities.

08 HYDROGEN↗
DOE OSTI2026

A hybrid Monte Carlo-deterministic second moment method with efficient variance reduction

In this work, we present a hybrid method that combines Monte Carlo with deterministic finite element methods to solve a linear Boltzmann transport equation. Our hybrid method runs orders of magnitude faster than Monte Carlo, without sacrificing accuracy, for a proxy problem from radiative transfer that contains both optically-thick and optically-thin material. We believe that this is the first demonstration of a hybrid Second Moment Method in more than one spatial dimension, the first to consider more than one material, and the first to use variance reduction. Our variance reduction approach arises from an asymptotic analysis in which we show that the magnitude of the scattering source grows without bound. We transform the problem to compute the deviation of the radiation intensity from isotropy. The magnitude of the source in the transformed problem is bounded, and the quality of the hybrid method solution is dramatically improved by a substantial reduction in the variance.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗