Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “shear correction”

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 163 records · Page 9

A computational model for the prediction of jet entrainment in the vicinity of nozzle boattails (the BOAT code)

The development of a computational model (BOAT) for calculating nearfield jet entrainment, and its incorporation in an existing methodology for the prediction of nozzle boattail pressures, is discussed. The model accounts for the detailed turbulence and thermochemical processes occurring in the mixing layer formed between a jet exhaust and surrounding external stream while interfacing with the inviscid exhaust and external flowfield regions in an overlaid, interactive manner. The ability of the BOAT model to analyze simple free shear flows is assessed by comparisons with fundamental laboratory data. The overlaid procedure for incorporating variable pressures into BOAT and the entrainment correction employed to yield an effective plume boundary for the inviscid external flow are demonstrated. This is accomplished via application of BOAT in conjunction with the codes comprising the NASA/LRC patched viscous/inviscid methodology for determining nozzle boattail drag for subsonic/transonic external flows.

Dash, S. M.↗

An alternating-direction-implicit algorithm for the unsteady potential equation in conservation form

An implicit finite difference scheme for an efficient computation of unsteady potential flow about airfoils is presented. The formulation uses density and velocity potential as dependent variables, and is cast in conservation form to assure the theoretically correct determination of shockwave location and speed. To enable boundary conditions to be imposed directly on the airfoil surface, a time varying sheared rectilinear coordinate transformation is employed. Calculated time history solutions on a pulsating airfoil are compared with the results of another unsteady transonic code. It is concluded that the method has excellent numerical stability and gives accurate solutions with sharply resolved shocks.

Chipman, R. R.↗

Second order modeling of boundary-free turbulent shear flows

A set of realizable second order models for boundary-free turbulent flows is presented. The constraints on second order models based on the realizability principle are re-examined. The rapid terms in the pressure correlations for both the Reynolds stress and the passive scalar flux equations are constructed to exactly satisfy the joint realizability. All other model terms (return-to-isotropy, third moments, and terms in the dissipation equations) already satisfy realizability. To correct the spreading rate of the axisymmetric jet, an extra term is added to the dissipation equation which accounts for the effect of mean vortex stretching on dissipation. The test flows used in this study are the mixing shear layer, plane jet, axisymmetric jet, and plane wake. The numerical solutions show that the unified model equations predict all these flows reasonably. It is expected that these models would be suitable for more complex and critical flows.

Shih, T.-H.↗

Second order modeling of boundary-free turbulent shear flows

A set of realizable second order models for boundary-free turbulent flows is presented. The constraints on second order models based on the realizability principle are re-examined. The rapid terms in the pressure correlations for both the Reynolds stress and the passive scalar flux equations are constructed to exactly satisfy the joint realizability. All other model terms (return-to-isotropy, third moments, and terms in the dissipation equations) already satisfy realizability. To correct the spreading rate of the axisymmetric jet, an extra term is added to the dissipation equation which accounts for the effect of mean vortex stretching on dissipation. The test flows used in this study are the mixing shear layer, plane jet, axisymmetric jet, and plane wake. The numerical solutions show that the unified model equations predict all these flows reasonably. It is expected that these models would be suitable for more complex and critical flows.

Shih, T.-H.↗

Hydrostatic Stress Effects Incorporated Into the Analysis of the High-Strain-Rate Deformation of Polymer Matrix Composites

Procedures for modeling the effect of high strain rate on composite materials are needed for designing reliable composite engine cases that are lighter than the metal cases in current use. The types of polymer matrix composites that are likely to be used in such an application have a deformation response that is nonlinear and that varies with strain rate. The nonlinearity and strain rate dependence of the composite response is primarily due to the matrix constituent. Therefore, in developing material models to be used in the design of impact-resistant composite engine cases, the deformation of the polymer matrix must be correctly analyzed. However, unlike in metals, the nonlinear response of polymers depends on the hydrostatic stresses, which must be accounted for within an analytical model. An experimental program has been carried out through a university grant with the Ohio State University to obtain tensile and shear deformation data for a representative polymer for strain rates ranging from quasi-static to high rates of several hundred per second. This information has been used at the NASA Glenn Research Center to develop, characterize, and correlate a material model in which the strain rate dependence and nonlinearity (including hydrostatic stress effects) of the polymer are correctly analyzed. To obtain the material data, Glenn s researchers designed and fabricated test specimens of a representative toughened epoxy resin. Quasi-static tests at low strain rates and split Hopkinson bar tests at high strain rates were then conducted at the Ohio State University. The experimental data confirmed the strong effects of strain rate on both the tensile and shear deformation of the polymer. For the analytical model, Glenn researchers modified state variable constitutive equations previously used for the viscoplastic analysis of metals to allow for the analysis of the nonlinear, strain-rate-dependent polymer deformation. Specifically, we accounted for the effects of hydrostatic stresses. An important discovery in the course of this work was that the hydrostatic stress effects varied during the loading process, which needed to be accounted for within the constitutive equations. The model is characterized primarily by shear data, with tensile data used to characterize the hydrostatic stress effects.

Goldberg, Robert K.↗

Density fluctuations and radiated noise for a high-temperature supersonic jet

Experimental data on density fluctuations were obtained by the laser Schlieren method in a supersonic jet which at the nozzle exit had a Mach number of 1.43 and a stagnation temperature of about 1090 K. The jet emerged into the ambient atmosphere in an anechoic chamber, correctly expanded from a nozzle which had an exit diameter of 10.8 cm. Using the information on the density fluctuations and the mean shear obtained by probes, the autocorrelation of the radiated noise was calculated by a theory that is suitable for Mach wave emission. This theory is a modification of that developed by Ffowcs Williams and Maidanik (1965). The calculated noise field agrees well with that obtained by using microphones outside the jet.

Parthasarathy, S. P.↗

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.↗

An alternating-direction implicit algorithm for unsteady potential flow

An implicit finite-difference scheme is presented for the efficient computation of unsteady potential flow about airfoils. The formulation uses density and the velocity potential as dependent variables and is cast in conservation form to assure the theoretically correct determination of shockwave location and speed. To enable boundary conditions to be imposed directly on the airfoil surface, a time-varying sheared-rectilinear coordinate transformation is employed. Calculated time-history solutions on a pulsating airfoil are compared with the results of other unsteady transonic codes, including a previous method of the authors. The present method is demonstrated to be unconditionally stable and to give accurate solutions with sharply resolved shocks.

Chipman, R.↗

Toward Relating Open- and Closed-Test Section Microphone Phased Array Aeroacoustic Measurements

Aeroacoustic measurements utilizing microphone phased arrays are becoming more commonplace in large-scale wind tunnels. However, for such costly tests it is difficult to assess the influence of a given facility's configuration on the acquired data. One such opportunity has recently arisen in the NASA Langley Research Center 14- by 22-Foot Subsonic Tunnel, where the High-Lift configuration of the Common Research Model was tested in both open-jet and closed wall configurations. Acoustic data were acquired on the same model with the same specific hardware at the same Mach number; and for similar model angles of attack, array apertures, and array angles relative to the model location. Initial deconvolution comparisons using DAMAS appear poor. The closed test section results show contamination likely due to background noise and reflections. The open-jet results show severe blurring due to the decorrelation of acoustic waves passing through a turbulent free shear layer. However, accounting for image sources in the DAMAS grid reduces visual artifacts in the closed test section results, and a preliminary coherence correction technique mitigates the blurring in the open-jet results. Qualitatively, both tests identify similar source behavior. Quantitatively, summed spectra of the slat noise agree in some frequency bands but differ in others. Proposals that may mitigate the difference are given for each configuration, such as correcting for the array protective cover and applying background subtraction in the closed wall configuration, and further developing the coherence loss correction for deconvolution methods in the open-jet setup.

microphone array↗

Airborne in situ computation of the wind shear hazard index

An algorithm for airborne in situ computation of the wind shear hazard index (F-factor) was developed and evaluated in simulation and verified in flight. The algorithm was implemented on NASA's B-737-100 airplane, and tested under severe maneuvering, nonhazardous wind conditions, and normal takeoffs and landings. The airplane was flown through actual microburst conditions in Orlando, FL, where the algorithm produced wind shear measurements which were confirmed by an independent, ground-based radar measurement. Flight test results indicated that the in situ F-factor algorithm correctly measured the effect of the wind environment on the airplane's performance, and produced no nuisance alerts.

Oseguera, Rosa M.↗

Recent Advancements in the PATO Material Response Code

Introduction: Predicting the complicated multiphysics phenomena during atmospheric entry requires high-fidelity modeling tools to refine estimates of mission risks during entry. To this end, new capabilities are being added to the Porous-material Analysis Toolbox based on OpenFOAM (PATO). PATO is an open-source software for Computational Material Response (CMR) of reactive porous materials submitted to high-temperature environments. The objective of this work is to highlight current efforts to add to and improve upon the modeling capabilities of PATO. These include efforts to loosely couple PATO with other discipline specialized codes including hypersonic Computational Fluid Dynamics (CFD), to assess the interaction effects between pyrolysis gas blowing and the boundary layer, and Computational Solid Mechanics (CSM), to address modeling of mechanical erosion. Other refinements include surface phenomena modeling capabilities to address the effects of silicone-based coatings applied to the TPS during flight preparation, and a unified multiphase solver for a mixed porous-material and plain-fluid domain. Coupling CMR with CFD (CMR/CFD): A loose coupling between PATO and the Data Parallel Line Relaxation (DPLR) CFD code has been achieved by making use of a blowing boundary condition at the heatshield surface available in DPLR. Starting with heat flux estimates with no pyrolysis gas blowing at the surface, blowing gases are computed by the CMR and passed to the CFD such that aerothermal properties of the environment can be recomputed for a new CMR computation. This leads to an iterative process which is supplemented with an estimate of the radiative heat flux using the Nonequilibrium air radiation (NEQAIR) program. The entire iterative process is illustrated in Figure 1. This coupling strategy has been utilized in computing the MSL material response. The goal is to compare the coupled CMR/CFD results with material response results obtained using traditional blowing corrections. Coupling CMS with CMR: A mechanical erosion model is currently being implemented in PATO to account for the additional mass removal induced by high shear conditions. The modeling process at each timestep consists of updating the mechanical properties as a function of temperature and computing the stress tensor and displacement fields of the material. Then, a failure criteria model determines the regions in which the stress exceeds the ultimate strength values resulting in mesh movement to account for mass removal. This model allows the material response simulation to compute the recession due to both oxidation and shear-induced erosion. The model is demonstrated by computing material response of sphere-cone arc jet samples. Surface Modeling Capabilities: NuSil, a silicone-based coating, was sprayed onto the MSL and Mars 2020 heatshields to mitigate shedding of phenolic dust. To better understand the effects of the NuSil coating on the material response, a novel model has been implemented in PATO. In this model, the equilibrium of the charred NuSil surface is modeled as pure silica, and a constant offset, inspired by the classical spallation model, is added to the the char blowing rate and wall enthalpy to reproduce HyMETS experimental results. The model has also been used to estimate the 3D material response of the MSL heatshield. Unified Solver: In addition to the iterative loose coupling approach mentioned above, a multiphase unified solver is being developed to couple the environment (plain-fluid phase) and the porous-material phase. The solver is based on the volume averaged conservation of mass, momentum, and energy for the macroscale with closure models which include microscale effects through effective physicochemical properties. The unified solver has been used to compute flow through a porous plug and solve the Beavers and Joseph problem. Since the strong coupling between phases is inherent to this solver, modeling assumptions present in other coupling methods of material response are mitigated. This strategy also makes it feasible to capture the competition between surface and volume ablation in the same computational domain, which is usually not possible with other coupling approaches.

Material Response↗

Recent Advancements in the PATO Material Response Code

Introduction: Predicting the complicated multiphysics phenomena during atmospheric entry requires high-fidelity modeling tools to refine estimates of mission risks during entry. To this end, new capabilities are being added to the Porous-material Analysis Toolbox based on OpenFOAM (PATO) [1,2,3]. PATO is an open-source software for Computational Material Response (CMR) of reactive porous materials submitted to high-temperature environments. The objective of this work is to highlight current efforts to add to and improve upon the modeling capabilities of PATO. These include efforts to loosely couple PATO with other discipline specialized codes including hypersonic Computational Fluid Dynamics (CFD), to assess the interaction effects between pyrolysis gas blowing and the boundary layer, and Computational Solid Mechanics (CSM), to address modeling of mechanical erosion. Other refinements include surface phenomena modeling capabilities to address the effects of silicone-based coatings applied to the TPS during flight preparation, and a unified multiphase solver for a mixed porous-material and plain-fluid domain. Coupling CMR with CFD (CMR/CFD): A loose coupling between PATO and the Data Parallel Line Relaxation (DPLR) [4] CFD code has been achieved by making use of a blowing boundary condition at the heatshield surface available in DPLR. Starting with heat flux estimates with no pyrolysis gas blowing at the surface, blowing gases are computed by the CMR and passed to the CFD such that aerothermal properties of the environment can be recomputed for a new CMR computation. This leads to an iterative process which is supplemented with an estimate of the radiative heat flux using the Nonequilibrium air radiation (NEQAIR) [5] program. The entire iterative process is illustrated in Figure 1. This coupling strategy has been utilized in computing the MSL material response. The goal is to compare the coupled CMR/CFD results with material response results obtained using traditional blowing corrections. Coupling CMS with CMR: A mechanical erosion model is currently being implemented in PATO to account for the additional mass removal induced by high shear conditions. The modeling process at each timestep consists of updating the mechanical properties as a function of temperature and computing the stress tensor and displacement fields of the material. Then, a failure criteria model determines the regions in which the stress exceeds the ultimate strength values resulting in mesh movement to account for mass removal. This model allows the material response simulation to compute the recession due to both oxidation and shear-induced erosion. The model is demonstrated by computing material response of sphere-cone arc jet samples. Surface Modeling Capabilities: NuSil, a silicone-based coating, was sprayed onto the MSL and Mars 2020 heatshields to mitigate shedding of phenolic dust. To better understand the effects of the NuSil coating on the material response, a novel model has been implemented in PATO. In this model, the equilibrium of the charred NuSil surface is modeled as pure silica, and a constant offset, inspired by the classical spallation model, is added to the the char blowing rate and wall enthalpy to reproduce HyMETS experimental results. The model has also been used to estimate the 3D material response of the MSL heatshield [6]. Unified Solver: In addition to the iterative loose coupling approach mentioned above, a multiphase unified solver is being developed to couple the environment (plain-fluid phase) and the porous-material phase. The solver is based on the volume averaged conservation of mass, momentum, and energy for the macroscale with closure models which include microscale effects through effective physicochemical properties. The unified solver has been used to compute flow through a porous plug and solve the Beavers and Joseph problem [7]. Since the strong coupling between phases is inherent to this solver, modeling assumptions present in other coupling methods of material response are mitigated. This strategy also makes it feasible to capture the competition between surface and volume ablation in the same computational domain, which is usually not possible with other coupling approaches.

Thermal Protection Systems↗

A Rotary Flow Channel for Shear Stress Sensor Calibration

A proposed shear sensor calibrator consists of a rotating wheel with the sensor mounted tangential to the rim and positioned in close proximity to the rim. The shear stress generated by the flow at the sensor position is simply tau(sub omega) = (mu)r(omega)/h, where mu is the viscosity of the ambient gas, r the wheel radius, omega the angular velocity of the wheel, and h the width of the gap between the wheel rim and the sensor. With numerical values of mu = 31 (mu)Pa s (neon at room temperature), r = 0.5 m, omega = 754 /s (7200 rpm), and h = 50.8 m, a shear stress of tau(sub omega) = 231 Pa can be generated. An analysis based on one-dimensional flow, with the flow velocity having only an angular component as a function of the axial and radial coordinates, yields corrections to the above simple formula for the curvature of the wheel, flatness of the sensor, and finite width of the wheel. It is assumed that the sensor mount contains a trough (sidewalls) to render a velocity release boundary condition at the edges of the rim. The Taylor number under maximum flow conditions is found to be 62.3, sufficiently low to obviate flow instability. The fact that the parameters entering into the evaluation of the shear stress can be measured to high accuracy with well-defined uncertainties makes the proposed calibrator suitable for a physical standard for shear stress calibration.

Zuckerwar, Allan J.↗

Correction of lunar seismograms for instrumental and near-surface effects and constraints on the velocity structure of the lunar interior

Long-period lunar seismograms were studied with the aim of identifying consistent sets of direct shear and secondary wave arrivals, thus constraining the velocities in the lunar mantle and the depths of the velocity discontinuities. Two velocity models were used to locate the natural impacts and the shallow moonquakes and to obtain the travel time residuals. Seismic sections were made of the radial, transverse, and vertical components of ground motion for impacts, shallow, and deep moonquakes in order to search for consistent sets of secondary wave arrivals. No conclusive set of secondary arrivals could be recognized on the seismic sections and thus the velocities and depths of the velocity discontinuities cannot be severely constrained by secondary arrivals. It is likely that the crust is thinner than 50 km and that a first-order discontinuity separates the upper and lower crust at a depth of between 20 and 30 km.

Horvath, P.↗

The aerodynamics of circulation control

Two dimensional subsonic wind tunnel tests were conducted on a 20% thickness: chord ratio circulation controlled elliptic aerofoil section equipped with forward and reverse blowing slots. Overall performance measurements were made over a range of trailing edge blowing momentum coefficients from 0 to 0.04; some included the effect of leading edge blowing. A detailed investigation of the trailing edge wall jet, using split film probes, hot wire probes and total head tubes, provided measurements of mean velocity components, Reynolds normal and shear stresses, and radial static pressure. The closure of the two dimensional angular momentum and continuity equations was examined using the measured data, with and without correction, and the difficulty of obtaining a satisfactory solution illustrated. Suggestions regarding the nature of the flow field which should aid the understanding of Coanda effect and the theoretical solution of highly curved wall jet flows are presented.

Wood, N. J.↗

Effects of space environment on space-based radar phased-array antenna: Status and preliminary observations

The overall objective is to evaluate the effect of the space environment on Kapton films considered for the Grumman space based radar (SBR) phased-array antenna. The most striking result is the overall good condition of the Kapton antenna planes and Kapton tensile specimens. This is largely attributable to the orientation of the Kapton (parallel and flush on the space end) and the stability of the LDEF in orbit. Results on elongation and mechanical properties of the plain and fiberglass reinforced Kapton will be described. Stress-dependent permanent deformation and some reductions in strain to failure were observed. Physical property testing of the materials is in progress. Electronic data acquisition and memory systems appeared to operate correctly, but functional tests were not yet performed. An evaluation of the high voltage-plasma interaction data is underway. Some minor systems anomalies (e.g., fastener sheared during removal, strong odor inside electronics container) were noted. Other observations such as radiation, contamination, impacts, and orientation features of atomic oxygen erosion are reported.

Whiteside, J. B.↗

A Note on the Wave Action Density of a Viscous Instability Mode on a Laminar Free-shear Flow

Using the assumptions of an incompressible and viscous flow at large Reynolds number, we derive the evolution equation for the wave action density of an instability wave traveling on top of a laminar free-shear flow. The instability is considered to be viscous; the purpose of the present work is to include the cumulative effect of the (locally) small viscous correction to the wave, over length and time scales on which the underlying base flow appears inhomogeneous owing to its viscous diffusion. As such, we generalize our previous work for inviscid waves. This generalization appears as an additional (but usually non-negligible) term in the equation for the wave action. The basic structure of the equation remains unaltered.

Balsa, Thomas F.↗

Crack tip deformation and fatigue crack growth

Recent research on fatigue crack growth is summarized. Topics discussed include the use of the differential stress intensity factor to characterize crack tip deformation, the use of the unzipping model to study the growth of microcracks and the fatigue crack growth in a ferritic-martensitic steel, and the development of a model of fatige crack growth threshold. It is shown that in the case of small yielding, the differential stress intensity factor provides an adequate description of cyclic plastic deformation at the crack tip and correlates well with the crack growth rate. The unzipping model based on crack tip shear decohesion process is found to be in good agreement with the measured crack growth and striation spacing measurements. The proposed model of crack growth threshold gives correct predictions of the crack growth behavior in the near-threshold region.

Liu, H.-W.↗