SEARCH · Engineering Papers

Results for “Closure model”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 901 records · Page 50

All sources NASA NTRS DOE OSTI

DOE OSTI2025

Yrast structures of 98 Nb and 99 Mo

Neutron-rich nuclei in the A ≈ 100 mass region are interesting due to a rapid shape transition, especially pronounced in the Zr isotopes, and more recently observed in Nb isotopes. 98 Nb, with only one proton and one neutron outside the subshell closure nucleus of 96 Zr, is amenable to shell model calculations. To further examine the rapid shape transition the yrast structure of 98 Nb was established in this work. This was the only yrast structure missing from all immediate neighbors to 96 Zr. The yrast structure of 98 Nb was studied in the fission of the compound systems formed in three heavy-ion induced reactions, 24 Mg (134.5 MeV) + 173 Yb, 23 Na (129 MeV) + 176 Yb, and 18 O (91 MeV) + 208 Pb. Prompt γ-ray spectroscopy was accomplished using the Gammasphere array. Excitation energies up to 3 MeV were observed for the first time in 98 Nb. The yrast structure above the previously known (5) + isomer was established. In the process of studying 98 Nb the yrast structure of positive-parity states in 99 Mo was extended to 3.7 MeV excitation energy, the previously-known 99 Nb level scheme was enriched, and two new levels were added in the level scheme of 97 Zr. The coupling of the odd proton occupying the g 9/2 orbital to the yrast states in the core nucleus of 97 Zr can account for all observed states in 98 Nb. The yrast structure for the positive-parity states of 99 Mo is compared to the deformed ground-state bands of the 101 Ru isotone and of 98,100 Mo.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2020

2020 LANL Contributions to the BATS Test in WIPP

The DOE Office of Nuclear Energy (DOE-NE) repository research and development program seeks to provide a sound technical basis for multiple viable disposal options, increase confidence in the robustness of generic disposal concepts, and develop science and engineering tools needed to support disposal concept implementation. Sandia, Los Alamos, and Lawrence Berkeley National labs are conducting research into salt, which includes the Brine Availability Test in Salt (BATS) field test at the Waste Isolation Pilot Plant (WIPP), a DOE Office of Environmental Management facility. BATS is leveraging the existing infrastructure associated with the WIPP to advance science informing generic disposal concepts. BATS is an ongoing heated borehole experiment being conducted in the WIPP underground. The goal of this experiment is to reduce the uncertainty associated with spent nuclear fuel disposition in geologic salt formations. The BATS experiments is designed to accomplish this by increasing our understanding of brine migration in salt, assessing damaged zones from mining and drilling, simulating a post-closure environment, confirming salt properties, and providing data for model validation. Phased testing for generic salt repository R&D was proposed in 2015. Subsequently, borehole thermal testing was developed as a first step in this process. This document summarizes the work conducted on the BATS experiment during fiscal year 2020 with an emphasis on the contributions by Los Alamos National Lab.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

DOE OSTI2024

Electron Influence on the Parallel Proton Firehose Instability in 10-moment, Multifluid Simulations

Instabilities driven by pressure anisotropy play a critical role in modulating the energy transfer in space and astrophysical plasmas. For the first time, we simulate the evolution and saturation of the parallel proton firehose instability using a multifluid model without adding artificial viscosity. These simulations are performed using a 10-moment, multifluid model with local and gradient relaxation heat-flux closures in high-β proton–electron plasmas. When these higher-order moments are included and pressure anisotropy is permitted to develop in all species, we find that the electrons have a significant impact on the saturation of the parallel proton firehose instability, modulating the proton pressure anisotropy as the instability saturates. Even for lower β's more relevant to heliospheric plasmas, we observe a pronounced electron energization in simulations using the gradient relaxation closure. Our results indicate that resolving the electron pressure anisotropy is important to correctly describe the behavior of multispecies plasma systems.

79 ASTRONOMY AND ASTROPHYSICS↗

NASA NTRS1978

An interaction solution algorithm for viscous aerodynamic performance prediction

A weak-interaction solution algorithm is established for aerodynamic flow field prediction about an isolated airfoil. It requires numerical solution of differential equations governing potential flow, viscous and turbulent boundary layer flow, and the turbulent wake flow downstream of a trailing edge. The algorithm accounts for computed viscous displacement effects on the potential flow. These in turn alter the viscous flow through imposed pressure gradients. Closure for turbulence is accomplished using a second order model. Numerical evaluations assess factors affecting solution accuracy, convergence and stability for the combined potential, boundary layer, and parabolic Navier-Stokes equation systems as solved using a finite element algorithm.

Baker, A. J.↗

NASA NTRS1989

A simplified approach to strain energy release rate computations for interlaminar fracture of composites

A simple approach for the strain energy release rate computations based on the finite element method and a singular fitting model is presented. The model uses the stress and displacement distributions at the delamination front. The method is applied to a mixed-mode double cracked-lap-shear composite configuration. The strain energy release rate components predicted by the model are compared with the finite element crack-closure method. The effect of the mesh size on the stress and displacement distribution is isolated. The strain energy release rates predicted by relatively coarse mesh sizes are in good agreement with the finite element crack closure method.

Armanios, Erian A.↗

NASA NTRS1991

Environment enhanced fatigue of advanced aluminum alloys and metal matrix composites

The environmental fatigue crack propagation behavior of advanced Al-Li-Cu based alloys and metal matrix composites is being characterized. Aqueous NaCl and water vapor, which produce atomic hydrogen by reactions on clean crack surfaces, are emphasized. The effects of environment sensitive crack closure, stress ratio, and precipitate microstructure are assessed. Mechanistic models are sought for intrinsic crack tip damage processes to enable predictions of cracking behavior outside of the data, metallurgical improvements in material cracking resistance, and insight on hydrogen compatibility.

Slavik, Donald C.↗

NASA NTRS1991

An efficient solution technique for shockwave-boundary layer interactions with flow separation and slot suction effects

An efficient method for computing two-dimensional compressible Navier-Stokes flow fields is presented. The solution algorithm is a fully-implicit approximate factorization technique based on an unsymmetric line Gauss-Seidel splitting of the equation system Jacobian matrix. Convergence characteristics are improved by the addition of acceleration techniques based on Shamanskii's method for nonlinear equations and Broyden's quasi-Newton update. Characteristic-based differencing of the equations is provided by means of Van Leer's flux vector splitting. In this investigation, emphasis is placed on the fast and accurate computation of shock-wave-boundary layer interactions with and without slot suction effects. In the latter context, a set of numerical boundary conditions for simulating the transpiration flow in an open slot is devised. Both laminar and turbulent cases are considered, with turbulent closure provided by a modified Cebeci-Smith algebraic model. Comparisons with computational and experimental data sets are presented for a variety of interactions, and a fully-coupled simulation of a plenum chamber/inlet flowfield with shock interaction and suction is also shown and discussed.

Edwards, Jack R.↗

NASA NTRS1993

A numerical study of confined turbulent jets

A numerical investigation is reported of turbulent incompressible jets confined in two ducts, one cylindrical and the other conical with a 5 degree divergence. In each case, three Craya-Curtet numbers are considered which correspond, respectively, to flow situations with no moderate and strong recirculation. Turbulence closure is achieved by using the k-epsilon model and a recently proposed realizable Reynolds stress algebraic equation model that relates the Reynolds stresses explicitly to the quadratic terms of the mean velocity gradients and ensures the positiveness of each component of the turbulent kinetic energy. Calculations are carried out with a finite-volume procedure using boundary-fitted curvilinear coordinates. A second-order accurate, bounded convection scheme and sufficiently fine grids are used to prevent the solutions from being contaminated by numerical diffusion. The calculated results are compared extensively with the available experimental data. It is shown that the numerical methods presented are capable of capturing the essential flow features observed in the experiments and that the realizable Reynolds stress algebraic equation model performs much better than the k-epsilon model for this class of flows of great practical importance.

Zhu, J.↗

NASA NTRS1993

2D axisymmetric analysis of SRM ignition transient

To analyze ignition transient of Space Shuttle solid rocket motor, a transient two-dimensional numerical model based on turbulent compressible Navier-Stokes equations in a generalized coordinate system was developed. One-dimensional numerical models (Peretz et al., 1973; Han, 1992; Pardue and Han, 1992) with empirical correlations data obtained from steady turbulent boundary layer flows agrees reasonably well with test rocket data by adjusting a few parameters. However, a 1D model can not provide a physical insight into the complex multidimensional thermal fields and flowfields in the chamber and the converging-diverging rocket nozzle. As an interim step, a 2D model was developed and compared with test data. A modified version of SIMPLE algorithm was used for the numerical model, and the standard k-epsilon model with a wall function was used for turbulence closure. Transient flowfields and thermal fields in the combustion chamber and the attached nozzle were obtained for a selected rocket geometry and propellant. Transient behaviors of the flow and thermal fields were analyzed, and were found to be in good agreement with physical expectations.

Bai, S. D.↗

NASA NTRS2005

Airborne Sunphotometer and Solar Spectral Flux Radiometer Measurements During INTEX/ITCT 2004

During the period 12 July - 8 August 2004, the NASA Ames 14-channel Airborne Tracking Sunphotometer (MTS-14) and Solar Spectral Flux Radiometer (SSFR) were operated aboard a Jetstream 31 (J31) aircraft and acquired measurements during nineteen science flights (approx. 53 flight hours) over the Gulf of Maine in support of the INTEX-NA (INtercontinental chemical Transport Experiment-North America) and ITCT (Intercontinental Transport and Chemical Transformation of anthropogenic pollution) field studies. In this paper, we will present results from analyses of those data sets. AATS-14 measures the direct solar beam transmission at fourteen discrete wavelengths (354-2138 nm), and provides instantaneous measurements of aerosol optical depth (AOD) spectra and water vapor column content, in addition to vertical profiles of aerosol extinction and water vapor density during suitable aircraft ascents and descents. SSFR consists of separate nadir and zenith viewing hemispheric FOV sensors that yield measurements of up- and downwelling solar irradiance at a spectral resolution of approx. 8-12 nm over the wavelength range 300-1700 nm. The objectives of the J31-based measurements during INTEX/ITCT were to provide AOD data for the evaluation of MODIS (MODerate-resolution Imaging Spectroradiometer) and MISR (Multi-angle Imaging Spectro-Radiometer) AOD retrievals, quantify sea surface spectral albedo (which can contribute the largest uncertainty to satellite aerosol retrievals for low aerosol loading), test closure (consistency) among suborbital results, test chemical-transport models using AOD profiles, and assess regional radiative forcing by combining satellite and suborbital results. Specific J31 flight patterns were designed to achieve these objectives, and they included a mixture of vertical profiles (spiral and ramped ascents and descents) and constant altitude horizontal transects at a variety of altitudes. Additional information is included in the original extended abstract.

Russell, P. B.↗

DOE OSTI2025

SCM overview & the EBR-II shutdown heat removal tests validation.

Pronghorn is an engineering-scale, coarse-mesh, thermal-hydraulics tool for supporting reactor-core simulations of advanced nuclear reactors. Most of the current efforts in Pronghorn have been devoted in developing porous finite-volume capabilities and adapting closure correlations for coarse-mesh thermal-hydraulics modeling. However, for liquid-metal reactors (LMRs) with wire-wrapped fuel pin assemblies, a pin-level thermal-hydraulic resolution is required for most safety case studies (pin rupture, channel blockage, etc.). For this purpose, a new Subchannel application is developed in MOOSE, which affords the required flow field resolution, while still preserving an engineering-scale approach. This new solver can be natively coupled to Pronghorn and other MOOSE objects to enable full-core, multi-physics, multi-scale engineering studies. This presentation presents the main features of the SCM code and demonstrates a validation case based on the EBR-II SHRT tests.

21 - SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLAN↗

NASA NTRS1991

High-pressure dielectric measurements of solid hydrogen to 170 GPa

Refractive-index measurements on solid hydrogen at visible frequencies at pressures up to 170 GPa are reported. No evidence is found for a divergence at 150 GPa, close to the low-temperature phase transition observed previously and suggested as being associated with metallization. The results are consistent with closure of the indirect gap. A fit of the data to a dielectric model indicates that the onset of visible absorption due to direct interband transition should occur above 200 GPa, consistent with previous direct observations. Pressure-induced molecular dissociation may occur before closure of the direct gap.

Hemley, R. J.↗

NASA NTRS2002

Wake Closure Characteristics and Afterbody Heating on a Mars Sample Return Orbiter

Aeroheating wind-tunnel tests were conducted on a 0.028 scale model of an orbiter concept considered for a possible Mars sample return mission. The primary experimental objectives were to characterize hypersonic near wake closure and determine if shear layer impingement would occur on the proposed orbiter afterbody at incidence angles necessary for a Martian aerocapture maneuver. Global heat transfer mappings, surface streamline patterns, and shock shapes were obtained in the NASA Langley 20-Inch Mach 6 Air and CF4 Tunnels for post-normal shock Reynolds numbers (based on forebody diameter) ranging from 1,400 to 415,000, angles of attack ranging from -5 to 10 degrees at 0, 3, and 6 degree sideslip, and normal-shock density ratios of 5 and 12. Laminar, transitional, and turbulent shear layer impingement on the cylindrical afterbody was inferred from the measurements and resulted in a localized heating maximum that ranged from 40 to 75 percent of the reference forebody stagnation point heating. Comparison of laminar heating prediction to experimental measurement along the orbiter afterbody highlight grid alignment challenges associated with numerical simulation of three- dimensional separated wake flows. Predicted values of a continuum breakdown parameter revealed significant regions of non-continuum flow downstream of the flow separation at the MSRO shoulder and in the region of the reattachment shock on the afterbody. The presence of these regions suggest that the Navier-Stokes predictions at the laminar wind-tunnel condition may encounter errors in the numerical calculation of the wake shear layer development and impingement due to non-continuum effects.

Horvath, Thomas J.↗

NASA NTRS1987

A transport model of the turbulent scalar-velocity

Performance tests of the third-order turbulence closure for predictions of separating and recirculating flows in backward-facing steps were studied. Computations of the momentum and temperature fields in the flow domain being considered entail the solution of time-averaged transport equations containing the second-order turbulent fluctuating products. The triple products, which are responsible for the diffusive transport of the second-order products, attain greater significance in separating and reattaching flows. The computations are compared with several algebraic models and with the experimental data. The prediction was improved considerably, particularly in the separated shear layer. Computations are further made for the temperature-velocity double products and triple products. Finally, several advantages were observed in the usage of the transport equations for the evaluation of the turbulence triple products; one of the most important features is that the transport model can always take the effects of convection and diffusion into account in strong convective shear flows such as reattaching separated layers while conventional algebraic models cannot account for these effects in the evaluation of turbulence variables.

Amano, R. S.↗

DOE OSTI2022

Modeling of a Generic Pebble Bed High-temperature Gas-cooled Reactor (PB-HTGR) with SAM

This report presents the modeling of the core of a generic pebble-bed reactor (PBR) at the system level using the System Analysis Module (SAM) code. This work is an extension of a previous work by the authors (Ooi et al. (2021)) that used the so-called 2-D ring model approach to model the PBMR-400. With the new approach, the pebble bed of the reactor is modeled with multiple PBCoreChannel components with spherical heat structures which allows the code to calculate thermal fluid parameters with built-in closure relations. The new core-channel approach is an improvement to the 2-D ring model approach as it does not introduce geometric distortions to the model and thus reduces the uncertainties of the predictions. In addition to thermal fluid simulations, point kinetics (PKE) are included to the model. Simulations are performed under a steady-state normal operation condition and a load-following transient scenario. This particular transient scenario is chosen as it tests both the thermal fluid and neutronics aspects of the model. The predicted results from both the steady-state and transient scenarios are compared with the results by Stew- art et al. (2021) who performed similar simulations with a Griffin-Pronghorn coupled tool. Despite the differences between the codes, with SAM being a system-analysis code and Pronghorn being a porous-medium code, both sets of results compare favorably. The over- all profiles and trends of the predicted temperatures and reactivities from the SAM and Griffin/Pronghorn simulations are similar, with some differences in their predicted values. The first part of the report covers the significance of a relatively fast-running approach that is capable of modeling the pebble bed reactor at the system-level while simultaneously capturing the radial thermal behavior of the core. Then, the modeling approach used in this work is discussed in details. Lastly, the results and comparisons with the Griffin/Pronghorn simulation are presented.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

NASA NTRS1989

A numerical and experimental study of confined swirling jets

A numerical and experimental study of a confined strong swirling flow is presented. Detailed velocity measurements are made using a two-component laser Doppler velocimeter (LDV) technique. Computations are performed using a differential second-moment (DSM) closure. The effect of inlet dissipation rate on calculated mean and turbulence fields is investigated. Various model constants are employed in the pressure-strain model to demonstrate their influences on the predicted results. Finally, comparison of the DSM calculations with the algebraic second-monent (ASM) closure results shows that the DSM is better suited for complex swirling flow analysis.

Nikjooy, M.↗

NASA NTRS2001

Toward a Unified Parameterization of the Boundary Layer and Moist Convection: Simulations of Clear and Cloudy Convection - Part 3

A model that employs a new form of mass-flux closure (described in Part I of this paper) is applied to a variety of clear and cloudy planetary boundary layers (PBLs) including dry convection from the Wangara Experiment, trade wind cumulus from the Barbados Oceanographic and Meteorological Experiment (BOMEX), and marine stratocumulus from the Atlantic Stratocumulus Experiment (ASTEX). For Wangara, the simulated variances and fluxes match that expected from similarity arguments, while the mean state is a little less mixed than the observations. In the BOMEX simulation. the shape and magnitude of the fluxes and the turbulence kinetic energy budget agree with LES results and observations. However, the liquid water mixing ratio is too large. This is attributed to an underprediction of the skewness. In agreement with observations from the ASTEX experiment, many of the model-simulated fields distinctly reflect a regime in transition between the trade wind cumulus and the classic stratocumulus-topped boundary layers. In general, the simulated entrainment rate tends to be a little underpredicted in regimes where there is little cloud-top radiative cooling (Wangara and BOMEX), while it is overpredicted in regimes where this process is more critical (e.g., ASTEX). Prior work suggests that this may be related to the manner in which the pressure terms are parameterized in the model. Overall, the model is able to capture some key physical features of these PBL regimes, and appears to have the potential to represent both cloud and boundary layer processes. Thus, this approach is a first step toward unifying these processes in large-scale models.

Lappen, Cara-Lyn↗

DOE OSTI2026

A New 1D Model for Thermal Mixing and Stratification in Advanced Reactor Transients

Thermal mixing and stratification in large pools and enclosures play a critical role in the safety and performance of pool-type nuclear reactors, particularly during transient scenarios involving significant temperature differences between incoming and bulk coolant. Accurate modeling of these phenomena is essential for predicting system behavior and supporting passive safety features such as natural circulation. Here, this paper presents a new 1D model for thermal mixing and stratification, developed and implemented in the SAM code. The model represents a large pool as 1D coolant jet channels and zero-dimensional bulk pool volumes, enabling the simulation of a wide range of flow configurations, including hot and cold jet interactions, stratified layers, and the influence of complex geometries such as ceilings, free surfaces, and internal obstacles. Heat exchange between jet and pool regions is governed by closure relations calibrated against 3D computational fluid dynamics (CFD) simulations. The model improves upon earlier approaches by incorporating time-dependent jet characteristics and capturing the associated delay effects more accurately. Code-to-code comparisons and validation against experimental data from the Thermal Stratification Test Facility demonstrate the model’s accuracy and flexibility. This work offers two key contributions: (1) an efficient and robust method for simulating thermal mixing and stratification at the system level, eliminating the need for external coupling between system analysis codes and CFD, and (2) a significant enhancement of SAM’s capabilities to analyze thermal stratification phenomena in advanced reactor systems.

SAM↗