Wind tunnel investigations of vortex breakdown on slender sharp edged wings Final report
Wind tunnel tests of vortex breakdown on sharp edged delta wings
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.
Wind tunnel tests of vortex breakdown on sharp edged delta wings
Axisymmetric flow of rotating stream analyzed to determine conditions allowing isolated vortex breakdown to develop
Analysis and computation of steady, compressible, quasi-axisymmetric flow of an isolated, slender vortex are considered. The compressible, Navier-Stokes equations are reduced to a simpler set by using the slenderness and quasi-axisymmetry assumptions. The resulting set along with a compatibility equation are transformed from the diverging physical domain to a rectangular computational domain. Solving for a compatible set of initial profiles and specifying a compatible set of boundary conditions, the equations are solved using a type-differencing scheme. Vortex breakdown locations are detected by the failure of the scheme to converge. Computational examples include isolated vortex flows at different Mach numbers, external axial-pressure gradients and swirl ratios.
In this paper, steady, axisymmetric inviscid, and viscous (laminar) swirling flows representing vortex breakdown phenomena are simulated using a stream function-vorticity-circulation formulation and two numerical methods. The first is based on an inverse iteration, where a norm of the solution is prescribed and the swirling parameter is calculated as a part of the output. The second is based on direct Newton iterations, where the linearized equations, for all the unknowns, are solved simultaneously by an efficient banded Gaussian elimination procedure. Several numerical solutions for inviscid and viscous flows are demonstrated, followed by a discussion of the results. Some improvements on previous work have been achieved: first order upwind differences are replaced by second order schemes, line relaxation procedure (with linear convergence rate) is replaced by Newton's iterations (which converge quadratically), and Reynolds numbers are extended from 200 up to 1000.
The incompressible axisymmetric steady Navier-Stokes equations in primitive variables are used to simulate vortex breakdown. The equations, discretized using a second-order, central-difference scheme, are linearized and then solved using an exact LU decomposition, Gaussian elimination, and Newton iteration. Solutions are presented for Reynolds numbers, based on vortex-core radius, as high as 1500. An attempt to study the stability of the axisymmetric solutions against three-dimensional perturbations is discussed.
The breakdown of a vortex flow in a tube is studied for a slightly diverging tube by means of a numerical solution of the Navier-Stokes equations for a quasi-stationary, incompressible, laminar flow. Numerical results are compared to experiments of Faler and Leibovich. The numerical data display features similar to the experimental results concerning breakdown conditions, location, and structure.
The breakdown of a vortex flow in a tube is studied for a slightly diverging tube by means of a numerical solution of the Navier-Stokes equations for a quasi-stationary, incompressible, laminar flow. Numerical results are compared to experiments of Faler and Leibovich. The numerical data display features similar to the experimental results concerning breakdown conditions, location, and structure.
A method of weighted residuals for the computation of rotationally symmetric quasi-cylindrical viscous incompressible vortex flow is presented and used to compute a wide variety of vortex flows. The method approximates the axial velocity and circulation profiles by series of exponentials having (N + 1) and N free parameters, respectively. Formal integration results in a set of (2N + 1) ordinary differential equations for the free parameters. The governing equations are shown to have an infinite number of discrete singularities corresponding to critical values of the swirl parameters. The computations point to the controlling influence of the inner core flow on vortex behavior. They also confirm the existence of two particular critical swirl parameter values: one separates vortex flow which decays smoothly from vortex flow which eventually breaks down, and the second is the first singularity of the quasi-cylindrical system, at which point physical vortex breakdown is thought to occur.
The breakdown of an isolated axisymmetric vortex embedded in an unbounded uniform flow is examined by numerical integration of the complete Navier-Stokes equations for unsteady axisymmetric flow. Results show that if the vortex strength is small, the solution approaches a steady flow and the vortex is stable. If the strength is large enough, the solution remains unsteady and a recirculating zone will appear near the axis, its form and internal structure resembling those of the axisymmetric breakdown bubbles with multi-cells observed by Faler and Leibovich (1978). For apppropriate combinations of flow parameters, the flow reveals quasi-periodicity. Parallel calculations with the quasi-cylindrical approximation indicate that so far as predicting of breakdown is concerned, its results coincide quite well with the results mentioned above. Both show that the vortex breakdown has little concern with the Reynolds number or with the critical classification of the upstream flow, at least for the lower range of Reynolds numbers.
The buffet response of the twin-tail configuration of the F/A-18 aircraft; a multidisciplinary problem, is investigated using three sets of equations on a multi-block grid structure. The first set is the unsteady, compressible, full Navier-Stokes equations. The second set is the coupled aeroelastic equations for bending and torsional twin-tail responses. The third set is the grid-displacement equations which are used to update the grid coordinates due to the tail deflections. The computational model consists of a 76 deg-swept back, sharp edged delta wing of aspect ratio of one and a swept-back F/A-18 twin-tails. The configuration is pitched at 32 deg angle of attack and the freestream Mach number and Reynolds number are 0.2 and 0.75 x 10(exp 6) respectively. The problem is solved for the initial flow conditions with the twin tail kept rigid. Next, the aeroelastic equations of the tails are turned on along with the grid-displacement equations to solve for the uncoupled bending and torsional tails response due to the unsteady loads produced by the vortex breakdown flow of the vortex cores of the delta wing. Two lateral locations of the twin tail are investigated. These locations are called the midspan and inboard locations.
The effects of freestream Mach number and angle of attack on the leading-edge vortex breakdown due to the terminating shock on a 65-degree, sharp-edged, cropped delta wing are investigated computationally, using the time-accurate solution of the laminar unsteady compressible full Navier-Stokes equations with the implicit upwind flux-difference splitting, finite-volume scheme. A fine O-H grid consisting of 125 x 85 x 84 points in the wrap-around, normal, and axial directions, respectively, is used for all the flow cases. Keeping the Reynolds number fixed at 3.23 x 10 exp 6, the Mach number is varied from 0.85 to 0.9 and the angle of attack is varied from 20 to 24 deg. The results show that, at 20-deg angle of attack, the increase of the Mach number from 0.85 to 0.9 results in moving the location of the terminating shock downstream. The results also show that, at 0.85 Mach number, the increase of the angle of attack from 20 to 24 deg results in moving the location of the terminating shock upstream. The results are in good agreement with the experimental data.
Flow-field measurements have been made to determine the effects of core blowing on vortex breakdown and control. The results of these proof-of-concept experiments clearly demonstrate the usefulness of water tunnels as test platforms for advanced flow-field simulation and measurement.
Flow-field measurements have been made to determine the effects of core blowing on vortex breakdown and control. The results of these proof-of-concept experiments clearly demonstrate the usefulness of water tunnels as test platforms for advanced flow-field simulation and measurement.
The results of an experimental study of spiraling flows in a slightly diverging, circular duct are reported. Seven types of flow disturbances were observed. In addition to the spiral and axisymmetric vortex breakdowns and the double helix mode, four other forms were identified and are reported. The type and axial location of the disturbance depended on the Reynolds and circulation numbers of the flow. Detailed velocity measurements were made by using a laser Doppler anemometer. Measurements made far upstream of any disturbance showed that the introduction of swirl resulted in the formation of a high axial velocity jet centered around the vortex center. A mapping of the velocity field of a so-called axisymmetric breakdown, formed at a Reynolds number of 2560, revealed that the recirculation zone is a two-celled structure, with four stagnation points on the vortex axis marking the axial extremes of the concentric cells. The dominant feature of the flow inside the bubble was the strong, periodic velocity fluctuations. Existing theoretical models do not predict the two-celled structure and the temporal velocity fluctuations that were observed.
Sinusoidal instability of aircraft trailing vortices was induced by differentially oscillating inboard and outboard flaps on a model wing in a ship towing tank. Measured amplification rates qualitatively agree with theoretical predictions. Axial velocities in the vortex core were directed towards the wing and measured to be approximately 25% of the tow speed. Vortex breakdown was observed to occur along trailers undergoing sinusoidal instability near but ahead (towards the wing) of positions of maximum trailer separation. It is shown that axisymmetric pressure gradients are imposed along the vortex core by the other sinusoidally deformed trailer and are responsible for the observed changes in core diameter. A theoretical model predicts the short-time behavior of the vortex core and shows that the response of the vortex depends crucially on the axial velocity in the vortex.
The integral form of the complete, unsteady, compressible, three-dimensional Navier-Stokes equations in the conservation form, cast in generalized coordinate system, are solved, numerically, to simulate the vortex breakdown phenomenon. The inviscid fluxes are discretized using Roe's upwind-biased flux-difference splitting scheme and the viscous fluxes are discretized using central differencing. Time integration is performed using a backward Euler ADI (alternating direction implicit) scheme. A full approximation multigrid is used to accelerate the convergence to steady state.
An inviscid model for a steady axisymmetric flow with swirl was studied. The governing equation is a nonlinear elliptic equation which has more than one solution for a certain range of the swirl parameter. The physically interesting solutions have closed streamlines that look like vortex breakdown (bubble-like solutions). A multigrid method is used to find these solutions. Using an FMG algorithm (nested iteration), the problem is solved in just a few multigrid cycles.
An extraordinary warming of the stratosphere in December-January 1976-77 was followed by tropospheric warming in the polar region and cooling in middle latitudes. During January 10-20, the associated polar anticyclone extended from the surface to 10 mb. Antecedents of the polar vortex breakdown are reviewed with the aid of results of zonal-harmonic analyses of planetary waves, for heights of the pressure surfaces (700-10 mb), temperature, and mean stratospheric temperature (the latter determined from satellite radiation measurements). Wave 1 in height and temperature played a dominant role in the stratosphere, attaining amplitudes of 1600 gpm and 25 C, respectively, at 10 mb. On the other hand, superposition of retrogressing wave 1 and quasi-stationary wave 2 in the height of the 300-mb surface, with individual amplitudes exceeding 300 gpm, is judged to have been an important factor in the overall development.