SEARCH · Engineering Papers
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.
Hypersonic laminar viscous interaction effects on aerodynamics of two-dimensional wedge and triangular planform wings
Fluid-structure interaction simulations have the potential to accelerate the certification process of parachute systems for interplanetary robotic exploration missions while helping to mitigate risks by covering a wider range of scenarios than can be flight tested. The Launch, Ascent, and Vehicle Aerodynamics team is working toward actualizing this potential by further developing the capability to perform such fluid-structure interaction simulations and validating it with the best supersonic parachute flight test data available. We perform several simulations of the Advanced Supersonic Parachute Inflation Research Experiments’ third flight test (SR03) where we independently vary several parameters to determine their sensitivity on the predicted drag. We investigate sensitivity of the drag to the added mass of the confluence fitting, to the effective porosity, to how the gore seams are modeled, to viscous effects, and to deceleration. Modeling the added mass of the confluence fitting has virtually no impact on the drag, but is worthwhile to include to better match the overall system’s mass, which is important for the simulations with deceleration. Parachute drag is only marginally sensitive to porosity and viscous effects. Maximum parachute drag is most sensitive to how the gore seams are modeled, causing a nearly 20% drop when they are modeled as stiffer radials, and to deceleration, causing an additional 10% reduction. Results show improved agreement with flight test measurements when the FSI model is closest to the as-built configuration and physics experienced in flight: with the confluence fitting, viscous effects, the gore seams modeled as stiff radials, and including deceleration.
Steady and unsteady measured pressures for a Rectangular Supercritical Wing (RSW) undergoing pitching oscillations have been presented. From the several hundred compiled data points, 27 static and 36 pitching oscillation cases have been proposed for computational Test Cases to illustrate the trends with Mach number, reduced frequency, and angle of attack. The wing was designed to be a simple configuration for Computational Fluid Dynamics (CFD) comparisons. The wing had an unswept rectangular planform plus a tip of revolution, a panel aspect ratio of 2.0, a twelve per cent thick supercritical airfoil section, and no twist. The model was tested over a wide range of Mach numbers, from 0.27 to 0.90, corresponding to low subsonic flows up to strong transonic flows. The higher Mach numbers are well beyond the design Mach number such as might be required for flutter verification beyond cruise conditions. The pitching oscillations covered a broad range of reduced frequencies. Some early calculations for this wing are given for lifting pressure as calculated from a linear lifting surface program and from a transonic small perturbation program. The unsteady results were given primarily for a mild transonic condition at M = 0.70. For these cases the agreement with the data was only fair, possibly resulting from the omission of viscous effects. Supercritical airfoil sections are known to be sensitive to viscous effects (for example, one case cited). Calculations using a higher level code with the full potential equations have been presented for one of the same cases, and with the Euler equations. The agreement around the leading edge was improved, but overall the agreement was not completely satisfactory. Typically for low-aspect-ratio rectangular wings, transonic shock waves on the wing tend to sweep forward from root to tip such that there are strong three-dimensional effects. It might also be noted that for most of the test, the model was tested with free transition, but a few points were taken with an added transition strip for comparison. Some unpublished results of a rigid wing of the same airfoil and planform that was tested on the pitch and plunge apparatus mount system (PAPA) showed effects of the lower surface transition Strip on flutter at the lower subsonic Mach numbers. Significant effects of a transition strip were also obtained on a wing with a thicker supercritical section on the PAPA mount system. Both of these flutter tests on the PAPA resulted in very low reduced frequencies that may be a factor in this influence of the transition strip. However, these results indicate that correlation studies for RSW may require some attention to the estimation of transition location to accurately treat viscous effects. In this report several Test Cases are selected to illustrate trends for a variety of different conditions with emphasis on transonic flow effects. An overview of the model and tests is given and the standard formulary for these data is listed. Sample data points are presented in both tabular and graphical form. A complete tabulation and plotting of all the Test Cases is given. Only the static pressures and the real and imaginary parts of the first harmonic of the unsteady pressures are available. All the data for the test are available in electronic file form. The Test Cases are also available as separate electronic files.
The thermal evaporation of a cold spherical cloud embedded in hot gas is analyzed, including the effects of viscous stresses as well as 'saturation' of the conductive heat flux, for the case of equal electron and ion temperatures. The relevant transport properties for an H-He plasma are discussed, including the modification of these transport coefficients when temperature or velocity gradients become large. The fluid-dynamical equations for steady spherical evaporation are given and cast into dimensinless form. Numerical integration of these equations and the results are discussed, emphasizing the novel features of the viscous solutions, astrophysical implications, and uncertainties in the physics.
Transonic flow solutions are obtained over a multielement airfoil (augmentor-wing) using the full-potential equation. Solutions obtained for a subcritical case and a strong shock case show good quantitative agreement with experiment in regions not dominated by viscous effects. In those regions where viscous effects are dominant, the results are still in good qualitative agreement. For the strong shock case, Mach number and angle-of-attack corrections were necessary to match experimental coefficient of lift. Typical results from the transonic augmentor-wing Potential Code on the Cray-1S computer require about 10 sec of CPU time for a three-order-of-magnitude drop in the maximum residual. The speed with which solutions can be generated, and the associated low cost, will make this code a practical tool for the design aerodynamicist.
Massive blowing effect on stagnation point radiative energy transfer of ablative heat shield during planetary entry at hypersonic velocities
An analytical model was developed to provide additional insight and understanding of the factors that influence the simulation and prediction of leak rates from small cracks in pressurized containers. Specifically, the analysis was aimed at developing an analytical model capable of predicting the hydrogen leak rates from a crack in the combustion chamber coolant discharge manifold on main engine 1 of the Orbiter Challenger that was discovered during flight readiness firings 1 and 2. This model was based on viscous pipe flow analyses and calibrated for the crack geometry by using helium leak-rate data obtained from both low- and high-pressure tests used to simulate the flight readiness firing test conditions. In addition, this model includes the effects of crack width changes caused by different working stresses associated with the different test conditions. Because of the combination of the small crack dimensions and the wide range of pressures used for the test conditions, either laminar or turbulent viscous effects dominated the flows at all test conditions. This model was used to illustrate the sensitivity of the predicted leak rates to considerations of test conditions, viscous flow effects, and geometric features of the crack. In addition, the model was certified by comparing the hydrogen leak-rate prediction for the flight readiness firing test condition to the actual measured leak rate. The prediction was within 9 percent of the measured value.
With the renewed interest in Cartesian gridding methodologies for the ease and speed of gridding complex geometries in addition to the simplicity of the control volumes used in the computations, it has become important to investigate ways of extending the existing Cartesian grid solver functionalities. This includes developing methods of modeling the viscous effects in order to utilize Cartesian grids solvers for accurate drag predictions and addressing the issues related to the distributed memory parallelization of Cartesian solvers. This research presents advances in two areas of interest in Cartesian grid solvers, viscous effects modeling and MPI parallelization. The development of viscous effects modeling using solely Cartesian grids has been hampered by the widely varying control volume sizes associated with the mesh refinement and the cut cells associated with the solid surface. This problem is being addressed by using physically based modeling techniques to update the state vectors of the cut cells and removing them from the finite volume integration scheme. This work is performed on a new Cartesian grid solver, NASCART-GT, with modifications to its cut cell functionality. The development of MPI parallelization addresses issues associated with utilizing Cartesian solvers on distributed memory parallel environments. This work is performed on an existing Cartesian grid solver, CART3D, with modifications to its parallelization methodology.
Viscous real gas and ideal gas supersonic flowfields over the forebody of a aeroassisted orbital transfer vehicle are determined using a unsteady factored implicit algorithm. Air in chemical equilibrium is considered and its local thermodynamic properties are computed by an equilibrium composition method. Numerical solutions are obtained for both real and ideal gases at a Mach number of 30 and at angles of attack up to 20 degrees. Shock stand-off distances and surface pressure distributions are presented for the gas models with and without viscous effects. For the freestream conditions considered, viscous effects dominate the flow.
Simplified theoretical approaches are shown, based on hypersonic similarity boundary-layer theory, which allow reasonably accurate estimates to be made of the surface pressures on plates on which viscous effects are important. The consideration of viscous effects includes the cases where curved surfaces, stream pressure gradients, and leadingedge bluntness are important factors.
An experimental aerodynamic investigation was conducted in the AEDC-VKF Hypervelocity Wind Tunnel (Tunnel F) at a nomial Mach number of 19 to determine hypersonic viscous interaction effects on the space shuttle orbiter. The tests were conducted at an angle of attack of 30 degrees over a free-stream Reynolds number (based on fuselage length) variation from 0.1 to 0.4 million. Viscous interaction parameter was varied from 0.02 to 0.06. Six component static stability force and moment data were measured by an internally compensated internal strain gage balance. Resulting data are presented.
Viscous drag model for plasma electrode boundary effect at ionic impact
We consider the effects of viscosity on the inviscid stability of the Batchelor vortex in a compressible flow. The problem is tackled asymptotically, in the limit of large (streamwise and azimuthal) wavenumbers, together with large Mach numbers. Previous studies, with viscous effects neglected, found that the nature of the solution passes through different regimes as the Mach number increases, relative to the wavenumber. This structure persists when viscous effects are included in the analysis. In the present study the mode present in the incompressible case ceases to be unstable at high Mach numbers and a center mode forms, whose stability characteristics are determined primarily by conditions close to the vortex axis. We find generally that viscosity has a stabilizing influence on the flow, while in the case of center modes, viscous effects become important at much larger Reynolds numbers than for the first class of disturbance.
Analytical and experimental investigation of the characteristics of a three-dimensional turbulent boundary layer in a rotating helical channel. Expressions are developed for the velocity profiles in the inner layer, where viscous effects dominate, and the outer layer, where viscous effects are small. The velocity profiles, wall shear stress, and limiting streamline angles are measured inside the passages of a flat-plate inducer at various radial and chordwise locations using rotating probes. Flow near the blade tip is found to be highly complex, due to interaction of blade boundary layers and annulus wall, resulting in appreciable radial inward flow, as well as a defect in mainstream velocity near the mid-passage. A wall shear stress correlation, which includes the effect of both Reynolds number and rotation parameter, is derived from the measured data.
A transonic flow past a boattailed afterbody under a small angle of attack was examined. It is known that the viscous effect offers significant modifications of the pressure distribution on the afterbody. Thus, the formulation for the inviscid flow was based on the consideration of a flow past a nonaxisymmetric body. The full three dimensional potential equation was solved through numerical relaxation, and quasi-axisymmetric boundary layer calculations were performed to estimate the displacement effect. It was observed again that the viscous effects were not negligible. The trend of the final results agreed well with the experimental data.
Rectangular-to-elliptical shape-transition (REST) inlets are of interest for use on scramjet engines because they are efficient and integrate well with the forebody of a planar vehicle. The classic design technique by Smart for these inlets produces an efficient inlet but the complex three-dimensional viscous effects are only approximately included. Certain undesirable viscous features often occur in these inlets. In the present work, a design toolset has been developed which allows for rapid design of REST-class inlet geometries and the subsequent Navier-Stokes analysis of the inlet performance. This gives the designer feedback on the complex viscous effects at each design iteration. This new tool is applied to design an inlet for on-design operation at Mach 8. The tool allows for rapid investigation of design features that was previously not possible. The outcome is that the inlet shape can be modified to affect aspects of the flow field in a positive way. In one particular example, the boundary layer build-up on the bodyside of the inlet was reduced by 20% of the thickness associated with the classically designed inlet shape.
Analytical and experimental investigations of the characteristics of three dimensional turbulent boundary layers in a rotating helical passage of an inducer rotor are reported. Expressions are developed for the velocity profiles in the inner layer, where the viscous effects dominate, in the outer layer, where the viscous effects are small, and in the interference layer, where the end walls influence the flow. The prediction of boundary layer growth is based on the momentum integral technique. The equations derived are general enough to be valid for all turbomachinery rotors with arbitrary pressure gradients. The experimental investigations are carried out in a flat plate inducer 3 feet in diameter. The mean velocity profiles, turbulence intensities and shear stresses, wall shear stress, and limiting streamline angles are measured at various radial and chordwise locations by using rotating probes. The measurements are in general agreement with the predictions. The radial flows are well represented by an expression which includes the effect of stagger angle and radial pressure gradient. The radial flows in the rotor channel are higher than those on a single blade. The collateral region exists only very near the blade surface. The radial component of turbulence intensity is higher than the streamwise component because of the effect of rotation.
This paper considers a disturbance evolving from a strictly linear finite-growth-rate instability wave, with nonlinear effects first becoming important in the critical layer. By incorporating viscous effects into the nonlinear critical-layer analysis of Goldstein and Leib (1988), it was possible to demonstrate how an initially linear instability wave evolves as it propagates downstream and how the viscous effects eventually become important, even when the viscosity is very small, due to continually decreasing scales generated by the nonlinear effects.