Engineering Papers⌕ Search

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.

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

At least 1,351 records · Page 75

Optical evidence for the metallization of Xenon at 132(5) GPa

Xe has been compressed in a diamond-anvil cell to approximately 200 GPa. The metallization of Xe by band-gap closure was investigated by obtaining optical data in both the metallic and insulating states. In the metallic state, the pressure dependence of the plasma frequency was determined from absorption data fitted with a free-electron model. In the insulating state, the pressure dependence of the band gap was determined from absorption data fitted with an indirect-band-gap model. The optical data indicate that the insulator-to-metal transition in Xe occurs at 132(5) GPa.

Goettel, Kenneth A.↗

Time Evolution of Modeled Reynolds Stresses in Planar Homogeneous Flows

The analytic expression of the time evolution of the Reynolds stress anisotropy tensor in all planar homogeneous flows is obtained by exact integration of the modeled differential Reynolds stress equations. The procedure is based on results of tensor representation theory, is applicable for general pressure-strain correlation tensors, and can account for any additional turbulence anisotropy effects included in the closure. An explicit solution of the resulting system of scalar ordinary differential equations is obtained for the case of a linear pressure-strain correlation tensor. The properties of this solution are discussed, and the dynamic behavior of the Reynolds stresses is studied, including limit cycles and sensitivity to initial anisotropies.

Jongen, T.↗

Overview of the Aerothermodynamics Analysis Conducted in Support of the STS-107 Accident Investigation

A graphic presentation of the aerothermodynamics analysis conducted in support of the STS-107 accident investigation. Investigation efforts were conducted as part of an integrated AATS team (Aero, Aerothermal, Thermal, Stress) directed by OVEWG. Graphics presented are: STS-107 Entry trajectory and timeline (1st off-nominal event to Post-LOS); Indications from OI telemetry data; Aero/aerothermo/thermal analysis process; Selected STS-107 side fuselage/OMS pod off-nominal temperatures; Leading edge structural subsystem; Relevant forensics evidence; External aerothermal environments; STS-107 Pre-entry EOM3 heating profile; Surface heating and temperatures; Orbiter wing leading edge damage survey; Internal aerothermal environments; Orbiter wing CAD model; Aerodynamic flight reconstruction; Chronology of aerodynamic/aerothermoydynamic contributions; Acreage TPS tile damage; Larger OML perturbations; Missing RCC panel(s); Localized damage to RCC panel/missing T-seal; RCC breach with flow ingestion; and Aero-aerothermal closure. NAIT served as the interface between the CAIB and NASA investigation teams; and CAIB requests for study were addressed.

Campbell, Charles H.↗

Analytics-at-scale of Sensor Data for Digital Monitoring in Nuclear Plants (4th Annual Report)

Nuclear plant sites collect and store large volumes of data collected from various equipment and systems. These datasets typically include plant process parameters, maintenance records, technical logs, online monitoring data, and equipment failure data. The collection of such data affords an opportunity to leverage data-driven machine learning and artificial intelligence technologies to provide diagnostic and prognostic capabilities within the nuclear power industry to reduce operating and maintenance costs. In this way, nuclear energy can become more economically competitive with other energy sources, and premature closures can be avoided. From a maintenance standpoint, savings can be achieved by leveraging machine learning and artificial intelligence technologies to develop data-driven algorithms to better diagnose and predict potential faults within the system. Improved model accuracy can lead to reductions in unnecessary maintenance and more efficient planning of future maintenance, thus lowering the costs associated with parts, labor, and unnecessary planned, forced, or extended outages. From an operations perspective, cost savings can be generated by shifting from route-based monitoring to wireless technologies for online monitoring, and by transitioning from onsite- to cloud-based computing and storage services. Wireless monitoring would reduce the operator manhours required for taking routine measurements, while cloud computing services would generate cost savings by reducing the amount of hardware needing to be purchased and maintained—all while scaling to both computational and storage demands. This report summarizes this project’s effort to shift from costly, labor-intensive preventative maintenance to cheaper predictive maintenance.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Dynamics of the middle atmosphere

It is demonstrated by means of model calculations that while the general configuration of middle atmosphere dynamics (i.e., winter westerlies and summer easterlies) is determined by differential solar heating, the deviations of the zonally-averaged temperature field from radiative equilibrium and the closure of the jet structures with increasing altitude result from the action of zonal mean momentum dissipation process. The apparent heating and acceleration of the mean zonal state by planetary wave heat and momentum fluxes are examined using an Eulerian framework, and it is demonstrated that these are overestimates of their net effect. It is argued that since decelerations of the mean zonal flow are required in both winter and summer, and planetary waves are known to be very weak in the summer middle atmosphere, gravity waves are probably responsible for most of the middle atmosphere momentum dissipation as a result of their attenuation with height above their wave breaking altitude. Radar observations of middle atmosphere dynamics together with some theoretical work indicates that the effects of breaking gravity waves should be included in the thermodynamic equation as well as the momentum equation.

Geller, M. A.↗

Closures and Simulation for Thermal Radiation Transport in Stochastic Media with Nonlinear Temperature Dependence

Because of the practical challenge of rendering very complex realistic spatial structures for numerical work, it is common practice to resort to characterizing such media as stochastic mixtures of materials, ideally parametrized with low order statistics such as the mean, variance, and correlation functions of the now random material properties. This enables realizations of the medium to be repeatedly generated and radiation transport computations to, in principle, be performed for a large ensemble of these realizations to obtain a statistically well-characterized radiation field. Statistical post-processing yields desired quantities such as conditional and unconditional mean radiation flux and probability distributions of transmitted radiation. However, such computations prove expensive for all but the simplest stochastic geometries and are most suited for benchmarking approximate models. The most common approximations lead to homogenized media so that transport computations are required only on a single medium realization but by construct provide only limited statistical information on the radiation field. Almost all approximate approaches to this problem attempt to develop equations for low order moments of the radiation intensity (mean, second moment, correlation function) but inevitably encounter a closure problem: the equation for any statistical moment will contain terms depending on unknown higher-order moments. Thus, the challenge shifts to one of developing closure relations that relate the unknown moments to the lower order moments. Under very special conditions, an exact closure can be derived but in general closures are heuristically stated constitutive relations. Also, closure approaches depend on whether the mixing statistics are spatially and/or temporally continuous as in fluctuating turbulent fields, or discontinuous as in randomly mixed solid chunks of material. Thus, unconditional averaging is generally applied in the former case but conditional averaging is more appropriate when the mixing is discontinuous. In this work, the emphasis is on binary statistical mixtures of immiscible fluids as as such the quantities of interest are averages (flux, temperature) conditioned on the material type.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

An implicit finite-difference code for a two-equation turbulence model for three-dimensional flows

An implicit finite difference code was developed which solves the transport equations for the turbulence kinetic energy and its dissipation rate in generalized coordinates in three dimensions. The finite difference equations are solved using the Beam-Warming algorithm. The kinetic energy-dissipation code, KEM, provides the closure; i.e., the turbulent viscosity for calculation of either compressible or incompressible flows. Turbulent internal flow over a backward-facing step has been calculated using the present code in conjunction with the Incompressible Navier-Stokes Code, INS3D. The results are in good agreement with experiments and two dimensional computations of other researchers.

Kaul, U. K.↗

Reduced order models for thermal radiative transfer problems based on moment equations and data-driven approximations of the Eddington tensor

Here a new group of structure and asymptotic preserving reduced-order models (ROMs) for multidimensional nonlinear thermal radiative transfer (TRT) problems is presented. They are formulated by means of the nonlinear projective approach and data compression techniques. The nonlinear projection is applied to the Boltzmann transport equation (BTE) to derive a hierarchy of low-order moment equations. Approximation of the Eddington tensor that provides exact closure for the system of moment equations is found with projection-based data-driven methodologies. These include the (i) proper orthogonal decomposition (POD), (ii) dynamic mode decomposition (DMD) and (iii) a variant of the DMD. A parameterization is derived for this ROM for the temperature of radiation incoming to the problem domain (the radiation drive temperature). This parameterization is informed from results of a dimensionless study of the TRT problem. Analysis of the ROMs is performed on the classical Fleck-Cummings TRT multigroup test problem in 2D geometry with a radiation-driven Marshak wave. Numerical results are presented to demonstrate the performance of these ROMs for the simulation of evolving radiation and heat waves. Results show these models to be sufficiently accurate for practical computations with rather low-rank representations of the Eddington tensor. As the rank of the approximation is increased, the errors of solutions generated by the ROMs gradually decreases.

42 ENGINEERING↗

Assessing Impacts on Pressure Stabilization and Leasing Acreage for CO 2 Storage Utilizing Oil Migration Concepts

Favorable geological storage for CO 2 has long been pictured as large anticlines with thick sandstones, similar to oil reservoirs in the petroleum system. Unlike oil, however, stored CO 2 does not need to be recoverable, which raises the possibility of using dissolution and residual trapping to augment the capacity of buoyant traps and tap more of the bulk rock volume. The work presented builds on that idea, asking the following question: If we inject CO 2 down to a syncline – analogous to the carrier bed in the petroleum system – how would this injection mechanism impact storage capacity and plume shape, migration, and stabilization? To address this question, we built a reservoir model, based on seismic interpretation of Middle Miocene strata, offshore Galveston, Texas. 3-D seismic and well logs were used to characterize key intervals. Reservoirs chosen for modeling are progradational-aggradational sands with mud intercalation. They have a higher degree of heterogeneity than the more conventional reservoirs commonly targeted for CO 2 storage. Modeling investigated how far the CO 2 plume would migrate under two scenarios: (1) injecting CO 2 at the base of the salt withdrawal basin (syncline scenario) and (2) injecting CO 2 at the base of the structural closure, similar to a common injection well location for EOR purposes (base scenario). For each scenario, we separately simulated injection of 30 MT of CO 2 and 60 MT of CO 2 continuously for 30 years and observe the plume and pressure evolution for 100 years after the injection stops. The simulation shows that injecting the CO 2 into a syncline limits the vertical migration of CO 2 , thus making synclinal injection more secure. In the syncline scenario, the geological layer around the injection point is more heterogeneous than the layer in the base scenario; thus, the CO 2 tends to migrate laterally. Additionally, in the syncline scenario, the plume does not even reach the upper part of the anticline, allowing us to safely store an additional amount of CO 2 into the reservoir. Furthermore, the simulation also shows that in the syncline scenario, the times needed for the reservoir to reach its stabilized pressure after the end of injections are faster. To summarize, CO 2 injection at the base of a syncline could provide additional storage, increase the safety of the project from the limited vertical plume migration, and expedite plume stabilization, which could result in the decrease of monitoring frequency as the project runs, thus lowering the operating cost of the project in the long run.

58 GEOSCIENCES↗

Interaction of a cumulus cloud ensemble with the large-scale environment. III - Semi-prognostic test of the Arakawa-Schubert cumulus parameterization

The verification of the Arakawa and Schubert (1974) cumulus parameterization is continued using a semiprognostic approach. Observed data from Phase III of GATE are used to provide estimates of the large-scale forcing of a cumulus ensemble at each observation time. Instantaneous values of the precipitation and the warming and drying due to cumulus convection are calculated using the parameterization. The results show that the calculated precipitation agrees very well with estimates from the observed large-scale moisture budget and from radar observations. The calculated vertical profiles of cumulus warming and drying also are quite similar to the observed. It is shown that the closure assumption adopted in the parameterization (the cloud-work function quasi-equilibrium) results in errors of generally less than 10% in the calculated precipitation. The sensitivity of the parameterization to some assumptions of the cloud ensemble model and the solution method for the cloud-base mass flux is investigated.

Lord, S. J.↗

Structural Sizing of a Horizontal Take-Off Launch Vehicle with an Air Collection and Enrichment System

In support of NASA's Next Generation Launch Technology (NGLT) program, the Andrews Gryphon booster was studied. The Andrews Gryphon concept is a horizontal lift-off, two-stage-to-orbit, reusable launch vehicle that uses an air collection and enrichment system (ACES). The purpose of the ACES is to collect atmospheric oxygen during a subsonic flight loiter phase and cool it to cryogenic temperature, ultimately resulting in a reduced initial take-off weight To study the performance and size of an air-collection based booster, an initial airplane like shape was established as a baseline and modeled in a vehicle sizing code. The code, SIZER, contains a general series of volume, surface area, and fuel fraction relationships that tie engine and ACES performance with propellant requirements and volumetric constraints in order to establish vehicle closure for the given mission. A key element of system level weight optimization is the use of the SIZER program that provides rapid convergence and a great deal of flexibility for different tank architectures and material suites in order to study their impact on gross lift-off weight. This paper discusses important elements of the sizing code architecture followed by highlights of the baseline booster study.

McCurdy, David R.↗

Development of a Unified SHOC-EDMF PBL Parameterization in the NASA GEOS Model

We present results from a novel extension of the boundary layer parameterization in the NASA Goddard Earth Observing System (GEOS) model. An Eddy Diffusivity-Mass Flux (EDMF) multi-plume updraft scheme was previously implemented in GEOS. Here we replace the K-profile eddy diffusivity scheme of Lock et. al. (2000) with the Simplified Higher Order Closure (SHOC). The scheme prognoses turbulent kinetic energy, and derives buoyancy flux, cloud fraction and condensate from a trivariate double-gaussian joint PDF. The original SHOC scheme is modified in two significant respects: (i) contributions from the EDMF updrafts are included in estimates of all higher order moments used to determine the PDF; and (ii) the second and third moments of total water are made prognostic variables, with source terms associated with EDMF updrafts. The predicted higher order moments as well as updraft and environmental properties are compared with large eddy simulations in several case studies.

SHOC-EDMF↗

A Unified SHOC-MF Boundary Layer Parameterization in the NASA Goddard Earth Observing System Model

We present results from a novel extension of the boundary layer parameterization in the NASA Goddard Earth Observing System (GEOS) model. An Eddy Diffusivity-Mass Flux (EDMF) multi-plume updraft scheme was previously implemented in GEOS. Here we replace the K-profile eddy diffusivity scheme of Lock et. al. (2000) with the Simplified Higher Order Closure (SHOC; Bogenschutz and Krueger, 2013). The scheme prognoses turbulent kinetic energy, and derives buoyancy flux, cloud fraction and condensate from a trivariate double-gaussian joint PDF. The original SHOC scheme is modified in two significant respects: (i) contributions from the EDMF updrafts are included in estimates of all higher order moments used to determine the PDF; and (ii) the second and third moments of total water are made prognostic variables, with source terms associated with EDMF updrafts. The predicted higher order moments as well as updraft and environmental properties are compared with large eddy simulations in several case studies. Results indicate that many higher order moments are underestimated when the mass flux contribution is neglected.

Nathan Arnold↗

A Unified SHOC-MF Boundary Layer Parameterization in the NASA Goddard Earth Observing System Model

We present results from a novel extension of the boundary layer parameterization in the NASA Goddard Earth Observing System (GEOS) model. An Eddy Diffusivity-Mass Flux (EDMF) multi-plume updraft scheme was previously implemented in GEOS. Here we replace the K-profile eddy diffusivity scheme of Lock et. al. (2000) with the Simplified Higher Order Closure (SHOC; Bogenschutz and Krueger, 2013). The scheme prognoses turbulent kinetic energy, and derives buoyancy flux, cloud fraction and condensate from a trivariate double-gaussian joint PDF. The original SHOC scheme is modified in two significant respects: (i) contributions from the EDMF updrafts are included in estimates of all higher order moments used to determine the PDF; and (ii) the second and third moments of total water are made prognostic variables, with source terms associated with EDMF updrafts. The predicted higher order moments as well as updraft and environmental properties are compared with large eddy simulations in several case studies. Results indicate that many higher order moments are underestimated when the mass flux contribution is neglected.

Nathan Arnold↗

Reversed-direction 2-point modelling applied to divertor conditions in DIII-D *

Abstract A predictive form of the extended 2-point model known as the ‘reverse 2-point model’, Rev2PM, is applied to a range of detachment levels in the open lower divertor of DIII-D, showing that the experimentally measured electron temperature ( T e ) and pressure ( p e ) at the divertor entrance can be calculated within 50% from target measurements, if and only if a posteriori corrections for convective heat flux are included in the model. Unlike the standard 2-point model, the Rev2PM calculates upstream scrape-off layer (SOL) quantities (such as separatrix T e and p e ) from target conditions (such as T e and parallel heat flux), with volumetric power and momentum losses depending solely on target T e . The Rev2PM is tested against a database of DIII-D inter-ELM divertor Thomson scattering measurements, built from a series of 6 MW, 1.3 MA, LSN H-mode discharges with varied main ion density, drift direction, and nitrogen puffing rate. Measured target T e ranged from 0.4–25 eV over this database, and upstream T e ranged from 5–60 eV. Poor agreement is found between upstream measurements and Rev2PM calculations that assume purely conductive parallel heat transport. However, introducing a posteriori corrections to account for convective heat transport brings the Rev2PM calculations within 50% of the measured upstream values across the dataset. These corrections imply that up to 99% of the parallel heat flux is carried by convection in detached conditions in the DIII-D open lower divertor, though further work is required to assess any potential dependencies on device size or divertor closure.

2-point model↗

Intruder band mixing in an ab initio description of 12 Be

The spectrum of 12 Be exhibits exotic features, e.g., an intruder ground state and shape coexistence, normally associated with the breakdown of a shell closure. While previous phenomenological treatments indicated the ground state has substantial contributions from intruder configurations, it is only with advances in computational abilities and improved interactions that this intruder mixing is observed in ab initio no-core shell model (NCSM) predictions. In this work, we extract electromagnetic observables and symmetry decompositions from the NCSM wave functions to demonstrate that the low-lying positive parity spectrum can be explained in terms of mixing of rotational bands with very different intrinsic structure coexisting within the low-lying spectrum. These observed bands exhibit an approximate SU(3) symmetry and are qualitatively consistent with Elliott model predictions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Crossing N = 28 Toward the Neutron Drip Line: First Measurement of Half-Lives at FRIB

Here, new half-lives for exotic isotopes approaching the neutron drip-line in the vicinity of N~28 for Z=12–15 were measured at the Facility for Rare Isotope Beams (FRIB) with the FRIB decay station initiator. The first experimental results are compared to the latest quasiparticle random phase approximation and shell-model calculations. Overall, the measured half-lives are consistent with the available theoretical descriptions and suggest a well-developed region of deformation below 48 Ca in the N=28 isotones. The erosion of the Z=14 subshell closure in Si is experimentally confirmed at N=28, and a reduction in the 38 Mg half-life is observed as compared with its isotopic neighbors, which does not seem to be predicted well based on the decay energy and deformation trends. This highlights the need for both additional data in this very exotic region, and for more advanced theoretical efforts.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

GDSA Framework Development and Process Model Integration FY2022

The Spent Fuel and Waste Science and Technology (SFWST) Campaign of the U.S. Department of Energy (DOE) Office of Nuclear Energy (NE), Office of Spent Fuel & Waste Disposition (SFWD) is conducting research and development (R&D) on geologic disposal of spent nuclear fuel (SNF) and high-level nuclear waste (HLW). A high priority for SFWST disposal R&D is disposal system modeling (Sassani et al. 2021). The SFWST Geologic Disposal Safety Assessment (GDSA) work package is charged with developing a disposal system modeling and analysis capability for evaluating generic disposal system performance for nuclear waste in geologic media. This report describes fiscal year (FY) 2022 advances of the Geologic Disposal Safety Assessment (GDSA) performance assessment (PA) development groups of the SFWST Campaign. The common mission of these groups is to develop a geologic disposal system modeling capability for nuclear waste that can be used to assess probabilistically the performance of generic disposal options and generic sites. The modeling capability under development is called GDSA Framework (pa.sandia.gov). GDSA Framework is a coordinated set of codes and databases designed for probabilistically simulating the release and transport of disposed radionuclides from a repository to the biosphere for post-closure performance assessment. Primary components of GDSA Framework include PFLOTRAN to simulate the major features, events, and processes (FEPs) over time, Dakota to propagate uncertainty and analyze sensitivities, meshing codes to define the domain, and various other software for rendering properties, processing data, and visualizing results.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗