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 415 records · Page 23

Research activities at the Center for Modeling of Turbulence and Transition

The main research activities at the Center for Modeling of Turbulence and Transition (CMOTT) are described. The research objective of CMOTT is to improve and/or develop turbulence and transition models for propulsion systems. The flows of interest in propulsion systems can be both compressible and incompressible, three dimensional, bounded by complex wall geometries, chemically reacting, and involve 'bypass' transition. The most relevant turbulence and transition models for the above flows are one- and two-equation eddy viscosity models, Reynolds stress algebraic- and transport-equation models, pdf models, and multiple-scale models. All these models are classified as one-point closure schemes since only one-point (in time and space) turbulent correlations, such as second moments (Reynolds stresses and turbulent heat fluxes) and third moments, are involved. In computational fluid dynamics, all turbulent quantities are one-point correlations. Therefore, the study of one-point turbulent closure schemes is the focus of our turbulence research. However, other research, such as the renormalization group theory, the direct interaction approximation method, and numerical simulations are also pursued to support the development of turbulence modeling.

Shih, Tsan-Hsing↗

Evaluation of Industry Standard Turbulence Models on an Axisymmetric Supersonic Compression Corner

Reynolds-averaged Navier-Stokes computations of a shock-wave/boundary-layer interaction (SWBLI) created by a Mach 2.85 flow over an axisymmetric 30-degree compression corner were carried out. The objectives were to evaluate four turbulence models commonly used in industry, for SWBLIs, and to evaluate the suitability of this test case for use in further turbulence model benchmarking. The Spalart-Allmaras model, Menter's Baseline and Shear Stress Transport models, and a low-Reynolds number k- model were evaluated. Results indicate that the models do not accurately predict the separation location; with the SST model predicting the separation onset too early and the other models predicting the onset too late. Overall the Spalart-Allmaras model did the best job in matching the experimental data. However there is significant room for improvement, most notably in the prediction of the turbulent shear stress. Density data showed that the simulations did not accurately predict the thermal boundary layer upstream of the SWBLI. The effect of turbulent Prandtl number and wall temperature were studied in an attempt to improve this prediction and understand their effects on the interaction. The data showed that both parameters can significantly affect the separation size and location, but did not improve the agreement with the experiment. This case proved challenging to compute and should provide a good test for future turbulence modeling work.

Turbulent Boundary Layer↗

Surface temperatures and glassy state investigations in tribology, part 3

Research related to the development of the limiting shear stress rheological model is reported. Techniques were developed for subjecting lubricants to isothermal compression in order to obtain relevant determinations of the limiting shear stress and elastic shear modulus. The isothermal compression limiting shear stress was found to predict very well the maximum traction for a given lubricant. Small amounts of side slip and twist incorporated in the model were shown to have great influence on the rising portion of the traction curve at low slide-roll ratio. The shear rheological model was also applied to a Grubin-like elastohydrodynamic inlet analysis for predicting film thicknesses when employing the limiting shear stress model material behavior.

Bair, S.↗

An electrochemical model for hot-salt stress-corrosion of titanium alloys

An electrochemical model of hot-salt stress-corrosion cracking of titanium alloys is proposed based on an oxygen-concentration cell. Hydrogen embrittlement is proposed as the direct cause of cracking, the hydrogen being generated as the results of the hydrolysis of complex halides formed at the shielded anode of the electrochemical cell. The model found to be consistent with the diverse observations made both in this study and by many investigators in this field.

Garfinkle, M.↗

Dynamics of coupled ice-ocean system in the marginal ice zone: Study of the mesoscale processes and of constitutive equations for sea ice

This study is aimed at the modelling of mesoscale processed such as up/downwelling and ice edge eddies in the marginal ice zones. A 2-dimensional coupled ice-ocean model is used for the study. The ice model is coupled to the reduced gravity ocean model (f-plane) through interfacial stresses. The constitutive equations of the sea ice are formulated on the basis of the Reiner-Rivlin theory. The internal ice stresses are important only at high ice concentrations (90-100%), otherwise the ice motion is essentially free drift, where the air-ice stress is balanced by the ice-water stress. The model was tested by studying the upwelling dynamics. Winds parallel to the ice edge with the ice on the right produce upwilling because the air-ice momentum flux is much greater that air-ocean momentum flux, and thus the Ekman transport is bigger under the ice than in the open water. The upwelling simulation was extended to include temporally varying forcing, which was chosen to vary sinusoidally with a 4 day period. This forcing resembles successive cyclone passings. In the model with a thin oceanic upper layer, ice bands were formed.

Hakkinen, S.↗

On the mathematical modeling of the Reynolds stress's equations

By considering the Reynolds stress equations as a possible descriptor of complex turbulent fields, pressure-velocity interaction and turbulence dissipation are studied as two of the main physical contributions to Reynolds stress balancing in turbulent flow fields. It is proven that the pressure interaction term contains turbulence generation elements. However, the usual 'return to isotropy' element appears more weakly than in the standard models. In addition, convection-like elements are discovered mathematically, but there is no mathematical evidence that the pressure fluctuations contribute to the turbulent transport mechanism. Calculations of some simple one-dimensional fields indicate that this extra convection, rather than the turbulent transport, is needed mathematically. Similarly, an expression for the turbulence dissipation is developed. The end result is a dynamic equation for the dissipation tensor which is based on the tensorial length scales.

Lin, Avi↗

Finite Element Modeling Used to Study Stress Distribution on the Foot

A method to study the stress distribution inside the forefoot during walking was developed at the Cleveland Clinic Foundation by a researcher from the NASA Glenn Research Center. In this method, a semiautomated process was outlined to create a three-dimensional, patient-specific, finite element model (FEM) of the forefoot using magnetic resonance images (MRI). The images were processed in Matlab using the k-nearest neighbor (k-NN) classification algorithm and Sobel edge detection to separate the different tissue types: bone, skin, fat, and muscle. This information was used to create curves and surfaces that were exported to an FEM preprocessor known as Truegrid. In Truegrid, eight-noded or brick elements were created by using surface mapping. The FEM was processed and postprocessed in Abaqus. Material properties of the models were obtained from past experiments such as fat pad confined compression, skin axial and biaxial tests, muscle in vivo compressive tests, and reference literature (bone properties). Nonlinear (hyperelastic) material models were used for the skin (epidermis and dermis), fat, and muscles; and a linear elastic model was used for the bones. Muscle activation during walking yielded uncertainties in the muscle material model since contracted muscles are stiffer than relaxed muscles. These uncertainties were resolved by performing a sensitivity analysis of the muscle material properties. The original properties were multiplied by arbitrary factors of 2, 3, 0.5, and 0.33. The strain and stress distributions, as well as the locations of peak values, were similar in all cases. The peak contact pressure P obtained for each case varied with respect to the applied factor f as follows:

Morales, Nelson↗

Three dimensional flow field inside compressor rotor, including blade boundary layers

The Reynolds stress equation, pressure strain correlation, and dissipative terms and diffusion are discussed in relation to turbulence modelling using the Reynolds stress model. Algebraic modeling of Reynolds stresses and calculation of the boundary layer over an axial cylinder are examined with regards to the kinetic energy model for turbulence modelling. The numerical analysis of blade and hub wall boundary layers, and an experimental study of rotor blade boundary layer in an axial flow compressor rotor are discussed. The Patankar-Spalding numerical method for two dimensional boundary layers is included.

Galmes, J. M.↗

Updates to PATO’s Thermo-Mechanical Model

The solid mechanics module within the Porous material Analysis Toolbox (PATO) material response code was updated to model the stress contribution from both internal pressure and pyrolysis shrinkage as function of material’s decomposition. Incorporating these physical phenomena, alongside the previously included temperature dependent mechanical properties, thermal expansion, and aerodynamic loads, allows for a more accurate thermo-mechanical behavior modeling of Thermal Protection System (TPS) materials, which is critical to assess spallation risks during atmospheric entry. This poster outlines the workflow for obtaining the necessary PATO inputs, derived from dilatometry experiments, which are used to model the stress field in TPS materials due to thermal expansion and pyrolysis shrinkage. Furthermore, to analyze the relevance of the stress contribution from the internal pressure, a study was conducted using TACOT with different permeability values.

TPS↗

Behavior of a turbulent boundary layer subjected to sudden transverse strain

Data from two experiments on the development of the components of the Reynolds stress tensor after a sudden application of transverse strain are compared. Computations were based on four different turbulence models: a first-order mixing length model, a second-order two-equation eddy viscosity model, and two second-order Reynolds stress models. The second-order models do not produce dramatic improvements over the simple mixing length model. The Reynolds stress models still need development to represent the physics of shear-strained turbulence well. Finally, it is demonstrated that the assumption of a scalar eddy viscosity, often used in engineering calculations of three-dimensional boundary layers, is quite reasonable.

Higuchi, H.↗

Calculation of 3D turbulent jets in crossflow with a multigrid method and a second-moment closure model

A multigrid method is presented for calculating turbulent jets in crossflow. Fairly rapid convergence is obtained with the k-epsilon turbulence model, but computations with a full Reynolds stress turbulence model (RSM) are not yet very efficient. Grid dependency tests show that there are slight differences between results obtained on the two finest grid levels. Computations using the RSM are significantly different from those with k-epsilon model and compare better to experimental data. Some work is still required to improve the efficiency of the computations with the RSM.

Demuren, A. O.↗

Transition Prediction for the 30P30N Airfoil and the CRM-HL using FUN3D

The ability to accurately model transition over high-lift configurations can lead to improved agreement between stationary Reynolds-Averaged Navier-Stokes (RANS) computations and experimental measurements. In this work, selected transport-equation-based RANS models for transition prediction in the NASA FUN3D flow solver are used to simulate flow over the 30P30N three-element airfoil and the High-Lift Common Research Model (CRM-HL). The primary RANS-based transition model is the Langtry-Menter (LM) 𝜸-Re𝜽𝒕 model coupled with the Menter shear-stress transport model. Additionally, the existing capability for transition prediction based on linear stability correlations is used to allow the modeling of transition over multiple elements of the 30P30N configuration. This approach pairs an algebraic transition model with the Spalart-Allmaras (SA) turbulence model that relies on an imposed transition location in FUN3D and the parabolized stability equations as implemented in the LASTRAC code. Predictions based on these different transition prediction approaches and models are compared with each other, experimental data, and previous simulations by evaluating the lift and drag coefficients, the surface-pressure/skin-friction distributions, and the locations of transition onset. Results are presented for multiple angles of attack for the 30P30N airfoil and CRM-HL. For the CRM-HL, a grid resolution study is performed for the LM transition model and the SA turbulence model.

High Lift↗

Thermoelastic stress in oceanic lithosphere due to hotspot reheating

The effect of hotspot reheating on the intraplate stress field is investigated by modeling the three-dimensional thermal stress field produced by nonuniform temperature changes in an elastic plate. Temperature perturbations are calculated assuming that the lithosphere is heated by a source in the lower part of the thermal lithosphere. A thermal stress model for the elastic lithosphere is calculated by superposing the stress fields resulting from temperature changes in small individual elements. The stress in an elastic plate resulting from a temperature change in each small element is expressed as an infinite series, wherein each term is a source or an image modified from a closed-from half-space solution. The thermal stress solution is applied to midplate swells in oceanic lithosphere with various thermal structures and plate velocities. The results predict a stress field with a maximum deviatoric stress on the order of 100 MPa covering a broad area around the hotspot plume. The predicted principal stress orientations show a complicated geographical pattern, with horizontal extension perpendicular to the hotspot track at shallow depths and compression along the track near the bottom of the elastic lithosphere.

Zhu, Anning↗

Anisotropic constitutive modeling for nickel base single crystal superalloy Rene N4 at 982 C

A back stress/drag stress constitutive model based on a crystallographic approach to model single crystal anisotropy is presented. Experimental results demonstrated the need for the back stress variable in the inelastic flow equations. Experimental findings suggested that back stress is orientation dependent and controls both strain hardening and recovery characteristics. Due to the observed stable fatigue loops at 1800 F, drag stress is considered constant for this temperature. The constitutive model operated with constraints determined only from tensile data was extensively tested from simple tensile and fatigue to complicated strain hold tests. The model predicted very well under those conditions.

Stouffer, D. C.↗

Comparison of Design Tools for Stress Analysis of Adhesively Bonded Joints

Analytical or semi-analytical models for stress analysis have long been a part of initial adhesive bond design and sizing. Even with the rise of general finite element software and methods, these design models have still remained a preferred method for fast and simple joint analysis. While these methods can yield fairly representative results, the models are constructed on the foundation of geometrical and material simplifications or assumptions that allow closed form or semi-closed form solutions. This study outlines the major differences in basic assumptions for three common design software packages under use at NASA, and shows the ramifications these assumptions in a few exemplar bonded joints.

Composites↗

A new approach to turbulence modeling

A new approach to Reynolds averaged turbulence modeling is proposed which has a computational cost comparable to two equation models but a predictive capability approaching that of Reynolds stress transport models. This approach isolates the crucial information contained within the Reynolds stress tensor, and solves transport equations only for a set of 'reduced' variables. In this work, Direct Numerical Simulation (DNS) data is used to analyze the nature of these newly proposed turbulence quantities and the source terms which appear in their respective transport equations. The physical relevance of these quantities is discussed and some initial modeling results for turbulent channel flow are presented.

Perot, B.↗

Progress Towards Computational Method for Circulation Control Airfoils

The compressible Reynolds-averaged Navier-Stokes equations are solved for circulation control airfoil flows. Numerical solutions are computed with both structured and unstructured grid solvers. Several turbulence models are considered, including the Spalart-Allmaras model with and without curvature corrections, the shear stress transport model of Menter, and the k-enstrophy model. Circulation control flows with jet momentum coefficients of 0.03, 0.10, and 0.226 are considered. Comparisons are made between computed and experimental pressure distributions, velocity profiles, Reynolds stress profiles, and streamline patterns. Including curvature effects yields the closest agreement with the measured data.

Swanson, R. C.↗