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.

SourceResearch connectionFrom yearThrough year

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

At least 289 records · Page 16

All sources NASA NTRS DOE OSTI

NASA NTRS1991

Investigations of turbulent scalar fields using probability density function approach

Scalar fields undergoing random advection have attracted much attention from researchers in both the theoretical and practical sectors. Research interest spans from the study of the small scale structures of turbulent scalar fields to the modeling and simulations of turbulent reacting flows. The probability density function (PDF) method is an effective tool in the study of turbulent scalar fields, especially for those which involve chemical reactions. It has been argued that a one-point, joint PDF approach is the one to choose from among many simulation and closure methods for turbulent combustion and chemically reacting flows based on its practical feasibility in the foreseeable future for multiple reactants. Instead of the multi-point PDF, the joint PDF of a scalar and its gradient which represents the roles of both scalar and scalar diffusion is introduced. A proper closure model for the molecular diffusion term in the PDF equation is investigated. Another direction in this research is to study the mapping closure method that has been recently proposed to deal with the PDF's in turbulent fields. This method seems to have captured the physics correctly when applied to diffusion problems. However, if the turbulent stretching is included, the amplitude mapping has to be supplemented by either adjusting the parameters representing turbulent stretching at each time step or by introducing the coordinate mapping. This technique is still under development and seems to be quite promising. The final objective of this project is to understand some fundamental properties of the turbulent scalar fields and to develop practical numerical schemes that are capable of handling turbulent reacting flows.

Gao, Feng↗

NASA NTRS1991

Turbulence modeling

The performance of existing two-equation eddy viscosity models was examined. An effort was made to develop better models for near-wall turbulence using direct numerical simulations of plane channel and boundary layer flows. The asymptotic near-wall behavior of turbulence was used to examine the problems of current second order closure models and develop new models with the correct near-wall behavior. Rapid Distortion Theory was used to analytically study the effects of mean deformation on turbulence, obtain analytical solutions for the spectrum tensor, Reynolds stress tensor, anisotropy tensor and its invariants, which can be used in the turbulence model development. The potential of the renormalization group theory in turbulence modeling was studied, as well as compressible turbulent flows, and modeling of bypass transition.

Shih, Tsan-Hsing↗

NASA NTRS1990

Proof test and fatigue crack growth modeling on 2024-T3 aluminum alloy

Pressure proof testing of aircraft fuselage structures has been suggested as a means of screening critical crack sizes and of extending their useful life. The objective of this paper is to study the proof-test concept and to model the crack-growth process on a ductile material. Simulated proof and operational fatigue life tests have been conducted on cracked panels made of 2024-T3 aluminum alloy sheet material. A fatigue crack-closure model was modified to simulate the proof test and operational fatigue cycling. Using crack-growth rate and resistance-curve data, the model was able to predict crack growth during and after the proof load. These tests and analyses indicate that the proof test increases fatigue life; but the beneficial life, after a 1.33 or 1.5 proof, was less than a few hundred cycles.

Newman, J. C., Jr.↗

DOE OSTI2023

High-Fidelity Modeling and Experiments to Inform Safety Analysis Codes for Heat Pipe Microreactors

Heat pipe microreactors are reactor designs that primarily use liquid-metal heat pipes to cool the core. The main interest in heat pipes is the fact that they can remove heat passively. This, along with the use of liquid metal, allows the reactor to operate at higher temperatures. Although the use of heat pipes in nuclear reactors is new, liquid-metal heat pipe technology is mature. Nevertheless, experimental data on heat pipes are scarce, and very little is known about their behavior during abnormal operations and close to their thermal limits. Therefore, new experiments and accurate heat pipe simulations are needed to develop reliable closure models. This work describes a joint experimental and numerical investigation into heat pipes that attempts an initial closure of this gap. Further, the numerical and experimental efforts are currently proceeding in parallel, aimed at different aspects of heat pipes. The numerical part is focused on gaps in local closures, and the experiments capture the overall heat pipe behavior.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2020

Baroclinic Residual Circulation and Mass Transport Due to Internal Tides

Baroclinic (BC) tidal residual circulation due to internal tides is investigated around islands over a shallow ridge using a numerical ocean model. Internal tides enhance vertical mixing over shallow slopes, leading to horizontal density gradients that drive BC residual circulation along the main thermocline. For a strongly stratified summer case, the vertical diffusivity estimated by the Mellor and Yamada turbulence closure model exceeds 1 × 10 -2 m 2 s -1 , and the velocity of BC residual circulations reaches 0.2 m s -1 . The magnitude of BC residual circulation is larger than that of barotropic residual tidal circulation, implying that BC residual circulation due to internal tides plays an important role in coastal ocean circulation. Furthermore, BC residual flow accounts for an equal percentage (5%) of the total tidal kinetic energy as the barotropic residual flow under summer stratification conditions. Results from a coupled sediment resuspension and transport model show the growth of intermediate nepheloid layers formed by strong bottom shear stress and BC residual flow. The residual component contributes largely (23% of the total in summer) to the total suspended sediment flux. Seasonal variability is explored, with weaker winter stratification leading to reduced mixing, and thus weaker BC residual circulation and sediment flux. Lastly, the magnitude of the BC residual circulation is also shown to be proportional to the square of the tidal amplitude.

54 ENVIRONMENTAL SCIENCES↗

NASA NTRS1986

Comparison of pressure-strain correlation models for the flow behind a disk

Attention is given to the behavior of Reynolds stresses in the separated wake region behind a disk that is attached in a normal fashion to a long cylinder of small diameter. Computations of the turbulent flow were made in a region beyond a disk by using the second-order closure model of turbulence. It is found that the models of Naot et al. (1970) and Launder et al. (1975) yield similar results and are reliable; the energy distribution may nevertheless be improved for the case of reattaching shear flows by taking the effects of mean strain into account.

Amano, R. S.↗

NASA NTRS2023

Fan Noise Predictions of the NASA Source Diagnostic Test Using Unsteady Simulations with LAVA Part I: Near-Field Aerodynamics and Turbulence

A sliding mesh technique within the Launch, Ascent, and Vehicle Aerodynamics (LAVA) computational framework is validated using the experimental dataset collected as part of the NASA Source Diagnostic Test (SDT) campaign. Two modeling approaches are explored: the unsteady Reynolds-Averaged Navier Stokes (URANS) with Spalart-Allmaras (SA) turbulence model closure, and a hybrid Reynolds-Averaged Navier Stokes/Large Eddy Simulation (RANS/LES) paradigm employing a Zonal Detached Eddy Simulation (ZDES) closure with enhanced shielding protection. Fan stage performance metrics, aerodynamic quantities and turbulent flow structures are analyzed in this work. Initial studies focusing on grid and time-step sensitivity are presented. Sensitivity to different variants of the SA turbulence model is analyzed, supporting the use of the baseline SA model in the production runs. Two conditions are analyzed in detail using URANS and hybrid RANS/LES (HRLES). Mean flow quantities are well-captured by both methods in the low-speed (approach) regime. While URANS misses all the upstream-propagating noise in the inlet due to the rotor-locked tones being evanescent in nature at subsonic fan tip speeds, HRLES captures this broadband component in its pressure field. At the high-speed (sideline) condition, URANS shows better agreement with the SDT data than HRLES in the interstage flow-field. In this regime, URANS captures the tonal content propagating through the inlet, since the tones are now cut-on. Both methods are suitable to capture fan stage performance metrics and mean flow quantities, but only HRLES is able to resolve the fine turbulent structures responsible for broadband noise. The results support the use of the sliding mesh technique implemented in this work for future turbomachinery applications within the LAVA solver framework.

AATT↗

NASA NTRS2023

Fan Noise Predictions of the NASA Source Diagnostic Test Using Unsteady Simulations with LAVA Part I: Near-Field Aerodynamics and Turbulence

AATT↗

DOE OSTI2023

Learning subgrid-scale models with neural ordinary differential equations

We propose a new approach to learning the subgrid-scale model when simulating partial differential equations (PDEs) solved by the method of lines and their representation in chaotic ordinary differential equations, based on neural ordinary differential equations (NODEs). Solving systems with fine temporal and spatial grid scales is an ongoing computational challenge, and closure models are generally difficult to tune. Machine learning approaches have increased the accuracy and efficiency of computational fluid dynamics solvers. In this approach neural networks are used to learn the coarse- to fine-grid map, which can be viewed as subgrid-scale parameterization. We propose a strategy that uses the NODE and partial knowledge to learn the source dynamics at a continuous level. Our method inherits the advantages of NODEs and can be used to parameterize subgrid scales, approximate coupling operators, and improve the efficiency of low-order solvers. Numerical results with the two-scale Lorenz 96 ODE, the convection-diffusion PDE, and the viscous Burgers' PDE are used to illustrate this approach.

97 MATHEMATICS AND COMPUTING↗

NASA NTRS2010

An Attempt to Derive the epsilon Equation from a Two-Point Closure

The goal of this paper is to derive the equation for the turbulence dissipation rate epsilon for a shear-driven flow. In 1961, Davydov used a one-point closure model to derive the epsilon equation from first principles but the final result contained undetermined terms and thus lacked predictive power. Both in 1987 and in 2001, attempts were made to derive the epsilon equation from first principles using a two-point closure, but their methods relied on a phenomenological assumption. The standard practice has thus been to employ a heuristic form of the equation that contains three empirical ingredients: two constants, c(sub 1 epsilon), and c(sub 2 epsilon), and a diffusion term D(sub epsilon) In this work, a two-point closure is employed, yielding the following results: 1) the empirical constants get replaced by c(sub 1), c(sub 2), which are now functions of Kappa and epsilon; 2) c(sub 1) and c(sub 2) are not independent because a general relation between the two that are valid for any Kappa and epsilon are derived; 3) c(sub 1), c(sub 2) become constant with values close to the empirical values c(sub 1 epsilon), c(sub epsilon 2), (i.e., homogenous flows); and 4) the empirical form of the diffusion term D(sub epsilon) is no longer needed because it gets substituted by the Kappa-epsilon dependence of c(sub 1), c(sub 2), which plays the role of the diffusion, together with the diffusion of the turbulent kinetic energy D(sub Kappa), which now enters the new equation (i.e., inhomogeneous flows). Thus, the three empirical ingredients c(sub 1 epsilon), c(sub epsilon 2), D (sub epsilon)are replaced by a single function c(sub 1)(Kappa, epsilon ) or c(sub 2)(Kappa, epsilon ), plus a D(sub Kappa)term. Three tests of the new equation for epsilon are presented: one concerning channel flow and two concerning the shear-driven planetary boundary layer (PBL).

Canuto, V. M.↗

DOE OSTI2021

Sapsan: Framework for Supernovae Turbulence Modeling with Machine Learning

Sapsan is a framework designed to make Machine Learning (ML) more accessible in the study of turbulence, with a focus on astrophysical applications. Sapsan includes modules to load, filter, subsample, batch, and split the data from hydrodynamic (HD) simulations for training and validation. Next, the framework includes built-in conventional and physically-motivated estimators that have been used for turbulence modeling. This ties into Sapsan’s custom estimator module, aimed at designing a custom ML model layer-by-layer, which is the core benefit of using the framework. To share your custom model, every new project created via Sapsan comes with pre-filled, ready-for-release Docker files. Furthermore, training and evaluation modules come with Sapsan as well. The latter, among other features, includes the construction of power spectra and comparison to established analytical turbulence closure models, such as a gradient model. Thus, Sapsan attempts to minimize the hard work required for data preparation and analysis, leaving one to focus on the ML model design itself.

79 ASTRONOMY AND ASTROPHYSICS↗

NASA NTRS1991

The effects of rotation on initially anisotropic homogeneous flows

Rotation of initially anisotropic homogeneous flows is studied using a model spectral tensor. It is shown that the anisotropy changes because of the influence of rotation through phase scrambling. Phase scrambling causes the Reynolds stresses to develop with damped oscillations. The final Reynolds stress anisotropy is found to be proportional to the initial structural tensor anisotropy. Closure models for the rapid pressure strain terms should reflect this change in anisotropy, and should drive the anisotropy to reach its final predicted state. Finally, it is shown that long-time integration using direct numerical simulations should be treated with care because phase scrambling effects on a discrete wave space can cause loss of resolution when time becomes large.

Mansour, Nagi N.↗

DOE OSTI2024

Data‐driven variational method for discrepancy modeling: Dynamics with small‐strain nonlinear elasticity and viscoelasticity

Abstract The effective inclusion of a priori knowledge when embedding known data in physics‐based models of dynamical systems can ensure that the reconstructed model respects physical principles, while simultaneously improving the accuracy of the solution in the previously unseen regions of state space. This paper presents a physics‐constrained data‐driven discrepancy modeling method that variationally embeds known data in the modeling framework. The hierarchical structure of the method yields fine scale variational equations that facilitate the derivation of residuals which are comprised of the first‐principles theory and sensor‐based data from the dynamical system. The embedding of the sensor data via residual terms leads to discrepancy‐informed closure models that yield a method which is driven not only by boundary and initial conditions, but also by measurements that are taken at only a few observation points in the target system. Specifically, the data‐embedding term serves as residual‐based least‐squares loss function, thus retaining variational consistency. Another important relation arises from the interpretation of the stabilization tensor as a kernel function, thereby incorporating a priori knowledge of the problem and adding computational intelligence to the modeling framework. Numerical test cases show that when known data is taken into account, the data driven variational (DDV) method can correctly predict the system response in the presence of several types of discrepancies. Specifically, the damped solution and correct energy time histories are recovered by including known data in the undamped situation. Morlet wavelet analyses reveal that the surrogate problem with embedded data recovers the fundamental frequency band of the target system. The enhanced stability and accuracy of the DDV method is manifested via reconstructed displacement and velocity fields that yield time histories of strain and kinetic energies which match the target systems. The proposed DDV method also serves as a procedure for restoring eigenvalues and eigenvectors of a deficient dynamical system when known data is taken into account, as shown in the numerical test cases presented here.

Masud, Arif↗

DOE OSTI2022

A whole-core steady-state thermal-hydraulic model for annular fuel type fluoride-salt-cooled reactors

A whole-core, steady-state thermal-hydraulic model is developed for the fluoride-salt-cooled small modular advanced high-temperature reactor (SmAHTR) that employs an annular fuel configuration. This pre-conceptual reactor design by Oak Ridge National Laboratory (ORNL) has the annular fuel and moderator pins arranged in a hexagonal layout. The FLiBe coolant flows from the bottom to the top of the core, parallel to the hexagonal bundle. The fuel and moderator pins in the core are discretized into finite volumes and the 3-D heat conduction equation is solved to obtain the temperature profile. Inter-fuel assembly conduction is also addressed. For this fuel assembly configuration, the coolant flows through two distinct regions – the hexagonal pin bundle and the annulus between the fuel pin and the tie rod. The fluid flow through the hexagonal bundles is modeled using the subchannel approach, in which the coolant region is discretized into corner, edge and interior subchannels and the resulting conservation equations are systematically solved. The 1-D mass, momentum and energy equations are solved for the annulus channels between the fuel pin and the tie rod. Pertinent closure models from the literature are employed to close the system of equations. We also performed a preliminary code-to-code comparison between the present model and a CFD model.. The resulting thermal-hydraulic model can provide temperature, flow rate and pressure drop profiles for the different solid and fluid regions throughout the entire core. Whole-core thermal-hydraulic results for a representative power profile are presented and discussed.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

NASA NTRS2009

Enhancement of USM3D Unstructured Flow Solver for High-Speed High-Temperature Shear Flows

Large temperature and pressure fluctuations have a profound effect on turbulence development in transonic and supersonic jets. For high-speed, high-temperature jet flows, standard turbulence models lack the ability to predict the observed mixing rate of a shear layer. Several proposals to address this deficiency have been advanced in the literature to modify the turbulence transport equations in a variety of ways. In the present study, some of the most proven and simple modifications to two-equation turbulence models have been selected and implemented in NASA's USM3D tetrahedral Navier-Stokes flow solver. The modifications include the addition of compressibility correction and pressure dilatation terms in the turbulence transport equations for high-speed flows, and the addition of a simple modification to the Boussinesq's closure model coefficient for high-temperature jets. The efficacy of the extended models is demonstrated by comparison with experimental data for two supersonic axisymmetric jet test cases at design pressure ratio.

Pandya, Mohagna J.↗

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↗

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.↗