Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “model code”

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 1,261 records · Page 70
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Cross sections and calculated yields of some radionuclides of yttrium, strontium and rubidium formed in proton-induced reactions on enriched strontium-86: possibility of production of 85g Sr, 83 Rb and 82m Rb in no-carrier-added form

Here, cross sections of the 86 Sr(p,3n) 84m Y, 86 Sr(p,αn) 82m Rb, and 86 Sr(p,x) 85g Sr reactions were measured from their respective thresholds up to 16.2 MeV and from 23.0 to 44.1 MeV at FZJ, and from 14.3 to 24.5 MeV at LBNL, using 96.4% enriched 86 SrCO 3 as target material. Thin targets prepared by sedimentation were irradiated with protons in a stacked-form, and the induced radioactivity was measured by high-resolution γ-ray spectrometry. Nuclear model calculations based on the code TALYS reproduced our experimental cross section data well. From the excitation functions, the integral yields of the above three radionuclides were calculated. The yield of 85g Sr via the natSr(n,γ) process was also measured using the TRIGA Mark-II reactor at AERE, Savar. A comparison of the reactor and cyclotron production of carrier-added 85g Sr is given. The production possibilities of the three investigated radionuclides in no-carrier-added forms at a 30 MeV cyclotron via new routes are discussed.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2023

Flexible, integrated modeling of tokamak stability, transport, equilibrium, and pedestal physics

The STEP (Stability, Transport, Equilibrium, and Pedestal) integrated-modeling tool has been developed in OMFIT to predict stable, tokamak equilibria self-consistently with core-transport and pedestal calculations. STEP couples theory-based codes to integrate a variety of physics, including magnetohydrodynamic stability, transport, equilibrium, pedestal formation, and current-drive, heating, and fueling. The input/output of each code is interfaced with a centralized ITER-Integrated Modelling & Analysis Suite data structure, allowing codes to be run in any order and enabling open-loop, feedback, and optimization workflows. This paradigm simplifies the integration of new codes, making STEP highly extensible. STEP has been verified against a published benchmark of six different integrated models. Core-pedestal calculations with STEP have been successfully validated against individual DIII-D H-mode discharges and across more than 500 discharges of the H98,y2 database, with a mean error in confinement time from experiment less than 19%. STEP has also reproduced results in less conventional DIII-D scenarios, including negative-central-shear and negative-triangularity plasmas. Predictive STEP modeling has been used to assess performance in several tokamak reactors. Simulations of a high-field, large-aspect-ratio reactor show significantly lower fusion power than predicted by a zero-dimensional study, demonstrating the limitations of scaling-law extrapolations. STEP predictions have found promising scenarios for an EXhaust and Confinement Integration Tokamak Experiment, including a high-pressure, 80%-bootstrap-fraction plasma. ITER modeling with STEP has shown that pellet fueling enhances fusion gain in both the baseline and advanced-inductive scenarios. Finally, STEP predictions for the SPARC baseline scenario are in good agreement with published results from the physics basis.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Nonlocal effects on thermal transport in hydrodynamic simulations of unmagnetized MagLIF-relevant gaspipes on NIF

We present simulations of heat flow relevant to gaspipe experiments on the National Ignition Facility to investigate kinetic effects on transport phenomena. D 2 and neopentane (C 5 H 12 ) filled targets are used to study the laser preheat stage of a MagLIF scheme where an axial magnetic field is sometimes applied to the target. Simulations were done with the radiation-MHD code HYDRA with a collision-dominated fluid model and the SNB nonlocal electron thermal conduction model. Using the SNB model to evolve the electron temperature increased the heat front propagation of neopentane gas targets compared to a local model by limiting radial heat flow. This increases electron temperature near the axis, which decreases laser absorption. We find that the effect of heat flow models on temperature profiles and laser propagation is modest. Beyond the SNB model, we utilize HYDRA to initialize plasma conditions for the Vlasov–Fokker–Planck K2 code. We run K2 until a quasi-steady state is reached and examine the impact of kinetic effects on heat transport. Although axial heat flow is well predicted by fluid models, the fluid model consistently overpredicts radial heat flow up to 150% in regions with the largest temperature gradient of D 2 filled gaspipes. On the other hand, the SNB nonlocal electron conduction model is found to be adequate for capturing kinetic heat flow in gaspipes.

Lau, Ryan Y. [Univ. of Colorado, Boulder, CO (Unit↗
DOE OSTI2026

Differential Privacy in Grid Kitchen: Implementation & Software Documentation

Sharing of power grid feeder models faces significant challenges due to the potential risk of exposing sensitive operational information. Traditional anonymization techniques have shown notable limitations in other sensitive domains, as evidenced by documented re-identification attacks that combine supposedly anonymized datasets with auxiliary information, raising concerns that similar vulnerabilities could affect power grid data. Consequently, there is a pressing need for a more rigorous privacy protection strategy that not only delivers formal mathematical guarantees but also preserves the analytical value of the shared models. To address this challenge, we have enhanced the Grid Kitchen framework by implementing differential privacy mechanisms within the distribution model dehydration pipeline. This implementation carefully calibrates and applies noise to sensitive attributes in feeder models according to configurable privacy levels—low, moderate, and high—each offering different balances between data utility and privacy protection. Our approach uses established noise functions (Gaussian for continuous data and Discrete Laplace for integer values) with parameters carefully calibrated so that the impact of individual data points is effectively masked in the final output. The integration leverages our Noise Catalog, which we developed to categorize feeder model properties by component type, data type, and sensitivity. This catalog guides the application of appropriate noise functions and privacy parameters ($\varepsilon$ and $\delta$) to each attribute, ensuring consistent privacy protection across the model while maintaining its structural integrity and analytical usefulness. This implementation also includes evaluation tools that allow model owners to assess the impact of privacy-preserving transformations before sharing data with external parties. This report provides documentation for the differential privacy capabilities added to the Grid Kitchen project. It includes a primer on differential privacy concepts and their importance in modern data sharing, details the architecture of our implementation, explains the privacy modes and parameter configurations, and offers practical guidance on using the code for applying differential privacy to grid feeder models. Through examples and code snippets, we demonstrate the effective application of these privacy-enhancing technologies, enabling utility operators and researchers to confidently share grid data while protecting sensitive information.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2020

Accurate determination of production data of the non-standard positron emitter 86 Y via the 86 Sr(p,n)-reaction

In view of several significant discrepancies in the excitation function of the 86 Sr(p,n) 86g+xm Y reaction which is the method of choice for the production of the non-standard positron emitter 86 Y for theranostic application, we carried out a careful measurement of the cross sections of this reaction from its threshold up to 16.2 MeV at Forschungszentrum Jülich (FZJ) and from 14.3 to 24.5 MeV at LBNL. Thin samples of 96.4% enriched 86 SrCO 3 were prepared by sedimentation and, after irradiation with protons in a stacked-form, the induced radioactivity was measured by high-resolution γ -ray spectrometry. The projectile flux was determined by using the monitor reactions nat Cu(p,xn) 62,63,65 Zn and nat Ti(p,x) 48 V, and the calculated proton energy for each sample was verified by considering the ratios of two reaction products of different thresholds. Additionally, the experimental cross section data obtained agreed well with the results of a nuclear model calculation based on the code TALYS. From the cross section data, the integral yield of 86 Y was calculated. Over the optimum production energy range E p = 14 → 7 MeV the yield of 86 Y amounts to 291 MBq/μA for 1 h irradiation time. This value is appreciably lower than the previous literature values calculated from measured and evaluated excitation functions. It is, however, more compatible with the experimental yields of 86 Y obtained in clinical scale production runs. The levels of the isotopic impurities 87m Y, 87g Y, and 88 Y were also estimated and found to be <2% in sum.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2026

Toward Energy-Efficient HPC: Insights from Power Profiling a Cloud-Resolving Earth System Model

Power is a fundamental constraint as supercomputing advances to exascale. Efficient operation within strict power budgets requires application-aware power management based on a detailed understanding of application-level power behavior. This work analyzes the Energy Exascale Earth System Model (E3SM) atmosphere component, SCREAM, on Perlmutter (NERSC) and Frontier (OLCF). We characterize power variation across inputs, concurrency levels, and power caps, evaluate the energy impact of code optimizations, and attribute energy within the code using a newly developed GPU energy model. Results show that SCREAM’s peak power remains stable during its core execution phase and decreases gradually as concurrency increases. Power capping experiments reveal a performance–energy "sweet spot". On Perlmutter, limiting GPU power to 50% of thermal design power (TDP) achieves up to 15% energy savings with a 7% performance penalty. On Frontier, a 40% TDP cap yields up to 10% energy savings with less than 10% performance loss. Code optimizations reduce SCREAM energy by shortening run time without increasing power. Modeling reveals a critical insight: data movement accounts for approximately 70% of SCREAM’s GPU energy. This fundamentally shifts the optimization focus from FLOPS to data transfer reduction for this class of applications, offering the most impactful strategy for improving energy efficiency. This work establishes a foundation for practical, application-aware power management at exascale.

Zhao, Zhengji [Lawrence Berkeley National Laborato↗
DOE OSTI2023

An extended Vlasov–Fokker–Planck approach for kinetic simulations of laser plasmas

Vlasov–Fokker–Planck simulation codes occupy an important niche in modeling laser-produced plasmas, since they are well suited to studying the effect of collisions on electron kinetic phenomena, especially energy transport. One of the most important elements of energy transport is the absorption of laser light by the plasma; however, simulating this in detail requires resolving oscillations of the laser light, whose characteristic timescale is orders of magnitude shorter than the simulation time needed to study transport physics. For this reason, most Vlasov–Fokker–Planck codes used to study electron transport in laser plasmas rely on simplified models of the laser–plasma coupling. Their underlying assumptions nominally preclude their use for modeling laser light having short-scale structure in space or time, such as broadband lasers. In this work, we derive a more general computational framework suitable for arbitrarily structured laser fields. Furthermore, our approach is based on an extended set of Vlasov–Fokker–Planck equations that separately solve for the low- and high-frequency plasma response. We implement these extended Vlasov–Fokker–Planck equations in the spherical harmonic code K2 and demonstrate the performance of the method on several laser absorption test problems, with particular attention to the judicious selection of time steps, time integrators, and spherical harmonic truncation, according to the intensity and spectrum of the laser light under consideration. Comparison with the widely used Langdon absorption operator shows the Langdon operator performs remarkably well for predicting laser heating in the simple cases considered here, even in situations that would seem to violate its underlying assumptions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Lifetimes of two- and four-neutron-hole states in 204 Pb and 206 Pb

The centroid-shift method was used to determine half-lives of excited levels in 204 Pb and 206 Pb. For the states with 𝐼 𝜋 = 12 + , 17 − , and 19 − in 204 Pb, values of 𝑇 1/2 = 4.9⁢(3), 3.9(3), and 0.9(3) ns, respectively, were determined. The half-life of the 𝐼 𝜋 = 9 − level in 206 Pb was measured to be 1.1(3) ns. In 204 Pb, the 12 + state arises from two neutron holes in the 𝑖 13/2 subshell, while the 17 − and 19 − levels are built on 𝜈⁡(𝑖$^{−3}_{13/2}$, 𝑝$^{−1}_{1/2}$) and 𝜈⁡(𝑖$^{−3}_{13/2}$, 𝑓$^{−1}_{5/2}$) four-neutron-hole configurations, respectively. The 9 − state in 206 Pb results from a 𝜈⁡(𝑖$^{−1}_{13/2}$, 𝑓${−1}_{5/2}$) excitation. As a result, the transition probabilities for the decay of these isomeric levels are reasonably consistent with the expected single-particle values, and these states are described well by shell-model calculations using the kshell code with the KHHE interaction.

190 ≤ A ≤ 219↗
DOE OSTI2024

The Santa Barbara Binary-disk Code Comparison

We have performed numerical calculations of a binary interacting with a gas disk, using 11 different numerical methods and a standard binary-disk setup. The goal of this study is to determine whether all codes agree on a numerically converged solution and to determine the necessary resolution for convergence and the number of binary orbits that must be computed to reach an agreed-upon relaxed state of the binary-disk system. We find that all codes can agree on a converged solution (depending on the diagnostic being measured). The zone spacing required for most codes to reach a converged measurement of the torques applied to the binary by the disk is roughly 1% of the binary separation in the vicinity of the binary components. For our disk model to reach a relaxed state, codes must be run for at least 200 binary orbits, corresponding to about a viscous time for our parameters, 0.2(a 2 Ω B /v) is the kinematic viscosity. The largest discrepancies between codes resulted from the dimensionality of the setup (3D vs. 2D disks). We find good agreement in the total torque on the binary between codes, although the partition of this torque between the gravitational torque, orbital accretion torque, and spin accretion torque depends sensitively on the sink prescriptions employed. In agreement with previous studies, we find a modest difference in torques and accretion variability between 2D and 3D disk models. We find cavity precession rates to be appreciably faster in 3D than in 2D.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

The development of a high-resolution Eulerian radiation-hydrodynamics simulation capability for laser-driven Hohlraums

Hohlraums are hollow cylindrical cavities with high-Z material walls used to convert laser energy into uniform x-ray radiation drives for inertial confinement fusion capsule implosions and high energy density physics experiments. Credible computational modeling of hohlraums requires detailed modeling and coupling of laser physics, hydrodynamics, radiation transport, heat transport, and atomic physics. We report on improvements to Los Alamos National Laboratory's xRAGE radiation-hydrodynamics code in order to enable hohlraum modeling. xRAGE's Eulerian hydrodynamics and adaptive mesh refinement make it uniquely well suited to study the impacts of multiscale features in hohlraums. In order to provide confidence in this new modeling capability, we demonstrate xRAGE's ability to produce reasonable agreement with data from several benchmark hohlraum experiments. We also use xRAGE to perform integrated simulations of a recent layered high density carbon capsule implosion on the National Ignition Facility in order to evaluate the potential impacts of the capsule support tent, mixed cell conductivity methodologies, plasma transport, and cross-beam energy transfer (XBT). We find that XBT, seeded by plasma flows in the laser entrance hole (LEH), causes a slight decrease in energy coupling to the capsule and that all of these impact the symmetry of the x-ray drive such that they have an appreciable impact on the capsule implosion shape.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

MSTAR2019 Code Description and User's Manual

This report describes an ideal cascade model of uranium enrichment that includes side feed and side product streams and flexible input options. It is based on the original MSTAR model developed by Ed Von Halle and implemented in a Visual Basic code. The current version allows the user to specify integer numbers of stages or the stage numbers of external flows instead of assays in those flows. The computational engine is written in FORTRAN-90 and is invoked by a user-friendly GUI written in C++. This version of MSTAR has been demonstrated to operate on Linux, Mac, and Windows platforms, and has undergone significant testing and quality analysis. A number of examples are presented to illustrate the operation of the code and guide the user. The code is extremely fast, and results are returned immediately. The input and output have been specially configured for analysis by the environmental sampling team of the International Atomic Energy Agency (IAEA). A number of additional output features are included to assist the user in visualizing computational results and downloading data to files for use in other analysis software.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2023

MACCS User Guide (V.4.2)

MACCS is used by the Nuclear Regulatory Commission (NRC) and various national and international organizations for probabilistic consequence analysis of nuclear power accidents. This User Guide is intended to assist analysts in understanding the MACCS/WinMACCS model and to provide information regarding the code. This user guide version describes MACCS Version 4.2. This User Guide provides a brief description of the model history, explains how to set up and execute a problem, and informs the user of the definition of various input parameters and any constraints placed on those parameters. This report is part of a series of reports documenting MACCS. Other reports include the MACCS Theory Manual, MACCS Verification Report, Technical Bases for Consequence Analyses Using MACCS, as well as documentation for preprocessor codes including SecPop, MelMACCS, and COMIDA2.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2021

Electron Dynamics Within a MITL Containing a Load

In this paper, we derive the vacuum electric fields within specific cylindrically symmetric magnetically insulated transmission lines (MITLs) in the limit of an infinite speed of light for an arbitrary time-dependent current. We focus our attention on two types of MITLs: the radial MITL and a spherically curved MITL. We then simulate the motion of charged particles, such as electrons, present in these MITLs due to the vacuum fields. In general, the motion of charged particles due to the vacuum fields is highly nonlinear since the fields are nonlinear functions of spatial coordinates and depend on an arbitrary time-dependent current drive. Using guiding center theory, however, one can describe the gross particle kinetics using a combination of $\textbf {E} \times \textbf {B}$ and $\nabla B$ drifts. In addition, we compare our approximate inner MITL field models and particle kinetics with those from a fully electromagnetic simulation code. We find that the agreement between the approximate model and the electromagnetic simulations is excellent.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Temporal evolution of the light emitted by a thin, laser-ionized plasma source

We present an experimental and simulation-based investigation of the temporal evolution of light emission from a thin, laser-ionized helium plasma source. We demonstrate an analytic model to calculate the approximate scaling of the time-integrated, on-axis light emission with the initial plasma density and temperature, supported by the experiment, which enhances the understanding of plasma light measurement for plasma wakefield accelerator (PWFA) plasma sources. Our model simulates the plasma density and temperature using a split-step Fourier code and a particle-in-cell code. A fluid simulation is then used to model the plasma and neutral density, and the electron temperature as a function of time and position. We then show the numerical results of the space-and-time-resolved light emission and that collisional excitation is the dominant source of light emission. We validate our model by measuring the light emitted by a laser-ionized plasma using a novel statistical method capable of resolving the nanosecond-scale temporal dynamics of the plasma light using a cost-effective camera with microsecond-scale timing jitter. This method is ideal for deployment in the high radiation environment of a particle accelerator that precludes the use of expensive nanosecond-gated cameras. Our results show that our models can effectively simulate the dynamics of a thin, laser-ionized plasma source. In addition, this work provides a detailed understanding of the plasma light measurement, which is one of the few diagnostic signals available for the direct measurement of PWFA plasma sources.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Feedback density and causal complexity of simulation model structure

Measures of simulation model complexity generally focus on outputs; we propose measuring the complexity of a model’s causal structure to gain insight into its fundamental character. This article introduces tools for measuring causal complexity. First, we introduce a method for developing a model’s causal structure diagram, which characterises the causal interactions present in the code. Causal structure diagrams facilitate comparison of simulation models, including those from different paradigms. Next, we develop metrics for evaluating a model’s causal complexity using its causal structure diagram. We discuss cyclomatic complexity as a measure of the intricacy of causal structure and introduce two new metrics that incorporate the concept of feedback, a fundamental component of causal structure. The first new metric introduced here is feedback density, a measure of the cycle-based interconnectedness of causal structure. The second metric combines cyclomatic complexity and feedback density into a comprehensive causal complexity measure. Finally, we demonstrate these complexity metrics on simulation models from multiple paradigms and discuss potential uses and interpretations. These tools enable direct comparison of models across paradigms and provide a mechanism for measuring and discussing complexity based on a model’s fundamental assumptions and design.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Interpretive modeling of tungsten divertor leakage during experiments with neon gas seeding

Abstract Many existing and future tokamaks with tungsten divertors operate, or will operate, with low- Z impurity seeding, but the direct effect of these seeded impurities on tungsten Scrape-off-Layer (SOL) transport has not been explored in detail. This paper reports on a DIII-D experiment designed to test how tungsten divertor leakage from the Small-Angle Slot V-Shaped, tungsten-coated divertor is impacted by neon seeding at a variety of injection rates and poloidal injection locations. Measurements from the experiment show an inverse relationship between the neon injection rate and the tungsten core penetration factor. Interpretive modeling is performed with a combination of the SOLPS-ITER and DIVIMP codes to assess the underlying tungsten behavior. The modeling results show that the reduction in tungsten divertor leakage is driven by both an increase in the divertor collisionality as well as a reduction in the ion temperature gradient near the divertor target. Collisions between low- Z impurities and tungsten impurities are found to have a significant impact on the tungsten SOL transport, such that ignoring the low- Z impurity collisional effects on the tungsten transport can result in an overestimate of the divertor leakage by an order-of-magnitude. Given the importance of these localized interactions, neon seeding from the closed, slot-like divertor has a clear advantage in being able to reduce tungsten divertor leakage without the high levels of neon core contamination that occur when seeding from other poloidal locations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Gapped lineon and fracton models on graphs

We introduce a $\mathbb{Z}_N$ stabilizer code that can be defined on any spatial lattice of the form $\Gamma\times C_{L_z}$, where $\Gamma$ is a general graph. We also present the low-energy limit of this stabilizer code as a Euclidean lattice action, which we refer to as the anisotropic $\mathbb{Z}_N$ Laplacian model. It is gapped, robust (i.e., stable under small deformations), and has lineons. Its ground state degeneracy (GSD) is expressed in terms of a "mod $N$-reduction" of the Jacobian group of the graph $\Gamma$. In the special case when space is an $L\times L\times L_z$ cubic lattice, the logarithm of the GSD depends on $L$ in an erratic way and grows no faster than $O(L)$. Furthermore, we also discuss another gapped model, the $\mathbb{Z}_N$ Laplacian model, which can be defined on any graph. It has fractons and a similarly strange GSD.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Supporting hierarchical soil biogeochemical modeling: version 2 of the Biogeochemical Transport and Reaction model (BeTR-v2)

Abstract. Reliable soil biogeochemical modeling is a prerequisite for credible projections of climate change and associated ecosystem feedbacks. This recognition has called for frameworks that can support flexible and efficient development and application of new or alternative soil biogeochemical modules in Earth system models (ESMs). The the Biogeochemical Transport and Reaction model version 1 (BeTR-v1) code (i.e., CLM4-BeTR) is one such framework designed to accelerate the development and integration of new soil biogeochemistry formulations into ESMs and to analyze structural uncertainty in ESM simulations. With a generic reactive transport capability, BeTR-v1 can represent multiphase (e.g., gaseous, aqueous, and solid), multi-tracer (e.g., nitrate and organic carbon), and multi-organism (e.g., plants, bacteria, and fungi) dynamics. Here, we describe the new version, Biogeochemical Transport and Reaction model version 2 (BeTR-v2), which adopts more robust numerical solvers for multiphase diffusion and advection and coupling between biogeochemical reactions and improves code modularization over BeTR-v1. BeTR-v2 better supports different mathematical formulations in a hierarchical manner by allowing the resultant model be run for a single topsoil layer or a vertically resolved soil column, and it allows the model to be fully coupled with the land component of the Energy Exascale Earth System Model (E3SM). We demonstrate the capability of BeTR-v2 with benchmark cases and example soil biogeochemical (BGC) implementations. By taking advantage of BeTR-v2's generic structure integrated in E3SM, we then found that calibration could not resolve biases introduced by different numerical coupling strategies of plant–soil biogeochemistry. These results highlight the importance of numerically robust implementation of soil biogeochemistry and coupling with hydrology, thermal dynamics, and plants – capabilities that the open-source BeTR-v2 provides. We contend that Earth system models should strive to minimize this uncertainty by applying better numerical solvers.

58 GEOSCIENCES↗