Search NASA⌕ Search

SEARCH · Search NASA

Results for “Simulation Use Case”

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 397 records · Page 22
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Determining Stability Margins in Adiabatic Superconducting Magnets with 3-D Finite Element Analysis

Superconducting magnets play a key role in the development of experiments at Fermilab; understanding the operating stability of these can allow us to utilize more potent magnets for future experiments (like the proposed Muon Collider), optimize the design of magnets in more immediate experiments (like Mu2e), and research the future use of more exotic materials (like high-temperature superconductors). This summer, I developed a 3-D parametric FEA program in ANSYS Mechanical APDL that simulates quench in superconducting magnets, and I also developed a parametric MATLAB program that predicts thermal behavior in magnet quench using the MIITS method. These programs can provide useful quenching parameters (like minimum quench energy and normal zone propagation velocity) for different cases of quench, leading to the previously mentioned objective of magnet design optimization. To test the programs, preliminary cases were run and the data produced was compared and analyzed. The results of these analyses, as well as the program operating methods, are discussed in this project.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

A Green's Function Wind Turbine Induction Model That Incorporates Complex Inflow Conditions

ABSTRACT In this work, we develop a new analytical turbine induction model that can incorporate complex inflow conditions including cases where the wind velocity and temperature profiles can vary as functions of height. This induction model is derived from the linearized Navier–Stokes and leads to a second‐order ODE that can be solved using a Green's function formulation. The corresponding Green's function for several configurations are found including the infinite domain, semi‐infinite domain with ground plane, and a power law velocity inflow profile. The results of this approach are then compared with simulations of the turbine induction field using the AMR‐Wind CFD solver with a uniformly loaded actuator disk model. These comparisons show that the Green's function approach captures the centerline blockage, three‐dimensional blockage flow field, and streamwise velocity slow down, with very good agreement for lower thrust conditions and at larger distances away from rotor disk. The effects of shear on the turbine blockage were also compared using a power law inflow profile, and we show that this approach matches the CFD predictions for the cases considered.

17 WIND ENERGY↗
DOE OSTI2025

Impact of nitrogen molecular breakup on divertor conditions in JET L-mode plasmas using SOLPS-ITER

SOLPS-ITER simulations of nitrogen-seeded, low-confinement mode plasmas in the Joint European Torus (JET) predict that the electron temperature in the low-field side (LFS) divertor leg is reduced locally by up to an order of magnitude when nitrogen is assumed to recycle as molecules (N 2 ) instead of atoms using a fixed nitrogen injection rate. The LFS divertor temperature reduction under the assumption of molecular recycling occurs due to a three-step mechanism: (1) the plasma penetration of nitrogen atoms is increased due to the strong triple bond of the N 2 molecule and the kinetic energy release in the dissociation event, both mechanisms contributing equally, (2) the abundance of (particularly multiply-charged) nitrogen ions in the divertor is increased and (3) the electron temperature is reduced due to the increase in radiation (by up to a factor of 4) from nitrogen ions. Setting the volume-integrated nitrogen radiated power to a constant value (0.6 MW) instead of the nitrogen injection rate, SOLPS-ITER predicts under the molecular nitrogen recycling assumption that the peak line-integrated N II, N III and N IV intensities in the LFS divertor are approximately within 15%, 35% and 5%, respectively, of the reference atomic nitrogen recycling case. The predicted peak N II, N III and N IV intensities under either assumption are within 30%, 65% and 5%, respectively, of measurements using the vertically viewing mirror-link divertor spectrometer (Meigs et al., 2010) in nitrogen-seeded JET L-mode plasmas (Lomanowski et al., 2019). ERO2.0 simulations using a constant nitrogen seeding rate on static background plasma solutions from EDGE2D-EIRENE (previously presented in Mäenpää et al., (2022), revised here to include fast reflections) predict that N II to N IV line emission is increased by 20% to 30% when nitrogen is assumed to recycle as molecules, demonstrating the importance of considering the effect of molecular dissociation reactions on the divertor plasma in a self-consistent manner.

Divertor↗
DOE OSTI2025

HydroEcoLSTM: A Python package with graphical user interface for hydro-ecological modeling with long short-term memory neural network

Machine learning (ML) is emerging as a promising tool for modeling hydro-ecological processes due to the increasing availability of large environmental data. However, the use of ML requires sufficient programming knowledge due to a lack of a graphical user interface (GUI). In this study, we introduced a GUI package, named HydroEcoLSTM, with the long short-term memory network (LSTM) as the core model, that allows non-ML experts to utilize their domain knowledge to construct complex ML models. We demonstrated the functionalities of HydroEcoLSTM with two practical examples, including (1) predictions of streamflow in both gauged and ungauged catchments and (2) predictions of multiple outputs (i.e., streamflow and isotope transport from two catchments). The simulation results obtained in both case experiments are satisfactory. In the first example, the average Nash–Sutcliffe Efficiency (NSE) for streamflow simulation during the testing period is 0.79 while the application of the trained model in two assumed ungauged catchments also achieves the average NSE of 0.68. In the second example, the average NSE for streamflow and instream isotope simulation during the testing period is 0.71. Ultimately, applications of HydroEcoLSTM with real-world examples demonstrate its potential use for practical applications and research without requiring extensive coding skills.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Accelerating magnonic simulations with the pseudospectral Landau-Lifshitz equation

The pseudospectral Landau-Lifshitz (PS-LL) model can describe atomic-scale magnetic exchange interactions within a continuum framework. This is achieved by employing a convolution kernel that models the nonlocal interaction in a grid-independent manner. Even though the PS-LL was originally introduced to address atomic exchange, any nonlocal kernel can be modeled. In the field of magnonics, the dipole field is fundamental to describe the dispersion relation of magnons, the quasiparticle representation of angular momentum. Because dipole-dipole interactions are long-range, numerical approaches typically rely on convolutions. Here, we demonstrate that the PS-LL model can be used to perform magnonic simulations with a single convolution kernel derived from analytical solutions. We demonstrate a twofold increase in computational speed compared with the full dipole calculation. This approach is valid insofar as the excitations are linear, which is typically the case for magnons. Our results have the potential to accelerate magnonic research, particularly for the inverse design method, where several simulations must be performed to achieve the desired outcome.

Mathematics and computing↗
DOE OSTI2025

Deciphering yield modification of hadron-triggered semi-inclusive recoil jets in heavy-ion collisions

Jet quenching is recognized as critical evidence for the existence of the quark-gluon plasma (QGP) and serves as an essential probe to study its transport properties. Measurements of hadron-triggered semi-inclusive recoil jets have gained popularity due to its capability to probe jets over an extended phase space at low transverse momenta (p T ) and large radii. Recent ALICE measurements showed that the I AA , yield ratio of recoil jets between heavy-ion and p+p collisions, rises with jet p T and exceeds unity at high p T , contradicting conventional expectations that jet quenching should result in I AA values less than one. In this contribution, we re-examine the surface bias and study the effects of energy losses for both trigger hadrons and recoil jets on I AA , employing the Linear Boltzmann Transport (LBT) model to simulate jet-medium interactions. Our findings suggest that a large portion of hadrons used for the triggers undergoes substantial energy loss, despite surface bias. In particular, the energy loss of the trigger hadrons elevates the I AA baseline, corresponding to the case of no energy loss for recoil jets, to be greatly larger than unity. This enhancement of the baseline implies that the measured I AA values being larger than unity could still signal jet quenching.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Electron temperature gradient instability and transport analysis in NSTX and NSTX-U plasmas

Extensive linear and nonlinear simulations to study electron temperature gradient (ETG) stability and thermal transport in National Spherical Torus Experiment (NSTX) and NSTX-U plasmas were performed using the fully electromagnetic gyrokinetic code CGYRO. Linear simulations were performed to determine ETG thresholds in different discharges, showing that ETG modes in spherical tokamaks can present different scalings compared to conventional aspect-ratio tokamaks. Nonlinear gyrokinetic simulations were conducted for selected cases to calculate electron thermal transport and compare to experimental values. Results are also compared with those of ETG modes in the multi-mode model and the Trapped-Gyro-Landau-Fluid reduced model codes, to better understand their applicability in spherical tokamaks.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Hidden-strangeness tetraquarks in the dynamical diquark model

The dynamical diquark model describes multiquark exotic hadrons in terms of diquark components nucleated by heavy quarks and successfully explains multiple features of hidden-charm and hidden-bottom exotics. Here, we apply the model to the marginally heavy case of hidden-strange states to probe whether mesons near 2 GeV with peculiar properties, such as 𝜙⁡(2170), 𝑓 2 ⁡(2340), and 𝑋⁡(2370), are possible tetraquark candidates. We calculate spin-multiplet average masses using potentials obtained through lattice simulations and quark models, and we also describe the detailed spectra of the expected multiplets as a diagnostic to discern the nature of future hadrons likely to be discovered in this mass region by experiments at facilities such as BESIII, JLab, and the EIC.

Bound states↗
DOE OSTI2025

The relative importance of building design parameters in reducing energy use and sensible heat release from buildings in light of forecasted future weather data and building coverage ratio

Buildings typically have a 60-to-75-year lifespan before they require significant maintenance or modifications. However, most builders evaluate the performance of their new buildings using whole-building energy simulation tools based on the current typical meteorological year (TMY) file or actual meteorological year. The energy use consumption and sensible heat release pattern observed from buildings could potentially change based on shifting global climates. Therefore, the recommended energy-efficiency design parameters might also change during these periods. In this study, we evaluate the role of different building design parameters, such as material reflectivity, HVAC COP, and insulation values, on building energy usage and sensible heat release from buildings with different building coverage ratios (BCR), based on the current and future weather file TMY (fTMY) for the middle of the century (2040–2060). The role of sensible heat release from buildings is not accounted for accurately while estimating building energy usage in most whole-building energy simulations. The study conducts a series of whole-building energy simulation analyses using EnergyPlus to evaluate the role of different design parameters based on TMY and fTMY weather conditions. The analysis is conducted for two hot desert climatic cities: Phoenix (USA) and Abu Dhabi (UAE). The results show that, for the base case in a future climate, the sensible heat release is reduced by an average of 30% due to the reduced delta T between the surface and ambient air. Further, the results show an increase in total energy consumption by 5% annually. The results also show that, for buildings with traditional coatings, shorter buildings release more heat than taller buildings. On the other hand, for buildings with reflective paints, shorter buildings release less heat than taller buildings. The findings from this study can be used by policymakers, utility companies, and builders to better understand the relative role of different building design parameters while constructing new and retrofitting existing buildings.

Alhazmi, Mansour [King Fahd University of Petroleu↗
DOE OSTI2026

PPPL Report on Reduced Modeling of Fusion Alpha Transport in ARC Burning Plasmas

We are reporting on the modeling of fusion alpha particle transport in the planned ARC fusion device being designed by the CFS (Commonwealth Fusion Systems: https://cfs.energy). The ARC tokamak is designed to operate in a burning-plasma regime characterized by a substantial population of fusion-born alpha particles. Alfvén eigenmode (AE) stability is assessed both analytically and numerically, incorporating alpha-particle drive, ion Landau, and radiative damping from thermal species and collisional damping from trapped electrons. Regions of unstable and near-threshold AE activity are mapped across ARC’s operational parameter space. Linear stability analysis with NOVA indicates multiple, often marginally unstable AEs, extending to toroidal mode numbers up to n= 30. The present report focuses on the ARC flat-top operating point prior to the sawtooth event. Alpha-particle transport on timescales exceeding the neoclassical slowing-down time is assessed using the NUBEAM module [1][2] of the TRANSP code [3], employing transport coefficients derived from the RBQ quasilinear modeling (cf. Appendix B). These global simulations identify favorable and unfavorable operating regimes with respect to alpha confinement, pressure redistribution, and overall alpha-heating efficiency. We also evaluate additional transport mechanisms—including neoclassical tearing mode (TM)–induced stochasticity, sawtooth-driven redistribution, and toroidal-field ripple using the kick model (cf. Appendix C) which makes use of the guiding-center code ORBIT, see Section 5. The kick model is integrated into TRANSP to enable self-consistent predictions of alpha-driven current formation and sustainment within the ARC scenario. Sensitivity scans are performed over the mode frequency, rational-surface alignment, island width, mode amplitude, and proximity of the limiter to the plasma. Our study provides an initial, physics-based guidance for machine design, operational planning, and equilibrium control, ensuring adequate alpha confinement and robust self-heating performance in ARC. Our simulations mostly targeted worst case scenarios, e.g. for TMs and sawteeth. Overall, we expect benign effects for the ARC scenario investigated in this work on fusion alpha confinement and losses in the presence of AEs, tearing modes and sawteeth. This report addresses three thrusts identified at the outset. The first thrust focuses on analytic estimates of the parametric dependencies of EP relaxation based on local AE stability simulations (Section 3). The second thrust involves global evaluations of AE stability using the NOVA, RBQ, and NUBEAM codes (Section 4). Finally, we investigate alpha-particle transport driven by low-frequency instabilities associated with sawteeth and tearing modes (Section 5).

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

A Data-driven approach to Core Power distribution reconstruction in a Nuclear Reactor

This report presents the initial development of a data-driven approach for reconstructing the core power distribution in a nuclear reactor (power shape synthesis) using ex-core sensors. Traditional techniques rely on deploying a large number of detectors throughout the reactor core. However, this approach is not feasible for innovative reactor concepts like Advanced Reactors and Microreactors. First, the tight lattice pitch, designed to maximize power density, limits the space available for sensors. Secondly, the harsh operating conditions are not compatible with commercially available detectors. The method proposed in this work integrates high-fidelity modeling with data-driven techniques to accurately reconstruct power distribution across various reactor types, thereby reducing the reliance on in-core sensors. Purdue University Reactor One (PUR-1) was selected as the test case. The CAD model representing the latest configuration of the PUR-1 core was imported into the OpenMC simulation framework, and the model was built. Additionally, the previously developed MCNP6 model was updated. The two models were assessed against the data collected during an experimental campaign conducted in July 2024. Thirty gold foils were placed in three Irradiation Assemblies in PUR-1 core. Using the measured activity of the irradiated foils, the neutron flux at different core locations was estimated.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2024

DEPLOYING FAST CHARGING INFRASTRUCTURE FOR ELECTRIC VEHICLES IN URBAN NETWORKS: AN ACTIVITY-BASED APPROACH

This paper explores an important problem under the domain of network modeling, the optimal configuration of charging infrastructure for electric vehicles (EVs) in urban networks considering EV users' daily activities and charging behavior. This study proposes a charging behavior simulation model considering different initial state of charge (SOC), travel distance, availability of home chargers, and the daily schedule of trips for each traveler. The proposed charging behavior simulation model examines the complete chain of trips for EV users as well as the interdependency of trips traveled by each driver. The problem of finding the optimum charging configuration is then formulated as a mixed-integer nonlinear programming problem that considers the dynamics of travel time and travel distance, the interdependency of trips made by each driver, limited range of EVs, remaining battery capacity for recharging, waiting time in queue, and detour to access a charging station. This problem is solved using a metaheuristic approach for a large-scale case network. A series of examples are presented to demonstrate the model efficacy and explore the impact of energy consumption on the final SOC and the optimum charging infrastructure.

Chain of Trips↗
DOE OSTI2024

Risk-Aware Framework Development for Disruption Prediction: Alcator C-Mod and DIII-D Survival Analysis

Abstract Survival regression models can achieve longer warning times at similar receiver operating characteristic performance than previously investigated models. Survival regression models are also shown to predict the time until a disruption will occur with lower error than other predictors. Time-to-event predictions from time-series data can be obtained with a survival analysis statistical framework, and there have been many tools developed for this task which we aim to apply to disruption prediction. Using the open-source Auton-Survival package we have implemented disruption predictors with the survival regression models Cox Proportional Hazards, Deep Cox Proportional Hazards, and Deep Survival Machines. To compare with previous work, we also include predictors using a Random Forest binary classifier, and a conditional Kaplan-Meier formalism. We benchmarked the performance of these five predictors using experimental data from the Alcator C-Mod and DIII-D tokamaks by simulating alarms on each individual shot. We find that developing machine-relevant metrics to evaluate models is an important area for future work. While this study finds cases where disruptive conditions are not predicted, there are instances where the desired outcome is produced. Giving the plasma control system the expected time-to-disruption will allow it to determine the optimal actuator response in real time to minimize risk of damage to the device.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Code Validation of SAM Using Forced and Natural Circulation Data from NACIE-UP Benchmark

Heavy liquid metals (HLMs) are promising candidates as coolants of Generation IV fast reactors due to their thermophysical properties. In the last decade, experimental work has been proposed as part of research and development efforts to develop such systems. In this context, researchers from the Brasimone Research Center have conducted many experiments using the Natural Circulation Experiment Upgrade (NACIE-UP) facility to study the thermofluid dynamic behavior of HLMs in rod bundle configurations with or without wire wrappers. This facility consists of a rectangular loop operated with lead-bismuth eutectic. Sensors across the loop monitor relevant parameters, i.e. temperatures, heat transfer, and flow conditions. Here, in the present work, we carefully select published data from NACIE-UP to validate the System Analysis Module (SAM), a modern system analysis code developed at Argonne National Laboratory. We developed one SAM model using specifications of the facility geometry and materials existing in relevant papers and reports. On top of that, these references provided the boundary conditions for simulating natural circulation and forced convection experiments in either steady or transient conditions. The SAM model simulates five test cases with diverse operating conditions. Ultimately, the code is proven to predict temperatures and mass flow rates that closely match the experiments. The discrepancies between numerical predictions and diverse transients are limited to a few degrees Celsius, showcasing that SAM is well suited for analyzing nuclear systems relying on HLM coolants.

advanced reactors↗
DOE OSTI2025

Using the optimal combined index weight ratio to improve the probability of anomaly detection in big area additive manufacturing

Big Area Additive Manufacturing (BAAM) of composites requires significant time, energy, and material, so it is critical to reduce production inefficiencies to make functional parts without multiple iterations. Statistical process control coupled with Principal Component Analysis (PCA) is a powerful technique that provides a quick, computationally inexpensive, and intuitive way for operators to detect defects that form in a manufacturing process without massive datasets. Recently, a combined index that is a weighted sum of the Hotelling's T 2 and squared residual error statistics has been proposed that can be monitored in one chart, improving interpretation accuracy and simplicity. However, the literature does not offer a formal method to optimise the weights. Here, we introduce two new approaches to the traditional weight selection approach using simulated and BAAM image data. Approach 1 uses a theoretically motivated optimum inspired by probabilistic principal component analysis. Approach 2 systematically varies the ratio of the weights to find the optimum. We show that approach 1 delivers optimal anomaly detection performance in select cases while approach 2 fares better in practice. Surprisingly, we also show that choosing a more complex PCA model has a minimal negative impact on anomaly detection performance compared to a more simplistic model.

3-dimensional printing↗
DOE OSTI2025

An Investigation of Thermal Properties of 2D Materials [Dissertation]

Studying the thermal conductivity of 2D materials is important due to the applications of 2D materials in fields such as thermal management, thermoelectricity, renewable energy, and sensors. As such, measurements of the thermal conductivity of these 2D materials become important to measure. Thermal conductivity is often difficult to measure for 2D materials due to their atomically thin nature and many experimental methods for doing so requiring contact with the sample, which can alter the thermal properties. A non-contact method for calculating the thermal conductivity of 2D materials supported on substrates in order to model the thermal conductivity of 2D materials for devices, is proposed and experimentally performed in this dissertation. The optothermal Raman technique is a useful non-contact diagnostic technique useful in determining the thermal conductivity of 2D materials. The optothermal Raman typically does not account for heat losses due to convection or radiation or substrate resistance, which are shown to be important factors to consider when developing an optothermal Raman model. Additionally, the calculation of the interfacial thermal conductance between the bottom surface of the sample and the top surface of the substrate, plays an important role in determining the final value of the thermal conductivity of a supported sample, and will yield differing results based on whether or not the conductance is calculated using an approach such as the Diffuse Mismatch Model (DMM) or calculated directly by varying the laser heating profile (usually done by changing the laser objective). This is shown to be the case for both graphene on Ni, graphene on Cu, and SnSe 2 on Cu. In addition to experimentally calculating the thermal conductivity of a 2D material with the optothermal Raman technique, the thermal conductivity of 2D materials can also be calculated using computational methods. The three-phonon method is a method which can be used to simulate phonon scattering processes and determine the thermal conductivity of semiconductors, wherein phonon scattering is the dominant mechanism which determines the thermal conductivity. The three-phonon method uses relaxation times for phonon scattering with other phonons, electrons, and other material system elements, such as isotopes or material defects, in order to create a single-mode relaxation time approximation (SMRTA), which is used to calculate the final value of the thermal conductivity. An important consideration when determining the thermal conductivity of a 2D material using this method is the device geometry, which is reflected in this work as the phonon-boundary scattering relaxation time. This inclusion is important along with the inclusion of phonon-electron scattering in accurately determining the thermal conductivity of a 2D material. In both the optothermal Raman experiments and the three-phonon method computations, strain is shown to have a demonstrable effect on the thermal conductivity of 2D materials. When a 1.1% strain was applied to the mechanical properties of SnSe, the three-phonon processes yielded a lower thermal conductivity than the no-strain case. For the optothermal Raman experiments, the strain induced in the Cu substrate and transferred to a single-layer graphene (SLG) sample yields a trend where the thermal conductivity of the SLG decreases with respect to strain applied. In the case where the interfacial thermal conductance was calculated directly, the conductance increased with respect to strain applied. This presents strain as a reliable and viable method for tuning the thermal properties of 2D materials for device applications.

36 MATERIALS SCIENCE↗
DOE OSTI2024

A new approach to the evaluation and solution of the relativistic kinetic dispersion relation and verification with continuum kinetic simulation

Here, the present work describes a new approach to evaluation and root finding for the kinetic dispersion relation of Langmuir waves, which is central to the analytical understanding of collisionless damping in plasmas. The plasma dispersion function is solved to machine precision using direct integration in the complex plane in combination with an analytic evaluation of the residue to account for the deformation along the Landau contour. To efficiently attain machine precision, the contour is displaced in the complex plane prior to integration, and numerical subtleties related to the placement of the contour are discussed. The approach is generic in that it applies to arbitrary distribution functions, with the present manuscript focused on relativistic cases. Detailed verification of results via direct kinetic simulation in a variety of configuration space dimensions is also presented. Finally, the technique is applied to the challenging case of highly relativistic (i.e. extremely hot) plasmas. Here we show both qualitative agreement with prior work, as well as the disappearance of the Landau root which would have significant implication for real-life observation or experiment.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

NSTXU Diagnostic Disruption Dynamic Loading Represented by Response Spectra

This article presents the results of transient dynamic simulations of loads due to disruption eddy currents on the NSTXU vacuum vessel. Dynamic loading at diagnostic mounting locations is expressed as response spectra derived from the time history results of the dynamic structural simulations of a variety of disruption scenarios. The disruption simulations draw on a history of the project assessments of worst case disruptions for specific components. Major efforts to assess disruption loading have included the vacuum vessel which is the major structural support for the machine, as well as the passive plates (PPs), high harmonic fast wave (HHFW) antenna, and centerstack casing. Each one of these efforts included transient electromagnetic simulations producing time-dependent eddy current Lorentz loads (and in some cases halo loads) which then were applied to time-dependent structural dynamic analyses intended to obtain the proper dynamic amplification factors. In some instances, the EM model and structural model were identical allowing direct transfer of EM forces to the structural model. In other cases, the EM and structural model were not identical and the vector potential (VP) transfer method was used. The results files from these analyses were available (or re-run) to post process in ANSYS Classic time history postprocessor. In conclusion, the ANSYS command is used to create response spectra from time history data at desired points on the vessel.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗