The problem of roughness drag at supersonic speeds
Methods of estimating roughness drag at supersonic speeds
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.
Methods of estimating roughness drag at supersonic speeds
Supersonic aircraft design is a computationally expensive optimization problem and multifidelity approaches over a significant opportunity to reduce design time and computational cost. This report presents tools developed to improve supersonic aircraft design capabilities including: aerodynamic tools for supersonic aircraft configurations; a systematic way to manage model uncertainty; and multifidelity model management concepts that incorporate uncertainty. The aerodynamic analysis tools developed are appropriate for use in a multifidelity optimization framework, and include four analysis routines to estimate the lift and drag of a supersonic airfoil, a multifidelity supersonic drag code that estimates the drag of aircraft configurations with three different methods: an area rule method, a panel method, and an Euler solver. In addition, five multifidelity optimization methods are developed, which include local and global methods as well as gradient-based and gradient-free techniques.
The NASA Robotic Conjunction Assessment Risk Analysis (CARA) team sends ephemeris data to the Joint Space Operations Center (JSpOC) for conjunction assessment screening against the JSpOC high accuracy catalog and then assesses risk posed to protected assets from predicted close approaches. Since most spacecraft supported by the CARA team are located in LEO orbits, atmospheric drag is the primary source of state estimate uncertainty. Drag magnitude and uncertainty is directly governed by atmospheric density and thus space weather. At present the actual effect of space weather on atmospheric density cannot be accurately predicted because most atmospheric density models are empirical in nature, which do not perform well in prediction. The Jacchia-Bowman-HASDM 2009 (JBH09) atmospheric density model used at the JSpOC employs a solar storm active compensation feature that predicts storm sizes and arrival times and thus the resulting neutral density alterations. With this feature, estimation errors can occur in either direction (i.e., over- or under-estimation of density and thus drag). Although the exact effect of a solar storm on atmospheric drag cannot be determined, one can explore the effects of JBH09 model error on conjuncting objects' trajectories to determine if a conjunction is likely to become riskier, less risky, or pass unaffected. The CARA team has constructed a Space Weather Trade-Space tool that systematically alters the drag situation for the conjuncting objects and recalculates the probability of collision for each case to determine the range of possible effects on the collision risk. In addition to a review of the theory and the particulars of the tool, the different types of observed output will be explained, along with statistics of their frequency.
A computational study of the wing for the distributed electric propulsion X-57 Maxwell airplane configuration at cruise and takeoff/landing conditions was completed. Two unstructured-mesh, Navier-Stokes computational fluid dynamics methods, FUN3D and USM3D, were used to predict the wing performance. The goal of the X-57 wing and distributed electric propulsion system design was to meet or exceed the required lift coefficient 3.95 for a stall speed of 58 knots, with a cruise speed of 150 knots at an altitude of 8,000 ft. The X-57 Maxwell airplane was designed with a small, high aspect ratio cruise wing that was designed for a high cruise lift coefficient (0.75) at angle of attack of 0deg. The cruise propulsors at the wingtip rotate counter to the wingtip vortex and reduce induced drag by 7.5 percent at an angle of attack of 0.6deg. The unblown maximum lift coefficient of the high-lift wing (with the 30deg flap setting) is 2.439. The stall speed goal performance metric was confirmed with a blown wing computed effective lift coefficient of 4.202. The lift augmentation from the high-lift, distributed electric propulsion system is 1.7. The predicted cruise wing drag coefficient of 0.02191 is 0.00076 above the drag allotted for the wing in the original estimate. However, the predicted drag overage for the wing would only use 10.1 percent of the original estimated drag margin, which is 0.00749.
This paper presents some practical methods for the aerodynamic design and evaluation of conventional aircraft. High-lift methodology which provides improved takeoff and landing and transonic maneuvering performance is discussed. Also, new techniques for estimating and minimizing cruise pressure drag are presented. These include a far-field theory to minimize trimmed induced drag, theories to estimate the spanwise variation of drag due to thickness and lift, and a far-field theory to estimate total pressure drag. In addition to the description of methods, aerodynamic design procedures are outlined and results from both the design and evaluation methods are presented.
Explore the source record for details and available documents.
Despite the 80-year history of classical wing theory, considerable research has recently been directed toward planform and wake effects on induced drag. Nonlinear interactions between the trailing wake and the wing offer the possibility of reducing drag. The nonlinear effect of compressibility on induced drag characteristics may also influence wing design. This thesis deals with the prediction of these nonlinear aspects of induced drag and ways to exploit them. A potential benefit of only a few percent of the drag represents a large fuel savings for the world's commercial transport fleet. Computational methods must be applied carefully to obtain accurate induced drag predictions. Trefftz-plane drag integration is far more reliable than surface pressure integration, but is very sensitive to the accuracy of the force-free wake model. The practical use of Trefftz plane drag integration was extended to transonic flow with the Tranair full-potential code. The induced drag characteristics of a typical transport wing were studied with Tranair, a full-potential method, and A502, a high-order linear panel method to investigate changes in lift distribution and span efficiency due to compressibility. Modeling the force-free wake is a nonlinear problem, even when the flow governing equation is linear. A novel method was developed for computing the force-free wake shape. This hybrid wake-relaxation scheme couples the well-behaved nature of the discrete vortex wake with viscous-core modeling and the high-accuracy velocity prediction of the high-order panel method. The hybrid scheme produced converged wake shapes that allowed accurate Trefftz-plane integration. An unusual split-tip wing concept was studied for exploiting nonlinear wake interaction to reduced induced drag. This design exhibits significant nonlinear interactions between the wing and wake that produced a 12% reduction in induced drag compared to an equivalent elliptical wing at a lift coefficient of 0.7. The performance of the split-tip wing was also investigated by wing tunnel experiments. Induced drag was determined from force measurements by subtracting the estimated viscous drag, and from an analytical drag-decomposition method using a wake survey. The experimental results confirm the computational prediction.
The results of calculations necessary to extrapolate performance data on an XB-70-1 wind tunnel model to full scale at Mach numbers from 0.76 to 2.53 are presented. The extrapolation was part of a joint program to evaluate performance prediction techniques for large flexible supersonic airplanes similar to a supersonic transport. The extrapolation procedure included: interpolation of the wind tunnel data at the specific conditions of the flight test points; determination of the drag increments to be applied to the wind tunnel data, such as spillage drag, boundary layer trip drag, and skin friction increments; and estimates of the drag items not represented on the wind tunnel model, such as bypass doors, roughness, protuberances, and leakage drag. In addition, estimates of the effects of flexibility of the airplane were determined.
Direct measurements have been made of the drag of a special test body and its stabilizing tail surfaces throughout free drops from high altitudes. The data obtained have been used to establish the relation between the drag coefficient and the Mach number for the body and for the tail surfaces over a range of Mach numbers from 0.85 to 1.15. For bodies of the form tested, the drag per square foot of frontal area increased abruptly from about 3 percent of atmospheric pressure at a Mach number of 0.95 to 17 percent of atmospheric pressure at a Mach number of 1.00, then linearly with Mach number to 28 percent of atmospheric pressure at a Mach number of approximately 1.15. Some doubt exists as to the applicability of the tail drag results to the estimation of wing drag at transonic speeds because of the possibility of appreciable interference effects between the vertical and the horizontal surfaces and between the body and the tail surfaces. Insofar as they are applicable, the tail drag results indicated that with symmetrical 6-percent-thick area may be expected to increase abruptly from 4 percent of atmospheric pressure at a Mach number of 0.88 to 36 percent of atmospheric pressure at a Mach number of 1.00, then linearly with Mach number to approximately 50 percent of atmospheric pressure at a Mach number of 1.15.
A detailed table of the standard equilibrium velocity and standard equilibrium time is presented for bodies falling in the standard atmosphere. This table gives the velocity at various altitudes and the time of fall from sea level to -4000 feet and from 80,000 feet to sea level. In addition to this standard table, there are given short tables and charts of an open-parachute descent and free-falls; the terminal velocity at sea level, and the variation of the weight-to-drag ratio (2w/cds)1/2 for various weight jumpers from 90 to 30 feet in open-parachute descent; and estimations of drag coefficients of silk and nylon parachutes. The table of standard equilibrium velocities and standard equilibrium times may be used directly for open-parachute descents, given the weight of the jumper, the diameter of the parachute, and the drag coefficient. For free-falls starting from horizontal flight, approximately 14 seconds must be added to the equilibrium time given in the table to obtain the total time to sea level. (author)
A coordinated experimental and numerical simulation effort is carried out to improve our understanding of the physics of acoustic liners in a grazing flow as well our computational aeroacoustics (CAA) method prediction capability. A numerical simulation code based on advanced CAA methods is developed. In a parallel effort, experiments are performed using the Grazing Flow Impedance Tube at the NASA Langley Research Center. In the experiment, a liner is installed in the upper wall of a rectangular flow duct with a 2 inch by 2.5 inch cross section. Spatial distribution of sound pressure levels and relative phases are measured on the wall opposite the liner in the presence of a Mach 0.3 grazing flow. The computer code is validated by comparing computed results with experimental measurements. Good agreements are found. The numerical simulation code is then used to investigate the physical properties of the acoustic liner. It is shown that an acoustic liner can produce self-noise in the presence of a grazing flow and that a feedback acoustic resonance mechanism is responsible for the generation of this liner self-noise. In addition, the same mechanism also creates additional liner drag. An estimate, based on numerical simulation data, indicates that for a resonant liner with a 10% open area ratio, the drag increase would be about 4% of the turbulent boundary layer drag over a flat wall.
The spacecraft of the Geopotential Research Mission (GRM) are cylindrical in form and designed to fly with their longitudinal axes parallel to their direction of flight. The ratio of length to diameter of these spacecraft is roughly equal to 5.0. Other spacecraft previously flown had corresponding ratios roughly equal to 1.0, and therefore the drag produced by impacts on the lateral surfaces of those spacecraft was not as large as it will be on the GRM spacecraft. Since the drag coefficient is essentially the drag force divided by the frontal area in flight, lateral impacts, when taken into account make the GRM drag coefficient significantly larger than the coefficients used before for shorter spacecraft. A simple formula is derived for the drag coefficient of a cylindrical body flying with its long axis along the direction of flight, and it is used to estimate the drag for the GRM. The formula shows that the drag due to lateral surface impacts depends on the ratio of length-to-diameter and on a coefficient C sub LS (lateral surface impact coefficient) which can be determined from previous cylindrical spacecraft flown with the same attitude, or can be obtained from laboratory measurements of momentum accommodation coefficients.
Within the first few seconds after launch of STS-124, debris traveling vertically near the vehicle was captured on two 16-mm film cameras surrounding the launch pad. One particular piece of debris caught the attention of engineers investigating the release of the flame trench fire bricks. The question to be answered was if the debris was a fire brick, and if it represented the first bricks that were ejected from the flame trench wall, or was the object one of the pieces of debris normally ejected from the vehicle during launch. If it was typical launch debris, such as SRB throat plug foam, why was it traveling vertically and parallel to the vehicle during launch, instead of following its normal trajectory, flying horizontally toward the north perimeter fence? By utilizing the Runge-Kutta integration method for velocity and the Verlet integration method for position, a method that suppresses trajectory computational instabilities due to noisy position data was obtained. This combination of integration methods provides a means to extract the best estimate of drag force and drag coefficient under the non-ideal conditions of limited position data. This integration strategy leads immediately to the best possible estimate of object density, within the constraints of unknown particle shape. These types of calculations do not exist in readily available off-the-shelf simulation software, especially where photogrammetry data is needed as an input.
Two methods for constructing performance functions for formation fight-for-drag-reduction suitable for use with an extreme-seeking control system are presented. The first method approximates an a prior measured or estimated drag-reduction performance function by combining real-time measurements of readily available parameters. The parameters are combined with weightings determined from a minimum squares optimization to form a blended performance function.
Tests were made in the NASA-Ames 40- by 80 Foot Wind Tunnel of a wing semispan with a nacelle (no propeller) from a typical, general aviation twin-engine aircraft. Measurements were made of the effect on drag of the cooling air flow through the nacelle. Internal and external nacelle pressures were measured. It was found that the cooling flow accounts for about 13% of the estimated airplane drag and about 42% of the cooling flow drag is associated with the internal flow. It was concluded that improvements could be made by relocating both the inlet and the outlet of the cooling air.
Tests were made in the Ames 40 by 80 ft Wind Tunnel of a semispan wing with a nacelle (no propeller) from a typical, general aviation twin-engine aircraft. Measurements were made of the effect on drag of the flow of cooling air through the nacelle. Internal and external nacelle pressures were measured. It was found that the cooling airflow accounts for about 13% of the total estimated airplane drag during both cruise and climb. The now of cooling air through the nacelle accounts for 30% of the airflow drag component during cruise and 42% during climb; the balance, in both cruise and climb, is attributed to [he external shape of the nacelle. It was suggested that improvements could possibly be made by relocating both the inlet and the outlet for the cooling air.
Included are results of an experimental investigation of low-speed turbulent flow over multiple two-dimensional transverse rigid wavy surfaces having a wavelength on the order of the boundary-layer thickness. Data include surface pressure and total drag measurements on symmetric and asymmetric wall waves under a low-speed turbulent boundary-layer flow. Several asymmetric wave configurations exhibited drag levels below the equivalent symmetric (sine) wave. The experimental results compare favorably with numerical predictions from a Reynolds-averaged Navier-Stokes spectral code. The reported results are of particular interest for the estimation of drag, the minimization of fabrication waviness effects, and the study of wind-wave interactions.
Estimation and control methods for a Drag-Free spacecraft are discussed. The functional and analytical synthesis of on-board estimators and controllers for an integrated attitude and translation control system is represented. The framework for detail definition and design of the baseline drag-free system is created. The techniques for solution of self-gravity and electrostatic charging problems are applicable generally, as is the control system development.