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.
A lifting surface theory has been developed for a helicopter rotor in forward flight for incompressible flow. The method utilized the concept of the linearized acceleration potential and make use of the vortex lattice procedures. Results in terms of lift coefficient slope for several forward flight conditions are given.
A procedure for solving the conical Euler equations on an adaptively refined mesh is presented, along with a method for determining which cells to refine. The solution procedure is a central-difference cell-vertex scheme. The adaptation procedure is made up of a parameter on which the refinement decision is based, and a method for choosing a threshold value of the parameter. The refinement parameter is a measure of mesh-convergence, constructed by comparison of locally coarse- and fine-grid solutions. The threshold for the refinement parameter is based on the curvature of the curve relating the number of cells flagged for refinement to the value of the refinement threshold. Results for three test cases are presented. The test problem is that of a delta wing at angle of attack in a supersonic free-stream. The resulting vortices and shocks are captured efficiently by the adaptive code.
A combined Eulerian/Lagrangian approach to calculating helicopter rotor flows with concentrated vortices is described. The method computes a general evolving vorticity distribution without any significant numerical diffusion. Concentrated vortices can be accurately propagated over long distances on relatively coarse grids with cores only several grid cells wide. The method is demonstrated for a blade/vortex impingement case in 2D and 3D where a vortex is cut by a rotor blade, and the results are compared to previous 2D calculations involving a fifth-order Navier-Stokes solver on a finer grid.
The linear stability of the trailing line vortex model of Batchelor (1964) is studied using a spectral collocation and matrix eigenvalue method. The entire unstable region in the swirl/axial wavenumber parameter space is mapped out for various azimuthal wavenumbers for both the inviscid and viscous stability problem. The results of the study provide a direct numerical validation of the large-azimuthal-wavenumber asymptotic analysis of Leibovich and Stewartson (1983). It is shown that accurate results are obtained up to azimuthal wavenumbers of 10,000 and greater, and the agreement with the asymptotic theory is excellent.
The simulation of the leading edge vortex flow about a series of conical delta wings through solution of the Navier-Stokes and Euler equations is studied. The occurrence, the validity, and the usefulness of separated flow solutions to the Euler equations of particular interest. Central and upwind difference solutions to the governing equations are compared for a series of cross sectional shapes, including both rounded and sharp tip geometries. For the rounded leading edge and the flight condition considered, viscous solutions obtained with either central or upwind difference methods predict the classic structure of vortical flow over a highly swept delta wing. Predicted features include the primary vortex due to leading edge separation and the secondary vortex due to crossflow separation. Central difference solutions to the Euler equations show a marked sensitivity to grid refinement. On a coarse grid, the flow separates due to numerical error and a primary vortex which resembles that of the viscous solution is predicted. In contrast, the upwind difference solutions to the Euler equations predict attached flow even for first-order solutions on coarse grids. On a sufficiently fine grid, both methods agree closely and correctly predict a shock-curvature-induced inviscid separation near the leeward plane of symmetry. Upwind difference solutions to the Navier-Stokes and Euler equations are presented for two sharp leading edge geometries. The viscous solutions are quite similar to the rounded leading edge results with vortices of similar shape and size. The upwind Euler solutions predict attached flow with no separation for both geometries. However, with sufficient grid refinement near the tip or through the use of more accurate spatial differencing, leading edge separation results. Once the leading edge separation is established, the upwind solution agrees with recently published central difference solutions to the Euler equations.
Counter-rotating vortices, formed by the interaction of film-cooling jets and the hot gas flow, adversely affect the performance of conventional film-cooling designs. Downstream vortex generators have been shown to improve cooling effectiveness by mitigating the effects of the counter-rotating vortices and by deflecting the cooling jet laterally. In this study, computational and experimental methods were used to examine how cylindrical film-cooling holes (D = 3.2 mm, L/D = 6, p/D = 3, α = 30 deg) with and without downstream vortex generators perform when the coolant supply channel is perpendicular to the direction of the hot gas. For this study, the hot gas had a temperature of 650 K and an average Mach number of 0.23. Here, the hot-gas-to-coolant temperature ratio was 1.9, and two blowing ratios (0.75 and 1.0) were studied. Results from the computational fluid dynamics study show how crossflow affects the interaction between the film-cooling jet and hot gas flow with and without downstream vortex generators. The experimental measurements were based on infrared thermography in a conjugate heat transfer environment. Results were obtained for film-cooling performance in terms of overall effectiveness, film effectiveness, and local heat transfer coefficients. The downstream vortex generators can increase the laterally averaged effectiveness by a factor of 1.5 relative to cylindrical holes, but this higher performance is restricted to low crossflow velocities and higher blowing ratios.
A three dimensional, forward-marching, viscous flow analysis is applied to the tip vortex generation problem. The equations include a streamwise momentum equation, a streamwise vorticity equation, a continuity equation, and a secondary flow stream function equation. The numerical method used combines a consistently split linearized scheme for parabolic equations with a scalar iterative ADI scheme for elliptic equations. The analysis is used to identify the source of the tip vortex generation process, as well as to obtain detailed flow results for a rectangular planform wing immersed in a high Reynolds number free stream at 6 degree incidence.
The problem of forced fluid sloshing in a partially filled spinning spherical tank is solved numerically using the finite element method. The governing equations include Coriolis acceleration, empirical fluid damping and spatially homogeneous vorticity first introduced by Pfeiffer. An exponential instability similar to flutter is detected in the present simulation for fill ratios below 50 percent. This instability appears in the model as a result of the homogeneous vortex assumption since the free slosh equations are neutrally stable in the Liapunov sense.
This paper will review the advances made recently in the Navier-Stokes CFD methods to simulate aerodynamics and aeroacoustics of helicopter rotors and rotor-body flows. Although a complete flowfield simulation of full helicopter is currently not feasible with these methods, impressive gains have been made in analyzing individual components of this complex problem in a very detailed manner. The use of the state-of-the-art numerical algorithms in solution methods, in conjunction with powerful supercomputers, like the Cray-2, have enabled noticeable progress to be made in modeling viscous-inviscid interactions, blade-vortex interactions, tip-vortex: simulation and wake effects, as well as high speed impulsive noise in hover and forward flight for isolated rotor blades. This paper will critically evaluate the presently available Euler and Navier-Stokes methods, both finite-difference and finite volume methods using structured and unstructured grids for helicopter applications for accuracy, suitability, and computational efficiency. The review will also include the recent progress made using overset grids to model rotor-body flows. All the material for this review will be drawn from the published material shown below.
A low-speed wind-tunnel study of the flow about a 76/40-deg double-delta wing is described for angles of attack ranging from -10 to 25 deg and Reynolds numbers ranging from 0.5 to 1.5 Million. The study was conducted to provide data for the purpose of understanding the vortical flow behavior and for validating Computational Fluid Dynamics methods. Flow visualization tests have provided insight into the effect of the angle of attack and Reynolds number of the vortex-dominated flow both on and off of the surface of the double-delta wing. Upper surface pressure recordings from pressure orifices and Pressure Sensitive Paint have provided data on the pressures induced by the vortices. Flowfield surveys were carried out at an angle of attack of 10 deg by using a thin 5-hole probe. Numerical solutions of the compressible thin-layer Navier-Stokes equations were conducted and compared to the experimental data.
Flow separation control on an adverse-pressure-gradient ramp model was investigated using various flow-control methods in the NASA Langley 15-Inch Wind Tunnel. The primary flow-control method studied used a sweeping jet actuator system to compare with more classic flow-control techniques such as micro-vortex generators, steady blowing, and steady- and unsteady-vortex generating jets. Surface pressure measurements and a new oilflow visualization technique were used to characterize the effects of these flow-control actuators. The sweeping jet actuators were run in three different modes to produce steady-straight, steady-angled, and unsteady-oscillating jets. It was observed that all of these flow-control methods are effective in controlling the separated flows on the ramp model. The steady-straight jet energizes the boundary layer by momentum addition and was found to be the least effective method for a fixed momentum coefficient. The steady-angled jets achieved better performance than the steady-straight jets because they generate streamwise vortices that energize the boundary layer by mixing high-momentum fluid with near wall low-momentum fluid. The unsteady-oscillating jets achieved the best performance by increasing the pressure recovery and reducing the downstream flow separation. Surface flow visualizations indicated that two out-of-phase counter-rotating vortices are generated per sweeping jet actuator, while one vortex is generated per vortex-generating jets. The extra vortex resulted in increased coverage, more pressure recovery, and reduced flow separation.
Small scale models of agricultural airplanes were tested and numerical methods were utilized to study interactions between the airplane wake and the dispersed spray and granular materials. Methods were developed to measure and predict dispersal transport and wake characteristics and dispersal techniques to obtain interactions more favorable to wide, uniform deposition patterns and reduced drift. In the full scale wind tunnel, full scale agricultural airplanes and dispersal systems for both liquid and solid applications were evaluated to improve aircraft aerodynamics and dispersal systems efficiency. The program status in these two facilities is reported with emphasis on wake interactions and dispersal systems research.
The goal of this current research is to establish a database that validate/calibrate wake encounter analysis methods for fleet-wide application; and measure/document atmospheric effects on wake decay. Two kinds of experiments, wind tunnel experiments and flight experiments, are performed. This paper discusses the different types of tests and compares their wake velocity measurement.
As the result of global, non-parallel flow stability analysis the single value of the disturbance growth-rate and respective frequency is obtained. This complex value characterizes the stability of the whole flow configuration and is not referred to any particular flow pattern. The global analysis assures that all the flow elements (wake, boundary and shear layer) are taken into account. The physical phenomena connected with the wake instability are properly reproduced by the global analysis. This enhances the investigations of instability of any 2-D flows, including ones in which the boundary layer instability effects are known to be of dominating importance. Assuming fully 2-D disturbance form, the global linear stability problem is formulated. The system of partial differential equations is solved for the eigenvalues and eigenvectors. The equations, written in the pure stream function formulation, are discretized via FDM using a curvilinear coordinate system. The complex eigenvalues and corresponding eigenvectors are evaluated by an iterative method. The investigations performed for various Reynolds numbers emphasize that the wake instability develops into the Karman vortex street. This phenomenon is shown to be connected with the first mode obtained from the non-parallel flow stability analysis. The higher modes are reflecting different physical phenomena as for example Tollmien-Schlichting waves, originating in the boundary layer and having the tendency to emerge as instabilities for the growing Reynolds number. The investigations are carried out for a circular cylinder, oblong ellipsis and airfoil. It is shown that the onset of the wake instability, the waves in the boundary layer, the shear layer instability are different solutions of the same eigenvalue problem, formulated using the non-parallel theory. The analysis offers large potential possibilities as the generalization of methods used till now for the stability analysis.
A discrete-vortex model was developed to account for the hazardous effects of the vortex trail issued from the edges of separation of a large leading wing on a small trailing wing. The model is divided into three main parts: the leading wing and its near wake, the near and far wakes of the leading wing, and the trailing wing and the portion of the far wake in its vicinity. The normal force, pitching moment, and rolling moment coefficients for the trailing wing are calculated. The circulation distribution in the vortex trail is calculated in the first part of the model where the leading wing is far upstream and hence is considered isolated. A numerical example is solved to demonstrate the feasibility of using this method to study interference between aircraft. The numerical results show the correct trends: The following wing experiences a loss in lift between the wing-tip vortex systems of the leading wing, a gain outside this region, and strong rolling moments which can change sign as the lateral relative position changes. All the results are strongly dependent on the vertical relative position.
The technique of Diode Laser Spectroscopy (DLS) with wavelength modulation is utilized to measure the concentration of methane in reacting vortex rings under microgravity conditions. From the measured concentration of methane, other major species such as water, carbon dioxide, nitrogen, and oxygen can be easily computed under the assumption of equilibrium chemistry with an iterative method called ITAC (Iterative Temperature with Assumed Chemistry). The conserved scalar approach in modelling the coupling between fluid dynamics and combustion is utilized to represent the unknown variables in terms of the mixture fraction and scalar dissipation rate in conjunction with ITAC. Post-processing of the DLS and the method used to compute the species concentration are discussed. From the flame luminosity results, ring circulation appears to increase the fuel consumption rate inside the reacting vortex ring and the flame height for cases with similar fuel volumes but different ring circulations. The concentrations of methane, water, and carbon dioxide agree well with available results from numerical simulations.
The author has identified the following significant results. LANDSAT MSS data of a three year time lapse indicate change of sea surface condition in Seto Inland Sea and coastal region. The red tide which formerly concentrated in the bay or inland sea now extends into an open sea. A small ocean vortex similar to mesoscale atmospheric vortex is revealed by the band 4 image of the satellite data. A manual photographic method applied to a single band image of MSS is effective in detecting sea surface pollution.
A computational code is developed for the integro-differential equations governing the motion of the centerlines of vortex filaments submerged in a background potential flow. These equations, which are derived from the method of matched asymptotic analysis, include the effect of the decaying large-magnitude circumferential and axial velocity components in the vortical cores. Numerical examples are presented to assess the effect of a large axial velocity and that of nonsimilar initial profiles in the vortical cores. The initial configurations of the filaments are chosen so as to fulfill the basic assumption of the asymptotic analysis, which is that the effective vortical core size is much smaller than all the other length scales in the flowfield, e.g., the radius of curvature and the interfilament distance. The computations are continued until the basic assumption is no longer valid, that is when the merging or intersection of filaments has begun. A classification of the various types of local or global merging or intersection of filaments is made and demonstrated by numerical examples. It is then shown that the asymptotic solution not only provides the initial data but also can be used to formulate the appropriate boundary conditions for the numerical solution of a merged region.