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 235 records · Page 13

Assessment of Turbulence Models for a Single-Injector Cooling Flow

This work assesses the ability of Reynolds-averaged Navier-Stokes turbulence models to predict a fundamental thermal cooling flow configuration. This geometry consists of a plate with a single cooling hole positioned downstream of a heated exhaust nozzle. The linear Shear Stress Transport (SST), non-linear Quadratic Constitutive Relation (QCR), and differential SSG/LRR Reynolds Stress Model (RSM) turbulence models were investigated for flow conditions with a blowing ratio of 1.0 and 2.0. Despite significant differences in turbulence model complexity, the results miss important details observed in the data and large eddy simulation (LES). Injection penetration is over-predicted. Mixing and centerline wake recovery is too slow. Temperature predictions are worse than velocity, suggesting that perhaps the simple gradient diffusion approximation used to model the turbulent heat flux may be to blame. These observations are put into context using other studies to identify shortcomings in the modeling of both the turbulent stress and turbulent heat flux.

Film Cooling↗

Assessment of Turbulence Models for a Single-Injector Cooling Flow

This work assesses the ability of Reynolds-averaged Navier-Stokes turbulence models to predict a fundamental thermal cooling flow configuration. This geometry consists of a plate with a single cooling hole positioned downstream of a heated exhaust nozzle. The linear Shear Stress Transport (SST), non-linear Quadratic Constitutive Relation (QCR), and differential SSG/LRR Reynolds Stress Model (RSM) turbulence models were investigated for flow conditions with a blowing ratio of 1.0 and 2.0. Despite significant differences in turbulence model complexity, the results miss important details observed in the data and large eddy simulation (LES). Injection penetration is over-predicted. Mixing and centerline wake recovery is too slow. Temperature predictions are worse than velocity, suggesting that perhaps the simple gradient diffusion approximation used to model the turbulent heat flux may be to blame. These observations are put into context using other studies to identify shortcomings in the modeling of both the turbulent stress and turbulent heat flux.

Film Cooling↗

Demonstration of Capability to Simulate Particle Irregular Shape and Poly-Disperse Mixtures Within Lunar Lander Plume-Surface Interaction

Plume-Surface Interaction (PSI) between lander engine plumes and regolith soil creates hazards in obscuration and contamination by particle clouds, high-energy ejecta streams, and landing area cratering damage. The MSFC Fluid Dynamics Branch is developing simulation tools to offer a predictive PSI capability to NASA customers such as the Human Lander System (HLS) and Commercial Lunar Payload Services (CLPS). The Gas-Granular Flow Solver (GGFS) is the main application tool for coupled gas-particle two-phase flow simulations to predict the range of PSI effects from onset of surface erosion to deep crater formation. GGFS features an Eulerian-Eulerian modeling approach, treating both gas and granular material as interacting continuum phases. Modeling the lunar regolith granular material fluidic characteristics poses special challenges due to complex particle shapes and mixture composition. The lunar regolith is poorly sorted with broad particle size distributions and large fines content. It has significant cohesion, due to interlocking jagged particle shapes. Eulerian granular material flow modeling requires closure formulations for the granular material constitutive models (stress, friction, collisional and kinetic energy dissipation, drag, etc.). While closure models for spherical particles are available from particle kinetic theory, closure models for realistic non-spherical particles must be extracted from unit physics Discrete Element Model (DEM) particle interaction simulations and provided in the form of tabular datasets. The effects of particle irregular shape (non-spherical shape factors, angular particle surface roughness, and interlocking features) are simulated by approximating the particle features in the form of grouped elemental spheres to form composite particles in the DEM simulations. The effects of the wide range of regolith mixture particle sizes and the strong effects of the presence of the small particle sizes results in high cohesion and low porosity of the regolith mixture. The range of particle sizes is simulated by binning the particle sizes into an appropriate finite number of particle-size species and solving the problem as a species mixture. Combining these two modeling approaches enables simulations to capture both, the contributions of the irregular particle shape and the particle size distribution. The integration and maturation of the DEM-based constitutive model database generation process and poly-disperse mixture binning approach into the GGFS simulation framework are proceeding under funding by the NASA Game Changing Development program. The status of current capabilities will be presented in comparisons of crater characteristics resulting for spherical and irregular shape particles, and for mono-, bi-, and tri-disperse mixture simulations of Apollo LM plume-surface interaction. The computational results confirm the significance of including the particle shape and mixture effects. Going forward plans for the full implementation of the general poly-disperse regolith modeling capability and maturation towards NASA project application readiness under the GCD program will be presented.

Peter A Liever↗

Demonstration of Capability to Simulate Particle Irregular Shape and Poly-Disperse Mixtures Within Lunar Lander Plume-Surface Interaction

Plume-Surface Interaction (PSI) between lander engine plumes and regolith soil creates hazards in obscuration and contamination by particle clouds, high-energy ejecta streams, and landing area cratering damage. The MSFC Fluid Dynamics Branch is developing simulation tools to offer a predictive PSI capability to NASA customers such as the Human Lander System (HLS) and Commercial Lunar Payload Services (CLPS). The Gas-Granular Flow Solver (GGFS) is the main application tool for coupled gas-particle two-phase flow simulations to predict the range of PSI effects from onset of surface erosion to deep crater formation. GGFS features an Eulerian-Eulerian modeling approach, treating both gas and granular material as interacting continuum phases. Modeling the lunar regolith granular material fluidic characteristics poses special challenges due to complex particle shapes and mixture composition. The lunar regolith is poorly sorted with broad particle size distributions and large fines content. It has significant cohesion, due to interlocking jagged particle shapes. Eulerian granular material flow modeling requires closure formulations for the granular material constitutive models (stress, friction, collisional and kinetic energy dissipation, drag, etc.). While closure models for spherical particles are available from particle kinetic theory, closure models for realistic non-spherical particles must be extracted from unit physics Discrete Element Model (DEM) particle interaction simulations and provided in the form of tabular datasets. The effects of particle irregular shape (non-spherical shape factors, angular particle surface roughness, and interlocking features) are simulated by approximating the particle features in the form of grouped elemental spheres to form composite particles in the DEM simulations. The effects of the wide range of regolith mixture particle sizes and the strong effects of the presence of the small particle sizes results in high cohesion and low porosity of the regolith mixture. The range of particle sizes is simulated by binning the particle sizes into an appropriate finite number of particle-size species and solving the problem as a species mixture. Combining these two modeling approaches enables simulations to capture both, the contributions of the irregular particle shape and the particle size distribution. The integration and maturation of the DEM-based constitutive model database generation process and poly-disperse mixture binning approach into the GGFS simulation framework are proceeding under funding by the NASA Game Changing Development program. The status of current capabilities will be presented in comparisons of crater characteristics resulting for spherical and irregular shape particles, and for mono-, bi-, and tri-disperse mixture simulations of Apollo LM plume-surface interaction. The computational results confirm the significance of including the particle shape and mixture effects. Going forward plans for the full implementation of the general poly-disperse regolith modeling capability and maturation towards NASA project application readiness under the GCD program will be presented.

Peter A. Liever↗

Bending effects of unsymmetric adhesively bonded composite repairs on cracked aluminum panels

The bending effects of unsymmetrically bonded composite repairs on cracked aluminum panels were quantified using a plate linear finite element model. Stress intensity factors and strain energy release rates were obtained from the model twice, once with out-of-plane displacement suppressed and another time without these restrictions. Several configurations were examined, crack growth stability was identified, and the effect of a debond was considered. The maximum stress intensity factor was also analyzed. Previous work by other authors was found to underpredict the bending effect.

Arendt, Cory↗

Model consistency in large eddy simulation of turbulent channel flows

Combinations of filters and subgrid scale stress models for large eddy simulation of the Navier-Stokes equations are examined by a priori tests and numerical simulations. The structure of the subgrid scales is found to depend strongly on the type of filter used, and consistency between model and filter is essential to ensure accurate results. The implementation of consistent combinations of filter and model gives more accurate turbulence statistics than those obtained in previous investigations in which the models were chosen independently from the filter. Results and limitations of the a priori test are discussed. The effect of grid refinement is also examined.

Piomelli, Ugo↗

Turbulent heat transfer in corrugated-wall channels with and without fins

A numerical study is performed examining flow and heat transfer characteristics in a channel with periodically corrugated walls. The complexity of the flow in this type of channel is demonstrated by such phenomena as flow impingement on the walls, separation at the bend corners, flow reattachment, and flow recirculation. Because of the strong nonisotropic nature of the turbulent flow in the channel, the full Reynolds-stress model was employed for the evaluation of turbulence quantities. Computations are made for several different corrugation periods and for different Reynolds numbers. The results computed by using the present model show excellent agreement with experimental data for mean velocities, the Reynolds stresses, and average Nusselt numbers. The study was further extended to a channel flow where fins are inserted at bends in the channel. It was observed that the insertion of fins in the flow passage has a visible effect on flow patterns and skin friction along the channel wall.

Amano, R. S.↗

Horizontal stresses induced by vertical processes in planetary lithospheres

Understanding the state of stress in the elastic lithosphere is of fundamental importance for planetary geophysics, as it is the link between the observed geologic structures on the surface and the processes which form and modify these structures. As such, it can provide valuable constraints for the difficult problem of determining interior structure and processes. On the Earth, most large scale, organized deformation can be related to lateral tectonics associated with plate dynamics; however, the tectonics on many extraterrestrial bodies (such as the Moon, Mars, and most of the outer-planet satellites) appears to be primarily vertical in nature, and the horizontal stresses induced by vertical motions and loads are expected to dominate the deformation of their lithospheres. The largest stress contributions from vertical loading come from the flexure of the lithosphere, which induces both bending moments and membrane stresses. We are concerned here only with nonflexural changes in the state of stress induced by processes such as sedimentary and volcanic deposition, erosional denudation, and changes in the thermal gradient that induce uplift or subsidence. This analysis is important both for evaluating stresses for specific regions in which the vertical stress history can be estimated, as well as for applying the proper loading conditions to global stress models. It is also of interest for providing a reference state of stress for interpreting stress measurements in the crust of the Earth.

Banerdt, W. B.↗

Comparison of Measured and Block Structured Simulations for the F-16XL Aircraft

This article presents a comparison of the predictions of three RANS codes for flight conditions of the F-16XL aircraft which feature vortical flow. The three codes, ENSOLV, PMB and PAB3D, solve on structured multi-block grids. Flight data for comparison was available in the form of surface pressures, skin friction, boundary layer data and photographs of tufts. The three codes provided predictions which were consistent with expectations based on the turbulence modelling used, which was k- , k- with vortex corrections and an Algebraic Stress Model. The agreement with flight data was good, with the exception of the outer wing primary vortex strength. The confidence in the application of the CFD codes to complex fighter configurations increased significantly through this study.

Boelens, O. J.↗

Model consistency in the large eddy simulation of turbulent channel flows

Various combinations of filters and subgrid scale stress models for large eddy simulation of the Navier-Stokes equations are studied by a priori tests and numerical simulations. Consistency between model and filter is found to be essential to ensure accurate results. Results and limitations of the a priori test are discussed. The effect of grid refinement is also examined.

Moin, Parviz↗

On sublayer streaks

Turbulent sublayer streaks are studied with the aid of a simplified theoretical model. In this the nonlinear activity is assumed to be intermittent and to act locally in space during a very short initial time so as to set up the initial conditions for the subsequent linear and inviscid evolution of the resulting three-dimensional flow disturbance. The mean shear flow is taken as a parallel one and a correction for the long-tern effects of viscosity is applied. A model for the initial nonlinear phase is chosen to represent the local Reynolds stresses that would be produced by a patch of local inflectional instability. The streamwise dimension of the resulting eddy is found to grow linearly with time in accordance with the algebraic instability mechanism (Landahl 1980). The associated Reynolds shear stress is expressible in a simple manner in terms of the liftup of the fluid elements and is suggestive of an algebraic-type Reynolds stress model similar to, but not identical to, that of Prandtl's (1925) mixing-length theory.

Landahl, M. T.↗

On prediction of equilibrium states in homogeneous compressible turbulence

Direct numerical simulations of compressible, homogeneous, turbulent shear flows are used to evaluate Reynolds stress models. Three pressure-strain models, which are either linear, quadratic, or cubic in the anisotropy tensor are considered. Dilatational dissipation and pressure-strain correlation models do not correctly capture the compressibility effects seen in the direct simulations. In particular, the increase in the anisotropy of normal stresses and the reduction in the shear stress are not reproduced by any of the models. Also, the use of the incompressible form of the dissipation-rate equation to determine the solenoidal part of the dissipation is found to be questionable.

Abid, Ridha↗

New results on the realizability of Reynolds stress turbulence closures

The realizability of Reynolds stress models in homogeneous turbulence is critically assessed from a theoretical standpoint. It is proven that a well known second-order closure formulated by Shih and Lumley using the strong realizability constraints of Schumann is, in fact, not a realizable model. The problem arises from the failure to properly satisfy the necessary positive second time derivative constraint when a principal Reynolds stress vanishes - a fatal flaw that becomes apparent when the non-analytic terms in their model are made single-valued as required on physical grounds. It is furthermore shown that the centrifugal acceleration generated by rotations of the principal axes of the Reynolds stress tensor can make the second derivative singular at the most extreme limits of realizable turbulence. This previously overlooked effect appears to make it impossible to identically satisfy the strong form of realizability in any version of the present generation of second-order closures. On the other hand, models properly formulated to satisfy the weak form of realizability - wherein states of one or two component turbulence are not accessible in finite time are found to be realizable. However, unlike the simpler and more commonly used second order closures, these models can be ill-behaved near the extreme limits of realizable turbulence due to the way that higher-degree nonlinearities are often unnecessarily introduced to satisfy realizability. Illustrative computations of homogeneous shear flows are presented to demonstrate these points which can have important implications for turbulence modeling.

Speziale, Charles G.↗

The Use of DNS in Turbulence Modeling

The use of Direct numerical simulations (DNS) data in developing and testing turbulence models is reviewed. The data is used to test turbulence models at all levels: algebraic, one-equation, two-equation and full Reynolds stress models were tested. Particular examples on the development of models for the dissipation rate equation are presented. Homogeneous flows are used to test new scaling arguments for the various terms in the dissipation rate equation. The channel flow data is used to develop modifications to the equation model that take into account near-wall effects. DNS of compressible flows under mean compression are used in testing new compressible modifications to the two-equation models.

Mansour, Nagi N.↗

Multiaxial and Thermomechanical Fatigue of Materials: A Historical Perspective and Some Future Challenges

Structural materials used in engineering applications routinely subjected to repetitive mechanical loads in multiple directions under non-isothermal conditions. Over past few decades, several multiaxial fatigue life estimation models (stress- and strain-based) developed for isothermal conditions. Historically, numerous fatigue life prediction models also developed for thermomechanical fatigue (TMF) life prediction, predominantly for uniaxial mechanical loading conditions. Realistic structural components encounter multiaxial loads and non-isothermal loading conditions, which increase potential for interaction of damage modes. A need exists for mechanical testing and development verification of life prediction models under such conditions.

Kalluri, Sreeramesh↗

An Active Model Split for Hybrid RANS/LES of Compressible Flows

Hybrid RANS/LES models have consistently shown superior accuracy over RANS in predictions of massively separated flows. However, these models have experienced difficulties with modeled stress depletion, smooth-wall separated flow, and shock-separated flows. The Active Model Split (AMS) hybridization was created to address some of the fundamental shortcomings of traditional hybrid RANS/LES models. While previous work considered incompressible flows, this paper demonstrates how the AMS framework can be extended to address compressible flows. AMS-based hybrid models for the turbulent heat flux and turbulent transport terms in the total energy equation are presented. This compressible hybrid RANS/LES framework is demonstrated on a subsonic channel, a supersonic channel, and a transonic shock-boundary layer interaction. Record includes video of presentation, to view download the mp4 video attached run time of 10 mins 18 secs.

hybrid RANS/LES↗

Turbulence modeling and experiments

The best way of verifying turbulence is to do a direct comparison between the various terms and their models. The success of this approach depends upon the availability of the data for the exact correlations (both experimental and DNS). The other approach involves numerically solving the differential equations and then comparing the results with the data. The results of such a computation will depend upon the accuracy of all the modeled terms and constants. Because of this it is sometimes difficult to find the cause of a poor performance by a model. However, such a calculation is still meaningful in other ways as it shows how a complete Reynolds stress model performs. Thirteen homogeneous flows are numerically computed using the second order closure models. We concentrate only on those models which use a linear (or quasi-linear) model for the rapid term. This, therefore, includes the Launder, Reece and Rodi (LRR) model; the isotropization of production (IP) model; and the Speziale, Sarkar, and Gatski (SSG) model. Which of the three models performs better is examined along with what are their weaknesses, if any. The other work reported deal with the experimental balances of the second moment equations for a buoyant plume. Despite the tremendous amount of activity toward the second order closure modeling of turbulence, very little experimental information is available about the budgets of the second moment equations. Part of the problem stems from our inability to measure the pressure correlations. However, if everything else appearing in these equations is known from the experiment, pressure correlations can be obtained as the closing terms. This is the closest we can come to in obtaining these terms from experiment, and despite the measurement errors which might be present in such balances, the resulting information will be extremely useful for the turbulence modelers. The purpose of this part of the work was to provide such balances of the Reynolds stress and heat flux equations for the buoyant plume.

Shabbir, Aamir↗

A 3-D CFD Analysis of the Space Shuttle RSRM With Propellant Fins @ 1 sec. Burn-Back

In this study 3-D Computational Fluid Dynamic (CFD) runs have been made for the Space Shuttle RSRM using 2 different grids and 4 different turbulent models, which were the Standard KE, the RNG KE, the Realizable KE, and the Reynolds stress model. The RSRM forward segment consists of 11 fins. By taking advantage of the forward fin symmetry only half of one fin along the axis had to be used in making the grid. This meant that the 3-D model consisted of a pie slice that encompassed 1/22nd of the motor circumference and went along the axis of the entire motor. The 3-D flow patterns in the forward fin region are of particular interest. Close inspection of these flow patterns indicate that 2 counter-rotating axial vortices emerge from each submerged solid propellant fin. Thus, the 3-D CFD analysis allows insight into complicated internal motor flow patterns that are not available from the simpler 2-D axi-symmetric studies. In addition, a comparison is made between the 3-D bore pressure drop and the 2-D axi-symmetric pressure drop.

Morstadt, Robert A.↗