Search NASA⌕ Search

SEARCH · Search NASA

Results for “Grid Stress”

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 253 records · Page 14

Revisiting Turbulence Model Validation for High-Mach Number Axisymmetric Compression Corner Flows

Two axisymmetric shock-wave/boundary-layer interaction (SWBLI) cases are used to benchmark one- and two-equation Reynolds-averaged Navier-Stokes (RANS) turbulence models. This validation exercise was executed in the philosophy of the NASA Turbulence Modeling Resource and the AIAA Turbulence Model Benchmarking Working Group. Both SWBLI cases are from the experiments of Kussoy and Horstman for axisymmetric compression corner geometries with SWBLI inducing flares of 20 and 30 degrees, respectively. The freestream Mach number was approximately 7. The RANS closures examined are the Spalart-Allmaras one-equation model and the Menter family of kappa − omega two equation models including the Baseline and Shear Stress Transport formulations. The Wind-US and CFL3D RANS solvers are employed to simulate the SWBLI cases. Comparisons of RANS solutions to experimental data are made for a boundary layer survey plane just upstream of the SWBLI region. In the SWBLI region, comparisons of surface pressure and heat transfer are made. The effects of inflow modeling strategy, grid resolution, grid orthogonality, turbulent Prandtl number, and code-to-code variations are also addressed.

Shock Wave↗

USM3D-ME Solutions for RANS Test Suite of High-Fidelity CFD Verification Workshop 2024

USM3D-ME solutions for the High-Fidelity CFD Verification Workshop 2024 are presented. Solutions are computed for the Reynolds-Averaged Navier-Stokes equations using the Spalart-Allmaras one equation turbulence model, SA-neg-QCR2000-R, that is enhanced with a nonlinear correction to the turbulence stresses and a correction for flow rotation. USM3D-ME solutions have been computed for the three verification cases, namely, a Joukowski airfoil, a subsonic three-dimensional flow around an extruded NACA 0012 wing in tunnel, and subsonic flow around a high-lift wing-body configuration. Iterative convergence of USM3D-ME solutions on several grids within various grid families is demonstrated. Grid convergence of integrated forces, pitching moment and sectional variations of surface pressure and skin friction are also presented.

CFD↗

A finite-difference outer layer and integral inner layer method for the solution of the turbulent boundary layer equations

A new turbulent boundary-layer method is developed which models the inner region with the law of the wall while the outer region uses Clauser's eddy viscosity in Matsuno's finite-difference method. The match point between the inner and outer regions as well as the wall shear stress are determined at each marching step during the computation. Results obtained for incompressible, two-dimensional flow over flat plates and ellipses are compared with solutions from a baseline method which uses a finite-difference method for the entire boundary layer. Since the present method used the finite-difference method in the outer region only, the number of grid points required was about half that needed for the baseline method. Accurate displacement and momentum thicknesses were predicted for all cases. Skin friction was predicted well for the flat plate, but the accuracy decreased significantly for the ellipses. Adding a wake functions to the law of the wall allows some of the pressure gradient effect to be taken into account thereby increasing the accuracy of the method.

Barnwell, R. W.↗

Orbiter Return-To-Flight Entry Aeroheating

The Columbia accident on February 1, 2003 began an unprecedented level of effort within the hypersonic aerothermodynamic community to support the Space Shuttle Program. During the approximately six month time frame of the primary Columbia Accident Investigation Board activity, many technical disciplines were involved in a concerted effort to reconstruct the last moments of the Columbia and her crew, and understand the critical events that led to that loss. Significant contributions to the CAIB activity were made by the hypersonic aerothermodynamic community(REF CAIB) in understanding the re-entry environments that led to the propagation of an ascent foam induced wing leading edge damage to a subsequent breech of the wing spar of Columbia, and the subsequent breakup of the vehicle. A core of the NASA hypersonic aerothermodynamics team that was involved in the CAIB investigation has been combined with the United Space Alliance and Boeing Orbiter engineering team in order to position the Space Shuttle Program with a process to perform in-flight Thermal Protection System damage assessments. This damage assessment process is now part of the baselined plan for Shuttle support, and is a direct out-growth of the Columbia accident and NASAs response. Multiple re-entry aeroheating tools are involved in this damage assessment process, many of which have been developed during the Return To Flight activity. In addition, because these aeroheating tools are part of an overall damage assessment process that also involves the thermal and stress analyses community, in addition to a much broader mission support team, an integrated process for performing the damage assessment activities has been developed by the Space Shuttle Program and the Orbiter engineering community. Several subsets of activity in the Orbiter aeroheating communities support to the Return To Flight effort have been described in previous publications (CFD?, Cavity Heating? Any BLT? Grid Generation?). This work will provide a description of the integrated process utilized to perform Orbiter tile damage assessment, and in particular will seek to provide a description of the integrated aeroheating tools utilized to perform these assessments. Individual aeroheating tools will be described which provide the nominal re-entry heating environment characterization for the Orbiter, the heating environments for tile damage, heating effects due to exposed Thermal Protection System substrates, the application of Computational Fluid Dynamics for the description of tile cavity heating, and boundary layer transition prediction. This paper is meant to provide an overall view of the integrated aeroheating assessment process for tile damage assessment as one of a sequence of papers on the development of the boundary layer transition prediction capability in support of Space Shuttle Return To Flight efforts.

Campbell, Charles H.↗

Measurement, data analysis and prediction of pyrotechnic shock from pin-pullers and separation joints

To evaluate the potential for spacecraft damage caused by activation of pyrotechnic mechanisms, pyrotechnic shock tests were conducted on three configurations: pin-pullers on an orthogonal double Hopkinson bar arrangement, pin-pullers on a mock-up of the HALOE (Halogen Occultation Experiment) structure, and a section of separation joint on a single Hopkinson bar. Strains and accelerations were measured. The strains were converted to output stresses, forces and moments. Acceleration shock response spectra are presented for both acceleration- and force-time signals. The devices were attached to the Hopkinson bars with adaptors, which are typical of attachments used in practice. In order to predict the effect of the adaptors, finite element models of the adaptor-bar combinations were analyzed and results are presented with regard to grid size, time interval, and frequency range required to predict the response to pulses of 10 and 100 microsecond duration.

Evans, Maria J.↗

Evolution of three-dimensionality in stable and unstable curved mixing layers

Results are presented of an experimental study tracing the effects of mild streamwise curvature (delta/R is less than 5 percent) on the 3D structure of a two-stream mixing layer at high Reynolds number (about 27,000). Measurements of the mean and fluctuating velocities were obtained on fine cross-plane grids at 11 streamwise locations using a rotatable cross-wide probe. Well-organized streamwise vorticity was generated in both cases through the braid instability. Although the vortex structures decayed with streamwise distance in both cases, the rate of decay for the unstable case was lower. As a result, the unstable layer exhibited noticeable spanwise variations in the mean velocity and Reynolds stresses. Both cases achieved linear growth, but the rate of growth for the unstable case was higher. The far-field spanwise-averaged peak Reynolds stresses were significantly higher for the destabilized case compared to the stabilized case, which exhibited levels comparable to those of a straight case.

Plesniak, Michael W.↗

An implicit multigrid scheme for hypersonic strong-interaction flowfields

An implicit multigrid method for the compressible Navier-Stokes equations is applied to a series of two-dimensional, steady, laminar hypersonic flows over a compression ramp, including both attached and separated flow conditions. The algorithm uses upwind-biased differencing for the convective and pressure terms and central differencing for the shear-stress and heat-transfer terms. An implicit spatially factored scheme is used to advance the solution in terms on a given mesh, and a V-cycle multigrid strategy is used to accelerate convergence. The multigrid algorithm provided an order of magnitude decrease in the computational time required to obtain a converged position of the separation location and enabled a substantially grid-independent result to be obtained for each of the flows considered.

Thomas, James L.↗

Fluid Structure Interaction in a Turbine Blade

An unsteady, three dimensional Navier-Stokes solution in rotating frame formulation for turbomachinery applications is presented. Casting the governing equations in a rotating frame enabled the freezing of grid motion and resulted in substantial savings in computer time. The turbine blade was computationally simulated and probabilistically evaluated in view of several uncertainties in the aerodynamic, structural, material and thermal variables that govern the turbine blade. The interconnection between the computational fluid dynamics code and finite element structural analysis code was necessary to couple the thermal profiles with the structural design. The stresses and their variations were evaluated at critical points on the Turbine blade. Cumulative distribution functions and sensitivity factors were computed for stress responses due to aerodynamic, geometric, mechanical and thermal random variables.

Gorla, Rama S. R.↗

Calculation of Thermally-Induced Displacements in Spherically Domed Ion Engine Grids

An analytical method for predicting the thermally-induced normal and tangential displacements of spherically domed ion optics grids under an axisymmetric thermal loading is presented. A fixed edge support that could be thermally expanded is used for this analysis. Equations for the displacements both normal and tangential to the surface of the spherical shell are derived. A simplified equation for the displacement at the center of the spherical dome is also derived. The effects of plate perforation on displacements and stresses are determined by modeling the perforated plate as an equivalent solid plate with modified, or effective, material properties. Analytical model results are compared to the results from a finite element model. For the solid shell, comparisons showed that the analytical model produces results that closely match the finite element model results. The simplified equation for the normal displacement of the spherical dome center is also found to accurately predict this displacement. For the perforated shells, the analytical solution and simplified equation produce accurate results for materials with low thermal expansion coefficients.

Soulas, George C.↗

High-Order Methods in NASA’s Next Generation of Computational Fluid Dynamics Tools

The missions of the National Aeronautics and Space Administration (NASA) routinely produce unique requirements and challenges for development and application of Computational Fluid Dynamics (CFD) methods. NASA presently embodies four distinct Mission Directorates: Aeronautics Research, Exploration Systems, Science, and Space Operations. These missions generate requirements for systems that operate in a wide variety of environments. They range from the high-speed flight of aerodynamically optimized vehicles operating in the earth’s atmosphere to spacecraft designed for missions that don’t favor aerodynamic optimization, some operating in the atmosphere of planets and planetary moons such as Mars and Venus or Saturn’s moon Titan. Systems supporting these vehicles, such as rocket and jet propulsion, reaction control systems, fluid and thermal transfer systems, etc. can also generate their own unique set of flow phenomena that challenge today’s CFD methodology. Through the NASA Engineering and Safety Center (NESC), NASA annually conducts state-of-the-discipline assessments in fifteen distinct engineering disciplines. These assessments are performed by the NASA Technical Fellows that lead Technical Discipline Teams (TDT) of recognized experts in these fifteen areas. In the Aerosciences discipline, three topics have been identified as the top challenges for the discipline: aero-plume interaction prediction, unsteady separated flows, and aerothermodynamic prediction. These challenge areas are defined by the Agency’s high-risk projects and problems on which the NESC is requested to perform independent testing, analysis, and assessments. When viewed as a whole, these tests, analyses, and assessments provide a clear view of the recurring technical challenges facing Agency engineers and researchers and can be used to guide future research and technology development. The present state-of-the-art in the application of CFD at NASA is the use of Reynolds-Averaged Navier- Stokes (RANS) solvers, primarily executed in a steady-state mode of operation. In isolated cases, Unsteady RANS (URANS) solvers have been employed when steady RANS solutions produce poorly converging or oscillating results or in cases, such as aeroelastic analysis, which require unsteady aerodynamic simulation. For most traditional external and internal aerodynamic flows, structured overset grids or unstructured grids are employed to minimize geometric modeling and grid generation times. Grid adaptation, primarily as a series of coarse-grain intermediate processing steps is also seeing use on particularly complex flow problems and configurations. In the case of aerothermodynamic flows, engineers have been forced to continue to employ structured grid techniques as the present unstructured grid methodology has proven inadequate in the prediction of surface heating. In the area of aero-plume interaction modeling, two-gas, frozen chemistry simulation is generally the state-of-the- art, with some production solvers capable of predicting flows with only a single gas component. Prediction of flows falling into the afore-mentioned top Aerosciences technical challenges have severely stressed the present state-of-the-art in CFD prediction and for some problems, such as unsteady separated flows and aero-plume interaction cases, engineers have begun employing Large Eddy Simulation (LES) and Hybrid RANS/LES techniques. In some isolated aero-propulsion interaction cases, chemically reacting flow simulations have been applied. These methods are highly evolutionary and engineers have little experience in their application, so they cannot be heavily relied upon in today’s application environment. Therefore, this leads one to muse over which numerical technologies will be included in the CFD tools that will be employed 30 years in the future. This presentation will describe specific technical problems that have stressed NASA’s traditional CFD methods to their breaking point and will link these issues to the Agency’s top Aerosciences technical challenges. The discussion will then shift to the characteristics of future CFD solvers that will be required to attack these challenges and how these characteristics differ from the present state-of-the- art. High-order methods certainly appear to have a place in the development of future CFD tools and some of the physical characteristics of our most challenging problems suggest that high-order methods are the only way to effectively solve them. But there are some relatively severe implementation issues that face these methods, particularly in the area of general applicability and robust operation as an engineering tool. Desired characteristics of next-generation CFD solvers will be discussed and the author’s view of which emerging numerical technologies might be employed to address these attributes will also be presented

David M Schuster↗

Freestream turbulence effects on the development of a rotor wake

The effect of freestream turbulence on the development of a three-dimensional wake of a compressor rotor blade was studied experimentally. The turbulence level at the inlet of a rotor was varied systematically using grids upstream of the rotor. The rotor wake was measured with inlet turbulence intensities of 0.5, 3, and 5%. The experimental results indicate that the maxium change in the mean velocity defect is 4% over the range of inlet turbulence levels employed, while the turbulence structure in the wake is altered more substantially. The freestream turbulence effect was also analyzed, numerically, using the modified Reynolds stress closure model. The comparison between numerical prediction and experimental data shows that the freestream turbulence effect can be represented successfully with the turbulence closure model employed in this paper.

Hah, C.↗

The gridded back contact and its effect on solar cell performance

The effect of using gridded back contacts on silicon solar cells made with Al baked surface field (BSF) was experimentally investigated and found to result in a significant increase in conversion efficiency due to higher short-circuit current. An increase of 5 to 8% was noted for thin cells while a smaller increase was observed for standard thickness cells. Several possible explanatory mechanisms are suggested and summarized, including stress relief, incomplete BSF formation, improved back surface recombination velocity, and improved optical reflection off the back silicon-air interface.

Giuliano, M.↗

Turbulent Simulations of Cooling Jets in Crossflow

Turbulent flow simulations were completed of the THX-III single hole film cooling experiment performed at NASA Glenn Research Center (GRC). Two different Large-Eddy Simulation (LES) approaches were completed along with a RANS based study. The first LES approach utilizes the Flux Reconstruction (FR) based finite element code GFR, while the second approach makes use of the compact difference solver FDL3DI. Both LES codes were run on grids larger than 1.5 Billion degrees of freedom (DoF). The RANS simulations used FUN3D and included a grid convergence as a turbulence closure study. A summary of the computational approaches and key findings are presented along with the experimental measurements. Significant progress has been made towards the completion of this work, as all but one simulation have been completed. The preliminary findings have shown the LES approaches have excellent agreement with each other, and closely follow the trends of the experiment. The Reynolds stress RANS closure, while the best of the RANS results, had difficulty matching the results of the LES and experimental measurements.

High-Order Methods↗

User's Manual and Final Report for Hot-SMAC GUI Development

A new software package called Higher Order Theory-Structural/Micro Analysis Code (HOT-SMAC) has been developed as an effective alternative to the finite element approach for Functionally Graded Material (FGM) modeling. HOT-SMAC is a self-contained package including pre- and post-processing through an intuitive graphical user interface, along with the well-established Higher Order Theory for Functionally Graded Materials (HOTFGM) thermomechanical analysis engine. This document represents a Getting Started/User's Manual for HOT-SMAC and a final report for its development. First, the features of the software are presented in a simple step-by-step example where a HOT-SMAC model representing a functionally graded material is created, mechanical and thermal boundary conditions are applied, the model is analyzed and results are reviewed. In a second step-by-step example, a HOT-SMAC model of an actively cooled metallic channel with ceramic thermal barrier coating is built and analyzed. HOT-SMAC results from this model are compared to recently published results (NASA/TM-2001-210702) for two grid densities. Finally, a prototype integration of HOTSMAC with the commercially available HyperSizer(R) structural analysis and sizing software is presented. In this integration, local strain results from HyperSizer's structural analysis are fed to a detailed HOT-SMAC model of the flange-to-facesheet bond region of a stiffened panel. HOT-SMAC is then used to determine the peak shear and peel (normal) stresses between the facesheet and bonded flange of the panel and determine the "free edge" effects.

Yarrington, Phil↗

Evaluation ofthe Middle East and North Africa Land Data Assimilation System

The Middle East and North Africa (MENA) region is dominated by dry, warm deserts, areas of dense population, and inefficient use of fresh water resources. Due to the scarcity, high intensity, and short duration of rainfall in the MENA, the region is prone to hydro climatic extremes that are realized by devastating floods and times of drought. However, given its widespread water stress and the considerable demand for water, the MENA remains relatively poorly monitored. This is due in part to the shortage of meteorological observations and the lack of data sharing between nations. As a result, the accurate monitoring of the dynamics of the water cycle in the MENA is difficult. The Land Data Assimilation System for the MENA region (MENA LDAS) has been developed to provide regional, gridded fields of hydrological states and fluxes relevant for water resources assessments. As an extension of the Global Land Data Assimilation System (GLDAS), the MENA LDAS was designed to aid in the identification and evaluation of regional hydrological anomalies by synergistically combining the physically-based Catchment Land Surface Model (CLSM) with observations from several independent data products including soil-water storage variations from the Gravity Recovery and Climate Experiment (GRACE) and irrigation intensity derived from the Moderate Resolution Imaging Spectroradiometer (MODIS). In this fashion, we estimate the mean and seasonal cycle of the water budget components across the MENA.

Bolten, John D.↗

Turbulence Modeling for Shock Wave/Turbulent Boundary Layer Interactions

Accurate aerodynamic computational predictions are essential for the safety of space vehicles, but these computations are of limited accuracy when large pressure gradients are present in the flow. The goal of the current project is to improve the state of compressible turbulence modeling for high speed flows with shock wave / turbulent boundary layer interactions (SWTBLI). Emphasis will be placed on models that can accurately predict the separated region caused by the SWTBLI. These flows are classified as nonequilibrium boundary layers because of the very large and variable adverse pressure gradients caused by the shock waves. The lag model was designed to model these nonequilibrium flows by incorporating history effects. Standard one- and two-equation models (Spalart Allmaras and SST) and the lag model will be run and compared to a new lag model. This new model, the Reynolds stress tensor lag model (lagRST), will be assessed against multiple wind tunnel tests and correlations. The basis of the lag and lagRST models are to preserve the accuracy of the standard turbulence models in equilibrium turbulence, when the Reynolds stresses are linearly related to the mean strain rates, but create a lag between mean strain rate effects and turbulence when nonequilibrium effects become important, such as in large pressure gradients. The affect this lag has on the results for SWBLI and massively separated flows will be determined. These computations will be done with a modified version of the OVERFLOW code. This code solves the RANS equations on overset grids. It was used for this study for its ability to input very complex geometries into the flow solver, such as the Space Shuttle in the full stack configuration. The model was successfully implemented within two versions of the OVERFLOW code. Results show a substantial improvement over the baseline models for transonic separated flows. The results are mixed for the SWBLI assessed. Separation predictions are not as good as the baseline models, but the over prediction of the peak heat flux downstream of the reattachment shock that plagues many models is reduced.

Lillard, Randolph P.↗

A Numerical Evaluation of Icing Effects on a Natural Laminar Flow Airfoil

As a part of CFD code validation efforts within the Icing Branch of NASA Glenn Research Center, computations were performed for natural laminar flow (NLF) airfoil, NLF-0414. with 6 and 22.5 minute ice accretions. Both 3-D ice castings and 2-D machine-generated ice shapes were used in wind tunnel tests to study the effects of natural ice is well as simulated ice. They were mounted in the test section of the Low Turbulence Pressure Tunnel (LTPT) at NASA Langley that the 2-dimensionality of the flow can be maintained. Aerodynamic properties predicted by computations were compared to data obtained through the experiment by the authors at the LTPT. Computations were performed only in 2-D and in the case of 3-D ice, the digitized ice shape obtained at one spanwise location was used. The comparisons were mainly concentrated on the lift characteristics over Reynolds numbers ranging from 3 to 10 million and Mach numbers ranging from 0.12 to 0.29. WIND code computations indicated that the predicted stall angles were in agreement with experiment within one or two degrees. The maximum lift values obtained by computations were in good agreement with those of the experiment for the 6 minute ice shapes and the minute 3-D ice, but were somewhat lower in the case of the 22.5 minute 2-D ice. In general, the Reynolds number variation did not cause much change in the lift values while the variation of Mach number showed more change in the lift. The Spalart-Allmaras (S-A) turbulence model was the best performing model for the airfoil with the 22.5 minute ice and the Shear Stress Turbulence (SST) turbulence model was the best for the airfoil with the 6 minute ice and also for the clean airfoil. The pressure distribution on the surface of the iced airfoil showed good agreement for the 6 minute ice. However, relatively poor agreement of the pressure distribution on the upper surface aft of the leading edge horn for the 22.5 minute ice suggests that improvements are needed in the grid or turbulence models.

Chung, James J.↗

PAB3D Simulations for the CAWAPI F-16XL

Numerical simulations of the flow around F-16XL are performed as a contribution to the Cranked Arrow Wing Aerodynamic Project International (CAWAPI) using the PAB3D CFD code. Two turbulence models are used in the calculations: a standard k-! model, and the Shih-Zhu-Lumley (SZL) algebraic stress model. Seven flight conditions are simulated for the flow around the F-16XL where the free stream Mach number varies from 0.242 to 0.97. The range of angles of attack varies from 0° to 20°. Computational results, surface static pressure, boundary layer velocity profiles, and skin friction are presented and compared with flight data. Numerical results are generally in good agreement with flight data, considering that only one grid resolution is utilized for the different flight conditions simulated in this study. The ASM results are closer to the flight data than the k-ε model results. The ASM predicted a stronger primary vortex, however, the origin of the vortex and footprint is approximately the same as in the k-ε predictions.

Flight Data↗