Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Closure model”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 271 records · Page 15

Fatigue-Life Prediction Methodology Using Small-Crack Theory

This paper reviews the capabilities of a plasticity-induced crack-closure model to predict fatigue lives of metallic materials using 'small-crack theory' for various materials and loading conditions. Crack-tip constraint factors, to account for three-dimensional state-of-stress effects, were selected to correlate large-crack growth rate data as a function of the effective-stress-intensity factor range (delta K(eff)) under constant-amplitude loading. Some modifications to the delta k(eff)-rate relations were needed in the near-threshold regime to fit measured small-crack growth rate behavior and fatigue endurance limits. The model was then used to calculate small- and large-crack growth rates, and to predict total fatigue lives, for notched and un-notched specimens made of two aluminum alloys and a steel under constant-amplitude and spectrum loading. Fatigue lives were calculated using the crack-growth relations and microstructural features like those that initiated cracks for the aluminum alloys and steel for edge-notched specimens. An equivalent-initial-flaw-size concept was used to calculate fatigue lives in other cases. Results from the tests and analyses agreed well.

Newmann, James C., Jr.↗

Numerical Solutions of the Complete Navier-Stokes Equations

This report details the development of a new two-equation turbulence closure model based on the exact turbulent kinetic energy k and the variance of vorticity, zeta. The model, which is applicable to three dimensional flowfields, employs one set of model constants and does not use damping or wall functions, or geometric factors.

Robinson, David F.↗

Fatigue Analyses Under Constant- and Variable-Amplitude Loading Using Small-Crack Theory

Studies on the growth of small cracks have led to the observation that fatigue life of many engineering materials is primarily "crack growth" from micro-structural features, such as inclusion particles, voids, slip-bands or from manufacturing defects. This paper reviews the capabilities of a plasticity-induced crack-closure model to predict fatigue lives of metallic materials using "small-crack theory" under various loading conditions. Constraint factors, to account for three-dimensional effects, were selected to correlate large-crack growth rate data as a function of the effective stress-intensity factor range (delta-Keff) under constant-amplitude loading. Modifications to the delta-Keff-rate relations in the near-threshold regime were needed to fit measured small-crack growth rate behavior. The model was then used to calculate small-and large-crack growth rates, and to predict total fatigue lives, for notched and un-notched specimens under constant-amplitude and spectrum loading. Fatigue lives were predicted using crack-growth relations and micro-structural features like those that initiated cracks in the fatigue specimens for most of the materials analyzed. Results from the tests and analyses agreed well.

Newman, J. C., Jr.↗

A Specimen Size Effect on the Fatigue Crack Growth Rate Threshold of IN 718

Fatigue crack growth rate (FCGR) tests were conducted on IN 718 in the solution annealed and aged condition at room temperature in accordance with E647-87. As part of each test, the FCGR threshold was measured using the decreasing Delta K method. A new heat of material was being tested and some of this material was sent to a different laboratory which wanted to use a specimen with a 127 mm width. Threshold data previously had been established on specimens with a width of 50.8 mm. As a check of the laboratory, tests were conducted at room temperature and R equal to 0.1 for comparison with the earlier data. The results were a threshold significantly higher than previously observed. Interchanging of specimen sizes and laboratories showed that the results were not due to a heat-to-heat or lab-to-lab variation. The results to be presented here are those obtained at the original laboratory. Growth rates were measured using the electric potential drop technique at R values of 0.1, 0.7, and 0.9. Compact tension specimen sizes with planer dimensions of 25.4 mm, 50.8 mm, and 127 mm were used. Crack growth rates at threshold were generally below 2.5 X 10(exp -8) mm / cycle. Closure measurements were made on some of the specimens by a manual procedure using a clip gage. When the crack growth rate data for the specimens tested at R equal to 0.1 were plotted as a function of applied Delta K, the thresholds varied with specimen width. The larger the width, the higher the threshold. The thresholds varied from 6.5 MPa-m(exp 1/2) for the 25.4 mm specimen to 15.4 MPa-m(exp 1/2) for the 127 mm specimen. At R equal to 0.7, the 25.4 mm and 50.8 mm specimens had essentially the same threshold, about 2.9 MPa-m(exp 1/2)while the 127 mm specimen had a threshold of 4.5 MPa-m(exp 1/2). When plotted as a function of effective Delta K, the R equal to 0.1 data are essentially normalized. Various aspects of the test procedure will be discussed as well as the results of analysis of the data using some different closure models.

Garr, K. R.↗

Low Speed Analysis of Mission Adaptive Flaps on a High Speed Civil Transport Configuration

Thin-layer Navier-Stokes analyses were done on a high speed civil transport configuration with mission adaptive leading-edge flaps. The flow conditions simulated were Mach = 0.22 and Reynolds number of 4.27 million for angles-of-attack ranging from 0 to 18 degrees. Two turbulence closure models were used. Analyses were done exclusively with the Baldwin-Lomax turbulence model at low angle-of-attack conditions. At high angles-of-attack where considerable flow separation and vortices occurred the Spalart-Allmaras turbulence model was also considered. The effects of flow transition were studied. Predicted aerodynamic forces, moment, and pressure are compared to experimental data obtained in the 14- by 22-Foot Subsonic Tunnel at NASA Langley. The forces and moments correlated well with experimental data in terms of trends. Drag and pitching moment were consistently underpredicted. Predicted surface pressures compared well with experiment at low angles-of-attack. Above 10 angle-of-attack the pressure comparisons were not as favorable. The two turbulent models affected the pressures on the flap considerably and neither produced correct results at the high angles-of-attack.

Lessard, Victor R.↗

Numerical Experiments on Homogeneous Strained Turbulence Subjected to Coriolis Force

Homogeneous turbulent flows with various combinations of strain-rate, rotation rate and coriolis force capture some important aspects of more complex flows with streamline curvature and rotation. Presently, a situation is considered in which as a box of turbulence rotates, strain axes rotate with it. This is to be contrasted with the elliptic streamline flow in which the strain axes are fixed in an inertial frame. The elliptic flow is known to exhibit (inviscid) growth of turbulent energy and one might expect even more rapid growth with the strain-axes following the box. Instead, it is found that the sign of the Reynolds shear stress is reversed leading to a negative production term for turbulent energy. Partial understanding of the phenomenon is obtained from a consideration of the rotation of inertial waves relative to the strain axes as well as the "pressure-less" RDT argument put forward by Cambon etal. [J. Fluid Mech, 278, 175]. Some comparisons with the predictions of second-order closure models will be presented.

Shariff, K.↗

Elliptic Relaxation of a Tensor Representation of the Pressure-Strain and Dissipation Rate

A formulation to include the effects of wall-proximity in a second moment closure model is presented that utilizes a tensor representation for the redistribution term in the Reynolds stress equations. The wall-proximity effects are modeled through an elliptic relaxation process of the tensor expansion coefficients that properly accounts for both correlation length and time scales as the wall is approached. DNS data and Reynolds stress solutions using a full differential approach at channel Reynolds number of 590 are compared to the new model.

Carlson, John R.↗

Advances in Structural Integrity Analysis Methods for Aging Metallic Airframe Structures with Local Damage

Analysis methodologies for predicting fatigue-crack growth from rivet holes in panels subjected to cyclic loads and for predicting the residual strength of aluminum fuselage structures with cracks and subjected to combined internal pressure and mechanical loads are described. The fatigue-crack growth analysis methodology is based on small-crack theory and a plasticity induced crack-closure model, and the effect of a corrosive environment on crack-growth rate is included. The residual strength analysis methodology is based on the critical crack-tip-opening-angle fracture criterion that characterizes the fracture behavior of a material of interest, and a geometric and material nonlinear finite element shell analysis code that performs the structural analysis of the fuselage structure of interest. The methodologies have been verified experimentally for structures ranging from laboratory coupons to full-scale structural components. Analytical and experimental results based on these methodologies are described and compared for laboratory coupons and flat panels, small-scale pressurized shells, and full-scale curved stiffened panels. The residual strength analysis methodology is sufficiently general to include the effects of multiple-site damage on structural behavior.

fatigue-crack growth↗

Analyses of Fatigue and Fatigue-Crack Growth under Constant- and Variable-Amplitude Loading

Studies on the growth of small cracks have led to the observation that fatigue life of many engineering materials is primarily crack growth from micro-structural features, such as inclusion particles, voids, slip-bands or from manufacturing defects. This paper reviews the capabilities of a plasticity-induced crack-closure model to predict fatigue lives of metallic materials using small-crack theory under various loading conditions. Constraint factors, to account for three-dimensional effects, were selected to correlate large-crack growth rate data as a function of the effective stress-intensity factor range (delta K(sub eff)) under constant-amplitude loading. Modifications to the delta K(sub eff)-rate relations in the near-threshold regime were needed to fit measured small-crack growth rate behavior. The model was then used to calculate small- and large-crack growth rates, and to predict total fatigue lives, for notched and un-notched specimens under constant-amplitude and spectrum loading. Fatigue lives were predicted using crack-growth relations and micro-structural features like those that initiated cracks in the fatigue specimens for most of the materials analyzed. Results from the tests and analyses agreed well.

Newman, J. C., Jr.↗

Large Eddy Simulation of Transonic Flow Field in NASA Rotor 37

The current paper reports on numerical investigations on the flow characteristics in a transonic axial compressor, NASA Rotor 37. The flow field was used previously as a CFD blind test case conducted by American Society of Mechanical Engineers in 1994. Since the CFD blind-test exercise, many numerical studies on the flow field in the NASA Rotor 37 have been reported. Although steady improvements have been reported in both numerical procedure and turbulence closure, it is believed that all the important aspects of the flow field have not been fully explained with numerical studies based on the Reynolds Averaged Navier-Stokes (RANS) solution. Experimental data show large dip in total pressure distribution near the hub at downstream of the rotor at 100% rotor speed. Most original numerical solutions from the blind test exercise did not predict this total pressure deficit correctly. This total pressure deficit at the rotor exit was attributed to a hub corner flow separation by the author. Several subsequent numerical studies with different turbulence closure model also calculated this dip in total pressure rise. Also, several studies attributed this total pressure deficit to a small leakage flow coming from the hub in the test article. As the experimental study cannot be repeated, either explanation cannot be validated. The primary purpose of the current investigation is to investigate the transonic flow field with both RANS and a Large Eddy Simulation (LES). The RANS approach gives similar results presented at the original blind test exercise. Although the RANS calculates higher overall total pressure rise, the total pressure deficit near the hub is calculated correctly. The numerical solution shows that the total pressure deficit is due to a hub corner flow separation. The calculated pressure rise from the LES agrees better with the measured total pressure rise especially near the casing area where the passage shock interacts with the tip clearance vortex and flow becomes unsteady due to this interaction. The LES simulation also calculates the total pressure rise deficit near the hub and it agrees well with the measured data.

Hah, Chunill↗

Towards a Comprehensive Model of Jet Noise Using an Acoustic Analogy and Steady RANS Solutions

An acoustic analogy is developed to predict the noise from jet flows. It contains two source models that independently predict the noise from turbulence and shock wave shear layer interactions. The acoustic analogy is based on the Euler equations and separates the sources from propagation. Propagation effects are taken into account by calculating the vector Green's function of the linearized Euler equations. The sources are modeled following the work of Tam and Auriault, Morris and Boluriaan, and Morris and Miller. A statistical model of the two-point cross-correlation of the velocity fluctuations is used to describe the turbulence. The acoustic analogy attempts to take into account the correct scaling of the sources for a wide range of nozzle pressure and temperature ratios. It does not make assumptions regarding fine- or large-scale turbulent noise sources, self- or shear-noise, or convective amplification. The acoustic analogy is partially informed by three-dimensional steady Reynolds-Averaged Navier-Stokes solutions that include the nozzle geometry. The predictions are compared with experiments of jets operating subsonically through supersonically and at unheated and heated temperatures. Predictions generally capture the scaling of both mixing noise and BBSAN for the conditions examined, but some discrepancies remain that are due to the accuracy of the steady RANS turbulence model closure, the equivalent sources, and the use of a simplified vector Green's function solver of the linearized Euler equations.

Miller, Steven A. E.↗

Assessment of Higher-Order RANS Closures in a Decelerated Planar Wall-Bounded Turbulent Flow

A reference DNS database is presented, which includes third- and fourth-order moment budgets for unstrained and strained planar channel flow. Existing RANS closure models for third- and fourth-order terms are surveyed, and new model ideas are introduced. The various models are then compared with the DNS data term by term using a priori testing of the higher-order budgets of turbulence transport, velocity-pressure-gradient, and dissipation for both the unstrained and strained databases. Generally, the models for the velocity-pressure-gradient terms are most in need of improvement.

Jeyapaul, Elbert↗

Detached Eddy Simulation for the F-16XL Aircraft Configuration

Numerical simulations for the flow around the F-16XL configuration as a contribution to the Cranked Arrow Wing Aerodynamic Project International 2 (CAWAPI-2) have been performed. The NASA Langley Tetrahedral Unstructured Software System (TetrUSS) with its USM3D solver was used to perform the unsteady flow field simulations for the subsonic high angle-of-attack case corresponding to flight condition (FC) 25. Two approaches were utilized to capture the unsteady vortex flow over the wing of the F-16XL. The first approach was to use Unsteady Reynolds-Averaged Navier-Stokes (URANS) coupled with standard turbulence closure models. The second approach was to use Detached Eddy Simulation (DES), which creates a hybrid model that attempts to combine the most favorable elements of URANS models and Large Eddy Simulation (LES). Computed surface static pressure profiles are presented and compared with flight data. Time-averaged and instantaneous results obtained on coarse, medium and fine grids are compared with the flight data. The intent of this study is to demonstrate that the DES module within the USM3D solver can be used to provide valuable data in predicting vortex-flow physics on a complex configuration.

Elmiligui, Alaa↗

Hybrid Transportation System Integrated Trajectory Design and Optimization for Mars Landing Site Accessibility

NASA’s Mars Study Capability Team continues the agency’s efforts to study and refine the nation’s plan to field a sustainable human Mars campaign. One of the primary open issues left unanswered during previous analysis cycles was the uncertainty of the impact to the vehicle performance requirement to deliver crew and cargo to the same landing site across multiple mission opportunities. The Mars Study Capability Team has recently developed an integrated trajectory optimization and system closure model to solve the complex interplanetary trajectory optimization using both low-thrust and high-thrust maneuvers. This paper demonstrates the capability of this new integrated trajectory design and optimization method as it applies to the landing site accessibility problem for the Hybrid transportation architecture. The results showed that the current vehicle is capable of reaching up to +20 degree latitude and down to -20 degree latitude in every mission opportunity from 2033 to 2052. However, reaching latitudes beyond +/- 20 required more propellant than the spacecraft is currently designed to carry for many of the mission opportunities.

Chai, Patrick R.↗

Pretest Simulations of a Supersonic Mixing and Combustion Validation Experiment to Assess Sensitivities

The reliance on CFD simulations to develop, design, and optimize scramjet systems (or components) has become commonplace. This reliance inevitably hinges on the ability of the computational analyst to quantify the level of confidence in their computational results. Unfortunately, nearly all the measured data available for this assessment comes from antiquated experimental datasets, or from tests that focused on the extraction of scramjet system (or component) performance. The objective of a CFD validation experiment is to quantify the predictive accuracy of one or more of the CFD physics submodels, implying that other uncertainties related to replicating the facility flow environment (e.g., knowledge of boundary conditions) must be minimized to the extent possible. This inevitably places stringent requirements on the quality and quantity of measurements taken to accurately specify inflow, outflow, and surface conditions for the CFD simulations; in addition to the measurements taken for the validation of physics submodels. This places additional demands on the experimental process above and beyond those for test article performance assessment. A recent high speed code credibility workshop series sponsored by AFRL identified a gap in existing validation data for fundamental assessments of turbulent mixing and combustion CFD closure models at scramjet engine relevant conditions. To address this gap, engineers at AFRL have designed a coaxial jet flame configuration that will be tested at two facilities (Research Cell 19 at the Air Force Research Lab, and at Purdue University). The effort described here documents pretest simulations of this validation experiment with the goal of fleshing out the extent of the facility flowpath that must be included to adequately reproduce the facility test section flow environment. The findings indicate that the flow around the support structure for the fuel injection centerbody upstream of the facility nozzle generates disturbances that persist throughout the nozzle expansion process; corrupting the azimuthal symmetry that was desired in the fuel/air mixing region of the test section. Simulations without this support structure maintained a high degree of azimuthal symmetry up until the fuel injection plane. However, even in this scenario the azimuthal symmetry was not maintained once the centerbody boundary layer transitioned to a wake flow downstream of the fuel injection plane.

CFD↗

Pretest Simulations of a Supersonic Mixing and Combustion Validation Experiment to Assess Sensitivities

The reliance on CFD simulations to develop, design, and optimize scramjet systems (or components) has become commonplace. This reliance inevitably hinges on the ability of the computational analyst to quantify the level of confidence in their computational results. Unfortunately, nearly all the measured data available for this assessment comes from antiquated experimental datasets, or from tests that focused on the extraction of scramjet system (or component) performance. The objective of a CFD validation experiment is to quantify the predictive accuracy of one or more of the CFD physics submodels, implying that other uncertainties related to replicating the facility flow environment (e.g., knowledge of boundary conditions) must be minimized to the extent possible. This inevitably places stringent requirements on the quality and quantity of measurements taken to accurately specify inflow, outflow, and surface conditions for the CFD simulations; in addition to the measurements taken for the validation of physics submodels. This places additional demands on the experimental process above and beyond those for test article performance assessment. A recent high speed code credibility workshop series sponsored by AFRL identified a gap in existing validation data for fundamental assessments of turbulent mixing and combustion CFD closure models at scramjet engine relevant conditions. To address this gap, engineers at AFRL have designed a coaxial jet flame configuration that will be tested at two facilities (Research Cell 19 at the Air Force Research Lab, and at Purdue University). The effort described here documents pretest simulations of this validation experiment with the goal of fleshing out the extent of the facility flowpath that must be included to adequately reproduce the facility test section flow environment. The findings indicate that the flow around the support structure for the fuel injection centerbody upstream of the facility nozzle generates disturbances that persist throughout the nozzle expansion process; corrupting the azimuthal symmetry that was desired in the fuel/air mixing region of the test section. Simulations without this support structure maintained a high degree of azimuthal symmetry up until the fuel injection plane. However, even in this scenario the azimuthal symmetry was not maintained once the centerbody boundary layer transitioned to a wake flow downstream of the fuel injection plane.

CFD↗

Fan Noise Predictions of the NASA Source Diagnostic Test using Unsteady Simulations with LAVA Part II - Tonal and Broadband Noise Assessment

The NASA Source Diagnostic Test (SDT) campaign experimental data is used for validation of a sliding mesh technique recently implemented within the Launch, Ascent, and Vehicle Aerodynamics computational framework for time-accurate simulation of rotating fans. The far-field acoustics are analyzed in this work, building upon the aerodynamic validation studies previously published in Part I. Two modeling approaches are explored: the unsteady Reynolds-averaged Navier Stokes (uRANS) with the negative Spalart-Allmaras (SA-neg) turbulence model closure, and a hybrid Reynolds-averaged Navier Stokes/large-eddy simulation (RANS/LES) paradigm employing a zonal detached-eddy simulation (ZDES) closure with enhanced shielding protection. Two convective flux scheme approaches with different dissipation properties are also explored with ZDES. The Ffowcs-Williams and Hawkings (FW-H) permeable surface approach is used for propagation of the near-field acoustics to the far-field microphone locations. This work analyzes the low-speed (approach) condition, characterized by a fan rotation speed of 7808 rotations-per-minute (RPM). Three different grid levels ranging between 200 million and 1.1 billion grid points are considered. Results show good prediction of the broadband noise levels at sideline angles ranging between 70° and 110°. The forward arc observers show an under-prediction of the overall sound pressure levels (OASPL) even at the fine grid level. A breakdown of the inlet and exhaust contributions reveals a steep drop-off in the broadband noise levels past a blade-passing frequency (BPF) of 1.5, potentially caused by a lack of resolved small-scale turbulent fluctuations in the interstage region. Past 110° the OASPL are over-predicted by up to 10 dB due to an over-prediction of the low-frequency broadband noise levels in the aft arc. The detuning of BPF2 caused by small deviations in the blade stagger angle around the wheel is captured, and a corresponding decrease in the sound power level for this cut-on tone is observed.

AATT↗

Combining Progressive Nodal Release with the Virtual Crack Closure Technique to Model Fatigue Delamination Growth Without Re-Meshing

The present work summarizes an approach to model mixed-mode 3D fatigue crack growth using the Virtual Crack Closure Technique (VCCT) without requiring re-meshing. It is demonstrated that the proposed approach can be used to simulate crack shapes that do not conform to the underlying mesh. The proposed approach relies solely on Paris Law characterization data to model delamination growth. Growth is determined as a post-processing step at the end of each increment, and hence no convergence issues associated with the progressive nodal release are encountered. This approach can be readily applied using standard solid element formulations and is implemented via an interface user element in Abaqus/Standard.

Carvalho, N. V. De↗