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 1,261 records · Page 70
All sourcesNASA NTRSDOE OSTI
NASA NTRS1975

Interaction between a shock wave and a turbulent boundary layer in transonic flow

Interaction between a shock wave and an unseparated turbulent boundary layer is considered. The method of matched asymptotic expansions is used, with solutions valid in the double limit as Reynolds number tends to infinity and Mach number tends to unity. The shock is weak enough that interaction effects can be considered as perturbations to the undisturbed flow; the case considered is that where the sonic line is near the outer edge of the boundary layer. It is shown that, with order estimates for Reynolds stress perturbations, the induced wall pressure distribution can be calculated using only the two outer interaction regions, independent of a specific closure condition and that this solution is in fact a turbulent free interaction solution. A detailed analysis of the inner regions, for which an eddy viscosity model for the Reynolds shear stress is used, provides a description of the variations in velocity, temperature and density near and at the wall.

Adamson, T. C., Jr.↗
DOE OSTI2021

Turbulence theories and statistical closure approaches

When discussing research in physics and in science more generally, it is common to ascribe equal importance to the three components of the scientific trinity: theoretical, experimental, and computational studies. This review will explore the future of modern turbulence theory by tracing its history, which began in earnest with Kolmogorov’s 1941 analysis of turbulence cascade and inertial range [A.N. Kolmogorov, Dokl. Akad. Nauk SSSR, 30, 299, (1941); 32, 19, (1941)]. The 80th Anniversary of Kolmogorov’s landmark study is a welcome opportunity to survey the achievements and evaluate the future of the theoretical approach of turbulence research. Over the years, turbulence theories have been critically important in laying the foundation of our understanding of the nature of turbulent flows. In particular, the Direct Interaction Approximation (DIA) [R.H. Kraichnan, J. Fluid Mech., 5, 497 (1959)] and its subsequent development, known as the statistical closure approach, can be identified as perhaps the most profound single advancement. The remarkable success of the statistical closure has furnished a platform to study such essential concepts as the energy transfer process and interacting scales, and the roles of the straining and sweeping motions. More recently, the quasi-Lagrangian formulation of V. L’vov & I. Procaccia and Kraichnan’s solvable passive scalar model provided powerful ways to explore another fundamental aspect of turbulent flows, the phenomena of intermittency, and the associated anomalous scaling exponents. In the meantime, the theory of fluid equilibria has been developed to describe the large-scale structures that can emerge from turbulent cascades of two-dimensional and geophysical flows at a later time. And yet, despite all these successes, analytical treatments suffer from mathematical complexities. As a result, the utility of theoretical approaches has been limited to relatively idealized flows. On the other hand, in recent decades, computational abilities and experimental facilities have reached an unprecedented scale. Looking beyond the horizon, the imminent deployment of exascale supercomputers will generate complete datasets of the entire flow field of key benchmark flows, allowing researchers to extract additional measurements concerning fully developed, complex turbulent flow fields far beyond those available from the statistical closure theories. Some other developments that could potentially influence the future course of turbulence theories include the advancement of machine learning, artificial intelligence, and data science; likely disruptions arising from the advent of quantum computation; and the increasingly prominent role of turbulence research in providing more accurate climate scientific data. Finally, turbulence theorists can leverage these developments by asking the right questions and developing advanced, sophisticated frameworks that will be able to predict and correlate vast amounts of data from the other two components of the trinity.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS2017

Global Energetics of Solar Flares. V. Energy Closure in Flares and Coronal Mass Ejections

In this study we synthesize the results of four previous studies on the global energetics of solar flares and associated coronal mass ejections (CMEs), which include magnetic, thermal, nonthermal, and CME energies in 399 solar M and X-class flare events observed during the first 3.5 years of the Solar Dynamics Observatory (SDO) mission. Our findings are: (1) The sum of the mean nonthermal energy of flare-accelerated particles (E(sub nt)), the energy of direct heating (E(sub dir)), and the energy in coronal mass ejections (E(sub CME)), which are the primary energy dissipation processes in a flare, is found to have a ratio of (E(sub nt)+E(sub dir)+E(sub CME))/E(sub mag) = 0.87±0.18, compared with the dissipated magnetic free energy E(sub mag), which confirms energy closure within the measurement uncertainties and corroborates the magnetic origin of flares and CMEs; (2) The energy partition of the dissipated magnetic free energy is: 0.51±0.17 in nonthermal energy of ≥ 6 keV electrons, 0.17± 0.17 in nonthermal ≥ 1 MeV ions, 0.07 ± 0.14 in CMEs, and 0.07 ± 0.17 in direct heating; (3) The thermal energy is almost always less than the nonthermal energy, which is consistent with the thick-target model; (4) The bolometric luminosity in white-light flares is comparable with the thermal energy in soft X-rays (SXR); (5) Solar Energetic Particle (SEP) events carry a fraction ≈ 0.03 of the CME energy, which is consistent with CME-driven shock acceleration; and (6) The warm-target model predicts a lower limit of the low-energy cutoff at ec ≈ 6 keV, based on the mean differential emission measure (DEM) peak temperature of T(sub e) = 8.6 MK during flares. This work represents the first statistical study that establishes energy closure in solar flare/CME events.

Solar activity↗
NASA NTRS1988

Advanced physical-chemical life support systems research

A proposed NASA space research and technology development program will provide adequate data for designing closed loop life support systems for long-duration manned space missions. This program, referred to as the Pathfinder Physical-Chemical Closed Loop Life Support Program, is to identify and develop critical chemical engineering technologies for the closure of air and water loops within the spacecraft, surface habitats or mobility devices. Computerized simulation can be used both as a research and management tool. Validated models will guide the selection of the best known applicable processes and in the development of new processes. For the integration of the habitat system, a biological subsystem would be introduced to provide food production and to enhance the physical-chemical life support functions on an ever-increasing basis.

Evanich, Peggy L.↗
NASA NTRS1995

Relaxation approximation in the theory of shear turbulence

Leslie's perturbative treatment of the direct interaction approximation for shear turbulence (Modern Developments in the Theory of Turbulence, 1972) is applied to derive a time dependent model for the Reynolds stresses. The stresses are decomposed into tensor components which satisfy coupled linear relaxation equations; the present theory therefore differs from phenomenological Reynolds stress closures in which the time derivatives of the stresses are expressed in terms of the stresses themselves. The theory accounts naturally for the time dependence of the Reynolds normal stress ratios in simple shear flow. The distortion of wavenumber space by the mean shear plays a crucial role in this theory.

Rubinstein, Robert↗
NASA NTRS1996

LES, DNS, and RANS for the Analysis of High-Speed Turbulent Reacting Flows

A filtered density function (FDF) method suitable for chemically reactive flows is developed in the context of large eddy simulation. The advantage of the FDF methodology is its inherent ability to resolve subgrid scales (SGS) scalar correlations that otherwise have to be modeled. Because of the lack of robust models to accurately predict these correlations in turbulent reactive flows, simulations involving turbulent combustion are often met with a degree of skepticism. The FDF methodology avoids the closure problem associated with these terms and treats the reaction in an exact manner. The scalar FDF approach is particularly attractive since it can be coupled with existing hydrodynamic computational fluid dynamics (CFD) codes.

Colucci, P. J.↗
NASA NTRS1998

Fracture Mechanics Analysis of Stitched Stiffener-Skin Debonding

An analysis based on plate finite elements and the virtual crack closure technique has been implemented to study the effect of stitching on mode I and mode II strain energy release rates for debond configurations. The stitches were modeled as discrete nonlinear fastener elements with a compliance determined by experiment. The axial and shear behavior of the stitches was considered, however, the two compliances and failure loads were assumed to be independent. Both a double cantilever beam (mode I) and a mixed mode skin-stiffener debond configuration were studied. In the double cantilever beam configurations, G(sub I) began to decrease once the debond had grown beyond the first row of stitches and was reduced to zero for long debonds. In the mixed-mode skin-stiffener configurations, G(sub I) showed a similar behavior as in the double cantilever beam configurations, however, G(sub u), continued to increase with increasing debond length.

Glaessgen, E. H.↗
NASA NTRS2023

Reduction in Computational Cost of Progressive Failure Analysis of Composite Structures

Designing aircraft structures requires efficient modeling approaches to iterate on multiple structural configurations to achieve an optimal design. Typically, damage tolerance is not considered at the design stage because of the high computational cost in its implementation within a finite element modeling approach. Therefore, analytical or empirical approaches are often used to size critical damage-tolerant structures once an optimal design is determined. In this study, the Progressive Release eXplicit Virtual Crack-Closure Technique (PRX-VCCT) is assessed for its capability to cost-effectively evaluate skin stringer separation of a blade-stiffened panel that is subjected to seven point bend loads. An initial verification study was performed to evaluate PRX VCCT to accurately simulate skin-stringer separation with respect to existing cohesive element approaches. Furthermore, the influence of element size, ranging from 0.10 in. to 0.40 in., on the total computational time using the PRX-VCCT is investigated. The results indicate that the PRX-VCCT can be used to accurately simulate skin stringer separation using large element lengths (0.40 in.). Additionally, a significant reduction in the computational time to simulate skin stringer separation is observed using the PRX-VCCT. Large-scale progressive damage analysis using PRX-VCCT can be implemented early in the design cycle of composite structures without requiring a global-to-local modeling approach.

Composites↗
DOE OSTI2024

Utah FORGE 5-2615: Determination and Analysis of Thermo-poromechanical Response of Fractured Rock - 2024 Annual Workshop Presentation

This is a presentation on the Determination and Modeling-Informed Analysis of Thermo-poromechanical Response of Fractured Rock for Application to FORGE by the University of Oklahoma, presented by Ahmad Ghassemi. This video presentation discusses how to improve understanding and control of coupled thermoporomechanical (or thermo-hydro-mechanical-THM) processes in reservoir development, and to demonstrate its role in interpretations of the fracture closure pressure. This presentation was featured in the Utah FORGE R&D Annual Workshop on August 13-15, 2024.

15 GEOTHERMAL ENERGY↗
NASA NTRS2017

Global Energetics of Solar Flares. V. Energy Closure in Flares and Coronal Mass Ejections

In this study we synthesize the results of four previous studies on the global energetics of solar flares and associated coronal mass ejections (CMEs), which include magnetic, thermal, nonthermal, and CME energies in 399 solar M and X-class flare events observed during the first 3.5 years of the Solar Dynamics Observatory (SDO) mission. Our findings are: (1) The sum of the mean nonthermal energy of flare-accelerated particles (E(sub nt)), the energy of direct heating (E(sub dir)), and the energy in coronal mass ejections (E(sub CME)), which are the primary energy dissipation processes in a flare, is found to have a ratio of (E(sub nt) + E(sub dir) + E(sub CME))/E(sub mag) = 0.87±0.18, compared with the dissipated magnetic free energy E(sub mag), which confirms energy closure within the measurement uncertainties and corroborates the magnetic origin of flares and CMEs; (2) The energy partition of the dissipated magnetic free energy is: 0.51±0.17 in nonthermal energy of ≥ 6 keV electrons, 0.17± 0.17 in nonthermal ≥ 1 MeV ions, 0.07 ± 0.14 in CMEs, and 0.07 ± 0.17 in direct heating; (3) The thermal energy is almost always less than the nonthermal energy, which is consistent with the thick-target model; (4) The bolometric luminosity in white-light flares is comparable with the thermal energy in soft X-rays (SXR); (5) Solar Energetic Particle (SEP) events carry a fraction ≈ 0.03 of the CME energy, which is consistent with CME-driven shock acceleration; and (6) The warm-target model predicts a lower limit of the low-energy cutoff at e(sub c) ≈ 6 keV, based on the mean differential emission measure (DEM) peak temperature of T(sub e) = 8.6 MK during flares. This work represents the first statistical study that establishes energy closure in solar flare/CME events.

Solar activity↗
NASA NTRS2017

Global Energetics of Solar Flares and CMEs: V. Energy Closure

In this study we synthesize the results of four previous studies on the global energetics of solar flares and associated coronal mass ejections (CMEs), which include magnetic, thermal, nonthermal, and CME energies in 399 solar M and X-class flare events observed during the first 3.5 years of the Solar Dynamics Observatory (SDO) mission. Our findings are: (1) The sum of the mean nonthermal energy of flare-accelerated particles (E(sub nt)), the energy of direct heating (E(sub dir)), and the energy in coronal mass ejections (E(sub CME)), which are the primary energy dissipation processes in a flare, is found to have a ratio of (E(sub nt)+E(sub dir)+E(sub CME))/E(sub mag) = 0.87±0.18, compared with the dissipated magnetic free energy E(sub mag), which confirms energy closure within the measurement uncertainties and corroborates the magnetic origin of flares and CMEs; (2) The energy partition of the dissipated magnetic free energy is: 0.51±0.17 in nonthermal energy of ≥ 6 keV electrons, 0.17± 0.17 in nonthermal ≥ 1 MeV ions, 0.07 ± 0.14 in CMEs, and 0.07 ± 0.17 in direct heating; (3) The thermal energy is almost always less than the nonthermal energy, which is consistent with the thick-target model; (4) The bolometric luminosity in white-light flares is comparable with the thermal energy in soft X-rays (SXR); (5) Solar Energetic Particle (SEP) events carry a fraction ≈ 0.03 of the CME energy, which is consistent with CME-driven shock acceleration; and (6) The warm-target model predicts a lower limit of the low-energy cutoff at e(sub c) ≈ 6 keV, based on the mean differential emission measure (DEM) peak temperature of T(sub e) = 8.6 MK during flares. This work represents the first statistical study that establishes energy closure in solar flare/CME events.

Solar activity↗
NASA NTRS2024

PALMO: An OVERFLOW Machine Learning Airfoil Performance Database

The OVERFLOW Machine Learning Airfoil Performance (PALMO) database has been created to enable robust modeling of airfoil performance in a variety of applications. The database uses OVERFLOW simulation data second-order accurate in time and fourth-order accurate in space with Spalart-Allmaras turbulence closure. The foundation of the in-development PALMO database is the airfoil base cube. Each base cube includes simulation data parametrized over a range of Mach numbers, Reynolds numbers, and angles-of-attack. This first release of the database includes the NACA 4-series airfoils, with parametrization in airfoil thickness and camber from an NACA 0006 to an NACA 4424. In total, 52,480 NACA 4-series calculations were run on the NASA High-End Compute Capability (HECC) supercomputer and the corresponding airfoil performance coefficients are embedded in the Appendix of this document for public distribution. This provides high-order-accurate simulation data covering a wide range of aerospace design applications, which enables users to develop OVERFLOW-quality airfoil performance look-up tables without additional high-performance computing. In addition to engineering design and analysis of aerospace vehicles, PALMO is well suited to be a benchmark dataset for the development and testing of machine learning methods in aerospace engineering. Downstream surrogate models enable OVERFLOW- quality airfoil performance predictions for any arbitrary combination of camber, thickness, Mach number, Reynolds number, and angle-of-attack within the bounds of the database.

Database↗
NASA NTRS1976

Streamline curvature effects on turbulent boundary layers

A theoretical tool has been developed for predicting, in a nonempirical manner, effects of streamline curvature and coordinate-system rotation on turbulent boundary layers. The second-order closure scheme developed by Wilcox and Traci has been generalized for curved streamline flow and for flow in a rotating coordinate system. A physically based straightforward argument shows that curvature/rotation primarily affects the turbulent mixing energy; the argument yields suitable curvature/rotation terms which are added to the mixing-energy equation. Singular-perturbation solutions valid in the wall layer of a curved-wall boundary layer and a fully developed rotating channel flow demonstrate that, with the curvature/rotation terms, the model predicts the curved-wall and the rotating coordinate system laws of the wall. Results of numerical computations of curved-wall boundary layers and of rotating channel flow show that curvature/rotation effects can be computed accurately with second-order closure.

Wilcox, D. C.↗
NASA NTRSDate not supplied

PALMO: An OVERFLOW Machine Learning Airfoil Performance Database

The OVERFLOW Machine Learning Airfoil Performance (PALMO) database has been created to enable robust modeling of airfoil performance in a variety of applications. The PALMO database uses OVERFLOW simulation data second-order accurate in time and fourth-order accurate in space with Spalart-Allmaras turbulence closure. The foundation of the in-development PALMO database is the airfoil base cube. Each base cube includes simulation data parametrized over a range of Mach numbers, Reynolds numbers, and angles-of-attack. This database includes the NACA 4-series airfoils, with parametrization in airfoil thickness and camber from an NACA 0006 to an NACA 4424. In total, 52,480 NACA 4-series OVERFLOW calculations were run on the NASA High-End Compute Capability (HECC) supercomputer. This provides high-order-accurate simulation data covering a wide range of aerospace design applications, which enables users to develop accurate airfoil performance look-up tables without additional high-performance computing. In addition to engineering design and analysis of aerospace vehicles, PALMO is well suited to be a benchmark dataset for the development and testing of machine learning methods in aerospace engineering. This work presents an example PALMO surrogate model that enables accurate airfoil performance predictions for any arbitrary combination of camber, thickness, Mach number, Reynolds number, and angle of attack within the bounds of the database. Airfoil performance tables predicted for an airfoil not used in training the model are used in three-dimensional OVERFLOW simulations to quantify the downstream accuracy on aggregate rotor performance metrics. For the NACA 3415 airfoil, which had no common thickness or camber with the training data, the surrogate predicted and CFD generated tables were within 2.1% of each other in the forward flight lift to drag metric. This suggests that performance tables generated for airfoils within the bounds of the PALMO database will yield aggregate rotor performance predictions on par with tables generated from directly running OVERFLOW airfoil calculations. The PALMO airfoil performance coefficients are available publicly.

Database↗
DOE OSTI2022

Coupling coarse-mesh CFD with fine-mesh CFD for modeling molten-salt reactors in the Virtual Test Bed (VTB)

The Nuclear Energy Advanced Modeling and Simulation (NEAMS) program aims at developing a simulation tool kit to accelerate the development and deployment of nuclear power technologies. NEAMS multiphysics tools have been designed to provide numerical simulation support for the design and licensing of GEN IV reactors. Pronghorn is NEAMS's coarse-mesh computational fluid dynamics (CFD) tool, which is designed to run 3D core transients in GEN IV reactors at a reduced computational cost. To increase their accuracy, coarse-mesh CFD simulations require calibrated closure coefficients. One way of computing these coefficients is via the Nek5000, NEAMS's high-fidelity CFD tool. This article discusses our current research lines in informing Pronghorn closure coefficients via Nek5000 to enable multiphysics simulations of the core cavity of the molten-salt fast reactor. We present an application in which Pronghorn mixing length turbulent viscosity has been calibrated from Nek5000 simulations. The resulting Pronghorn thermal-hydraulics model is then coupled to Griffin, the NEAMS neutron transport solver, to solve for the steady-state configuration of the molten-salt fast reactor. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2022

Coupling coarse mesh CFD with fine mesh CFD for modeling for Molten Salt Reactors in the Virtual Test Bed VTB

The Nuclear Energy Advanced Modeling and Simulation (NEAMS) program aims at developing a simulation tool kit to accelerate the development and deployment of nuclear power technologies. NEAMS multiphysics tools have been designed to provide numerical simulation support for the design and licensing of GEN IV reactors. Pronghorn is NEAMS’s coarse-mesh computational fluid dynamics (CFD) tool, which is designed to run 3D core transients in GEN IV reactors at a reduced computational cost. To increase their accuracy, coarse-mesh CFD simulations require calibrated closure coefficients. One way of computing these coefficients is via the Nek5000, NEAMS’s high-fidelity CFD tool. This article discusses our current research lines in informing Pronghorn closure coefficients via Nek5000 to enable multiphysics simulations of the core cavity of the molten-salt fast reactor. We present an application in which Pronghorn mixing length turbulent viscosity has been calibrated from Nek5000 simulations. The resulting Pronghorn thermal-hydraulics model is then coupled to Griffin, the NEAMS neutron transport solver, to solve for the steady-state configuration of the molten-salt fast reactor.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
NASA NTRS1999

A Shell/3D Modeling Technique for Delaminations in Composite Laminates

A shell/3D modeling technique was developed for which a local solid finite element model is used only in the immediate vicinity of the delamination front. The goal was to combine the accuracy of the full three-dimensional solution with the computational efficiency of a plate or shell finite element model. Multi-point constraints provide a kinematically compatible interface between the local 3D model and the global structural model which has been meshed with plate or shell finite elements. For simple double cantilever beam (DCB), end notched flexure (ENF), and single leg bending (SLB) specimens, mixed mode energy release rate distributions were computed across the width from nonlinear finite element analyses using the virtual crack closure technique. The analyses served to test the accuracy of the shell/3D technique for the pure mode I case (DCB), mode II case (ENF) and a mixed mode I/II case (SLB). Specimens with a unidirectional layup where the delamination is located between two 0 plies, as well as a multidirectional layup where the delamination is located between two non-zero degree plies, were simulated. For a local 3D model extending to a minimum of about three specimen thicknesses in front of and behind the delamination front, the results were in good agreement with mixed mode strain energy release rates obtained from computations where the entire specimen had been modeled with solid elements. For large built-up composite structures modeled with plate elements, the shell/3D modeling technique offers a great potential, since only a relatively small section in the vicinity of the delamination front needs to be modeled with solid elements.

Krueger, Ronald↗