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 199 records · Page 11

All sources NASA NTRS DOE OSTI

NASA NTRS1997

Prediction of Crack Growth under Variable-Amplitude Loading in Thin-Sheet 2024-T3 Aluminum Alloys

The present paper is concerned with the application of a "plasticity-induced" crack closure model to study fatigue crack growth under various load histories. The model was based on the Dugdale model but modified to leave plastically deformed material in the wake of the advancing crack. The model was used to correlate crack growth rates under constant-amplitude loading and then used to predict crack growth under variable-amplitude and spectrum loading on thin-sheet 2024- T3 aluminum alloys. Predicted crack-opening stresses agreed well with test data from the literature. The crack-growth lives agreed within a factor of two for single and repeated spike overloads/underloads and within 20 percent for spectrum loading. Differences were attributed to fretting-product-debris-induced closure and three-dimensional affects not included in the model.

Newman, J. C., Jr.↗

NASA NTRS1995

On the consistency of Reynolds stress turbulence closures with hydrodynamic stability theory

The consistency of second-order closure models with results from hydrodynamic stability theory is analyzed for the simplified case of homogeneous turbulence. In a recent study, Speziale, Gatski, and MacGiolla Mhuiris showed that second-order closures are capable of yielding results that are consistent with hydrodynamic stability theory for the case of homogeneous shear flow in a rotating frame. It is demonstrated in this paper that this success is due to the fact that the stability boundaries for rotating homogeneous shear flow are not dependent on the details of the spatial structure of the disturbances. For those instances where they are -- such as in the case of elliptical flows where the instability mechanism is more subtle -- the results are not so favorable. The origins and extent of this modeling problem are examined in detail along with a possible resolution based on rapid distortion theory (RDT) and its implications for turbulence modeling.

Speziale, Charles G.↗

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

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↗

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

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

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

NASA NTRSDate not supplied

Recent Advancements in the PATO Material Response Code

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

Material Response↗

NASA NTRS1996

Fatigue response of perforated titanium for application in laminar flow control

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

Johnson, W. Steven↗

NASA NTRSDate not supplied

Recent Advancements in the PATO Material Response Code

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

Thermal Protection Systems↗

NASA NTRS1983

Shock noise features using the SCIPVIS code

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

Seiner, J. M.↗