Search NASA⌕ Search

SEARCH · Search NASA

Results for “numerical validation”

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

Anti-Ultralocality and Plateau Models of Inflation

Anti-ultralocality refers to the growth of spatial gradient terms relative to velocity terms in the coupled Einstein--scalar field equations. It is a characteristic feature of decelerated expansion before the onset of inflation. Previous numerical relativity studies have shown that anti-ultralocality prevents the onset of inflation in models with power-law inflaton potentials. In this paper, we show that models with plateau-shaped inflaton potentials, which are considered to be the simplest way to generate a tensor-to-scalar ratio below current observational upper limits, are especially vulnerable to anti-ultralocality effects. The reasons are the flatness of the plateau and the energy density gap of $\sim 10$ orders of magnitude between the Planck density and the plateau potential energy. To study the problem, we develop a protocol for assessing the viability of inflationary models in general, and we apply it to a plateau potential using a previously validated numerical relativity code. We find that, starting from generic initial conditions, the growth of gradient terms in the Einstein equations relative to non-gradient terms either prevents inflation from lasting for enough $e$-folds or triggers a phase of quantum runaway. We show that the fine-tuning of initial conditions necessary to avoid these issues becomes more severe as the energy scale of inflation is made smaller, disfavoring common approaches for reducing the tensor-to-scalar ratio.

FOS: Physical sciences↗
DOE OSTI2025

Modeling the formation of Sedan Crater using the FLAG and HOSS codes

Numerical modeling of explosion crater formation requires accounting for complex physical processes. Numerical validation of explosion cratering is an important step in modeling and requires experimental data for comparison. Models using discrete elements and continuum models have both benefits and drawbacks to their approaches. In this work, we consider both an arbitrary Lagrangian–Eulerian (ALE) hydrocode and a finite discrete element method (FDEM) approach to modeling the formation of the Sedan crater, the largest human-made crater in the United States. The Sedan crater formed from an underground nuclear detonation in the Nevada desert as part of Project Plowshare. Our models show that the continuum approach of the hydrocode matched well compared to early test time prior to the mound rupture and subsequent fireball venting, when most of the alluvium exhibited fluid behavior. Our FDEM approach matched the final crater dimensions well, after material had settled back into the crater, when material strength and solid mechanics play key roles. Our work shows how leveraging the benefits of multiple numerical approaches can lead to better understanding of complex physical problems, especially problems with limited experimental data. By using a continuum approach to early-time hydrodynamics and an FDEM approach to later-time solid mechanics, we can better understand the different physical regimes of explosion crater formation.

36 MATERIALS SCIENCE↗
NASA NTRS2014

An Experimental Investigation of the 30P30N Multi-Element High-Lift Airfoil

High-lift devices often generate an unsteady flow field producing both broadband and tonal noise which radiates from the aircraft. In particular, the leading edge slat is often a dominant contributor to the noise signature. An experimental study of a simplified unswept high-lift configuration, the 30P30N, has been conducted to understand and identify the various flow-induced noise sources around the slat. Closed-wall wind tunnel tests are performed in the Florida State Aeroacoustic Tunnel (FSAT) to characterize the slat cove flow field using a combination of surface and off-body measurements. Mean surface pressures compare well with numerical predictions for the free-air configuration. Consistent with previous measurements and computations for 2D high-lift configurations, the frequency spectra of unsteady surface pressures on the slat surface display several narrowband peaks that decrease in strength as the angle of attack is increased. At positive angles of attack, there are four prominent peaks. The three higher frequency peaks correspond, approximately, to a harmonic sequence related to a feedback resonance involving unstable disturbances in the slat cove shear layer. The Strouhal numbers associated with these three peaks are nearly insensitive to the range of flow speeds (41-58 m/s) and the angles of attack tested (3-8.5 degrees). The first narrow-band peak has an order of magnitude lower frequency than the remaining peaks and displays noticeable sensitivity to the angle of attack. Stereoscopic particle image velocimetry (SPIV) measurements provide supplementary information about the shear layer characteristics and turbulence statistics that may be used for validating numerical simulations.

Pascioni, Kyle A.↗
DOE OSTI2024

(Doublon) Benchmarking of Different Inverse Point Kinetics Implementations for an Autocorrected Reactimeter Algorithm

In November 2017, the Transient Reactor Test Facility returned to operation. Since that time, many transient test series have been completed, such as the Transient Heatsink Overpower Response capsule (THOR), the Transient Water Irradiation System for TREAT (TWIST), and Sirius. Each has provided valuable data for materials performance and reactor safety that can be applied in future designs. During each experimental series, detector count rates provided important information on the core behavior during transients. However, a limitation of these data is that variations in the neutron distribution during experiments can cause errors when attempting to infer reactivity evolution from detector signals. Neutron physics codes can be used to compute the flux shape variations. However, this is a poor solution when the experimental data is used for code verification, validation and uncertainty quantification. Indeed, if the output of the code is used both as a reference and to correct what the reference is compared to, the circular dependency limits the quality of the verification, validation and uncertainty quantification approach. To overcome this problem, the autocorrected reactimeter algorithm (ACRA) has been developed. This approach infers a time-dependent reactivity evolution by testing different spatial corrections and selecting the one that minimizes reactivity variations when the core is in a frozen configuration (i.e., when there is no variation in parameters affecting reactivity). However, the scope of this method was limited to transients where there were negligible thermal feedback. Indeed, the core is never in a frozen configuration when the fuel temperature varies during the whole transient. This is our motivation for developing an improved version of the ACRA that does not require frozen configurations. To develop this new algorithm, we need a precise and unbiased implementation of the inverse point kinetic equations (IPKEs) as any error in the reactivity evaluation will be propagated into the choice of the optimal spatial correction. Indeed, the previous reactimeter algorithm would use approximations, such as a negligible flux amplitude derivative, to focus on rapidity. For the numerical validation of ACRA, we aim at absolute error under for reactivity derived from signals similar to the one of this study. In this summary, we test eight different IPKE implementations. Each will process a mockup signal built for this study, similar to those that the future ACRA will process. Each reactivity output will be compared to the reference reactivity that has been used to generate the mockup signal. The implementation minimizing the difference with the reference reactivity will be used in the development of a new ACRA formulation.

73 - NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS1994

Interface Shape and Convection During Solidification and Melting of Succinonitrile

An experimental study was conducted of the crystal growth of succinonitrile during solidification, melting, and no-growth conditions using a horizontal Bridgman furnace and square glass ampoule. For use as input boundary conditions to numerical codes, thermal profiles on the outside of the ampoule at five locations around its periphery were measured along the ampoule's length. Temperatures inside the ampoule were also measured. The shapes of the s/l interface in various two dimensional planes were quantitatively determined. Though interfaces were nondendritic and noncellular, they were not flat, but were highly curved and symmetric in only one unique longitudinal y-z plane (at x=O). The shapes of the interface were dominated by the primary longitudinal flow cell characteristic of shallow cavity flow in horizontal Bridgman; this flow cell was driven by the imposed furnace temperature gradient and caused a 'radical' thermal gradient such that the upper half of the ampoule was hotter than the bottom half. We believe that due to the strong convection, the release of latent heat does not significantly influence the thermal conditions near the interface. We hope that the interface shape and thermal data presented in this paper can be used to optimize crystal growth processes and validate numerical models.

Degroh, Henry C., III↗
DOE OSTI2020

NUMERICAL INVESTIGATION OF FUEL PROPERTY EFFECTS ON MIXED-MODE COMBUSTION IN A SPARK-IGNITION ENGINE

In the present study, mixed-mode combustion of an E30 fuel in a direct-injection spark-ignition engine is numerically investigated at a fuel-lean operating condition using multidimensional computational fluid dynamics (CFD). A fuel surrogate matching Research Octane Number (RON) and Motor Octane Number (MON) of E30 is first developed using neural network based non-linear regression model. To enable efficient 3D engine simulations, a 164-species skeletal reaction mechanism incorporating NO, chemistry is reducedfrom a detailed chemical kinetic model. A hybrid approach that incorporates the G-equation model for tracking turbulentflame front, and the multi-zone well-stirred reactor model for predicting auto-ignition in the end gas, is employed to account for turbulent combustion interactions in the engine cylinder. Predicted in-cylinder pressure and heat release rate traces agree well with experimental measurements. The proposed modelling approach also captures moderated cyclic variability. Two different types of combustion cycles, corresponding to purely deflagrative and mixed-mode combustion, are observed. In contrast to the purely deflagrative cycles, mixed-mode combustion cycles feature early flame propagation followed by end-gas auto-ignition, leading to two distinctive peaks in heat release rate traces. The positive correlation between mixed-mode combustion cycles and early flame propagation is well captured by simulations. With the validated numerical setup, effects of NOx chemistry on mixed-mode combustion predictions are investigated. NOx chemistry is found to promote auto-ignition through residual gas recirculation, while the deflagrative flame propagation phase remains largely unaffected. Local sensitivity analysis is then performed to understand effects of physical and chemical properties of the fuel, i.e., heat of evaporation (HoV) and laminar flame speed (S-L). An increased HoV tends to suppress end-gas auto-ignition due to increased vaporization cooling, while the impact of HoV on flame propagation is insignificant. In contrast, an increased SL is found to significantly promote both flame propagation and auto-ignition. The promoting effect of SL on auto-ignition is not a direct chemical effect; it is rather caused by an advancement of the combustion phasing, which increases compression heating of the end gas.

detailed chemical kinetics↗
DOE OSTI2024

Bench Testing Data and Report for an Early Prototype Pitch Resonator WEC

This dataset encompasses data and documentation from bench tests conducted on an early prototype of a "pitch resonator" wave energy converter (WEC). The testing aimed to validate numerical models and reduce risks associated with the pitch resonator concept, which is designed to convert the pitching and rolling motions of a buoy into electrical power. The project's goal is to provide supplementary power, in the range of 10-100 watts, to the National Science Foundation's Ocean Observatories Initiative Pioneer Array. Two distinct testing phases are documented: one using a single degree of freedom (1DOF) test rig, and another employing a six degree of freedom (6DOF) Stewart platform, known as the Large Amplitude Motion Platform (LAMP). These tests assessed various factors, such as system performance in different motion scenarios, the torque exerted by wave forces, and the impact of mounting configurations. The dataset includes raw test data in MATLAB (.mat) format, detailed metadata, and a report describing the experimental procedures and preliminary findings.

16 TIDAL AND WAVE POWER↗
DOE OSTI2022

Measuring binary fluidization of nonspherical and spherical particles using machine learning aided image processing

Abstract The binary fluidization of Geldart D type nonspherical wood particles and spherical low density polyethylene (LDPE) particles was investigated in a laboratory‐scale bed. The experiment was performed for varying static bed height, wood particles count, as well as superficial gas velocity. The LDPE velocity field were quantified using particle image velocimetry (PIV). The wood particles orientation and velocity are measured using particle tracking velocimetry (PTV). A machine learning pixel‐wise classification model was trained and applied to acquire wood and LDPE particle masks for PIV and PTV processing, respectively. The results show significant differences in the fluidization behavior between LDPE only case and binary fluidization case. The effects of wood particles on the slugging frequency, mean, and variation of bed height, and characteristics of the particle velocities/orientations were quantified and compared. This comprehensive experimental dataset serves as a benchmark for validating numerical models.

Li, Cheng↗
DOE OSTI2023

Characterizing the seismic response of a molten salt nuclear reactor

Abstract One Generation IV nuclear reactor, which uses a fluoride salt as a coolant, graphite reflector blocks as a moderator, and circulating buoyant TRISO pebbles as fuel is at an advanced stage of development. To characterize the seismic behavior of components of this reactor, validate numerical models for analysis, and develop recommendations for design, a set of earthquake‐simulator experiments on a scaled model of the reactor vessel and its internals was executed on a six‐degree‐of‐freedom earthquake simulator. The model was seismically isolated at its base using two types of spherical sliding bearings. The scaled model involved representations of the prototype reactor vessel, core barrel, reflector blocks, coolant, and spherical fuel pebbles. The material and geometric properties of different test components were selected based on a dynamic similitude scaling analysis and an approximate length scale of 0.4. Four sets of three‐component earthquake motions were used as inputs for testing. Instrumentation on the test specimen recorded the dynamic responses of the outer vessel, core barrel, and reflector‐block assembly, the hydrodynamic responses (sloshing and hydrodynamic pressure) of the liquid coolant, pebble consolidation under earthquake shaking, and the behavior of the isolation systems. This paper describes the design of the experiments and presents key results from the tests. The dynamic responses of the outer vessel, core barrel, and the reflector blocks revealed that the components responded as a unit for the intense shaking used in the experiments. The sloshing response of the fluid in a thin annulus near the perimeter of the vessel was heavily damped. The change in the packing fraction of the pebble bed under repeated, intense 3D earthquake shaking was less than 3%. Seismically isolating the vessel substantially reduced demands on its internal components.

Engineering↗
DOE OSTI2021

A duality-based coupling of Cosserat crystal plasticity and phase field theories for modeling grain refinement

High-rate deformation processes of metals entail intense grain refinement and special attention needs to be paid to capture the evolution of microstructure. In this article, a new formulation for coupling Cosserat crystal plasticity and phase field is developed. A common approach is to penalize kinematic incompatibility between lattice orientation and displacement-based elastic rotation. However, this can lead to significant solution sensitivity to the penalty parameter, resulting in low accuracy and convergence rates. To address these issues, a duality-based formulation is developed which directly imposes the rotational kinematic compatibility. A weak inf-sup-based skew-symmetric stress projection is introduced to suppress instabilities present in the dual formulation. An additional least squares stabilization is introduced to suppress the spurious lattice rotation with a suitable parameter range derived analytically and validated numerically. The required high-order continuity is attained by the reproducing kernel approximation. It is observed that equal order displacement-rotation-phase field approximations are stable, which allows efficient employment of the same set of shape functions for all independent variables. The proposed formulation is shown to yield superior accuracy and convergence with marginal parameter sensitivity compared to the penalty-based approach and successfully captures the dominant rotational recrystallization mechanism including block dislocation structures and grain boundary migration.

42 ENGINEERING↗
DOE OSTI2026

Quantum dynamics of cosmological particle production: interacting quantum field theories with matrix product states

Understanding real-time dynamics of interacting quantum fields in curved space-time remains a major theoretical challenge. We employ tensor network methods to study such dynamics using interacting scalar and gauge theories in 1+1 spacetime dimensions, subject to a quench modeling a homogeneously expanding gravitational background. The models considered are the scalar λϕ 4 theory and the Schwinger model, i.e. a Dirac fermion coupled to a U(1) gauge field which is equivalent via bosonization to a scalar field with a cosine self-interaction. In the free scalar limit, both theories reproduce known analytical results, providing a nontrivial numerical validation of bosonization in curved spacetime for the Schwinger model. Our central finding is that self-interactions lead to a suppression of gravitational particle production compared to the free-field case, as evidenced by two-point functions and the spectra of produced particles. We further examine the behavior of entanglement generation and find that interactions suppress entanglement growth in the λϕ 4 theory, while in the Schwinger model, the interplay between suppressed particle production and enhanced inter-particle correlations leads to more complex entanglement behavior. Our results pave the way for further explorations of nonperturbative quantum real-time dynamics of interacting scalar and gauge theories in arbitrary gravitational backgrounds.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Origin of nontopological soliton dark matter: solitosynthesis or phase transition

This work demonstrates that nontopological solitons with large global charges and masses, even above the Planck scale, can form in the early universe and dominate the dark matter abundance. In solitosynthesis, solitons prefer to grow as large as possible under equilibrium dynamics when an initial global charge asymmetry is present. Their abundance is set by when soliton formation via particle fusion freezes out, and their charges are set by the time it takes to accumulate free particles. This work improves the estimation of both quantities, and in particular shows that much larger-charged solitons form than previously thought. The results are estimated analytically and validated numerically by solving the coupled Boltzmann equations. Without solitosynthesis, phase transitions can still form solitons from particles left inside false-vacuum pockets and determine their present-day abundance and properties. Even with zero charge asymmetry, solitons formed in this way can have very large charges on account of statistical fluctuations in the numbers of (anti)particles inside each pocket.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

A novel framework for elastoplastic behaviour of anisotropic solids

A general framework for developing nonlinear hyperelastic/plastic constitutive laws for anisotropic solids experiencing large strains and strain rates has been developed. The proposed framework does not rely on the “a priori” known strain energy function, but instead introduces a physical decomposition of the material element into seven physically independent stress bearing mechanisms, each of which has a constitutive law in terms of internal moments described by a scalar function of a single variable. Furthermore, the model has been encoded into a combined finite-discrete element method and tested against static geomechanical test data. The numerical validation experiments show the model can reproduce plastic anisotropic behavior in both biaxial and uniaxial loading of a geomaterial.

58 GEOSCIENCES↗
DOE OSTI2021

Design and performance evaluation of a dual-circuit thermal energy storage module for air conditioners

We present experimental results and a validated numerical model of a dual-circuit phase-change thermal energy storage module for air conditioners. The module incorporates a phase-change material encapsulated in compressed expanded natural graphite foam. We used n-tetradecane as the PCM with a transition temperature (~4.5 °C) suitable for air-conditioning applications. Heat exchange to and from the module is accomplished through two fluid loops operating as a heat source and sink embedded inside multiple slabs of the composite material. This dual-circuit design enables easier integration with air-conditioning equipment and provides enhanced flexibility in system operation as compared to the state-of-the-art thermal storage systems. When integrated with an air-conditioner, this design will enable peak-load shaving and enhances operational efficiency. The thermal storage device was designed for a nominal storage capacity of ~ 3.5 kWh. We evaluated the heat transfer and energy storage performance of this device using standalone heat transfer experiments to estimate key thermal resistances and identify design improvements before integration with an air conditioner. The numerical model of the heat exchanger uses a combination of discretized and lumped parameter approaches to maintain a balance between accuracy and computational expense. Our analyses show that the geometric features and integration of fluid tubes are key contributors to the thermal contact resistance between the fluid and the thermal storage material, and consequently, to the overall performance of the thermal storage module. Our standalone experiments also identified important operating scenarios in which this thermal storage module can be used for air-conditioning in buildings.

25 ENERGY STORAGE↗
DOE OSTI2022

Rethinking the ill-posedness of the spectral function reconstruction — Why is it fundamentally hard and how Artificial Neural Networks can help

Reconstructing hadron spectral functions through Euclidean correlation functions are of the important missions in lattice QCD calculations. However, in a Källen–Lehmann (KL) spectral representation, the reconstruction is observed to be ill-posed in practice. It is usually ascribed to the fewer observation points compared to the number of points in the spectral function. In this paper, by solving the eigenvalue problem of continuous KL convolution, we show analytically that the ill-posedness of the inversion is fundamental and it exists even for continuous correlation functions. We discussed how to introduce regulators to alleviate the predicament, in which include the Artificial Neural Networks (ANNs) representations recently proposed by the Authors in another study. The uniqueness of solutions using ANNs representations is manifested analytically and validated numerically. Reconstructed spectral functions using different regularization schemes are also demonstrated, together with their eigen-mode decomposition. We observe that components with large eigenvalues can be reliably reconstructed by all methods, whereas those with low eigenvalues need to be constrained by regulators.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Performance evaluation of an air conditioning system based on quasi isothermal dehumidifcation

This work presents a first study of novel air conditioning system based on quasi isothermal dehumidification combined with dew-point evaporative cooling through the Maisotsenko cycle (M-Cycle). The dehumidification process is realized in desiccant coated heat exchanger with heat rejection system. The key idea behind the solution is to maximize the cooling effectiveness by creating favourable conditions for dehumidification and evaporative cooling process. This is achieved by rejecting the heat of sorption in the desiccant system and providing cold and dry air for M-Cycle air cooler, which allows the system to be powered by using low-grade heat energy at 50 °C. The analysis is performed on the basis of experimentally validated numerical models of each component combined in a unified model of the system. The goal of the study is to determine the performance of the solution under different operational conditions. In addition two types of M-Cycle air coolers (i.e. cross-flow and regenerative unit) are compared to establish which unit is more suitable as the source of cooling energy for the system. It was established that proposed system is able to provide comfortable room parameters (supply air temperature equal 17 °C and humidity ratio equal 0.010 kg/kg-da) for outdoor air temperatures ranging from 28 to 36 °C and outdoor air humidity ranging from 0.010 to 0.018 kg/kg-da. This means that the proposed solution is able to provide comfortable conditions in temperate and humid climates. Results also indicate that the system consisting of the regenerative M-Cycle HMX allows better cooling load management, offers higher thermal Coefficient of Performance (COP th ) and uses less water compared to the system with the cross-flow HMX. However, the system with cross-flow heat and mass exchanger (HMX) has a relatively greater electrical Coefficient of Performance (COP el ) in most cases.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2024

Experimental evaluation of thermal-fluids performance of Helium-Cooled Flat Plate (HCFP) divertor

Various helium (He)-cooled solid-tungsten (W) divertor concepts have been proposed for long-pulse magnetic fusion energy reactors. Among these concepts, the He-cooled flat plate divertor (HCFP) modules have the largest plasma-facing surface area of ~0.2 m 2 . Simulations have shown that the most recent version of the design can withstand heat fluxes as great as 8 MW/m 2 . Earlier experimental studies of a single shortened HCFP cooling unit with a slot length of 7.6 cm used air at ambient temperature and pressures below 0.6 MPa. Here, we present initial experimental studies of a copper alloy and steel test section modeling a shortened HCFP cooling unit using He at prototypical pressure of 10 MPa, inlet temperatures T i ≤ 200 °C and steady-state incident heat fluxes q" ≤ 1.2 MW/m 2 . Results for Nusselt number Nu as a function of Reynolds number Re were obtained for Re = 1.2×10 4 –3.4×10 4 , and used to develop a Nu(Re) correlation and validate numerical models of the test section using commercial computational fluid dynamics (CFD) software. Simulations with this model are performed to evaluate the effect of the shortened slot. Furthermore, these analyses are used to estimate the thermal-fluids performance of the HCFP under prototypical conditions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Experimental measurements of fluid flow in an 84-pin hexagonal rod bundle with spacer grid for a gas-cooled fast modular reactor

A 50 MW e helium-cooled fast modular reactor (FMR) is under development, and as part of the Department of Energy Integrated Research Project (IRP), Texas A&M University is conducting the thermal-hydraulic characterization of its baseline core configuration. Here, we experimentally investigated the axial and cross flow fields characteristics in the hexagonal fuel rod bundle composed of 84 rods, a central rod, and spacer grids at a Reynolds number of 12,000. A fully transparent experimental facility resembling of one unit of the fuel assembly was constructed. Time-resolved particle image velocimetry (TR-PIV) measurements were performed to characterize hydraulic behavior downstream the spacer grid. Velocity measurements were conducted in the axial and radial direction of the rod bundle. From the PIV velocity vector fields, the full-field flow statistics were computed for the mean velocity, vorticity, and Reynolds stresses. An analysis of the energy decay downstream of the spacer grid was conducted with calculations of secondary-flow intensities, turbulent kinetic energy, power spectrum, and spatial-temporal velocity cross correlations. The vorticity field was obtained and pairs of counter-rotating vortexes were identified using a 3D reconstruction of several measurement planes. The data from this experimental campaign will be used to validate numerical models based on Computational Fluid Dynamics and other codes.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗