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 721 records · Page 40

All sources NASA NTRS DOE OSTI

NASA NTRS1988

A re-evaluation of finite-element models and stress-intensity factors for surface cracks emanating from stress concentrations

A re-evaluation of the 3-D finite-element models and methods used to analyze surface crack at stress concentrations is presented. Previous finite-element models used by Raju and Newman for surface and corner cracks at holes were shown to have ill-shaped elements at the intersection of the hole and crack boundaries. These ill-shaped elements tended to make the model too stiff and, hence, gave lower stress-intensity factors near the hole-crack intersection than models without these elements. Improved models, without these ill-shaped elements, were developed for a surface crack at a circular hole and at a semi-circular edge notch. Stress-intensity factors were calculated by both the nodal-force and virtual-crack-closure methods. Both methods and different models gave essentially the same results. Comparisons made between the previously developed stress-intensity factor equations and the results from the improved models agreed well except for configurations with large notch-radii-to-plate-thickness ratios. Stress-intensity factors for a semi-elliptical surface crack located at the center of a semi-circular edge notch in a plate subjected to remote tensile loadings were calculated using the improved models. The ratio of crack depth to crack length ranged form 0.4 to 2; the ratio of crack depth to plate thickness ranged from 0.2 to 0.8; and the ratio of notch radius to the plate thickness ranged from 1 to 3. The models had about 15,000 degrees-of-freedom. Stress-intensity factors were calculated by using the nodal-force method.

Tan, P. W.↗

DOE OSTI2026

elm-diagnostics

elm-diagnostics is a Python package for computing diagnostic analyses and visualizations for the E3SM Land Model (ELM) component and is meant to support new feature development in ELM. The tool reads model history files and performs quantitative analyses including budget-closure checking, variable transformations, temporal aggregations, and statistical summaries to support model evaluation, validation, and scientific interpretation. The framework is designed for extensibility, with modular architecture enabling straightforward addition of new diagnostic methods, derived variables, analysis types, visualization approaches, and model-specific adaptations

Hoffman, Matt [Los Alamos National Laboratory]↗

NASA NTRS1989

Prediction of separated transonic wing flows with a non-equilibrium algebraic model

A nonequilibrium algebraic turbulence model, which is based on the turbulence closure scheme of Johnson and King (1985), is proposed to predict separated transonic wing flows. The influence of history effects are modeled by solving a partial differential equation for the maximum total Reynolds shear stress, which is then used to scale the eddy viscosity of an algebraic model. The turbulence model is implemented in a three-dimensional, Reynolds-averaged Navier-Stokes code. Comparisons with experimental data are presented which show clearly that the nonequilibrium type of turbulence model is essential for accurate prediction of transonic separated flows.

Abid, Ridha↗

NASA NTRS2008

Hybrid LES/RANS Simulation of Transverse Sonic Injection into a Mach 2 Flow

A computational study of transverse sonic injection of air and helium into a Mach 1.98 cross-flow is presented. A hybrid large-eddy simulation / Reynolds-averaged Navier-Stokes (LES/RANS) turbulence model is used, with the two-equation Menter baseline (Menter-BSL) closure for the RANS part of the flow and a Smagorinsky-type model for the LES part of the flow. A time-dependent blending function, dependent on modeled turbulence variables, is used to shift the closure from RANS to LES. Turbulent structures are initiated and sustained through the use of a recycling / rescaling technique. Two higher-order discretizations, the Piecewise Parabolic Method (PPM) of Colella and Woodward, and the SONIC-A ENO scheme of Suresh and Huyhn are used in the study. The results using the hybrid model show reasonably good agreement with time-averaged Mie scattering data and with experimental surface pressure distributions, even though the penetration of the jet into the cross-flow is slightly over-predicted. The LES/RANS results are used to examine the validity of commonly-used assumptions of constant Schmidt and Prandtl numbers in the intense mixing zone downstream of the injection location.

Boles, John A.↗

NASA NTRS1991

A review of near-wall Reynolds-stress

The advances made in second-order near-wall turbulence closures are summarized. All closures examined are based on some form of high Reynolds number models for the Reynolds stress and the turbulent kinetic energy dissipation rate equations. Consequently, most near-wall closures proposed to data attempt to modify the high Reynolds number models for the dissipation rate equation so that the resultant models are applicable all the way to the wall. The near-wall closures are examined for their asymptotic behavior so that they can be compared with the proper near-wall behavior of the exact equations. A comparison of the closure's performance in the calculation of a low Reynolds number plane channel flow is carried out. In addition, the closures are evaluated for their ability to predict the turbulence statistics and the limiting behavior of the structure parameters compared to direct simulation data.

So, R. M. C.↗

DOE OSTI2021

Benchmarking NIMROD continuum kinetic formulations through the steady-state poloidal flow

In this work, continuum kinetic formulations are employed as a mechanism to include closure physics in an extended magnetohydrodynamics model. Two continuum kinetic approaches have been implemented in the plasma fluid code NIMROD [Sovinec et al., “Nonlinear magnetohydrodynamics with high-order finite elements,” J. Comput. Phys. 195, 355 (2004)] including a Chapman–Enskog-like (CEL) formulation and a more conventional df approach. Ion kinetic closure schemes are employed to describe the neoclassical flow properties in axisymmetric toroidal geometry. In particular, predictions for steady-state values of poloidal flow profiles in tokamak geometry are provided using both the df formulation and two different solution techniques for the CEL approach. These results are benchmarked against analytic theory predictions as well as results from the drift kinetic code DK4D. The continuum kinetic formulations employed here show agreement with both the analytic theory and DK4D results, and offer a novel velocity space representation involving higher-order finite elements in pitch angle.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRS1973

Solutions of the chemical kinetic equations for initially inhomogeneous mixtures.

Following the recent discussions by O'Brien (1971) and Donaldson and Hilst (1972) of the effects of inhomogeneous mixing and turbulent diffusion on simple chemical reaction rates, the present report provides a more extensive analysis of when inhomogeneous mixing has a significant effect on chemical reaction rates. The analysis is then extended to the development of an approximate chemical sub-model which provides much improved predictions of chemical reaction rates over a wide range of inhomogeneities and pathological distributions of the concentrations of the reacting chemical species. In particular, the development of an approximate representation of the third-order correlations of the joint concentration fluctuations permits closure of the chemical sub-model at the level of the second-order moments of these fluctuations and the mean concentrations.

Hilst, G. R.↗

NASA NTRS1973

High-performance space shuttle auxiliary propellant valve system

Several potential valve closures for the space shuttle auxiliary propulsion system (SS/APS) were investigated analytically and experimentally in a modeling program. The most promising of these were analyzed and experimentally evaluated in a full-size functional valve test fixture of novel design. The engineering investigations conducted for both model and scale evaluations of the SS/APS valve closures and functional valve fixture are described. Preliminary designs, laboratory tests, and overall valve test fixture designs are presented, and a final recommended flightweight SS/APS valve design is presented.

Smith, G. M.↗

NASA NTRS1996

Parallel Simulation of Unsteady Turbulent Flames

Time-accurate simulation of turbulent flames in high Reynolds number flows is a challenging task since both fluid dynamics and combustion must be modeled accurately. To numerically simulate this phenomenon, very large computer resources (both time and memory) are required. Although current vector supercomputers are capable of providing adequate resources for simulations of this nature, the high cost and their limited availability, makes practical use of such machines less than satisfactory. At the same time, the explicit time integration algorithms used in unsteady flow simulations often possess a very high degree of parallelism, making them very amenable to efficient implementation on large-scale parallel computers. Under these circumstances, distributed memory parallel computers offer an excellent near-term solution for greatly increased computational speed and memory, at a cost that may render the unsteady simulations of the type discussed above more feasible and affordable.This paper discusses the study of unsteady turbulent flames using a simulation algorithm that is capable of retaining high parallel efficiency on distributed memory parallel architectures. Numerical studies are carried out using large-eddy simulation (LES). In LES, the scales larger than the grid are computed using a time- and space-accurate scheme, while the unresolved small scales are modeled using eddy viscosity based subgrid models. This is acceptable for the moment/energy closure since the small scales primarily provide a dissipative mechanism for the energy transferred from the large scales. However, for combustion to occur, the species must first undergo mixing at the small scales and then come into molecular contact. Therefore, global models cannot be used. Recently, a new model for turbulent combustion was developed, in which the combustion is modeled, within the subgrid (small-scales) using a methodology that simulates the mixing and the molecular transport and the chemical kinetics within each LES grid cell. Finite-rate kinetics can be included without any closure and this approach actually provides a means to predict the turbulent rates and the turbulent flame speed. The subgrid combustion model requires resolution of the local time scales associated with small-scale mixing, molecular diffusion and chemical kinetics and, therefore, within each grid cell, a significant amount of computations must be carried out before the large-scale (LES resolved) effects are incorporated. Therefore, this approach is uniquely suited for parallel processing and has been implemented on various systems such as: Intel Paragon, IBM SP-2, Cray T3D and SGI Power Challenge (PC) using the system independent Message Passing Interface (MPI) compiler. In this paper, timing data on these machines is reported along with some characteristic results.

Menon, Suresh↗

NASA NTRS1973

A study of crack closure in fatigue

Crack closure phenomenon in fatigue was studied by using a Ti-6Al-4V titanium alloy. The occurrence of crack closure was directly measured by an electrical-potential method, and indirectly by load-strain measurement. The experimental results showed that the onset of crack closure depends on both the stress ratio, and the maximum stress intensity factor. No crack closure was observed for stress ratio, greater than 0.3 in this alloy. A two-dimensional elastic model was used to explain the behavior of the recorded load-strain curves. Closure force was estimated by using this model. Yield level stress was found near the crack tip. Based on this estimated closure force, the crack opening displacement was calculated. This result showed that onset of crack closure detected by electrical-potential measurement and crack-opening-displacement measurement is the same. The implications of crack closure on fatigue crack are considered. The experimental results show that crack closure cannot fully account for the effect of stress ratio, on crack growth, and that it cannot be regarded as the sole cause for delay.

Shih, T. T.↗

NASA NTRS1977

A dual-loop model of the human controller in single-axis tracking tasks

A dual loop model of the human controller in single axis compensatory tracking tasks is introduced. This model possesses an inner-loop closure which involves feeding back that portion of the controlled element output rate which is due to control activity. The sensory inputs to the human controller are assumed to be system error and control force. The former is assumed to be sensed via visual, aural, or tactile displays while the latter is assumed to be sensed in kinesthetic fashion. A nonlinear form of the model is briefly discussed. This model is then linearized and parameterized. A set of general adaptive characteristics for the parameterized model is hypothesized. These characteristics describe the manner in which the parameters in the linearized model will vary with such things as display quality. It is demonstrated that the parameterized model can produce controller describing functions which closely approximate those measured in laboratory tracking tasks for a wide variety of controlled elements.

Hess, R. A.↗

NASA NTRS2011

An Overview of NASA's Integrated Design and Engineering Analysis (IDEA) Environment

Historically, the design of subsonic and supersonic aircraft has been divided into separate technical disciplines (such as propulsion, aerodynamics and structures), each of which performs design and analysis in relative isolation from others. This is possible, in most cases, either because the amount of interdisciplinary coupling is minimal, or because the interactions can be treated as linear. The design of hypersonic airbreathing vehicles, like NASA's X-43, is quite the opposite. Such systems are dominated by strong non-linear interactions between disciplines. The design of these systems demands that a multi-disciplinary approach be taken. Furthermore, increased analytical fidelity at the conceptual design phase is highly desirable, as many of the non-linearities are not captured by lower fidelity tools. Only when these systems are designed from a true multi-disciplinary perspective, can the real performance benefits be achieved and complete vehicle systems be fielded. Toward this end, the Vehicle Analysis Branch at NASA Langley Research Center has been developing the Integrated Design and Engineering Analysis (IDEA) Environment. IDEA is a collaborative environment for parametrically modeling conceptual and preliminary designs for launch vehicle and high speed atmospheric flight configurations using the Adaptive Modeling Language (AML) as the underlying framework. The environment integrates geometry, packaging, propulsion, trajectory, aerodynamics, aerothermodynamics, engine and airframe subsystem design, thermal and structural analysis, and vehicle closure into a generative, parametric, unified computational model where data is shared seamlessly between the different disciplines. Plans are also in place to incorporate life cycle analysis tools into the environment which will estimate vehicle operability, reliability and cost. IDEA is currently being funded by NASA?s Hypersonics Project, a part of the Fundamental Aeronautics Program within the Aeronautics Research Mission Directorate. The environment is currently focused around a two-stage-to-orbit configuration with a turbine-based combined cycle (TBCC) first stage and a reusable rocket second stage. IDEA will be rolled out in generations, with each successive generation providing a significant increase in capability, either through increased analytic fidelity, expansion of vehicle classes considered, or by the inclusion of advanced modeling techniques. This paper provides the motivation behind the current effort, an overview of the development of the IDEA environment (including the contents and capabilities to be included in Generation 1 and Generation 2), and a description of the current status and detail of future plans.

Robinson, Jeffrey S.↗

NASA NTRS2015

TPSAS-NF1676L-22800-DND

This lecture will provide an overview of an advanced turbulence closure scheme used in cloud-resolving models (CRMs) and climate modeling effort with the superparameterization approach. The CRM is embedded into each atmospheric column of a climate model as a super parameterization. It replaces the traditional subgrid-scale physical parameterizations in a conventional climate model. Two sensitivity simulations with doubling of CO2 and an increase of sea surface temperature, respectively, are performed and compared to similar simulations without this advanced turbulence closure in CRM to better understand the role of low-level clouds in climate sensitivity.

Kuan-Man Xu↗

NASA NTRS1995

Classical closure theory and Lam's interpretation of epsilon-RNG

Lam's phenomenological epsilon-renormalization group (RNG) model is quite different from the other members of that group. It does not make use of the correspondence principle and the epsilon-expansion procedure. We demonstrate that Lam's epsilon-RNG model is essentially the physical space version of the classical closure theory in spectral space and consider the corresponding treatment of the eddy viscosity and energy backscatter.

Zhou, YE↗

NASA NTRS1995

Applications of direct numerical simulation of turbulence in second order closures

This paper discusses two methods of developing models for the rapid pressure-strain correlation term in the Reynolds stress transport equation using direct numerical simulation (DNS) data. One is a perturbation about isotropic turbulence, the other is a perturbation about two-component turbulence -- an extremely anisotropic turbulence. A model based on the latter method is proposed and is found to be very promising when compared with DNS data and other models.

Shih, Tsan-Hsing↗

NASA NTRS1992

The evolution equation for the flame surface density in turbulent premixed combustion

One central ingredient in flamelet models for turbulent premixed combustion is the flame surface density. This quantity conveys most of the effects of the turbulence on the rate of energy release and is obtained via a modeled transport equation, called the Sigma-equation. Past theoretical work has produced a rigorous approach that leads to an exact, but unclosed, formulation for the turbulent Sigma-equation. In this exact Sigma-equation, it appears that the dynamical properties of the flame surface density are determined by a single parameter, namely the turbulent flame stretch. Unfortunately, the flame surface density and the turbulent flame stretch are not available from experiments and, in the absence of experimental data, little is known on the validity of the closure assumptions used in current flamelet models. Direct Numerical Simulation (DNS) is the obvious, complementary approach to get basic information on these fundamental quantities. Three-dimensional DNS of premixed flames in isotropic turbulent flow is used to estimate the different terms appearing in the Sigma-equation. A new methodology is proposed to provide the source and sink terms for the flame surface density, resolved both temporally and spatially throughout the turbulent flame brush. Using this methodology, the effects of the Lewis number on the rate of production of flame surface area are described in great detail and meaningful comparisons with flamelet models can be performed. The analysis reveals in particular the tendency of the models to overpredict flame surface dissipation as well as their inability to reproduce variations due to thermo-diffusive phenomena. Thanks to the detailed information produced by a DNS-based analysis, this type of comparison not only underscores the shortcomings of current models but also suggests ways to improve them.

Trouve, A.↗

NASA NTRS1993

The evolution equation for the flame surface density in turbulent premixed combustion

The mean reaction rate in flamelet models for turbulent premixed combustion depends on two basic quantities: a mean chemical rate, called the flamelet speed, and the flame surface density. Our previous work had been primarily focused on the problem of the structure and topology of turbulent premixed flames, and it was then determined that the flamelet speed, when space-averaged, is only weakly sensitive to the turbulent flow field. Consequently, the flame surface density is the key quantity that conveys most of the effects of the turbulence on the rate of energy release. In flamelet models, this quantity is obtained via a modeled transport equation called the Sigma-equation. Past theoretical work has produced a rigorous approach that leads to an exact but unclosed formulation for the turbulent Sigma-equation. In the exact Sigma-equation, it appears that the dynamical properties of the flame surface density are determined by a single parameter, namely the turbulent flame stretch. Unfortunately, the turbulent flame stretch as well as the flame surface density is not available from experiments, and, in the absence of experimental data, little is known on the validity of the closure assumptions used in current flamelet models. Direct Numerical Simulation (DNS) is the alternative approach to get basic information on these fundamental quantities. In the present work, three-dimensional DNS of premixed flames in isotropic turbulent flow is used to estimate the different terms appearing in the Sigma-equation. A new methodology is proposed to provide the source and sink terms for the flame surface density, resolved both temporally and spatially throughout the turbulent flame brush. Using this methodology, our objective is to extract the turbulent flame stretch from the DNS data base and then perform extensive comparisons with flamelet models. Thanks to the detailed information produced by the DNS-based analysis, it is expected that this type of comparison will not only underscore the shortcomings of current models, but also suggest ways to improve them.

Trouve, Arnaud↗

NASA NTRS1997

Testing and Implementation of Advanced Reynolds Stress Models

A research program was proposed for the testing and implementation of advanced turbulence models for non-equilibrium turbulent flows of aerodynamic importance that are of interest to NASA. Turbulence models that are being developed in connection with the Office of Naval Research ARI in Non-equilibrium are provided for implementation and testing in aerodynamic flows at NASA Langley Research Center. Close interactions were established with researchers at Nasa Langley RC and refinements to the models were made based on the results of these tests. The models that have been considered include two-equation models with an anisotropic eddy viscosity as well as full second-order closures. Three types of non-equilibrium corrections to the models have been considered in connection with the ARI on Nonequilibrium Turbulence: conducted for ONR.

Speziale, Charles G.↗