Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “upwind”

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.

At least 757 records · Page 42

Calculations of separated 3-D flows with a pressure-staggered Navier-Stokes equations solver

A Navier-Stokes equations solver based on a pressure correction method with a pressure-staggered mesh and calculations of separated three-dimensional flows are presented. It is shown that the velocity pressure decoupling, which occurs when various pressure correction algorithms are used for pressure-staggered meshes, is caused by the ill-conditioned discrete pressure correction equation. The use of a partial differential equation for the incremental pressure eliminates the velocity pressure decoupling mechanism by itself and yields accurate numerical results. Example flows considered are a three-dimensional lid driven cavity flow and a laminar flow through a 90 degree bend square duct. For the lid driven cavity flow, the present numerical results compare more favorably with the measured data than those obtained using a formally third order accurate quadratic upwind interpolation scheme. For the curved duct flow, the present numerical method yields a grid independent solution with a very small number of grid points. The calculated velocity profiles are in good agreement with the measured data.

Kim, S.-W.↗

Euler solutions to nonlinear acoustics of non-lifting hovering rotor blades

For the first time a computational fluid dynamics (CFD) method is used to calculate directly the high-speed impulsive (HSI) noise of a non-lifting hovering rotor blade out to a distance of over three rotor radii. In order to accurately propagate the acoustic wave in a stable and efficient manner, an implicit upwind-biased Euler method is solved on a grid with points clustered along the line of propagation. A detailed validation of the code is performed for a rectangular rotor blade at tip Mach numbers ranging from 0.88 to 0.92. The agreement with experiment is excellent at both the sonic cylinder and at 2.18 rotor radii. The agreement at 3.09 rotor radii is still very good, showing improvements over the results from the best previous method. Grid sensitivity studies indicate that with special attention to the location of the boundaries a grid with approximately 60,000 points is adequate. This results in a computational time of approximately 40 minutes on a Cray-XMP. The practicality of the method to calculate HSI noise is demonstrated by expanding the scope of the investigation to examine the rectangular blade as well as a highly swept and tapered blade over a tip Mach number range of 0.80 to 0.95. Comparisons with experimental data are excellent and the advantages of planform modifications are clearly evident. New insight is gained into the mechanisms of nonlinear propagation and the minimum distance at which a valid comparison of different rotors can be made: approximately two rotor radii from the center of rotation.

Baeder, J. D.↗

Spatial adaption procedures on unstructured meshes for accurate unsteady aerodynamic flow computation

Spatial adaption procedures for the accurate and efficient solution of steady and unsteady inviscid flow problems are described. The adaption procedures were developed and implemented within a two-dimensional unstructured-grid upwind-type Euler code. These procedures involve mesh enrichment and mesh coarsening to either add points in a high gradient region or the flow or remove points where they are not needed, respectively, to produce solutions of high spatial accuracy at minimal computational costs. A detailed description is given of the enrichment and coarsening procedures and comparisons with alternative results and experimental data are presented to provide an assessment of the accuracy and efficiency of the capability. Steady and unsteady transonic results, obtained using spatial adaption for the NACA 0012 airfoil, are shown to be of high spatial accuracy, primarily in that the shock waves are very sharply captured. The results were obtained with a computational savings of a factor of approximately fifty-three for a steady case and as much as twenty-five for the unsteady cases.

Rausch, Russ D.↗

Numerical solutions of 2-D multi-stage rotor/stator unsteady flow interactions

The Rai method of single-stage rotor/stator flow interaction is extended to handle multistage configurations. In this study, a two-dimensional Navier-Stokes multi-zone approach was used to investigate unsteady flow interactions within two multistage axial turbines. The governing equations are solved by an iterative, factored, implicit finite-difference, upwind algorithm. Numerical accuracy is checked by investigating the effect of time step size, the effect of subiteration in the Newton-Raphson technique, and the effect of full viscous versus thin-layer approximation. Computer results compared well with experimental data. Unsteady flow interactions, wake cutting, and the associated evolution of vortical entities are discussed.

Yang, R.-J.↗

Cost-effective accurate coarse-grid method for highly convective multidimensional unsteady flows

A fundamentally multidimensional convection scheme is described based on vector transient interpolation modeling rewritten in conservative control-volume form. Vector third-order upwinding is used as the basis of the algorithm; this automatically introduces important cross-difference terms that are absent from schemes using component-wise one-dimensional formulas. Third-order phase accuracy is good; this is important for coarse-grid large-eddy or full simulation. Potential overshoots or undershoots are avoided by using a recently developed universal limiter. Higher order accuracy is obtained locally, where needed, by the cost-effective strategy of adaptive stencil expansion in a direction normal to each control-volume face; this is controlled by monitoring the absolute normal gradient and curvature across the face. Higher (than third) order cross-terms do not appear to be needed. Since the wider stencil is used only in isolated narrow regions (near discontinuities), extremely high (in this case, seventh) order accuracy can be achieved for little more than the cost of a globally third-order scheme.

Leonard, B. P.↗

Numerical flux formulas for the Euler and Navier-Stokes equations. 2: Progress in flux-vector splitting

The accuracy of various numerical flux functions for the inviscid fluxes when used for Navier-Stokes computations is studied. The flux functions are benchmarked for solutions of the viscous, hypersonic flow past a 10 degree cone at zero angle of attack using first order, upwind spatial differencing. The Harten-Lax/Roe flux is found to give a good boundary layer representation, although its robustness is an issue. Some hybrid flux formulas, where the concepts of flux-vector and flux-difference splitting are combined, are shown to give unsatisfactory pressure distributions; there is still room for improvement. Investigations of low diffusion, pure flux-vector splittings indicate that a pure flux-vector splitting can be developed that eliminates spurious diffusion across the boundary layer. The resulting first-order scheme is marginally stable and not monotone.

Coirier, William J.↗

A new flux splitting scheme

A new flux splitting scheme is proposed. The scheme is remarkably simple and yet its accuracy rivals and in some cases surpasses that of Roe's solver in the Euler and Navier-Stokes solutions performed in this study. The scheme is robust and converges as fast as the Roe splitting. An approximately defined cell-face advection Mach number is proposed using values from the two straddling cells via associated characteristic speeds. This interface Mach number is then used to determine the upwind extrapolation for the convective quantities. Accordingly, the name of the scheme is coined as Advection Upstream Splitting Method (AUSM). A new pressure splitting is introduced which is shown to behave successfully, yielding much smoother results than other existing pressure splittings. Of particular interest is the supersonic blunt body problem in which the Roe scheme gives anomalous solutions. The AUSM produces correct solutions without difficulty for a wide range of flow conditions as well as grids.

Liou, Meng-Sing↗

Computation of compressible quasi-axisymmetric slender vortex flow and breakdown

The unsteady, compressible Navier-Stokes equations are used to compute and analyze compressible quasi-axisymmetric isolated vortices. The Navier-Stokes equations are solved using an implicit, upwind, flux difference splitting finite volume scheme. The developed three dimensional solver was verified by comparing its solution profiles with those of a slender, quasi-axisymmetric vortex solver for a subsonic, quasi-axisymmetric vortex in an unbounded domain. The Navier-Stokes solver is then used to solve for a supersonic, quasi-axisymmetric vortex flow in a configured circular duct. Steady and unsteady vortex-shock interactions and breakdown were captured. The problem was also calculated using the Euler solver of the same code; the results were compared with those of the Navier-Stokes solver. The effect of the initial swirl was investigated.

Kandil, Osama A.↗

Supersonic quasi-axisymmetric vortex breakdown

An extensive computational study of supersonic quasi-axisymmetric vortex breakdown in a configured circular duct is presented. The unsteady, compressible, full Navier-Stokes (NS) equations are used. The NS equations are solved for the quasi-axisymmetric flows using an implicit, upwind, flux difference splitting, finite volume scheme. The quasi-axisymmetric solutions are time accurate and are obtained by forcing the components of the flowfield vector to be equal on two axial planes, which are in close proximity of each other. The effect of Reynolds number, for laminar flows, on the evolution and persistence of vortex breakdown, is studied. Finally, the effect of swirl ration at the duct inlet is investigated.

Kandil, Osama A.↗

Unstructured-grid methods development: Lessons le arned

The development is summarized of unstructured grid methods for the solution of the equations of fluid flow and some of the lessons learned are shared. The 3-D Euler equations are solved, including spatial discretizations, temporal discretizations, and boundary conditions. An example calculation with an upwind implicit method using a CFL (Courant Friedricks Lewy) number of infinity is presented for the Boeing 747 aircraft. The results obtained in less than one hour of CPU time on a Cray-2 computer, thus demonstrating the speed and robustness of the present capability.

Batina, John T.↗

A patched-grid algorithm for complex aircraft configurations

A patched-grid algorithm for the analysis of complex configurations with an implicit, upwind-biased Navier-Stokes solver is presented. Through the use of a generalized coordinate transformation at the zonal interface between two or more blocks, the algorithm can be applied to highly stretched viscous grids and to arbitrarily-shaped patch boundaries. Applications are made to the SR71 reconnaissance aircraft in a high-altitude environment at a supersonic speed and to the F/A-18 forebody-strake configuration at subsonic, high-alpha conditions, in support of the NASA High-Alpha Research Program.

Walters, Robert W.↗

Scatterometer azimuthal response and wind wave directionality

Azimuthal response of a scatterometer to radiation scattered by the sea surface was studied in a wind-wave tank. The variation of the normalized radar cross section with the azimuth angle is fitted by a three-term series. Results show that the upwind-downwind asymmetry decreases as the wind speed increases. The crosswind modulation depends on the wind velocity. The results show that the evolution of the long-wind-crosswind ratio evolves with wind speed in a manner similar to the evolution of the isotropy of short capillary-gravity waves. The maximum of the isotropy of the short wind waves is obtained for wind velocities close to 4 m/s. For the same value of the velocity, the variations of radar response between long-wind and crosswind directions is minimum. For lower or higher values of wind velocities the directional accuracy of the radar increases, since the wind-wave field tends to align in the wind direction.

Giovanangeli, J. P.↗

Scatterometer directional response during rain

Rainfall modification of directional scatterometer response from the sea surface was simulated in wind-wave tank experiments. Data show that for the range of conditions in laboratory experiments, rain enhances radar cross section for all azimuthal angles relative to wind direction. This result broadens previous measurements, which showed that scatterometer response increases with increasing rainfall for radars pointing upwind. But more to the point, the data also show that the directional dynamic-range of scatterometry diminishes rapidly as rainfall rate increases. Thus, while it may be possible to determine wind speed and direction during rain, it will require adequate system sensitivity.

Bliven, L.↗

Experimental investigation of the dependence of radar backscattering on wind speed, wind stress and wave height

During summer 1988, radar measurements were performed in conjunction with detailed environmental observations on Lake Washington at the University of Washington Sand Point field station. Radar data were collected at 5.3 and 10 GHz for incidence angles between 30 and 60 deg with VV-polarization. The environmental measurements included wind speed and direction, large-wave heights, the high-frequency portion of the wave spectrum, humidity, and air and water temperatures. The small-scale wave spectrum was measured using a resistance wire gauge. The results show that backscatter increased with wind speed as expected. However, little difference was observed in the scattering coefficient for upwind and crosswind directions. The results also indicated an increase in the amplitude of small waves with friction velocity.

Gogineni, S. P.↗

Navier-Stokes solutions about the F/A-18 forebody-leading-edge extension configuration

Three-dimensional viscous flow computations are presented for the F/A-18 forebody including the leading-edge extension geometry. Solutions are obtained from an algorithm for the compressible Navier-Stokes equations that incorporates an upwind-biased, flux-difference-splitting approach along with longitudinally patched grids. Results are presented for both laminar and fully turbulent flow assumptions and include correlations with wind-tunnel as well as flight-test results. A good quantitative agreement for the forebody surface pressure distribution is achieved between the turbulent computations and wind-tunnel measurements at M ∞ = 0.6. The computed turbulent surface flow patterns on the forebody qualitatively agree well with in-flight surface flow patterns obtained on an F/A-18 aircraft at M ∞ = 0.34

Leading edge extension↗

Application of a patched-grid algorithm to the F/A-18 forebody-leading-edge extension configuration

A patched-grid algorithm for the analysis of complex configurations with an implicit, upwind-biased Navier-Stokes solver is presented. Results from both a spatial-flux and a time-flux conservation approach to patching across zonal boundaries are presented. A generalized coordinate transformation with a biquadratic geometric element is used at the zonal interface in order to treat highly stretched viscous grids and arbitrarily shaped zonal boundaries. Applications are made to the F/A-18 forebody-Ieading-edge extension configuration at subsonic, high-alpha conditions. Computed surface flow patterns compare well with ground-based and flight-test results; the large effect of Reynolds number on the forebody flowfield is shown.

Leading edge extension↗

Time-marching methods for three-dimensional steady and unsteady viscous imcompressible flows

An implicit algorithm for the solution of three-dimensional, steady and unsteady, viscous, incompressible flows is presented. The algorithm is based on an upwind-relaxation finite-difference method. Steady-state solutions are carried out using a time-marching solution technique in combination with a local time-stepping strategy. To obtain time-accurate solutions, a subiterative procedure is employed at each physical time step using a global time step to ensure the divergence-free condition. Steady-state flows in several straight ducts and in a square duct with a 90-degree bend are computed and compared with analytical and experimental results. The classical problem of starting flow in a circular pipe is chosen to verify the time accuracy of the present scheme. Finally, the three-dimensional bubble-type vortex breakdown of a slender cylindrical vortex in an unbounded flow is investigated.

Hsu, C.-H.↗

A grid-independent approximate Riemann solver with applications to the Euler and Navier-Stokes equations

A new two-dimensional approximate Riemann solver has been developed that obtains fluxes on grid faces via wave decomposition. By utilizing information propagation in the velocity-difference directions rather than in the grid-normal directions, this flux function more appropriately interprets and hence more sharply resolves shock and shear waves when they lie oblique to the grid. The model uses five waves to describe the difference in states at a grid face. Two acoustic waves, one shear wave, and one entropy wave propagate in the direction defined by the local velocity difference vector, while the fifth wave is a shear wave that propagates at a right angle to the other four. Test cases presented include a shock reflecting off a wall, a pure shear wave, supersonic flow over an airfoil, and viscous separated airfoil flow. Results using the new model give significantly sharper shock and shear contours than a grid-aligned solver. Navier-Stokes computations over an aifoil show reduced pressure distortions in the separated region as a result of the grid-independent upwinding.

Rumsey, Christopher L.↗