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 73 records · Page 4

A comparison of turbulence closure models for transonic flows about airfoils

Navier-Stokes transonic airfoil calculations using three different turbulence closure models are presented and compared with available experimental data. Two of the models, Cebeci-Smith and Baldwin-Lomax, are equilibrium models in the sense that the turbulent shear stress is assumed to depend only on the local properties of the mean flow. The third model, Johnson and King, is a nonequilibrium model in which an ordinary differential equation is employed to account for the streamwise development of the maximum Reynolds shear stress. An eddy viscosity distribution across the boundary layer is assumed which is functionally dependent on this stress. For attached flows with little inviscid-viscous interaction, little difference in the results with the different models is noted. With stronger interactions and with separation occurring on the airfoils, the nonequilibrium model is shown to perform better.

King, Lyndell S.↗

The present state and future direction of second order closure models for compressible flows

The topics are presented in viewgraph form and include: (1) Reynolds stress closure models; (2) Favre averages and governing equations; (3) the model for the deviatoric part of the pressure-strain rate correlation; (4) the SSG pressure-strain correlation model; (5) a compressible turbulent dissipation rate model; (6) variable viscosity effects; (7) near-wall stiffness problems; (8) models of the Reynolds mass and heat flux; and (9) a numerical solution of the compressible turbulent transport equation.

Gatski, Thomas B.↗

ASSESSMENT OF MULTIPHASE CFD WITH ZERO CLOSURE MODEL FOR BOILING WATER REACTOR FUEL ASSEMBLIES

Extension of existing boiling closures to the high void fraction conditions of BWRs is one of challenges for M-CFD applications. This paper presents the assessment of an Eulerian M-CFD model with Zero closure model that has been specifically designed for two-phase flow and heat transfer applications in BWR fuel assemblies. The baseline Zero closure aims at including the simplest and most robust representation of the key physical mechanisms, rather than aiming at absolute accuracy of the predictions. The capabilities of M-CFD with Zero were evaluated against the international OECD/NRC BFBT benchmark data through the comparative analyses of local void fraction, sub-channel void fraction error and exit quality. INL high performance computing system has been employed to demonstrate the scalability of the M-CFD application to nuclear reactor fuel assemblies. The M-CFD simulation results show a good agreement with the experimental data for both quantitative and qualitative results, with the highest uncertainty being observed in the prediction of the local void fraction profile at low void fraction condition. This paper presents a recent view of the best practices for M-CFD of BWR fuel assemblies to provide a useful reference to identify and overcome common challenges in the applications of M-CFD tools. We provide recommendations for geometric modeling, mesh generation and solver configurations to optimize the methodology application.

97 MATHEMATICS AND COMPUTING↗

Assessment and application of Reynolds stress closure models to high-speed compressible flows

The paper presents results from the development of higher order closure models for the phenomological modeling of high-speed compressible flows. The work presented includes the introduction of an improved pressure-strain correlationi model applicable in both the low- and high-speed regime as well as modifications to the isotropic dissipation rate to account for dilatational effects. Finally, the question of stiffness commonly associated with the solution of two-equation and Reynolds stress transport equations in wall-bounded flows is examined and ways of relaxing these restrictions are discussed.

Gatski, T. B.↗

Design of a nonsingular Level 2.5 second-order closure model for the prediction of atmospheric turbulence

The suitability of applying the Mellor and Yamada (1974, 1982) Level 2.5 second-order turbulence closure model to general circulation models is investigated by examining not only the scheme's simulation of fully (or nearly fully) developed turbulence, but also its simulation of rapidly growing or strongly decaying turbulence. The behavior of the model is presented over its entire domain of definition, with special consideration given to the pathologies of the model. The model is then modified for the case of growing turbulence to rectify some of its physical shortcomings for that case, and to remove the pathologies that prohibit its use in a general circulation model. The performance of the modified Level 2.5 model is compared to the performance of various other modified versions through the numerical simulation for a growing convective PBL. The results show that the modified Level 2.5 model is a viable candidate for the prediction of turbulence and the simulation of the PBL in general circulation models.

Helfand, H. M.↗

Sensible-Enthalpy-Based Conditional Moment Closure Model for Homogeneous Charge Compression Ignition with Temperature Inhomogeneity

Abstract This paper presents a sensible-enthalpy-based conditional moment closure (CMC) model for homogeneous charge compression ignition (HCCI) with temperature inhomogeneity. The focus of study is placed on the effects of a conditioning variable on the prediction of heat release rates and developing the mixing models for sensible enthalpy, a reacting scalar. The sensible-enthalpy-based CMC model is validated a priori with 2-D direct numerical simulation (DNS) data for HCCI, and compared with the total-enthalpy-based formulation. Sensible enthalpy is shown to be strongly correlated with temperature, which makes sensible enthalpy well suited for a conditioning variable in CMC for HCCI. Results show that sensible-enthalpy-based CMC performs well for the prediction of heat release rates for all temperature stratification levels investigated in this study, while total-enthalpy-based CMC leads to the overprediction of heat release rates when the temperature stratification level is high. While conditional fluctuations of species mass fractions are substantial, sensible-enthalpy-based CMC is found to accurately estimate heat release rates during ignition, as well as reaction rates of key species, due to substantially suppressed conditional temperature fluctuations. The mapping closure models for the probability density function (PDF) and conditional dissipation rates of sensible enthalpy are proposed. The proposed models are found to capture the characteristics of reactive scalar mixing.

42 ENGINEERING↗

On Master-Length Scale Formulations for Stable Conditions in Turbulence Closure Models

Three formulations of the turbulence-length scales used in numerical modelling of atmospheric flows are compared. The comparison is made using the Mellor–Yamada–Nakanishi–Niino turbulence closure model within the stable boundary layer local similarity framework. With an appropriate choice of model constants, the model predictions are barely discernible and compare well with the empirical data obtained from the SHEBA campaign.

54 ENVIRONMENTAL SCIENCES↗

Roles of drizzle in a one-dimensional third-order turbulence closure model of the nocturnal stratus-topped marine boundary layer

This study focuses on the effects of drizzle in a one-dimensional third-order turbulence closure model of the nocturnal stratus-topped marine boundary layer. When the simulated drizzle rate is relatively small (maximum approximately equal to 0.6 mm/day), steady-state solutions are obtained. The boundary layer stabilizes essentially because drizzle causes evaporative cooling of the subcloud layer. This stabilization considerably reduces the buoyancy flux and turbulence kinetic energy below the stratus cloud. Thus, drizzle tends to decouple the cloud from the subcloud layer in the model, as suggested by many observational studies. In addition, the evaporation of drizzle in the subcloud layer creates small scattered clouds, which are likely to represent cumulus clouds, below the solid stratus cloud in the model. The sensitivity experiments show that these scattered clouds help maintain a coupled boundary layer. When the drizzle rate is relatively large (maximum approximately equal to 0.9 mm/day), the response of the model becomes transient with bursts in turbulent fluxes. This phenomenon is related to the formation of the scattered cloud layer below the solid stratus cloud. It appears that the model is inadequate to represent the heat and moisture transport by strong updrafts covering a small fractional area in cumulus convection.

Wang, Shouping↗

Higher order turbulence closure models

Theoretical models are developed and numerical studies conducted on various types of flows including both elliptic and parabolic. The purpose of this study is to find better higher order closure models for the computations of complex flows. This report summarizes three new achievements: (1) completion of the Reynolds-stress closure by developing a new pressure-strain correlation; (2) development of a parabolic code to compute jets and wakes; and, (3) application to a flow through a 180 deg turnaround duct by adopting a boundary fitted coordinate system. In the above mentioned models near-wall models are developed for pressure-strain correlation and third-moment, and incorporated into the transport equations. This addition improved the results considerably and is recommended for future computations. A new parabolic code to solve shear flows without coordinate tranformations is developed and incorporated in this study. This code uses the structure of the finite volume method to solve the governing equations implicitly. The code was validated with the experimental results available in the literature.

Amano, Ryoichi S.↗

Turbulence closure modeling near rigid boundaries

The near-wall region plays an essential role in turbulent boundary layers: it is a region of high shear; the peak rate of production and peak intensity of turbulence occurs there; and the peak rate of dissipation occurs right at the wall. Nevertheless, this region has received less attention from modelers than have more nearly homogeneous flows. One reason for this is that when the boundary layer is near equilibrium, experimental data can be used to prescribe the flow in the wall layer. Another reason is that most turbulence models are developed under assumptions of near homogeneity. This is a poor approximation in the wall region. A single-point moment closure model for the strongly non-homogeneous A turbulent flow near a rigid boundary is developed.

Durbin, Paul A.↗

Deep learning closure models for large-eddy simulation of flows around bluff bodies

Near-wall flow simulation remains a central challenge in aerodynamics modelling: Reynolds-averaged Navier–Stokes predictions of separated flows are often inaccurate, and large-eddy simulation (LES) can require prohibitively small near-wall mesh sizes. A deep learning (DL) closure model for LES is developed by introducing untrained neural networks into the governing equations and training in situ for incompressible flows around rectangular prisms at moderate Reynolds numbers. The DL-LES models are trained using adjoint partial differential equation (PDE) optimization methods to match, as closely as possible, direct numerical simulation (DNS) data. They are then evaluated out-of-sample – for aspect ratios, Reynolds numbers and bluff-body geometries not included in the training data – and compared with standard LES models. The DL-LES models outperform these models and are able to achieve accurate LES predictions on a relatively coarse mesh (downsampled from the DNS mesh by factors of four or eight in each Cartesian direction). We study the accuracy of the DL-LES model for predicting the drag coefficient, near-wall and far-field mean flow, and resolved Reynolds stress. A crucial challenge is that the LES quantities of interest are the steady-state flow statistics; for example, a time-averaged velocity component $\langle {u}_i\rangle (x) = \lim _{t \rightarrow \infty } ({1}/{t}) \int _0^t u_i(s,x)\, {\rm d}s$ . Calculating the steady-state flow statistics therefore requires simulating the DL-LES equations over a large number of flow times through the domain. It is a non-trivial question whether an unsteady PDE model with a functional form defined by a deep neural network can remain stable and accurate on $t \in [0, \infty )$ , especially when trained over comparatively short time intervals. Our results demonstrate that the DL-LES models are accurate and stable over long time horizons, which enables the estimation of the steady-state mean velocity, fluctuations and drag coefficient of turbulent flows around bluff bodies relevant to aerodynamics applications.

Mechanics↗

Analysis of Extensive Cross-Flow Separation using Higher-Order RANS Closure Models

The turbulent flow fields associated with the incompressible flow over a 6:1 prolate spheroid at high angle of attack, and the supersonic flow over an ogive cylinder are studied. Both these flows are characterized by large separation and vortical flow regions and therefore provide a challenging database for comparison of turbulent closure models. Of interest is the ability to predict the effects of separation and associated vortical motion common to both flows. Two turbulent models are investigated that each represent the class of linear eddy-viscosity models (LEVMs) and explicit algebraic stress models (EASMs). Since the EASM accounts for anisotropic effects, the in- fluence of these effects on flow field predictions can be assessed. The EASM model is shown to both improve the separation location prediction and pressure trough under the secondary vortex on the 6:1 prolate spheroid at high angle of attack and high Reynolds number, and improve the prediction of the separation location on a supersonic ogive cylinder.

Morison, J. H.↗

Analysis of Extensive Cross-Flow Separation using Higher-Order RANS Closure Models

The turbulent flow fields associated with the incompressible flow over a 6:1 prolate spheroid at high angle of attack, and the supersonic flow over an ogive cylinder are studied. Both these flows are characterized by large separation and vortical flow regions and therefore provide a challenging database for comparison of turbulent closure models. Of interest is the ability to predict the effects of separation and associated vortical motion common to both flows. Two turbulent models are investigated that each represent the class of linear eddy-viscosity models (LEVMs) and explicit algebraic stress models (EASMs). Since the EASM accounts for anisotropic effects, the influence of these effects on flow field predictions can be assessed. The EASM model is shown to both improve the separation location prediction and pressure trough under the secondary vortex on the 6:l prolate spheroid at high angle of attack and high Reynolds number, and improve the prediction of the separation location on a supersonic ogive cylinder.

Morrison, J. H.↗

Scientific machine learning for closure models in multiscale problems: A review

Here, closure problems are omnipresent when simulating multiscale systems, where some quantities and processes cannot be fully prescribed despite their effects on the simulation's accuracy. Recently, scientific machine learning approaches have been proposed as a way to tackle the closure problem, combining traditional (physics-based) modeling with data-driven (machine-learned) techniques, typically through enriching differential equations with neural networks. This paper reviews the different reduced model forms, distinguished by the degree to which they include known physics, and the different objectives of a priori and a posteriori learning. The importance of adhering to physical laws (such as symmetries and conservation laws) in choosing the reduced model form and choosing the learning method is discussed. The effect of spatial and temporal discretization and recent trends toward discretization-invariant models are reviewed. In addition, we make the connections between closure problems and several other research disciplines: inverse problems, Mori-Zwanzig theory, and multi-fidelity methods. In conclusion, much progress has been made with scientific machine learning approaches for solving closure problems, but many challenges remain. In particular, the generalizability and interpretability of learned models is a major issue that needs to be addressed further.

97 MATHEMATICS AND COMPUTING↗

The relationship between eddy-transport and second-order closure models for stratified media and for vortices

The question is considered of how complex a model should be used for the calculation of turbulent shear flows. At the present time there are models varying in complexity from very simple eddy-transport models to models in which all the equations for the nonzero second-order correlations are solved simultaneously with the equations for the mean variables. A discussion is presented of the relationship between these two models of turbulent shear flow. Two types of motion are discussed: first, turbulent shear flow in a stratified medium and, second, the motion in a turbulent line vortex. These two cases are instructive because in the first example eddy-transport methods have proven reasonably effective, whereas in the second, they have led to erroneous conclusions. It is not generally appreciated that the simplest form of eddy-transport theory can be derived from second-order closure models of turbulent flow by a suitably limiting process. This limiting process and the suitability of eddy-transport modeling for stratified media and line vortices are discussed.

Donaldson, C. D.↗

Neural network representability of fully ionized plasma fluid model closures

The closure problem in fluid modeling is a well-known challenge to modelers aiming to accurately describe their systems of interest. Over many years, analytic formulations in a wide range of regimes have been presented but a practical, generalized fluid closure for magnetized plasmas remains an elusive goal. Here, as a first step toward constructing a novel data-based approach to this problem, we apply ever-maturing machine learning methods to assess the capability of neural network architectures to reproduce crucial physics inherent in known magnetized plasma closures. We find encouraging results, indicating the applicability of neural networks to closure physics, but also arrive at recommendations on how one should choose appropriate network architectures for the given locality properties dictated by the underlying physics of the plasma.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Improvement of two-phase closure models in CTF using Bayesian inference

Under the Consortium for Advanced Simulation of Light Water Reactors (CASL) program, extensive capabilities have been developed in CTF to analyze light-water reactors (LWRs) for normal operating conditions, departure from nucleate boiling (DNB), and system transients. However, further improvements are required in the modeling and simulation of boiling water reactors (BWRs), which is a focus of the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program. In this work, CTF validation results were used to optimize selected modeling coefficients by calibrating to experimental data using a Bayesian inference approach. Here, calibration studies were conducted to improve (vapor) void fraction prediction without worsening the two-phase pressure drop prediction, as well as to improve the two-phase pressure drop prediction. Calibration was performed for interfacial drag and wall shear models. Surrogates were developed to alleviate the computational expense required for sampling the parameter space using Markov chain Monte Carlo (MCMC). An assessment performed with calibrated models demonstrated an improvement of CTF in its prediction of key parameters such as void fraction and two-phase pressure drop.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Study of the Alfven eigenmodes stability in CFQS plasma using a Landau closure model

The aim of this study is to analyze the stability of the Alfven eigenmodes (AE) in the Chinese First Quasi-axisymmetric Stellarator (CFQS). The AE stability is calculated using the code FAR3d that solves the reduced MHD equations to describe the linear evolution of the poloidal flux and the toroidal component of the vorticity in a full 3D system, coupled with equations of density and parallel velocity moment for the energetic particles (EP) species including the effect of the helical couplings and acoustic modes. The Landau damping and resonant destabilization effects are added in the model by a given closure relation. The simulation results indicate the destabilization of n = 1 to 4 AEs by EP during the slowing down process, particularly n = 1 and n = 2 toroidal AEs (TAE), n = 3 elliptical AE (EAE) and n = 4 non circular AE (NAE). If the resonance is caused by EPs with an energy above 17 keV (weakly thermalized EP), n = 2 EAEs and n = 3 NAEs are unstable. On the other hand, EPs with an energy below 17 keV (late thermalization stage) lead to the destabilization of n = 3 and n = 4 TAEs. The simulations for an off-axis NBI injection indicate the further destabilization of n = 2 to 4 AEs although the growth rate of the n = 1 AEs slightly decreases, so no clear optimization trend with respect to the NBI deposition region is identified. In addition, n = 2, 4 helical AE (HAE) are unstable above an EP β threshold. Also, if the thermal β of the simulation increases (higher thermal plasma density) the AE stability of the plasma improves. Overall, the simulations including the effect of the finite Larmor radius and electron-ion Landau damping show the stabilization of the n = 1 to 4 EAE/NAEs as well as a decrease of the growth rate and frequency of the n = 1 to 4 BAE/TAEs.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗