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 505 records · Page 28
All sourcesNASA NTRSDOE OSTI
NASA NTRS2000

Troughs on Martian Ice Sheets: Analysis of Their Closure and Mass Balance

At the Copenhagen workshop on Martian polar processes, Ralf Greve commented that the flow regime surrounding scarps and troughs of the Martian polar ice sheets cannot be modeled using traditional "plan view" ice-sheet models. Such models are inadequate because they typically use reduced equations that embody certain simplifications applicable only to terrestrial ice sheets where the upper ice sheet surface is smooth. In response to this suggestion, we have constructed a 2-dimensional, time dependent "side view" (two spatial dimensions: one horizontal, one vertical) model of scarp closure that is designed to overcome the difficulties described by Greve. The purpose of the model is to evaluate the scales of stress variation and styles of flow closure so as to estimate errors that may be encountered by "plan view" models. We show that there may be avenues whereby the complications associated with scarp closure can be overcome in "plan view" models through appropriate parameterizations of 3-dimensional effects. Following this, we apply the flow model to simulate the evolution of a typical scarp on the North Polar Cap of Mars. Our simulations investigate: (a) the role of "radiation trapping" (see our companion abstract) in creating and maintaining "spiral-like" scarps on the ice sheet, (b) the consequences of different flowlaws and ice compositions on scarp evolution and, in particular, scarp age, and (c) the role of dust and debris in scarp evolution.

Fountain, A.↗
DOE OSTI2020

Interface learning in fluid dynamics: Statistical inference of closures within micro–macro-coupling models

Many complex multiphysics systems in fluid dynamics involve using solvers with varied levels of approximations in different regions of the computational domain to resolve multiple spatiotemporal scales present in the flow. The accuracy of the solution is governed by how the information is exchanged between these solvers at the interface and several methods have been devised for such coupling problems. In this article, we construct a data-driven model by spatially coupling a microscale lattice Boltzmann method (LBM) solver and macroscale finite difference method (FDM) solver for reaction-diffusion systems. The coupling between the micro-macro solvers has one to many mapping at the interface leading to the interface closure problem, and we propose a statistical inference method based on neural networks to learn this closure relation. The performance of the proposed framework in a bifidelity setting partitioned between the FDM and LBM domain shows its promise for complex systems where analytical relations between micro-macro solvers are not available.

42 ENGINEERING↗
DOE OSTI2020

PORFLOW modeling of Vadose zone flow and transport for the E-area intermediate level vault

In support of the E-Area Performance Assessment, a two-dimensional model of water flow and radionuclide transport through the E-Area Intermediate Level Vault (ILV) and local vadose zone has been developed using the PORFLOWTM software. The purpose of the model is to calculate flux to the water table for radionuclides eluted from the ILV during its operational life, the period of institutional control, and times following site closure. Results of model calculations will be used by a three-dimensional PORFLOW model of transport through the aquifer to determine radionuclide concentrations at a hypothetical 100 meter well and at the site boundary where contaminated groundwater could be assessed by members of the public following site closure. While newly developed, the model structure closely follows that used in the 2008 PA while incorporating a refined computational mesh, updated material properties for the vadose zone soil and vault concrete, revised infiltration rates, and new solid-liquid distribution coefficients. The model also addresses degradation of concrete hydraulic properties by blending soil and concrete water retention curves over a 500 year period. That is, the hydraulic properties of ILV concrete start out as fresh concrete at the time of site closure and degrade to soil properties over a 500 year time period. This approach has not been used previously in E-Area PA’s. The model has been used to calculate water flow through the ILV vadose zone and, on a trial basis, transport of a limited number of radionuclides from the waste region inside the ILV to the water table. Figure ES-1 provides a view of the model and an example of water flow after the concrete has substantially degraded. Figure dimensions are in centimeters. The red area in the left hand side figure is the waste disposal region within the ILV. Colored bands in the figures indicate different soil and concrete regions. Figure ES-2 shows a calculation of flux to the water table (mol/yr) for U-234 and its one-year half-life decay chain (U-234 → Th-230 → Ra-226 → Pb-210). Pb-210 flux is too small to show on the same scale used for Th-230 and Ra-226 in Figure ES-2. The model has also been used to make a preliminary evaluation of the release of tritium from TPBAR disposal containers that will be place in the ILV. A separately calculated tritium release rate as a function of time was used as a source in the PORFLOW model. The TPBAR disposal schedule and tritium release rate differ from the values used in the 2008 PA. Results from these initial trial runs are compared to calculations made for the 2008 PA ILV analysis. In brief, a new PORFLOW model of flow and transport through the ILV vadose zone has been developed and tested. The model is ready to be used for PA calculations.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
NASA NTRS1992

Mechanisms of fatigue crack retardation following single tensile overloads in powder metallurgy aluminum alloys

In ingot metallurgy (IM) alloys, the number of delay cycles following a single tensile overload typically increases from a minimum at an intermediate baseline stress intensity range, Delta-K(B), with decreasing Delta-K(B) approaching threshold and increasing Delta-K(B) approaching unstable fracture to produce a characteristic 'U' shaped curve. Two models have been proposed to explain this behavior. One model is based on the interaction between roughness and plasticity-induced closure, while the other model only utilizes plasticity-induced closure. This article examines these models, using experimental results from constant amplitude and single overload fatigue tests performed on two powder metallurgy (PM) aluminum alloys, AL-905XL and AA 8009. The results indicate that the 'U'-shaped curve is primarily due to plasticity-induced closure, and that the plasticity-induced retardation effect is through-thickness in nature, occurring in both the surface and interior regions. However, the retardation effect is greater at the surface, because the increase in plastic strain at the crack tip and overload plastic zone size are larger in the plane-stress surface regions than in the plane-strain interior regions. These results are not entirely consistent with either of the proposed models.

Bray, G. H.↗
DOE OSTI2020

ZAM Modeling Study to Support the Tank Closure Cesium Removal (TCCR) 1A Unit

Currently at the Savannah River Site (SRS), the Tank Closure Cesium Removal (TCCR) is an “at-tank” process designed to remove cesium from aqueous tank waste. Cesium will be removed by ion exchange using the engineered IONSIV® R9120 form of Crystalline Silicotitanate (CST) media. The current TCCR design has two columns online in a lead-lag configuration to optimize media usage and achieve the target decontamination. Once the lead column is saturated with cesium, it will be removed from service, the lag column will rotate into the lead position, and a new column with fresh ion-exchange media will be placed into the lag position. The TCCR process for cesium removal from Tank 10H is detailed in X-SOW-H- 00002. Demonstration of the system began in early calendar year 2019 with two batches of salt solution generated by dissolving saltcake in Tank 10H, followed by processing of these batches through the TCCR system. A third TCCR Tank 10H dissolved saltcake batch is scheduled for processing soon. Upon completion of the demonstration with Tank 10H dissolved saltcake, Tank 9H salt solution will be transferred to Tank 10H and subsequently processed through the TCCR unit with new CST media (referred to as R9120-B 30x60) added to new IX columns. The TCCR processing campaign of Tank 9H salt solution is referred to as TCCR-1A .

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2025

Implications of a weakening N = 126 shell closure away from stability for r -process astrophysical conditions

The formation of the third r-process abundance peak near A ∼ 195 is highly sensitive to both nuclear structure far from stability and the astrophysical conditions that produce the heaviest elements. In particular, the N = 126 shell closure plays a crucial role in shaping this peak. Experimental data hints that the shell weakens as proton number departs from Z = 82, a trend largely missed by global mass models. To investigate its impact on r-process nucleosynthesis, we employ both standard global models with strong closures and modified Duflo-Zuker (DZ) models that reproduce the weakening, combined with three sets of β − -decay rates. Strong shell closures generate sharply peaked abundances, whereas weakened closures consistent with the experimental trend produce broader, flatter patterns. Accurately reproducing the solar third peak under weakened shell strength requires sufficiently neutron-rich conditions that significant fission occurs, and slower decay rates. These results demonstrate that a weakening N = 126 shell closure away from stability imposes significant constraints on the astrophysical environments of the r-process and underscores the need for precise mass measurements and improved characterization of β − -decay properties in this region.

N = 126 closed shells↗
NASA NTRS2006

Second Moment Closure Near the Two-component Limit

The purpose of this paper is to explore some wider implications of the two-component limit for both single point turbulence models and spectral closure theories. Although the two-component limit arises most naturally in inhomogeneous problems like wall-bounded turbulence, the analysis will be restricted to homogeneous turbulence. But since homogeneous turbulence is the crucial case for realizability, the conclusions will nevertheless be applicable to modeling. Th essential point of our argument is that whereas the evolution of the stochastic velocity field is Markovian because it is governed by the Navier-Stokes equations, the exact stress evolution equation is not Markovian because it is unclosed. This property of moment evolution has been stressed by Kraichnan (1959). We will show that modeling stress evolution at the two-component limit with a closure that is Markovian in the stresses alone leads to basic inconsistencies in single-point modeling and, perhaps surprisingly, in spectral modes as well.

Rubinstein, Robert↗
DOE OSTI2022

A coupled ground heat flux–surface energy balance model of evaporation using thermal remote sensing observations

Abstract. One of the major undetermined problems in evaporation (ET) retrieval using thermal infrared remote sensing is the lack of a physically based ground heat flux (G) model and its integration within the surface energy balance (SEB) equation. Here, we present a novel approach based on coupling a thermal inertia (TI)-based mechanistic G model with an analytical surface energy balance model, Surface Temperature Initiated Closure (STIC, version STIC1.2). The coupled model is named STIC-TI. The model is driven by noon–night (13:30 and 01:30 local time) land surface temperature, surface albedo, and a vegetation index from MODIS Aqua in conjunction with a clear-sky net radiation sub-model and ancillary meteorological information. SEB flux estimates from STIC-TI were evaluated with respect to the in situ fluxes from eddy covariance measurements in diverse ecosystems of contrasting aridity in both the Northern Hemisphere and Southern Hemisphere. Sensitivity analysis revealed substantial sensitivity of STIC-TI-derived fluxes due to the land surface temperature uncertainty. An evaluation of noontime G (Gi) estimates showed 12 %–21 % error across six flux tower sites, and a comparison between STIC-TI versus empirical G models also revealed the substantially better performance of the former. While the instantaneous noontime net radiation (RNi) and latent heat flux (LEi) were overestimated (15 % and 25 %), sensible heat flux (Hi) was underestimated (22 %). Overestimation (underestimation) of LEi (Hi) was associated with the overestimation of net available energy (RNi−Gi) and use of unclosed surface energy balance flux measurements in LEi (Hi) validation. The mean percent deviations in Gi and Hi estimates were found to be strongly correlated with satellite day–night view angle difference in parabolic and linear pattern, and a relatively weak correlation was found between day–night view angle difference versus LEi deviation. Findings from this parameter-sparse coupled G–ET model can make a valuable contribution to mapping and monitoring the spatiotemporal variability of ecosystem water stress and evaporation using noon–night thermal infrared observations from future Earth observation satellite missions such as TRISHNA, LSTM, and SBG.

Bhattacharya, Bimal K.↗
NASA NTRS1997

Multistage Schemes with Multigrid for Euler and Navier-Strokes Equations: Components and Analysis

A class of explicit multistage time-stepping schemes with centered spatial differencing and multigrids are considered for the compressible Euler and Navier-Stokes equations. These schemes are the basis for a family of computer programs (flow codes with multigrid (FLOMG) series) currently used to solve a wide range of fluid dynamics problems, including internal and external flows. In this paper, the components of these multistage time-stepping schemes are defined, discussed, and in many cases analyzed to provide additional insight into their behavior. Special emphasis is given to numerical dissipation, stability of Runge-Kutta schemes, and the convergence acceleration techniques of multigrid and implicit residual smoothing. Both the Baldwin and Lomax algebraic equilibrium model and the Johnson and King one-half equation nonequilibrium model are used to establish turbulence closure. Implementation of these models is described.

Swanson, R. C.↗
NASA NTRS1980

Numerical solutions of the Navier-Stokes equations for transonic afterbody flows

The time dependent Navier-Stokes equations in mass averaged variables are solved for transonic flow over axisymmetric boattail plume simulator configurations. Numerical solution of these equations is accomplished with the unsplit explict finite difference algorithm of MacCormack. A grid subcycling procedure and computer code vectorization are used to improve computational efficiency. The two layer algebraic turbulence models of Cebeci-Smith and Baldwin-Lomax are employed for investigating turbulence closure. Two relaxation models based on these baseline models are also considered. Results in the form of surface pressure distribution for three different circular arc boattails at two free stream Mach numbers are compared with experimental data. The pressures in the recirculating flow region for all separated cases are poorly predicted with the baseline turbulence models. Significant improvements in the predictions are usually obtained by using the relaxation models.

Swanson, R. C., Jr.↗
NASA NTRS1980

A numerical study of a marine subtropical stratus cloud layer and its stability

A model for numerical simulation of stratus cloud layers is constructed by combining a second-order closure, turbulent transfer model with a thermal radiative transfer model. The turbulent transfer model allows water vapor saturation. The combined turbulence-radiation model is applied to both a horizontally uniform one-dimensional case and a horizontally nonuniform two-dimensional case. In the latter, the dynamics of mesoscale circulations are also incorporated. Results of the two-dimensional simulation show that the layer cloud instability occurs where the sea surface temperature is high and the large-scale subsidence is weak. The simulated instability is analyzed in view of an instability criterion, the eddy kinetic energy budget, and evaporative cooling near the cloud top.

Moeng, C.-H.↗
DOE OSTI2024

Estimation of General Area Dose Rates During Closure Weld Operations on the DOE Standard Canister Pacific Basin Nuclear Conference 2024

The Department of Energy (DOE) Spent Nuclear Fuel (SNF) Road Ready Demonstration Project is an effort to package DOE-managed SNF for road-ready dry storage. A key operation in this demonstration is performing the closure weld on the DOE Standard Canister after it has been loaded with the SNF. Planning the closure weld operations is ongoing and must consider what parts of the operations (if any) can be done manually and what parts must be done remotely. A key parameter to inform the evolving planning of these operations is the dose rate in the vicinity of the canister and the dose hazard posed to the closure weld operators and collocated workers. This study developed a Monte Carlo N-Particle (MCNP) model of the closure weld operation, which was used to estimate general area dose rates as well as dose rates to the welding location, weld operators, and collocated workers. These estimates provide some insight into the potential radiation exposure hazard for workers during the operations and are a factor to be considered in the ongoing development of designs, technology, processes, and administrative controls for packaging DOE-managed SNF for road-ready dry storage.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

Analysis of Compressibility Effects and Nonlinear Property Variations in a Supercritical CO2 Mixing Layer

Calculations are performed using the large eddy simulation technique to qualitatively assess the effects of fluid compressibility and thermodynamic nonlinearities on the dynamics of a supercritical flow field. Resulting data are also used to interrogate instantaneous subfilter-scale fields, as modeled by the mixed-dynamic Smagorinsky closure and gradient diffusion model. A three-dimensional wall-resolved calculation of a spatially evolving mixing layer composed solely of carbon dioxide is performed. This is followed by a coarse grid calculation to examine model performance as a function of resolution. Theoretical analysis indicates that the partial derivatives of density with respect to pressure and temperature are important modulators of the temperature and pressure field evolution, respectively. Results pertaining to subfilter velocity-velocity and velocity-temperature fields indicate that the models are most stressed in the regions where compressibility and thermo-fluid nonlinearities dominate the physics.

42 ENGINEERING↗
NASA NTRS1989

The generalized drift flux approach: Identification of the void-drift closure law

The main characteristics and the potential advantages of generalized drift flux models are presented. In particular it is stressed that the issue on the propagation properties and on the mathematical nature (hyperbolic or not) of the model and the problem of closure are easier to tackle than in two fluid models. The problem of identifying the differential void-drift closure law inherent to generalized drift flux models is then addressed. Such a void-drift closure, based on wave properties, is proposed for bubbly flows. It involves a drift relaxation time which is of the order of 0.25 s. It is observed that, although wave properties provide essential closure validity tests, they do not represent an easily usable source of quantitative information on the closure laws.

Boure, J. A.↗
NASA NTRS1991

Large-eddy simulation of compressible turbulence

The increase in the range of length scales with increasing Reynolds number limits the direct simulation of turbulent flows to relatively simple geometries and low Reynolds numbers. However, since most flows of engineering interest occur at much higher Reynolds number than is currently within the capabilities of full simulation, prediction of these flow fields can only be obtained by solving some suitably-averaged set of governing equations. In the traditional Reynolds-averaged approach, the Navier-Stokes equations are averaged over time. This in turn yields correlations between various turbulence fluctuations. It is these terms, e.g. the Reynolds stresses, for which a turbulence model must be derived. Turbulence modeling of incompressible flows has received a great amount of attention in the literature. An area of research that has received comparatively less attention is the modeling of compressible turbulent flows. An approach to simulating compressible turbulence at high Reynolds numbers is through the use of Large-Eddy Simulation (LES). In LES the dependent variables are decomposed into a large-scale (resolved) component and a sub-grid scale component. It is the small-scale components of the velocity field which are presumably more homogeneous than the large scales and, therefore, more easily modeled. Thus, it seems plausible that simpler models, which should be more universal in character than those employed in second-order closure schemes, may be developed for LES of compressible turbulence. The objective of the present research, therefore, is to explore models for the Large-Eddy Simulation of compressible turbulent flows. Given the recent successes of Zeman in second order closure modeling of compressible turbulence, model development was guided by principals employed in second-order closures.

Squires, Kyle D.↗
NASA NTRS1995

Dynamics and structure of turbulent premixed flames

In earlier work (Mantel & Bilger, 1994) the structure of the turbulent premixed flame was investigated using statistics based on conditional averaging with the reaction progress variable as the conditioning variable. The DNS data base of Trouve and Poinsot (1994) was used in this investigation. Attention was focused on the conditional dissipation and conditional axial velocity in the flame with a view to modeling these quantities for use in the conditional moment closure (CMC) approach to analysis of kinetics in premixed flames (Bilger, 1993). Two remarkable findings were made: there was almost no acceleration of the axial velocity in the flame front itself; and the conditional scalar dissipation remained as high, or higher, than that found in laminar premixed flames. The first finding was surprising since in laminar flames all the fluid acceleration occurs through the flame front, and this could be expected also for turbulent premixed flames at the flamelet limit. The finding gave hope of inventing a new approach to the dynamics of turbulent premixed flames through use of rapid distortion theory or an unsteady Bernoulli equation. This could lead to a new second order closure for turbulent premixed flames. The second finding was contrary to our measurements with laser diagnostics in lean hydrocarbon flames where it is found that conditional scalar dissipation drops dramatically below that for laminar flamelets when the turbulence intensity becomes high. Such behavior was not explainable with a one-step kinetic model, even at non-unity Lewis number. It could be due to depletion of H2 from the reaction zone by preferential diffusion. The capacity of the flame to generate radicals is critically dependent on the levels of H2 present (Bilger, et al., 1991). It seemed that a DNS computation with a multistep reduced mechanism would be worthwhile if a way could be found to make this feasible. Truly innovative approaches to complex problems often come only when there is the opportunity to work close at hand with the (in this case numerical) experimental data. Not only can one spot patterns and relationships in the data which could be important, but one can also get to know the limitations of the technique being used, so that when the next experiment is being designed it will address resolvable questions. A three-year grant from the Australian Research Council has enabled us to develop a small capability at the University of Sydney to work on DNS of turbulent reacting flow, and to analyze data bases generated at CTR. Collaboration between the University of Sydney and CTR is essential to this project and finding a workable modus operandum for this collaboration, given the constraints involved, has been a major objective of the past year's effort. The overall objectives of the project are: (1) to obtain a quantitative understanding of the dynamics of turbulent premixed flames at high turbulence levels with a view to developing improved second order closure models; and (2) to carry out new DNS experiments on turbulent premixed flames using a carefully chosen multistep reduced mechanism for the chemical kinetics, with a view to elucidating the laser diagnostic findings that are contrary to the findings for DNS using one-step kinetics. In this first year the objectives have been to make the existing CTR data base more accessible to coworkers at the University of Sydney, to make progress on understanding the dynamics of the flame in this existing CTR data base, and to carefully construct a suitable multistep reduced mechanism for use in a new set of DNS experiments on turbulent premixed flames.

Bilger, R. W.↗