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 379 records · Page 21
All sourcesNASA NTRSDOE OSTI
NASA NTRS1976

Calculation methods for compressible turbulent boundary layers

Calculation procedures for non-reacting compressible two- and three-dimensional turbulent boundary layers were reviewed. Integral, transformation and correlation methods, as well as finite difference solutions of the complete boundary layer equations summarized. Alternative numerical solution procedures were examined, and both mean field and mean turbulence field closure models were considered. Physics and related calculation problems peculiar to compressible turbulent boundary layers are described. A catalog of available solution procedures of the finite difference, finite element, and method of weighted residuals genre is included. Influence of compressibility, low Reynolds number, wall blowing, and pressure gradient upon mean field closure constants are reported.

Bushnell, D. M.↗
NASA NTRS2001

An Improved Heat Budget Estimation Including Bottom Effects for General Ocean Circulation Models

This paper studies the effects of the underwater light field on heat-budget calculations of general ocean circulation models for shallow waters. The presence of a bottom significantly alters the estimated heat budget in shallow waters, which affects the corresponding thermal stratification and hence modifies the circulation. Based on the data collected during the COBOP field experiment near the Bahamas, we have used a one-dimensional turbulence closure model to show the influence of the bottom reflection and absorption on the sea surface temperature field. The water depth has an almost one-to-one correlation with the temperature rise. Effects of varying the bottom albedo by replacing the sea grass bed with a coral sand bottom, also has an appreciable effect on the heat budget of the shallow regions. We believe that the differences in the heat budget for the shallow areas will have an influence on the local circulation processes and especially on the evaporative and long-wave heat losses for these areas. The ultimate effects on humidity and cloudiness of the region are expected to be significant as well.

Carder, Kendall↗
DOE OSTI2019

Explicit Filtering and Reconstruction to Reduce Grid Dependence in Convective Boundary Layer Simulations Using WRF-LES

As model grid resolutions move from the mesoscale to the microscale, turbulent structures represented in atmospheric boundary layer simulations change dramatically. At intermediate resolutions, the so-called gray zone, turbulent motions are not resolved accurately, posing a challenge to numerical simulations. The representation of turbulence is also highly sensitive to the choice of closure model. Here, we examine explicit filtering and reconstruction in the gray zone as a technique to better represent atmospheric turbulence. The convective boundary layer is simulated using the Weather Research and Forecasting (WRF) Model with horizontal resolutions ranging from 25 m to 1 km. Four large-eddy simulation (LES) turbulence models are considered: the Smagorinsky model, the TKE-1.5 model, and two versions of the dynamic reconstruction model (DRM). The models are evaluated by their ability to produce consistent mean potential temperature profiles, heat and momentum fluxes, velocity fields, and turbulent kinetic energy spectra as the grids become coarser. The DRM, a mixed model that uses an explicit filtering and reconstruction technique to account for resolvable subfilter-scale (RSFS) stresses, performs very well at resolutions of 500 m and 1 km without any special tuning, whereas the Smagorinsky and TKE-1.5 models produce heavily grid-dependent results.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1985

Tone excited jets. V - A theoretical model and comparison with experiment

A mathematical model of tone-excited jets is developed which consists of two major components: a mathematical description of the process by which the intrinsic instability waves of the jet are excited by the upstream tones; and the modeling of the nonlinear interaction between the mean flow of the jet, the excited large-scale instability waves or turbulence structure, and the fine scale turbulence. It is assumed that each of these jet flow components can be characterized by a few parameters, which are then related by conservation equations that are supplemented by closure models. This quasi-linear model's results are compared with experimental measurements, and good agreement is obtained over a wide range of excitation frequencies and excitation levels.

Tam, C. K. W.↗
DOE OSTI2022

Operator inference for non-intrusive model reduction with quadratic manifolds

This paper proposes a novel approach for learning a data-driven quadratic manifold from high-dimensional data, then employing this quadratic manifold to derive efficient physics-based reduced-order models. The key ingredient of the approach is a polynomial mapping between high-dimensional states and a low-dimensional embedding. This mapping consists of two parts: a representation in a linear subspace (computed in this work using the proper orthogonal decomposition) and a quadratic component. The approach can be viewed as a form of data-driven closure modeling, since the quadratic component introduces directions into the approximation that lie in the orthogonal complement of the linear subspace, but without introducing any additional degrees of freedom to the low-dimensional representation. Combining the quadratic manifold approximation with the operator inference method for projection-based model reduction leads to a scalable non-intrusive approach for learning reduced-order models of dynamical systems. Further, applying the new approach to transport-dominated systems of partial differential equations illustrates the gains in efficiency that can be achieved over approximation in a linear subspace.

42 ENGINEERING↗
DOE OSTI2022

CTF Improved Drag Model and Flow Regime Transition Criteria

The demand for accurate prediction of two-phase flow behavior in a boiling water reactor (BWR) requires a comprehensive understanding of flow regime, void fraction, heat transfer, and pressure drop. The CTF subchannel code, which is used for the Thermal/Hydraulic (T/H) solution in the Consortium for Advanced Simulation of Light Water Reactors (CASL)-developed Virtual Environment for Reactor Application (VERA) core simulator, is being further developed for BWR applications. In support of this goal, the present work highlights some of the two-phase closure model developments towards improving the CTF void fraction prediction, especially for subcooled boiling. The drift-flux approach has been well-developed for upward dispersed two-phase flows and proven to be accurate in predicting void fraction in bubbly and slug flow regimes. In this work, these kinematic constitutive relations for the drift-flux velocity have been implemented into CTF to describe the interfacial drag of bubbly flow as an alternative to the existing model for better void fraction prediction. The success of these constitutive relations also relies on a good flow regime map that accounts for flow conditions and channel geometry. A more reliable flow regime transition criteria that account for the flow condition has also been implemented in this study for modeling the flow regime transition criteria. The newly implemented models are shown to give improved void fraction predictions in comparison to experimental data.

Hizoum, Belgacem↗
DOE OSTI2025

2025 Advances in NekRS: Supporting improved performance for nuclear applications

This report presents several 2025 advancements in NekRS, a high-fidelity spectral element CFD code developed at Argonne National Laboratory to support the NEAMS thermal-hydraulics program. The forthcoming v25 release consolidates several of these advances, adding new features for portability across heterogeneous GPU architectures, real-time in situ visualization, improved turbulence modeling, and conjugate heat transfer coupling. Over the past year, NekRS has demonstrated strong scalability and performance on DOE’s leading exascale platforms, including Aurora and Frontier, confirming its readiness for some of the largest and most complex simulations attempted to date. These achievements provide a powerful new platform for high-fidelity data generation, which in turn supports the development and validation of advanced closure models critical for reactor safety and design. Significant algorithmic innovations have also been introduced. A new global runtime h-refinement capability simplifies workflows by reducing mesh preparation burdens and enabling coarse-to-fine restarts. Building on this, a novel multigrid strategy was implemented to accelerate pressure and transport solves at scale, addressing long-standing bottlenecks in exascale CFD. Together, these developments improve both the efficiency and accessibility of high-fidelity simulations for reactor-relevant problems. Collectively, these enhancements represent a major step forward in simulation technology, positioning NekRS as a cornerstone of NEAMS efforts to enable accurate, efficient, and scalable high-fidelity analysis of advanced nuclear systems.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
NASA NTRSDate not supplied

Recent Advancements in the PATO Material Response Code

Introduction: Predicting the complicated multiphysics phenomena during atmospheric entry requires high-fidelity modeling tools to refine estimates of mission risks during entry. To this end, new capabilities are being added to the Porous-material Analysis Toolbox based on OpenFOAM (PATO). PATO is an open-source software for Computational Material Response (CMR) of reactive porous materials submitted to high-temperature environments. The objective of this work is to highlight current efforts to add to and improve upon the modeling capabilities of PATO. These include efforts to loosely couple PATO with other discipline specialized codes including hypersonic Computational Fluid Dynamics (CFD), to assess the interaction effects between pyrolysis gas blowing and the boundary layer, and Computational Solid Mechanics (CSM), to address modeling of mechanical erosion. Other refinements include surface phenomena modeling capabilities to address the effects of silicone-based coatings applied to the TPS during flight preparation, and a unified multiphase solver for a mixed porous-material and plain-fluid domain. Coupling CMR with CFD (CMR/CFD): A loose coupling between PATO and the Data Parallel Line Relaxation (DPLR) CFD code has been achieved by making use of a blowing boundary condition at the heatshield surface available in DPLR. Starting with heat flux estimates with no pyrolysis gas blowing at the surface, blowing gases are computed by the CMR and passed to the CFD such that aerothermal properties of the environment can be recomputed for a new CMR computation. This leads to an iterative process which is supplemented with an estimate of the radiative heat flux using the Nonequilibrium air radiation (NEQAIR) program. The entire iterative process is illustrated in Figure 1. This coupling strategy has been utilized in computing the MSL material response. The goal is to compare the coupled CMR/CFD results with material response results obtained using traditional blowing corrections. Coupling CMS with CMR: A mechanical erosion model is currently being implemented in PATO to account for the additional mass removal induced by high shear conditions. The modeling process at each timestep consists of updating the mechanical properties as a function of temperature and computing the stress tensor and displacement fields of the material. Then, a failure criteria model determines the regions in which the stress exceeds the ultimate strength values resulting in mesh movement to account for mass removal. This model allows the material response simulation to compute the recession due to both oxidation and shear-induced erosion. The model is demonstrated by computing material response of sphere-cone arc jet samples. Surface Modeling Capabilities: NuSil, a silicone-based coating, was sprayed onto the MSL and Mars 2020 heatshields to mitigate shedding of phenolic dust. To better understand the effects of the NuSil coating on the material response, a novel model has been implemented in PATO. In this model, the equilibrium of the charred NuSil surface is modeled as pure silica, and a constant offset, inspired by the classical spallation model, is added to the the char blowing rate and wall enthalpy to reproduce HyMETS experimental results. The model has also been used to estimate the 3D material response of the MSL heatshield. Unified Solver: In addition to the iterative loose coupling approach mentioned above, a multiphase unified solver is being developed to couple the environment (plain-fluid phase) and the porous-material phase. The solver is based on the volume averaged conservation of mass, momentum, and energy for the macroscale with closure models which include microscale effects through effective physicochemical properties. The unified solver has been used to compute flow through a porous plug and solve the Beavers and Joseph problem. Since the strong coupling between phases is inherent to this solver, modeling assumptions present in other coupling methods of material response are mitigated. This strategy also makes it feasible to capture the competition between surface and volume ablation in the same computational domain, which is usually not possible with other coupling approaches.

Material Response↗
NASA NTRS1996

Fatigue response of perforated titanium for application in laminar flow control

The room temperature tensile and fatigue response of non-perforated and perforated titanium for laminar flow control application was investigated both experimentally and analytically. Results showed that multiple perforations did not affect the tensile response, but did reduce the fatigue life. A two dimensional finite element stress analysis was used to determine that the stress fields from adjacent perforations did not influence one another. The stress fields around the holes did not overlap one another, allowing the materials to be modeled as a plate with a center hole. Fatigue life was predicted using an equivalent MW flow size approach to relate the experimental results to microstructural features of the titanium. Predictions using flaw sizes ranging from 1 to 15 microns correlated within a factor of 2 with the experimental results by using a flow stress of 260 MPa. By using two different flow stresses in the crack closure model and correcting for plasticity, the experimental results were bounded by the predictions for high applied stresses. Further analysis of the complex geometry of the perforations and the local material chemistry is needed to further understand the fatigue behavior of the perforated titanium.

Johnson, W. Steven↗
NASA NTRSDate not supplied

Recent Advancements in the PATO Material Response Code

Introduction: Predicting the complicated multiphysics phenomena during atmospheric entry requires high-fidelity modeling tools to refine estimates of mission risks during entry. To this end, new capabilities are being added to the Porous-material Analysis Toolbox based on OpenFOAM (PATO) [1,2,3]. PATO is an open-source software for Computational Material Response (CMR) of reactive porous materials submitted to high-temperature environments. The objective of this work is to highlight current efforts to add to and improve upon the modeling capabilities of PATO. These include efforts to loosely couple PATO with other discipline specialized codes including hypersonic Computational Fluid Dynamics (CFD), to assess the interaction effects between pyrolysis gas blowing and the boundary layer, and Computational Solid Mechanics (CSM), to address modeling of mechanical erosion. Other refinements include surface phenomena modeling capabilities to address the effects of silicone-based coatings applied to the TPS during flight preparation, and a unified multiphase solver for a mixed porous-material and plain-fluid domain. Coupling CMR with CFD (CMR/CFD): A loose coupling between PATO and the Data Parallel Line Relaxation (DPLR) [4] CFD code has been achieved by making use of a blowing boundary condition at the heatshield surface available in DPLR. Starting with heat flux estimates with no pyrolysis gas blowing at the surface, blowing gases are computed by the CMR and passed to the CFD such that aerothermal properties of the environment can be recomputed for a new CMR computation. This leads to an iterative process which is supplemented with an estimate of the radiative heat flux using the Nonequilibrium air radiation (NEQAIR) [5] program. The entire iterative process is illustrated in Figure 1. This coupling strategy has been utilized in computing the MSL material response. The goal is to compare the coupled CMR/CFD results with material response results obtained using traditional blowing corrections. Coupling CMS with CMR: A mechanical erosion model is currently being implemented in PATO to account for the additional mass removal induced by high shear conditions. The modeling process at each timestep consists of updating the mechanical properties as a function of temperature and computing the stress tensor and displacement fields of the material. Then, a failure criteria model determines the regions in which the stress exceeds the ultimate strength values resulting in mesh movement to account for mass removal. This model allows the material response simulation to compute the recession due to both oxidation and shear-induced erosion. The model is demonstrated by computing material response of sphere-cone arc jet samples. Surface Modeling Capabilities: NuSil, a silicone-based coating, was sprayed onto the MSL and Mars 2020 heatshields to mitigate shedding of phenolic dust. To better understand the effects of the NuSil coating on the material response, a novel model has been implemented in PATO. In this model, the equilibrium of the charred NuSil surface is modeled as pure silica, and a constant offset, inspired by the classical spallation model, is added to the the char blowing rate and wall enthalpy to reproduce HyMETS experimental results. The model has also been used to estimate the 3D material response of the MSL heatshield [6]. Unified Solver: In addition to the iterative loose coupling approach mentioned above, a multiphase unified solver is being developed to couple the environment (plain-fluid phase) and the porous-material phase. The solver is based on the volume averaged conservation of mass, momentum, and energy for the macroscale with closure models which include microscale effects through effective physicochemical properties. The unified solver has been used to compute flow through a porous plug and solve the Beavers and Joseph problem [7]. Since the strong coupling between phases is inherent to this solver, modeling assumptions present in other coupling methods of material response are mitigated. This strategy also makes it feasible to capture the competition between surface and volume ablation in the same computational domain, which is usually not possible with other coupling approaches.

Thermal Protection Systems↗
NASA NTRS1983

Shock noise features using the SCIPVIS code

The numerical computational supersonic jet plume model SCIPVIS was analyzed in order to determine its effectiveness for predicting the complicated flow structure associated with supersonic shock containing jet plumes. The analysis focused on assessing the applicability of using this code for predicting the broadband shock noise radiation. The analysis was carried out for several exit Mach number cases (including 1.00, 1.41, and 1.99) covering both underexpanded and overexpanded flow conditions. The three primary turbulence closure models incorporated in the SCIPVIS code were utilized for comparisons to measured plume static pressure distributions. Among other results, the prediction of plume static pressures for all three nozzle exit Mach numbers was found to perform exceptionally well when compared with experimental data for both underexpanded and overexpanded cases. Computations performed using this code for calculating the average shock wavelength for the 1.99 exit Mach number case were shown to provide an excellent prediction of the peak frequency for broadband shock noise. The spatial variation of the turbulent kinetic energy was found to be in qualitative agreement with wedge hot-film measurements.

Seiner, J. M.↗
NASA NTRS1993

Toward modeling wingtip vortices

Wingtip vortices are generated by lifting airfoils; their salient features are compactness and relatively slow rate of decay. The principal motivation for studying the far field evolution of wingtip vortices is the need to understand and predict the extent of the vortex influence during aircraft take-off or landing. On submarines a wingtip vortex ingested into a propeller can be a source of undesirable noise. The main objectives of this research are (1) to establish theoretical understanding of the principal mechanisms that govern the later (diffusive) stages of a turbulent vortex, (2) to develop a turbulence closure model representing the basic physical mechanisms that control the vortex diffusive stage, and further (3) to investigate coupling between the near and far field evolutions; in other words, to study the effect of initial conditions on the vortex lifetime and the ultimate state.

Zeman, O.↗
DOE OSTI2016

STAR-CCM+ Verification and Validation Plan

The commercial Computational Fluid Dynamics (CFD) code STAR-CCM+ provides general purpose finite volume method solutions for fluid dynamics and energy transport. This document defines plans for verification and validation (V&V) of the base code and models implemented within the code by the Consortium for Advanced Simulation of Light water reactors (CASL). The software quality assurance activities described herein are port of the overall software life cycle defined in the CASL Software Quality Assurance (SQA) Plan [Sieger, 2015]. STAR-CCM+ serves as the principal foundation for development of an advanced predictive multi-phase boiling simulation capability within CASL. The CASL Thermal Hydraulics Methods (THM) team develops advanced closure models required to describe the subgrid-resolution behavior of secondary fluids or fluid phases in multiphase boiling flows within the Eulerian-Eulerian framework of the code. These include wall heat partitioning models that describe the formation of vapor on the surface and the forces the define bubble/droplet dynamic motion. The CASL models are implemented as user coding or field functions within the general framework of the code. This report defines procedures and requirements for V&V of the multi-phase CFD capability developed by CASL THM. Results of V&V evaluations will be documented in a separate STAR-CCM+ V&V assessment report. This report is expected to be a living document and will be updated as additional validation cases are identified and adopted as part of the CASL THM V&V suite.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2016

Initial conditions and modeling for simulations of shock driven turbulent material mixing

Here, we focus on the simulation of shock-driven material mixing driven by flow instabilities and initial conditions (IC). Beyond complex multi-scale resolution issues of shocks and variable density turbulence, me must address the equally difficult problem of predicting flow transition promoted by energy deposited at the material interfacial layer during the shock interface interactions. Transition involves unsteady large-scale coherent-structure dynamics capturable by a large eddy simulation (LES) strategy, but not by an unsteady Reynolds-Averaged Navier–Stokes (URANS) approach based on developed equilibrium turbulence assumptions and single-point-closure modeling. On the engineering end of computations, such URANS with reduced 1D/2D dimensionality and coarser grids, tend to be preferred for faster turnaround in full-scale configurations.

36 MATERIALS SCIENCE↗
NASA NTRS1985

Analysis of turbulent underexpanded jets. II - Shock noise features using SCIPVIS

SCIPVIS, the computational model discussed by Dash et al. (1985), is assessed in predicting the complicated flow structure associated with shock-containing plumes. In addition, the analysis in this study examines this code's applicability as a basic part of a program for estimating broadband shock noise radiation. The results of this study show that excellent agreement exists between predicted and measured static pressure distributions for both underexpanded and overexpanded flow cases considered. Of the three turbulence closure models incorporated in the SCIPVIS code, the kW model of Spalding produces the most uniform agreement with measurement. The k-epsilon-2 model of Launder consistently overestimates plume spreading for supersonic jets with exit Mach numbers in the 1-2 range. Dash's (1983) k-epsilon-2-cc, compressibility-corrected version of Launder's model underestimates plume spreading. Good qualitative agreement was also obtained between the measured longitudinal turbulence intensity and that predicted by the code for the same trial case. Comparison of measured and predicted broadband shock noise spectrum peak values were found to be in excellent agreement. This utilized a variant of the Harper-Bourne and Fisher (1973) phase-array model: the effective shock spacing was reinterpreted as the value of the end of the plume potential core, determined herein by the SCIPVIS code.

Seiner, J. M.↗
NASA NTRS1981

Methods and models for predicting fatigue crack growth under random loading

Papers are presented in the volume summarizing the baseline data, methodology, procedures, and results of a round-robin analysis which was conducted to predict the fatigue crack growth in 2219-T851 aluminum center-cracked specimens subjected to flight loading in random cycle-by-cycle format. The objective of the analysis was to assess whether data from constant-amplitude fatigue crack growth tests on center-cracked specimens can be used to predict fatigue crack growth lives under random loading. The following approaches are discussed in detail: a root-mean-square approach, a crack-closure model, a multi-parameter yield zone model, and a load-interaction model.

Chang, J. B.↗
DOE OSTI2016

Investigation on numerical method of one-dimensional drift flux model for boiling two-phase flow

A numerical method for the solution of one-dimension drift-flux model has been developed. Correlations for heat and mass transfer and kinematic constitutive equation were used to close the set of the partial differential equations. The finite volume method based on the staggered grid and fully implicit backward Euler method were applied to discretize the spatial and time differential items, respectively. The discretization equations were solved iteratively in a same manner as in SIMPLE method for a single-phase flow. The verification and validation of the numerical method and closure model for one-dimension drift-flux model have been performed by comparing the code solution with the experimental data. Furthermore, the axial distribution of void fraction and pressure also have been investigated under different heated case. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
NASA NTRS1991

Progress in the development of PDF turbulence models for combustion

A combined Monte Carlo-computational fluid dynamic (CFD) algorithm was developed recently at Lewis Research Center (LeRC) for turbulent reacting flows. In this algorithm, conventional CFD schemes are employed to obtain the velocity field and other velocity related turbulent quantities, and a Monte Carlo scheme is used to solve the evolution equation for the probability density function (pdf) of species mass fraction and temperature. In combustion computations, the predictions of chemical reaction rates (the source terms in the species conservation equation) are poor if conventional turbulence modles are used. The main difficulty lies in the fact that the reaction rate is highly nonlinear, and the use of averaged temperature produces excessively large errors. Moment closure models for the source terms have attained only limited success. The probability density function (pdf) method seems to be the only alternative at the present time that uses local instantaneous values of the temperature, density, etc., in predicting chemical reaction rates, and thus may be the only viable approach for more accurate turbulent combustion calculations. Assumed pdf's are useful in simple problems; however, for more general combustion problems, the solution of an evolution equation for the pdf is necessary.

Hsu, Andrew T.↗