SEARCH · Search NASA
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.
Method of axisymmetric irrotational characteristics used for analyzing supersonic and hypersonic flow of calorically perfect gas through conical nozzles
The calculation of the outer inviscid flow about a rectangular wing moving at supersonic speeds is reported. The inviscid equations of motion governing the flow generated by the wing form a set of hyperbolic differential equations. The flow field about the rectangular wing is separated into three regions consisting of the forebody, the afterbody, and the wing wake. Solutions for the forebody are obtained using conical flow techniques while the afterbody and the wing wake regions are treated as initial value problems. The numerical solutions are compared in the two dimensional regions with known exact solutions.
The unsteady, compressible, thin-layer, Navier-Stokes equations are solved to simulate steady and unsteady, asymmetric, vortical laminar flow around cones at high incidences and supersonic Mach numbers. The equations are solved by using an implicit, upwind, flux-difference splitting, finite-volume scheme. Locally conical flows are assumed and the solutions are obtained by forcing the conserved components of the flowfield vector to be equal at two axial stations located at 0.95 and 1.0. Computational examples cover steady and unsteady asymmetric flows around a circular cone and its control using side strakes. Unsteady asymmetric flows are also presented for elliptic- and diamond-section cones, which model asymmetric vortex shedding around round and sharp-edged delta wing flows.
Supersonic and hypersonic flow of inviscid ideal gas over conical delta wings using three dimensional method of characteristics
Solutions of the time dependent, Reynolds-averaged, Navier-Stokes equations are presented and are compared with a family of experimental results for the three-dimensional interaction of a shock wave with a turbulent boundary layer. The solutions correctly predict the major features of the flow field independent of the shock strength and the extent of separation when using a two-equation turbulence model with wall functions. The experimentally observed boundary between cylindrical and conical flow regimes is also predicted. However, for the flow fields with large separated zones, the details of the measured pressure distributions are not accurately computed, which indicates a need for improved turbulence modeling and/or grid resolution.
The unsteady, compressible, thin-layer, Navier-Stokes (NS) equations are solved to simulate steady and unsteady, asymmetric, vortical laminar flow around cones at high incidences and supersonic Mach numbers. The equations are solved by using an implicit, upwind, flux-difference splitting (FDS), finite-volume scheme. The locally conical flow assumption is used and the solutions are obtained by forcing the conserved components of the flowfield vector to be equal at two axial stations located at 0.95 and 1.0. Computational examples cover steady and unsteady asymmetric flows around a circular cone and its control using side strakes. The unsteady asymmetric flow solution around the circular cone has also been validated using the upwind, flux-vector splitting (FVS) scheme with the thin-layer NS equations and the upwind FDS with the full NS equations. The results are in excellent agreement with each other. Unsteady asymmetric flows are also presented for elliptic- and diamond-section cones, which model asymmetric vortex shedding around round- and sharp-edged delta winds.
This progress report documents the accomplishments achieved in the period from December 1, 1992 until November 30, 1993. These accomplishments include publications, national and international presentations, NASA presentations, and the research group supported under this grant. Topics covered by documents incorporated into this progress report include: active control of asymmetric conical flow using spinning and rotary oscillation; supersonic vortex breakdown over a delta wing in transonic flow; shock-vortex interaction over a 65-degree delta wing in transonic flow; three dimensional supersonic vortex breakdown; numerical simulation and physical aspects of supersonic vortex breakdown; and prediction of asymmetric vortical flows around slender bodies using Navier-Stokes equations.
Solutions of the Euler equations are presented for M = 1.5 flow past a 70-degree-swept delta wing. At an angle of attack of 10 degrees, strong leading-edge vortices are produced. Two computational approaches are taken, based upon fully three-dimensional and conical flow theory. Both methods utilize a finite-volume discretization solved by a pseudounsteady multistage scheme. Results from the two approaches are in good agreement. Computations have been done on a 16-million-word CYBER 205 using 196 x 56 x 96 and 128 x 128 cells for the two methods. A sizable data base is generated, and some of the practical aspects of manipulating it are mentioned. The results reveal many interesting physical features of the compressible vortical flow field and also suggest new areas needing research.
Measurements of the vibrational temperature in a shock layer produced by a cylindrical blunt body in an arc jet nitrogen flow are made for both the neutral molecule and the molecular ion. Spectra were obtained at two locations in the shock layer produced by a 20-cm long cylindrical blunt body inserted in the arc jet conical flow. The technique used for determining the temperatures from the spectral radiation involves obtaining ratios of integrated intensities for certain spectral regions of the measured spectrum and making a comparison with computed ratios as a function of temperature for the same spectral regions. The technique is difficult to apply due to radiation overlapping bands from other species. However, using a trial and error procedure, the vibrational and rotational temperatures at the two locations are inferred. The results indicate that the shock layer is significantly out of thermal equilibrium.
An aerodynamic prediction technique based on the full potential equation in conservation form is developed for the treatment of supersonic flows. This technique bridges the gap between simplistic linear theory methods and complex Euler solvers. A novel local density linearization concept and a second order accurate retarded density scheme, both producing the correct artificial viscosity, are introduced in developing an implicit marching scheme for solving the scalar potential. Results for conical flows over delta wings and a wing-body combination and for non-conical flows over bodies of revolution at angles of attack are compared with Euler and nonconservative full potential calculations and experimental data. The present formulation requires an order of magnitude less computer time and significantly less computer memory over Euler codes and exhibits a considerable improvement in computational efficiency and generality over an existing nonconservative full potential code.
The effectiveness of certain active-control methods for asymmetric flows around circular cones is investigated by using computational solution of the unsteady, compressible full Navier-Stokes equations. Two main methods of active control which include flow injection and surface heating are used. For the flow-injection-control method, flow injection is used either in the normal direction to the surface or in the tangential direction to the surface. For the surface-heating-control method, the temperature of the cone surface is increased. The effectiveness of a hybrid method of flow control which combines normal injection with surface heating has also been studied. The Navier-Stokes equations, subjected to various surface boundary conditions, are solved by using an implicit, upwind, flux-difference splitting, finite-volume scheme for locally-conical flow solutions.
The purpose is to describe the development of unstructured grid methods which have several advantages when compared to methods which make use of structured grids. Unstructured grids, for example, easily allow the treatment of complex geometries, allow for general mesh movement for realistic motions and structural deformations of complete aircraft configurations which is important for aeroelastic analysis, and enable adaptive mesh refinement to more accurately resolve the physics of the flow. Steady Euler calculations for a supersonic fighter configuration to demonstrate the complex geometry capability; unsteady Euler calculations for the supersonic fighter undergoing harmonic oscillations in a complete-vehicle bending mode to demonstrate the general mesh movement capability; and vortex-dominated conical-flow calculations for highly-swept delta wings to demonstrate the adaptive mesh refinement capability are discussed. The basic solution algorithm is a multi-stage Runge-Kutta time-stepping scheme with a finite-volume spatial discretization based on an unstructured grid of triangles in 2D or tetrahedra in 3D. The moving mesh capability is a general procedure which models each edge of each triangle (2D) or tetrahedra (3D) with a spring. The resulting static equilibrium equations which result from a summation of forces are then used to move the mesh to allow it to continuously conform to the instantaneous position or shape of the aircraft. The adaptive mesh refinement procedure enriches the unstructured mesh locally to more accurately resolve the vortical flow features. These capabilities are described in detail along with representative results which demonstrate several advantages of unstructured grid methods. The applicability of the unstructured grid methodology to steady and unsteady aerodynamic problems and directions for future work are discussed.
The component integration of a class of hypersonic high-lift configurations known as waveriders into hypersonic cruise vehicles was evaluated. A wind-tunnel model was developed which integrates realistic vehicle components with two waverider shapes, referred to as the straight-wing and cranked-wing shapes. Both shapes were conical-flow-derived waveriders for a design Mach number of 4.0. Experimental data and limited computational fluid dynamics (CFD) predictions were obtained over a Mach number range of 1.6 to 4.63 at a Reynolds number of 2.0 x 10(exp 6) per foot. The CFD predictions and flow visualization data confirmed the shock attachment characteristics of the baseline waverider shapes and illustrated the waverider flow-field properties. Experimental data showed that no significant performance degradations, in terms of maximum lift-to-drag ratios, occur at off-design Mach numbers for the waverider shapes and the integrated configurations. A comparison of the fully-integrated waverider vehicles to the baseline shapes showed that the performance was significantly degraded when all of the components were added to the waveriders, with the most significant degradation resulting from aftbody closure and the addition of control surfaces. Both fully-integrated configurations were longitudinally unstable over the Mach number range studied with the selected center of gravity location and for unpowered conditions. The cranked-wing configuration provided better lateral-directional stability characteristics than the straight-wing configuration.
A multigrid method using different smoothers has been developed to solve the Euler equations discretized by a nonoscillatory scheme up to fourth order accuracy. The best smoothing property is provided by a five-stage Runge-Kutta technique with optimized coefficients, yet the most efficient smoother is a backward Euler technique in factored and diagonalized form. The singlegrid solution for a hypersonic, viscous conic flow is in excellent agreement with the solution obtained by the third order MUSCL and Roe's method. Mach 8 inviscid flow computations for a complete entry probe have shown that the accuracy is at least as good as the symmetric TVD scheme of Yee and Harten. The implicit multigrid method is four times more efficient than the explicit multigrid technique and 3.5 times faster than the single-grid implicit technique. For a Mach 8.7 inviscid flow over a blunt delta wing at 30 deg incidence, the CPU reduction factor from the three-level multigrid computation is 2.2 on a grid of 37 x 41 x 73 nodes.
A new inverse boundary layer method is developed and applied to incompressible flows with laminar separation and reattachment. Test cases for two dimensional flows are computed and the results are compared with those of other inverse methods. One advantage of the present method is that the calculation of the inviscid velocities may be determined at each marching step without having to iterate. The inverse method was incorporated with the direct method to calculate the incompressible, conical flow over a slender delta wing at incidence. The location of the secondary separation line on the leeward surface of the wing is determined and compared with experiment for a unit aspect ratio wing at 20.5 deg incidence. The viscous flow in the separated region was calculated using prescribed skin friction coefficients.
The unsteady supersonic flow around a rigid sharp-edged delta wing is solved using the unsteady Euler and thin-layer Navier-Stokes equations. The problem is formulated relative to a moving frame of reference so that computation of the grid motion is not required. A three-dimensional, implicit, aproximately factored, central-differencing, finite-volume scheme is used to obtain the time-accurate, unsteady, locally conical flow. The inviscid and viscous results are compared during the periodic response of a wing undergoing forced rolling-oscillation motion. The results show substantial differences in the surface pressure results and other distributed flow characteristics. These differences are attributed to the inviscid and viscous flow predictions in the vortex-dominated regions. The results of the total loads do not show any significant differences.
A review is presented of recent progress in a research program directed towards the development of an improved vortex-flow technology base. It is pointed out that separation induced vortex-flows from the leading and side edges play an important role in the high angle-of-attack aerodynamic characteristics of a wide range of modern aircraft. In the analysis and design of high-speed aircraft, a detailed knowledge of this type of separation is required, particularly with regard to critical wind loads and the stability and performance at various off-design conditions. A description of analytical methods is presented. The theoretical methods employed are divided into two classes which are dependent upon the underlying aerodynamic assumptions. One conical flow method is considered along with three different nonconical flow methods. Comparisons are conducted between the described methods and available aerodynamic data. Attention is also given to a vortex flow drag study and a vortex flow wing design using suction analogy.
This paper reports the results of a comprehensive linear stability analysis of the boundary layer flow over the common research model with natural laminar flow (CRM-NLF) aircraft configuration. The flow conditions match selected test conditions from a recent wind tunnel experiment in the National Transonic Facility at the NASA Langley Research Center. Previous work has shown that the measured onset of laminar-turbulent transition during the experiments can be correlated with the linear amplification of Tollmien-Schlichting (TS) and stationary crossflow (CF) instabilities in the swept wing boundary layer. However, a significant scatter ( N ∈ (4,9)) was observed in the values of the logarithmic amplification factor along the measured transition front. This previous analysis was based on an approximate basic state (based on a boundary layer code with conical flow approximation) and parallel stability computations without surface curvature effects. Here, we examine the effects of these various approximations with the goal of quantifying the resulting changes in the N-factor correlations. Specifically, both linear stability theory (LST) and the parabolized stability equations (PSE)are used in conjunction with an accurate definition of the laminar boundary layer flow as computed with a Navier-Stokes solver with a Reynolds-Averaged-Navier-Stokes (RANS) based turbulence model within the turbulent parts of the flow. Furthermore, the effects of instability wave propagation within a fully three-dimensional boundary layer are also evaluated by integrating the disturbance growth rates along suitably chosen, curvilinear (i.e., nonplanar) propagation trajectories. The results of this analysis are also used in an accompanying paper by Venkatachari et al. to develop improved, physics based transition predictions for the same CRM-NLF configuration.