Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Stress Modeling”

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 289 records · Page 16

Analytical methods for the development of Reynolds-stress closures in turbulence

The derivation of Reynolds-stress models for viscous incompressible turbulent flow on the basis of the Navier-Stokes and continuity equations is explored in an analytical review. The formulation of the basic equations is outlined, and particular attention is given to zero-equation and one-equation models based on eddy viscosity, two-equation (k-l, k-epsilon, and k-omega) models, and second-order closure models. The superior performance of the second-order models is demonstrated by numerical results for (1) the return to isotropy from anisotropic homogeneous turbulence, (2) homogeneous turbulent shear flow in a rotating frame, (3) and fully developed inhomogeneous turbulent channel flow in a rotating frame. The inherent limitations of the Reynolds-stress approach and areas for further improvement are discussed.

Speziale, Charles G.↗

Computation and turbulence modeling for three-dimensional boundary layers including turbomachinery rotor flows

A method is developed for predicting the behavior of three-dimensional, turbulent boundary layers occurring in internal flows, including those on turbomachinery rotor blades. These boundary layers are complex, turbulent, and subject to Coriolis and centrifugal forces. The major thrust of this paper is the development and use of an algebraic Reynolds stress model (ARSM) that captures the changes in turbulent flow structure arising from curvature, rotation, and three dimensionality. The prediction of pressure-driven secondary flow agrees well with the measured data, and all three turbulence models (k-epsilon, algebraic eddy viscosity, and ARSM) show the same level of agreement. The prediction of boundary-layer development on rotor blades shows much better agreement with measurements with the ARSM. It is essential to employ higher-order turbulence models to capture the effects of rotation, curvature, and three dimensionality on boundary layers in turbomachinery.

Zhang, J.↗

Calculation of three-dimensional compressible turbulent boundary layers

From comparisons of high speed data with low speed closure procedures using variable mean density, there does not appear to be any appreciable influence of compressibility upon turbulent shear stress modeling in compressible turbulent boundary layers, even for extreme cases such as Mach 14 to 20 with a change in density across the layer of up to a factor of 100. Other evidence of apparent lack of compressibility caused new physics which may alter the shear stress for the compressible boundary layer cases including: (1) fluctuation Mach number was generally less than 1; (2) the shear stress distribution through the boundary layer was not a function of Mach number for zero pressure gradient flows; (3) the Morkovin hypothesis was valid up to Mach 5 (based on fluctuation data); (4) profile N power was not a function of Mach number, at least up to Mach 10; and (5) the nondimensional burst period was approximately the same as that for low speed.

Source record↗

Numerical computation of aerodynamics and heat transfer in a turbine cascade and a turn-around duct using advanced turbulence models

The objective of this research is to develop turbulence models to predict the flow and heat transfer fields dominated by the curvature effect such as those encountered in turbine cascades and turn-around ducts. A Navier-Stokes code has been developed using an explicit Runge-Kutta method with a two layer k-epsilon/ARSM (Algebraic Reynolds Stress Model), Chien's Low Reynolds Number (LRN) k-epsilon model and Coakley's LRN q-omega model. The near wall pressure strain correlation term was included in the ARSM. The formulation is applied to Favre-averaged N-S equations and no thin-layer approximations are made in either the mean flow or turbulence transport equations. Anisotropic scaling of artificial dissipation terms was used. Locally variable timestep was also used to improve convergence. Detailed comparisons were made between computations and data measured in a turbine cascade by Arts et al. at Von Karman Institute. The surface pressure distributions and wake profiles were predicted well by all the models. The blade heat transfer is predicted well by k-epsilon/ARSM model, as well as the k-epsilon model. It's found that the onset of boundary layer transition on both surfaces is highly dependent upon the level of local freestream turbulence intensity, which is strongly influenced by the streamline curvature. Detailed computation of the flow in the turn around duct has been carried out and validated against the data by Monson as well as Sandborn. The computed results at various streamwise locations both on the concave and convex sides are compared with flow and turbulence data including the separation zone on the inner well. The k-epsilon/ARSM model yielded relatively better results than the two-equation turbulence models. A detailed assessment of the turbulence models has been made with regard to their applicability to curved flows.

Lakshminarayana, B.↗

Geophysical models for the formation and evolution of coronae on Venus

The proposition that Venusian coronae form over sites of mantle upwelling and are modified by subsequent gravitational relaxation is examined using two geophysical models to determine whether and under what conditions these mechanisms can produce the topography and tectonics exhibited by coronae in the Magellan altimetry data and radar images. It is shown that mantle diapirism can produce the domical topography of novae, which may be coronae in the earliest stage of formation. The model stresses induced at the surface by a mantle diapir imply the formation of radially oriented extensional fracturing as observed in novae. The novae dimensions indicate that the diapirs responsible for them are smaller than about 100 km in radius and that the elastic lithosphere is less than 32 km thick. A flattened diapir at the top of the mantle is modeled and shown to result in plateaulike uplift. The volume of the flattened model diapir is similar to that of the spherical diapirs derived for novae.

Janes, Daniel M.↗

Juncture Flow Computations using kL-Based Turbulence Models

The development and implementation of kL-based Reynolds-Average Navier-Stokes (RANS) two-equation turbulence models are reported in this paper. The kL model is based on Abdol-Hamid's closure and Menter's modification to Rotta's two-equation model. Rotta showed that a reliable transport equation can be formed from the turbulence length scale L, and the turbulence kinetic energy k. Rotta's kL equation is well suited for a term-by-term modeling and displays useful features compared to other scale formulations. One of the important differences is the inclusion of higher-order velocity derivatives in the source terms of the scale equation. This can enhance the ability of RANS solvers to simulate unsteady flows in URANS mode. The k-kL scheme has been modified to include the Algebraic Reynolds Stress Model (ARSM) and the Quadratic Constitutive Relation (QCR) as nonlinear models. These models have gone through extensive validations using two-dimensional benchmark problems. The present study documents the application of the k-kL schemes to simulate flow around the Juncture Flow Model (JFM). The k-kL prediction results show generally good comparisons with measurements. The results from this formulation are similar to, or better than results using the SA-RC-QCR2000 two-equation turbulence model. The ARSM and the QCR formulations of k-kL show promise with a similar level of computational resources as basic two-equation turbulence models.

K S Abdol-Hamid↗

Developments in the computation of turbulent boundary layers

Computational techniques applicable to turbulent boundary layers are classified into solutions of Reynolds-averaged equations, in which all the effects of the turbulence are modelled, and solutions of three-dimensional, time dependent Navier-Stokes equations, in which the large eddies are calculated and only the turbulence at scales smaller than the computational mesh spacings has to be modelled. Current computation costs place engineering computations in the first of these categories; large eddy simulations are appropriate currently for special studies of the dynamical processes of turbulence in idealized flow fields. It is shown that the two methods are interrelated and that each can gain from advances in the other. The degree of success of a pair of increasingly complex Reynolds stress models to broaden their range of applicability is examined through comparisons with experimental data for a variety of flow conditions. An example of a large-eddy simulation is presented, compared with experimental results, and used to evaluate the models for pressure rate-of-strain correlations and dissipation in the Reynolds-averaged equations.

Rubesin, M. W.↗

Modelling the transitional boundary layer

Recent developments in the modelling of the transition zone in the boundary layer are reviewed (the zone being defined as extending from the station where intermittency begins to depart from zero to that where it is nearly unity). The value of using a new non-dimensional spot formation rate parameter, and the importance of allowing for so-called subtransitions within the transition zone, are both stressed. Models do reasonably well in constant pressure 2-dimensional flows, but in the presence of strong pressure gradients further improvements are needed. The linear combination approach works surprisingly well in most cases, but would not be so successful in situations where a purely laminar boundary layer would separate but a transitional one would not. Intermittency-weighted eddy viscosity methods do not predict peak surface parameters well without the introduction of an overshooting transition function whose connection with the spot theory of transition is obscure. Suggestions are made for further work that now appears necessary for developing improved models of the transition zone.

Narasimha, R.↗

Irrigation in the Earth System

Irrigation accounts for ~70% of global freshwater withdrawals and ~90% of consumptive water use, driving myriad Earth system impacts. In this Review, we summarize how irrigation currently impacts key components of the Earth system. Estimates suggest that more than 3.6 million km 2 of currently irrigated land, with hot spots in the intensively cultivated US High Plains, California Central Valley, Indo-Gangetic Basin and northern China. Process-based models estimate that ~2,700 ± 540 km 3 irrigation water is withdrawn globally each year, broadly consistent with country-reported values despite these estimates embedding substantial uncertainties. Expansive irrigation has modified surface energy balance and biogeochemical cycling. A shift from sensible to latent heat fluxes, and resulting land–atmosphere feedbacks, generally reduce regional growing season surface temperatures by ~1–3 °C. Irrigation can ameliorate temperature extremes in some regions, but conversely exacerbates moist heat stress. Modelled precipitation responses are more varied, with some intensive cropping regions exhibiting suppressed local precipitation but enhanced precipitation downstream owing to atmospheric circulation interactions. Additionally, irrigation could enhance cropland carbon uptake; however, it can also contribute to elevated methane fluxes in rice systems and mobilize nitrogen loading to groundwater. Cross-disciplinary, integrative research efforts can help advance understanding of these irrigation–Earth system interactions, and identify and reduce uncertainties, biases and limitations.

Climate sciences↗

Computational Analysis of a Chevron Nozzle Uniquely Tailored for Propulsion Airframe Aeroacoustics

A computational flow field and predicted jet noise source analysis is presented for asymmetrical fan chevrons on a modern separate flow nozzle at take off conditions. The propulsion airframe aeroacoustic asymmetric fan nozzle is designed with an azimuthally varying chevron pattern with longer chevrons close to the pylon. A baseline round nozzle without chevrons and a reference nozzle with azimuthally uniform chevrons are also studied. The intent of the asymmetric fan chevron nozzle was to improve the noise reduction potential by creating a favorable propulsion airframe aeroacoustic interaction effect between the pylon and chevron nozzle. This favorable interaction and improved noise reduction was observed in model scale tests and flight test data and has been reported in other studies. The goal of this study was to identify the fundamental flow and noise source mechanisms. The flow simulation uses the asymptotically steady, compressible Reynolds averaged Navier-Stokes equations on a structured grid. Flow computations are performed using the parallel, multi-block, structured grid code PAB3D. Local noise sources were mapped and integrated computationally using the Jet3D code based upon the Lighthill Acoustic Analogy with anisotropic Reynolds stress modeling. In this study, trends of noise reduction were correctly predicted. Jet3D was also utilized to produce noise source maps that were then correlated to local flow features. The flow studies show that asymmetry of the longer fan chevrons near the pylon work to reduce the strength of the secondary flow induced by the pylon itself, such that the asymmetric merging of the fan and core shear layers is significantly delayed. The effect is to reduce the peak turbulence kinetic energy and shift it downstream, reducing overall noise production. This combined flow and noise prediction approach has yielded considerable understanding of the physics of a fan chevron nozzle designed to include propulsion airframe aeroacoustic interaction effects.

Massey, Steven J.↗

Unsteady CFD Simulations of a Compression Corner Geometry Using Wall-Modeled LES Methods in Loci/CHEM

Several wall-modeled large eddy simulation (WMLES) methods are tested by simulating an unsteady Mach 2.0 compression corner geometry in Mississippi State’s Loci/CHEM solver. This study is conducted to evaluate the usage and requirements of these WMLES methods for applications regarding fluctuating pressure environments on launch vehicles with computational fluid dynamics (CFD). Two hybrid Reynolds-averaged Navier-Stokes (RANS)-large eddy simulation (LES) methods, Dynamic Hybrid RANS-LES (DHRL) and Improved Delayed Detached Eddy Simulation (IDDES), and one wall-stress-model, the Algebraic Wall Model for Wall-Modeled LES (AWMLES), are tested on varying grid and timestep refinement levels. These grid and timestep sizes are chosen to test the minimum requirements for successfully running these WMLES methods. The simulations are evaluated based on turbulent boundary layer properties in the developed boundary layer as well as unsteady quantities relating to fluctuating pressure environments in the region of the compression corner. The DHRL method shows good agreement with the comparison wind tunnel data and shows good grid and timestep convergence. The results from the IDDES and WMLES simulations show good agreement for several quantities with some discrepancies regarding others. The results presented in this paper will be used to inform further studies in predicting unsteady environments on higher-complexity geometries.

Bryson Frank↗

The NASA/Industry Design Analysis Methods for Vibrations (DAMVIBS) Program - Boeing helicopters airframe finite element modeling

Finite-element modeling of the airframe vibration of the Army/Boeing CH-47D helicopter is conducted with comparisons to experimental data in an effort to improve the design process. A NASTRAN FEM is developed that is fully representative of the test configuration and includes the support fixture, shakers, and the aircraft/shaker suspension system. The analysis is conducted with specific attention given to the prediction of reasonable forced amplitudes throughout the airframe. Reasonable correlation is noted between the FEM and experimental results, although improved correlation can be obtained by including more accurate damping values and secondary effects such as stringer shear loading. It is shown that the general stress model does not provide an adequate dynamic analysis on which to base design improvements. A more detailed model is required that emphasizes highly detailed helicopter elements and employs a finer mesh particularly in the description of the mass distribution.

Gabel, R.↗

Static structural tests of a 1.5-meter diameter fabric attached inflatable decelerator

Meridional tape forces, permeability, and change in model contours were measured on the aft half of a 1.5-m-diameter attached inflatable decelerator (AID). Inflation pressures up to 103 kN per square meter and external pressures of 0.14 to 102 kN per square meter were used. The results indicated that the model stresses were near the desired isotensoid condition. Future AID designs should consider both stiffness and strength of the meridional tapes to obtain the optimum division of load between tapes and fabric. Permeability of pressurized fabric structures is a few orders of magnitude higher than that obtained in the standard low-pressure test on material specimens.

Willis, C. M.↗

Computation of turbulent flows

The essential ingredients and the effectiveness of several levels of turbulent-flow partial differential equation models are considered. Zero-equation models use only the partial differential equation (pde) for the mean velocity field and do not employ turbulence pde's. One-equation models make use of an additional pde relating to the turbulence velocity scale. Attention is also given to two-equation models, stress-equation models, and large eddy simulations. Large-eddy simulations are concerned with a three-dimensional time-dependent numerical computation of the large-scale turbulence.

Reynolds, W. C.↗

Some TEM observations of Al2O3 scales formed on NiCrAl alloys

The microstructural development of Al2O3 scales on NiCrAl alloys has been examined by transmission electron microscopy. Voids were observed within grains in scales formed on a pure NiCrAl alloy. Both voids and oxide grains grew measurably with oxidation time at 1100 C. The size and amount of porosity decreased towards the oxide-metal growth interface. The voids resulted from an excess number of oxygen vacancies near the oxidemetal interface. Short-circuit diffusion paths were discussed in reference to current growth stress models for oxide scales. Transient oxidation of pure, Y-doped, and Zr-doped NiCrAl was also examined. Oriented alpha-(Al, Cr)2O3 and Ni(Al, Cr)2O4 scales often coexisted in layered structures on all three alloys. Close-packed oxygen planes and directions in the corundum and spinel layers were parallel. The close relationship between oxide layers provided a gradual transition from initial transient scales to steady state Al2O3 growth.

Smialek, J.↗

Developments in the computation of turbulent boundary layers

Two methods of turbulence computation are discussed in terms of their basic simularities. It is shown that the two methods are interrelated and that each can gain from advances in the other. The degree of success of a pair of increasingly complex Reynolds stress models to broaden their range of applicability is examined through comparison with experimental data for a variety of flow conditions. An example of a large eddy simulation is presented, compared with experimental results, and used to evaluate the models for pressure rate of strain correlation and dissipation in the Reynolds averaged equations.

Rubesin, M. W.↗

Some TEM observations of Al2O3 scales formed on NiCrAl alloys

The microstructural development of Al2O3 scales on NiCrAl alloys has been examined by transmission electron microscopy. Voids have been observed within grains in scales formed on a pure NiCrAl alloy. Both voids and oxide grains grew measurably with oxidation time at 1100 C. The size and amount of porosity decreased towards the oxide-metal growth interface. It was postulated that the voids resulted from an excess number of oxygen vacancies near the oxide-metal interface. Short-circuit diffusion paths were discussed in reference to current growth stress models for oxide scales. Transient oxidations of pure, Y-doped, and Zr-doped NiCrAl was also examined. Oriented alpha-(Al,Cr)2O3 and Ni(Al,Cr)2O4 scales often coexisted in layered structures on all three alloys. Close-packed oxygen planes and directions in the corundum and spinel layers were parallel. The close relationships between oxide layers provided a gradual transition from initial transient scales to steady state Al2O3 growth.

Smialek, J.↗