Search NASA⌕ Search

SEARCH · Search NASA

Results for “Numerical”

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

Estimating Residual Gas Beam Loss in the BNL AGS Booster:Model, Numerical Tools, and Implications for a Polarimeter Chamber

Beamloss from residual gas scattering is calculated for a storage ring of circumference C = 200m, covering a linear magnetic rigidity ramp from B ρ inj = 1.25Tm to B ρ ext = 9.50Tm over tramp = 1.00s. Three loss mechanisms are considered: nuclear inelastic reactions, single Coulomb scattering beyond the acceptance angle of the machine, and charge-state exchange for heavy ions. Calculations are performed for protons and Au 32+ ions, the latter being the charge state delivered by the EBIS pre-injector into the BNL AGS Booster. At a baseline vacuum of P 0 = 5 × 10 −11 mbar the integrated fractional beam loss over the full ramp is negligible for protons (< 0.001%) and 33.1% for Au 32+ . The impact of a 1m polarimeter section at elevated pressure P s = 5 × 10 −9 mbar is evaluated: the proton loss remains negligible, while the Au 32+ loss rises to 45.1%, establishing the vacuum requirement for the polarimeter chamber. In the present compact approximate model, Au 32+ losses are dominated by charge exchange at injection, while proton losses are dominated by single Coulomb scattering. The design, outgassing budget, pumping strategy, and commissioning sequence for a UHV polarimeter chamber housing silicon detectors and ASICs are described in detail.

43 PARTICLE ACCELERATORS↗
DOE OSTI2024

Experimental and Numerical Characterization of High Temperature Deformation Behavior of 347H Stainless Steel

This study investigates how temperature affects the plasticity and thermal creep behavior of 347H stainless steel under uniaxial tension. The research combined experimental testing with advanced computational modeling. Two types of experiments were conducted: uniaxial tensile tests at temperatures from 100°C to 750°C using strain rates of ~10⁻⁴ s⁻¹, and creep tests at temperatures between 600°C and 750°C under various stress levels. These experimental results were used to develop and validate a new integrated mechanistic model that can predict material behavior under any loading condition while accounting for both stress and temperature effects. The model was implemented using a polycrystalline microstructure simulation framework based on elasto-viscoplastic Fast Fourier Transform (EVPFFT). It incorporates three key deformation mechanisms: thermally activated dislocation glide, dislocation climb, and vacancy diffusional creep. The model accounts for internal stress distribution within single crystals and considers how precipitates and solute atoms (both interstitial and substitutional) affect dislocation movement. After validation against experimental data, the model was used to generate Ashby-Weertman deformation mechanism maps for 347H steel, providing new insights into how microstructure influences the activation of different creep mechanisms.

Joy, Jobin K.↗
DOE OSTI2025

Numerical simulations of liquid jetting with solid inclusions

The dynamics of finite-sized particles in fluids, and their influence on the overall flow, are of great interest across several industrial, environmental, and medical fields. In the context of inkjet printing, the presence of solid inclusions can be either intentional, as in additive manufacturing, or unintentional, as in standard printing processes. These inclusions can strongly impact the jetting process, causing effects such as jet asymmetry, bubble entrapment, and the formation of satellite droplets. Understanding and controlling particle behavior is therefore essential, particularly to predict how and when particles are ejected over multiple jetting cycles. It is therefore critical to develop reliable models that allow for a deeper understanding of the complex interplay between particle and fluid during the whole printing process. To address this, we present a tailored implementation of the Color-Gradient multicomponent Lattice Boltzmann Method for fully resolved three-dimensional (3D) simulations of multicycle liquid jetting with particles. Our method supports realistic parameter settings aligned with industrial inkjet systems, and we provide both qualitative and quantitative validation against experimental data. Additionally, we introduce a simplified model based on the Stokes drag law, in which solid particles are represented as point particles and do not influence the fluid flow. Despite this limitation, the model offers a computationally efficient means to explore the vast parameter space typically encountered in industrial applications, allowing, e.g., identifying critical ejection regions and estimating the number of cycles required for particle release. These qualitative insights are valuable for guiding and complement fully two-way coupled simulations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Numerical Investigation of Butterfly Valve Performance in Variable Valve Sizes, Positions and Flow Regimes

Reliability and efficiency of valves are necessary for precise control and sufficient heat-flow to heat application plants for the integrated energy systems of nuclear power plants (NPPs). Strategic Management Analysis Requirement and Technology (SMART) valves’ ability to control flow and assess environmental parameters stands out for these requirements. Their ability to sustain the downstream flow rate, prevent reverse flow, and maintain pressure in the heat transport loop is much more efficient with the integration of sensors and intelligent algorithms. For assessing valve performance and monitoring, mechanical design and operating conditions are two important parameters. In this study, the butterfly valves of three different sizes are simulated with water and steam using STAR-CCM+ in various flow regimes and positions to analyze performance parameters to strategize an automated control system for efficiently balancing the heat–transport network. Also, flow behavior is studied using velocity and pressure fields for valve–body geometry optimization. It can be observed, through performance parameters, that the valves are suitable for operation between 30° and 90° positions with significantly low loss coefficients and high flow coefficients, and the performance parameters follow a certain pattern in both water and steam flow in each scenario.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2025

A Numerical Model for the Dynamics of Pickup Ions Outside the Heliopause and IBEX “Ribbon” Observation

The leading mechanism for the origin of the Interstellar Boundary Explorer (IBEX) “ribbon” of enhanced energetic neutral atoms (ENAs) from the outer heliosphere is the so-called secondary ENA process. In this mechanism, neutralized ions escape the heliosphere, and after several times of charge exchange collision,s some propagate back toward Earth in directions nearly perpendicular to the local interstellar magnetic field. However, the physical processes governing the distribution of the pickup ions (PUIs) outside the heliopause are still unclear. In this study, we build a new global model of the IBEX ribbon, where a key component is to calculate the dynamics of the PUIs outside the heliopause by solving the gyrophase-averaged focused transport equations on top of the background magnetohydrodynamic-kinetic model for the heliosphere-interstellar medium interaction. We discuss how the properties of the simulated ribbon change with different scattering parameters and show simulation results of the ENA sky maps, ribbon centers, ribbon profiles at several azimuthal slices, the ribbon ENA source region, and the velocity distribution of PUIs in the ENA source region for different scattering frequencies. Our results show that the model can provide reasonable ribbon sky maps comparable to the observed ribbon when the PUI mean free path is on the order of several thousand astronomical units.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Numerical Treatment of Shock-induced Nuclear Burning in Double Detonation Type Ia Supernovae

We present a benchmark problem to assess the treatment of shock-induced nuclear burning in the context of double detonation Type Ia supernovae. In a stratified white dwarf model, we implement a shock-detection criterion that suppresses burning in zones characterized by compression and significant pressure gradients, controlled by a tunable parameter, f shock . One-dimensional simulations, using the open-source Castro suite, were conducted across three treatments—burning fully enabled, and burning suppressed with f shock = 2/3 and f shock = 1—across three spatial resolutions (5.0, 2.5, and 0.3125 km). At the finest resolution, the burning-enabled and f shock = 1 models converge, while the f shock = 2/3 front continues to show slight offset behavior. Since most simulations are carried out at much lower resolutions, our tests support the idea that burning in shocks should always be disabled in practice. We also observe that the behavior of lower-resolution simulations remains extremely sensitive to the choice of f shock .

GPU computing↗
DOE OSTI2026

Wake-Resolving Acoustic Tomography: Advances through Numerical Covariance Methods

Acoustic tomography offers path-integrated measurements of atmospheric velocity and temperature fluctuations with high spatial resolution. Classical implementations of time-dependent stochastic inversion rely on homogeneous, isotropic covariance models that are poorly suited to the anisotropic structure of wind turbine wakes. By directly estimating heterogeneous covariances from large-eddy simulations (LESs) into the time-dependent stochastic inversion operator, we relax implicit assumptions in the analytical models used historically. Retrievals using these LES-informed models improve agreement with true fields in variance, turbulent kinetic energy, and spectral content compared to analytical and precursor-based covariance models. The results indicate that LES-informed covariance models can enhance the accuracy of acoustic tomography retrievals in complex, anisotropic flows such as wind turbine wakes in some cases and highlight instances where analytical models still offer competitive performance, despite their simplifying assumptions.

17 WIND ENERGY↗
DOE OSTI2026

Numerical Modeling of a Two-Stage Ocean Current Turbine

The Equinox Ocean Turbines (EQOT) current energy converter has a unique design with power generation in two small-diameter turbines attached to the tips of a large-diameter passive rotor. This configuration offers some key advantages for capturing ocean currents. With no centrally placed generator, almost no reaction torque is required at the nacelle of the main large-diameter rotor, and the small-diameter tip turbine generators operate at a higher speed and lower torque. The physics that determine the performance and loads on the turbine are also unique. The interactions of the flow field between the two stages and the general architecture of the system cannot be captured with traditional mid-fidelity modeling tools. For design iterations and large sets of load cases, it is important to have mid-fidelity models that can capture the important phenomenon with enough accuracy to identify global trends. This work uses a limited set of high-fidelity computational fluid dynamics (CFD) simulations to help inform the selection of and construction of a custom mid-fidelity model. Mid-fidelity modeling approaches were verified by comparing key turbine performance quantities to those found with the CFD model. Hydrodynamic interactions of the two-stage rotor were identified through high-fidelity CFD modeling. This highlighted the impact of the main rotor tip vortex and wake on the secondary rotor apparent inflow. This results in a relative flow rotation and sharp deficit, that change the optimal secondary rotor rotation speed and adds unsteadiness to the blade loading respectively. Multiple mid-fidelity approaches were evaluated for their ability to capture these effects. A simple approximation of the combined-stage performance based on single-stage BEM provides a reasonable rough prediction, especially near the peak TSR values, with some larger discrepancy at higher TSRs. Predicting the combined-stage performance based on single-stage CFD data improves this prediction across the TSR range. Although the combined-stage modeling in OLAF was not successful in this stage of the project, it showed promise as a mid-fidelity method, assuming the parameters can be tuned to account for the significant differences in time and length scales between the main and secondary rotors. This may be addressed through code changes in future work. A significant finding from the OLAF work was the agreement between the vortex core radius values found independently via a parameter space search and via CFD. The technique of using single-stage secondary rotor BEM, with a custom inflow taken from single-stage main rotor CFD or OLAF, provides an efficient method to capture one-way coupled flow interactions. This method provided generally good predictions of the impact of the flow rotation on the secondary rotor but struggled to accurately predict the peaks of the unsteady load progression. Future work could include some superposition of a tuned main rotor trailing edge viscous wake into the custom inflow to better predict this interaction.

16 TIDAL AND WAVE POWER↗