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 559 records · Page 31
All sourcesNASA NTRSDOE OSTI
NASA NTRS1992

The pressure-dilatation correlation in compressible flows

Simulations of simple compressible flows have been performed to enable the direct estimation of the pressure-dilatation correlation. The generally accepted belief that this correlation may be important in high-speed flows has been verified by the simulations. The pressure-dilatation correlation is theoretically investigated by considering the equation for fluctuating pressure in an arbitrary compressible flow. This leads to the isolation of a component of the pressure-dilatation that exhibits temporal oscillations on a fast time scale. Direct numerical simulations of homogeneous shear turbulence and isotropic turbulence show that this fast component has a negligible contribution to the evolution of turbulent kinetic energy. Then, an analysis for the case of homogeneous turbulence is performed to obtain a formal solution for the nonoscillatory pressure-dilatation. Simplifications lead to a model that algebraically relates the pressure-dilatation to quantities traditionally obtained in incompressible turbulence closures. The model is validated by direct comparison with the simulations.

Sarkar, S.↗
DOE OSTI2014

Theory and Implementation of Nuclear Safety System Codes - Part II: System Code Closure Relations, Validation, and Limitations

This is Part II of two articles describing the details of thermal-hydraulic sys- tem codes. In this second part of the article series, the system code closure relationships (used to model thermal and mechanical non-equilibrium and the coupling of the phases) for the governing equations are discussed and evaluated. These include several thermal and hydraulic models, such as heat transfer coefficients for various flow regimes, two phase pressure correlations, two phase friction correlations, drag coefficients and interfacial models be- tween the fields. These models are often developed from experimental data. The experiment conditions should be understood to evaluate the efficacy of the closure models. Code verification and validation, including Separate Effects Tests (SETs) and Integral effects tests (IETs) is also assessed. It can be shown from the assessments that the test cases cover a significant section of the system code capabilities, but some of the more advanced reactor designs will push the limits of validation for the codes. Lastly, the limitations of the codes are discussed by considering next generation power plants, such as Small Modular Reactors (SMRs), analyz- ing not only existing nuclear power plants, but also next generation nuclear power plants. The nuclear industry is developing new, innovative reactor designs, such as Small Modular Reactors (SMRs), High-Temperature Gas-cooled Reactors (HTGRs) and others. Sub-types of these reactor designs utilize pebbles, prismatic graphite moderators, helical steam generators, in- novative fuel types, and many other design features that may not be fully analyzed by current system codes. This second part completes the series on the comparison and evaluation of the selected reactor system codes by discussing the closure relations, val- idation and limitations. These two articles indicate areas where the models can be improved to adequately address issues with new reactor design and development.

99 GENERAL AND MISCELLANEOUS↗
NASA NTRS2011

Two-Point Turbulence Closure Applied to Variable Resolution Modeling

Variable resolution methods have become frontline CFD tools, but in order to take full advantage of this promising new technology, more formal theoretical development is desirable. Two general classes of variable resolution methods can be identified: hybrid or zonal methods in which RANS and LES models are solved in different flow regions, and bridging or seamless models which interpolate smoothly between RANS and LES. This paper considers the formulation of bridging methods using methods of two-point closure theory. The fundamental problem is to derive a subgrid two-equation model. We compare and reconcile two different approaches to this goal: the Partially Integrated Transport Model, and the Partially Averaged Navier-Stokes method.

Girimaji, Sharath S.↗
DOE OSTI2020

An improved ten-moment closure for reconnection and instabilities

The integration of kinetic effects in fluid models is important for global simulations of Earth's magnetosphere. The use of the two-fluid ten-moment model, which includes the pressure tensor for both species, has had some success in simulating Ganymede and Mercury with a simple closure model. Here, we discuss a heat flux closure which accounts for some limitations of the earlier work while remaining computationally tractable. Comparisons with kinetic simulations for magnetic reconnection and lower-hybrid drift instabilities show good agreement with kinetic results and improvements on previous closure models.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1992

Assessment of closure coefficients for compressible-flow turbulence models

A critical assessment is made of the closure coefficients used for turbulence length scale in existing models of the transport equation, with reference to the extension of these models to compressible flow. It is shown that to satisfy the compressible 'law of the wall', the model coefficients must actually be functions of density gradients. The magnitude of the errors that result from neglecting this dependence on density varies with the variable used to specify the length scale. Among the models investigated, the k-omega model yields the best performance, although it is not completely free from errors associated with density terms. Models designed to reduce the density-gradient effect to an insignificant level are proposed.

Huang, P. G.↗
DOE OSTI2021

A practical extension of the recursive multi-fidelity model for the emulation of hole closure experiments

We report in regimes of high strain rate, the strength of materials often cannot be measured directly in experiments. Instead, the strength is inferred based on an experimental observable, such as a change in shape, that is matched by simulations supported by a known strength model. In hole closure experiments, the rate and degree to which a central hole in a plate of material closes during a dynamic loading event are used to infer material strength parameters. Due to the complexity of the experiment, many computationally expensive, three-dimensional simulations are necessary to train an emulator for calibration or other analyses. These simulations can be run at multiple grid resolutions, where dense grids are slower but more accurate. In an effort to reduce the computational cost, a combination of simulations with different resolutions can be combined to develop an accurate emulator within a limited training time. We explore the novel design and construction of an appropriate functional recursive multi-fidelity emulator of a strength model for tantalum in hole closure experiments that can be applied to arbitrarily large training data. Hence, by formulating a multi-fidelity model to employ low-fidelity simulations, we were able to reduce the error of our emulator by approximately 81% with only an approximately 1.6% increase in computing resource utilization.

36 MATERIALS SCIENCE↗
DOE OSTI2018

CTF_

CTF is a thermal-hydraulic simulation tool designed for modeling the reactor vessel region of light water reactors. The code is capable of modeling both coolant and solid regions and includes a nuclear fuel rod model. Coolant is represented with a 9-equation (two phases plus droplets and non-condensable gas) model with closure relationships being provided for normal and accident conditions. Meshing is done using the subchannel approach, which models the coolant region as a collection of parallel 1-D flow paths that communicate laterally through "gaps". CTF also includes a collection of supporting software tools that include: Xml2ctf - A utility for creating native CTF input decks from a simplified input file; pyCTF - A collection of Python utilities used for post-processing CTF results, manipulating CTF input models, and testing differences between output files; CTFFuel - A separate front-end to the CTF fuel rod solver package; Multistate - A separate driver for CTF that communicates through the code coupling interface; and veraAPI - A c++ interface to the code.

Salko, Robert K↗
DOE OSTI2026

A Hybrid Biophysical‐Machine Learning Framework for Diurnal Surface Energy Flux Estimation Using Proximal Sensing

Thermal infrared-based remote sensing of surface energy fluxes has traditionally relied on high spatial resolution satellite data with revisit frequencies on the order of weeks. In this study, we evaluate a biophysics-based analytical surface energy balance model for predicting latent energy ( LE ) and sensible heat ( H ) fluxes using proximal sensing observations. The Surface Temperature Initiated Closure (STIC1.2) model has been extensively validated across a wide range of spatial and temporal scales using various satellite-derived thermal infrared data sets. Here we extend this validation by applying STIC at sub-hourly temporal resolution over multiple growing seasons for four distinct agricultural systems. We further develop and evaluate novel STIC variants that incorporate machine learning (ML) techniques to eliminate the need for surface energy balance observations, specifically net radiation and soil heat flux, thereby enhancing model applicability in data-sparse settings. The integration of a ML component to estimate surface available energy is shown to have strong predictive performance for both LE (R 2 = 0.81–0.94) and H (R 2 = 0.46–0.72) across all agricultural systems examined here, demonstrating the potential of hybrid biophysical-machine learning approaches for surface energy balance modeling with minimal data requirements. This study concludes with a novel application of explainable machine learning (exML) to diagnose sources of model error. This exML framework attributes residual prediction errors to both model input variables and environmental drivers not explicitly included in the simulation experiments. This approach provides a new pathway for improving model design and integrating previously overlooked yet influential variables into future model iterations.

evapotranspiration↗
DOE OSTI2026

Data for: A hybrid biophysical-machine learning framework for diurnal surface energy flux estimation using proximal sensing

Thermal-based remote sensing of surface energy fluxes has traditionally relied on high spatial resolution satellite data with revisit frequencies on the order of weeks. In this study, we evaluate a biophysics-based analytical surface energy balance model for predicting latent energy (LE) and sensible heat (H) fluxes using proximal sensing observations. The Surface Temperature Initiated Closure (STIC1.2) model has been extensively validated across a wide range of spatial and temporal scales using various satellite-derived thermal datasets. Here we extend this validation by applying STIC at sub-hourly temporal resolution over multiple growing seasons for four distinct agricultural systems. We further develop and evaluate novel STIC variants that incorporate machine learning (ML) techniques to eliminate the need for specific surface energy balance observations, specifically net radiation and soil heat flux, thereby enhancing model applicability in data-sparse settings. The integration of an ML component to estimate surface available energy is shown to have strong predictive performance for both LE (R2 = 0.81-0.94) and H (R2 = 0.46-0.72) across all agricultural systems examined here, demonstrating the potential of hybrid biophysical – machine learning approaches for surface energy balance modeling with minimal data requirements. This study concludes with a novel application of explainable machine learning (exML) to diagnose sources of model error. This exML framework attributes residual prediction errors to both model input variables and environmental drivers not explicitly included in the simulation experiments. This approach provides a new pathway for improving model design and integrating previously overlooked yet influential variables into future model iterations.

Agricultural Sciences↗
NASA NTRS1994

Single point modeling of rotating turbulent flows

A model for the effects of rotation on turbulence is proposed and tested. These effects which influence mainly the rate of turbulence decay are modeled in a modified turbulent energy dissipation rate equation that has explicit dependence on the mean rotation rate. An appropriate definition of the rotation rate derived from critical point theory and based on the invariants of the deformation tensor is proposed. The modeled dissipation rate equation is numerically well behaved and can be used in conjunction with any level of turbulence closure. The model is applied to the two-equation kappa-epsilon turbulence model and is used to compute separated flows in a backward-facing step and an axisymmetric swirling coaxial jets into a sudden expansion. In general, the rotation modified dissipation rate model shows some improvements over the standard kappa-epsilon model.

Hadid, A. H.↗
DOE OSTI2018

Nonadiabatic Flamelet Formulation for Predicting Wall Heat Transfer in Rocket Engines

A flamelet-based combustion model is introduced for the prediction of wall heat transfer in rocket engines and confined combustion systems. To account for convective heat loss due to the interaction of the flame with the wall, a permeable thermal boundary condition is introduced in the counterflow diffusion flame configuration. The solution of the resulting nonadiabatic flame structure forms a three-dimensional manifold, which is parameterized in terms of mixture fraction, progress variable, and temperature. The performance of the model is first assessed through a direct numerical simulation analysis of an H 2 /O 2 diffusion flame that is stabilized at an inert isothermal wall. The developed nonadiabatic flamelet model is shown to accurately predict the temperature, chemical composition, and wall heat transfer. Combined with a presumed probability density function closure, the model is then applied to large-eddy simulation of a single-injector rocket combustor to examine effects of heat transfer on the turbulent flame structure in rocket engines.

42 ENGINEERING↗
DOE OSTI2024

A robust spectral element implementation of the $k - τ$ RANS model in Nek5000/NekRS

The $k - ω$ Reynolds Averaged Navier Stokes (RANS) model is one of the industry standard approaches for modeling of turbulent flows. It performs better than the $k - ϵ$ model for low Reynolds number flows and is also more suitable for boundary layers with adverse pressure gradients. Major drawback of the model, however, is that the asymptotic value of $ω$ at the walls is singular, necessitating the use of a contrived “sufficiently” large value for $ω$ as the boundary condition for its transport equation. Here, this invariably leads to the solution being sensitive to near wall grid spacing. While an acceptable solution for low order (finite volume) methods, the excessive near wall gradients lead to persistent numerical stability issues in high order codes. To alleviate the problem, specifically in the context of the high order spectral element code Nek5000, a regularized $k - ω$ approach was formulated in our prior work (Tomboulides et al., 2018). The formulation, however, relies on the use of wall distance and its gradients for modeling the closure terms and can pose problems for simulations in complex geometries. This work presents a novel implementation of the $k - τ$ RANS model in Nek5000, where $τ = 1/ω$, eliminating the need for regularization, owing to the asymptotically bounded behavior of the source terms in the $τ$ transport equation, and also eliminating dependence on wall distance. Robustness and stability of the $k - τ$ model is ensured through implicit treatment of the source terms and their careful numerical implementation and demonstrated through several cases aimed at verification and validation. Studies include both canonical and engineering relevant problems, viz., turbulent channel flow, pipe flow, backward facing step, flow over NACA 0012 airfoil and flow in a T-junction. Results from the $k - τ$ model are shown to be consistent with regularized $k - ω$ model and also with the $k - ω$ SST model in OpenFOAM (for select studies). Comparison with experimental data is also shown, where available, to bolster validation efforts for the $k - τ$ model implementation through prediction of key turbulent quantities of interest.

Nek5000↗
DOE OSTI2020

Assessing Global and Local Radiative Feedbacks Based on AGCM Simulations for 1980–2014/2017

In order to avoid contamination of unobservable and uncertain effective radiative forcing (ERF) on the diagnosis of radiative feedbacks based on short-term climate variability, we verify the Kernel-Gregory feedback calculation method using atmospheric model experiments with prescribed ERF. We show that both clear-sky radiative fluxes and all-sky radiative feedbacks have a closure between model simulations and kernel derivations. A near-zero global mean net cloud feedback is found in the simulation with prescribed ERF, which results in a more negative global net climate feedback, -2 W m -2 K -1 . Consistent with AMIP6 ensemble mean results, the lapse rate feedback is the largest contributor among all feedbacks to temperature amplification over the three poles (Arctic, Antarctic and Tibetan Plateau), followed by surface albedo feedback and Planck feedback deviation from its global mean. Except for the higher surface albedo feedback in the Antarctic, other feedbacks are almost same between the Arctic and Antarctic.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1991

Modeling of crack bridging in a unidirectional metal matrix composite

The effective fatigue crack driving force and crack opening profiles were determined analytically for fatigue tested unidirectional composite specimens exhibiting fiber bridging. The crack closure pressure due to bridging was modeled using two approaches; the fiber pressure model and the shear lag model. For both closure models, the Bueckner weight function method and the finite element method were used to calculate crack opening displacements and the crack driving force. The predicted near crack tip opening profile agreed well with the experimentally measured profiles for single edge notch SCS-6/Ti-15-3 metal matrix composite specimens. The numerically determined effective crack driving force, Delta K(sup eff), was calculated using both models to correlate the measure crack growth rate in the composite. The calculated Delta K(sup eff) from both models accounted for the crack bridging by showing a good agreement between the measured fatigue crack growth rates of the bridged composite and that of unreinforced, unbridged titanium matrix alloy specimens.

Ghosn, Louis J.↗
NASA NTRS1992

Modeling of crack bridging in a unidirectional metal matrix composite

The effective fatigue crack driving force and crack opening profiles were determined analytically for fatigue tested unidirectional composite specimens exhibiting fiber bridging. The crack closure pressure due to bridging was modeled using two approaches: the fiber pressure model and the shear lag model. For both closure models, the Bueckner weight function method and the finite element method were used to calculate crack opening displacements and the crack driving force. The predicted near crack tip opening profile agreed well with the experimentally measured profiles for single edge notch SCS-6/Ti-15-3 metal matrix composite specimens. The numerically determined effective crack driving force, Delta K(eff), was calculated using both models to correlate the measure crack growth rate in the composite. The calculated Delta K(eff) from both models accounted for the crack bridging by showing a good agreement between the measured fatigue crack growth rates of the bridged composite and that of unreinforced, unbridged titanium matrix alloy specimens.

Ghosn, Louis J.↗
DOE OSTI2015

Multivariate moment closure techniques for stochastic kinetic models

Stochastic effects dominate many chemical and biochemical processes. Their analysis, however, can be computationally prohibitively expensive and a range of approximation schemes have been proposed to lighten the computational burden. These, notably the increasingly popular linear noise approximation and the more general moment expansion methods, perform well for many dynamical regimes, especially linear systems. At higher levels of nonlinearity, it comes to an interplay between the nonlinearities and the stochastic dynamics, which is much harder to capture correctly by such approximations to the true stochastic processes. Moment-closure approaches promise to address this problem by capturing higher-order terms of the temporally evolving probability distribution. Here, we develop a set of multivariate moment-closures that allows us to describe the stochastic dynamics of nonlinear systems. Multivariate closure captures the way that correlations between different molecular species, induced by the reaction dynamics, interact with stochastic effects. We use multivariate Gaussian, gamma, and lognormal closure and illustrate their use in the context of two models that have proved challenging to the previous attempts at approximating stochastic dynamics: oscillations in p53 and Hes1. In addition, we consider a larger system, Erk-mediated mitogen-activated protein kinases signalling, where conventional stochastic simulation approaches incur unacceptably high computational costs.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1992

The present state of DIA models and their impact on one point closures

An outline is presented of the DIA model. Another outline is presented of TSDIA (two scale DIA) and some suggestions for turbulence modeling. Two proposals are made: helicity for the study of the effects of swirling flow and cross flow; and density variance for the study of highly compressed flows.

Yoshizawa, Akira↗
DOE OSTI2019

Dynamics of potential vorticity staircase evolution and step mergers in a reduced model of beta-plane turbulence

A two-field model of potential vorticity (PV) staircase structure and dynamics relevant to both beta-plane and drift-wave plasma turbulence is studied numerically and analytically. The model evolves averaged PV whose flux is both driven by and regulates a potential enstrophy field, ε. The model employs a closure using a mixing length model. Its link to bistability, vital to staircase generation, is analyzed and verified by integrating the equations numerically. Long-time staircase evolution consistently manifests a pattern of metastable quasiperiodic configurations, lasting for hundreds of time units, yet interspersed with abrupt ( Δt << 1 ) mergers of adjacent steps in the staircase. The mergers occur at the staircase lattice defects where the pattern has not completely relaxed to a strictly periodic solution that can be obtained analytically. Other types of stationary solutions are solitons and kinks in the PV gradient and ε-profiles. The waiting time between mergers increases strongly as the number of steps in the staircase decreases. This is because of an exponential decrease in interstep coupling strength with growing spacing. The long-time staircase dynamics is shown numerically be determined by local interaction with adjacent steps. In conclusion, mergers reveal themselves through the explosive growth of the turbulent PV flux, which, however, abruptly drops to its global constant value once the merger is completed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗