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 991 records · Page 55

A k-Omega Turbulence Model for Quasi-Three-Dimensional Turbomachinery Flows

A two-equation k-omega turbulence model has been developed and applied to a quasi-three-dimensional viscous analysis code for blade-to-blade flows in turbomachinery. the code includes the effects of rotation, radius change, and variable stream sheet thickness. The flow equations are given and the explicit runge-Kutta solution scheme is described. the k-omega model equations are also given and the upwind implicit approximate-factorization solution scheme is described. Three cases were calculated: transitional flow over a flat plate, a transonic compressor rotor, and transonic turbine vane with heat transfer. Results were compared to theory, experimental data, and to results using the Baldwin-Lomax turbulence model. The two models compared reasonably well with the data and surprisingly well with each other. Although the k-omega model behaves well numerically and simulates effects of transition, freestream turbulence, and wall roughness, it was not decisively better than the Baldwin-Lomax model for the cases considered here.

Chima, Rodrick V.↗

Flow simulations about steady-complex and unsteady moving configurations using structured-overlapped and unstructured grids

The limiting factor in simulating flows past realistic configurations of interest has been the discretization of the physical domain on which the governing equations of fluid flow may be solved. In an attempt to circumvent this problem, many Computational Fluid Dynamic (CFD) methodologies that are based on different grid generation and domain decomposition techniques have been developed. However, due to the costs involved and expertise required, very few comparative studies between these methods have been performed. In the present work, the two CFD methodologies which show the most promise for treating complex three-dimensional configurations as well as unsteady moving boundary problems are evaluated. These are namely the structured-overlapped and the unstructured grid schemes. Both methods use a cell centered, finite volume, upwind approach. The structured-overlapped algorithm uses an approximately factored, alternating direction implicit scheme to perform the time integration, whereas, the unstructured algorithm uses an explicit Runge-Kutta method. To examine the accuracy, efficiency, and limitations of each scheme, they are applied to the same steady complex multicomponent configurations and unsteady moving boundary problems. The steady complex cases consist of computing the subsonic flow about a two-dimensional high-lift multielement airfoil and the transonic flow about a three-dimensional wing/pylon/finned store assembly. The unsteady moving boundary problems are a forced pitching oscillation of an airfoil in a transonic freestream and a two-dimensional, subsonic airfoil/store separation sequence. Accuracy was accessed through the comparison of computed and experimentally measured pressure coefficient data on several of the wing/pylon/finned store assembly's components and at numerous angles-of-attack for the pitching airfoil. From this study, it was found that both the structured-overlapped and the unstructured grid schemes yielded flow solutions of comparable accuracy for these simulations. This study also indicated that, overall, the structured-overlapped scheme was slightly more CPU efficient than the unstructured approach.

Newman, James C., III↗

Three-dimensional unstructured method for flows past bodies in 6-DOF relative motion

A three dimensional, unstructured-mesh methodology was developed to simulate unsteady flows past bodies in relative motion, where the trajectory was determined from the instantaneous aerodynamics. The method coupled the equations of fluid flow and those of rigid-body dynamics, and captured the time-dependent interference between stationary and moving boundaries. The unsteady, compressible Euler equations were solved on dynamic, unstructured meshes by an explicit, finite-volume, upwind method. The grid adaptation was performed within a window placed around the moving body. The Euler equations of dynamics were solved by a Runge-Kutta integration scheme. The flow solver and the adaptation scheme were validated by simulating the transonic, unsteady flow around a wing undergoing a forced, periodic pitching motion, then comparing the results with the experimental data. To validate the trajectory code, the six-degrees-of-freedom (DOF) motion of a store separating from a wing was computed using the experimentally determined force and moment fields, then comparing with an independently generated trajectory. Finally, the overall methodology was demonstrated by simulating the unsteady flowfield and the trajectory of a store dropped from a wing. The methodology, its computational cost notwithstanding, has proven to be accurate, automated, easy for dynamic gridding, and relatively efficient for the required man-hours.

Singh, K. P.↗

Design of 3-D Nacelle near Flat-Plate Wing Using Multiblock Sensitivity Analysis (ADOS)

One of the major design tasks involved in reducing aircraft drag is the integration of the engine nacelles and airframe. With this impetus, nacelle shapes with and without the presence of a flat-plate wing nearby were optimized. This also served as a demonstration of the 3-D version of the recently developed aerodynamic design optimization methodology using sensitivity analysis, ADOS. The required flow analyses were obtained by solving the three-dimensional, compressible, thin-layer Navier-Stokes equations using an implicit, upwind-biased, finite volume scheme. The sensitivity analyses were performed using the preconditioned version of the SADD scheme (sensitivity analysis on domain decomposition). In addition to demonstrating the present method's capability for automatic optimization, the results offered some insight into two important issues related to optimizing the shapes of multicomponent configurations in close proximity. First, inclusion of the mutual interference between the components resulted in a different shape as opposed to shaping an isolated component. Secondly, exclusion of the viscous effects compromised not only the flow physics but also the optimized shapes even for isolated components.

Eleshaky, Mohamed E.↗

Shape optimization of single- and two-element airfoils on multiblock grids

A multiblock, discrete sensitivity analysis method is used to couple a direct optimization method and a flow analysis method. The domain is divided into smaller subdomains for which the sensitivities are obtained separately. Then, an effective sensitivity equation is solved to complete the coupling of all the sensitivity information. The flow analysis is based on the thin-layer Navier-Stokes equations solved by an implicit, upwind-biased, finite-volume method. The method of feasible directions is used for the present gradient-based optimization approach. First, a transonic airfoil is optimized to investigate the behavior of the method in highly nonlinear flows as well as the effect of different blocking strategies on the procedure. A supercritical airfoil is produced from an initially symmetric airfoil with multiblocking affecting the path but not the final shape. Secondly, a two-element airfoil is shape optimized in subsonic flow to demonstrate the present method's capability of shaping aerodynamically interfering elements simultaneously. For a very low and a very high Reynolds number cases, the shape of the main airfoil and the flap are optimized to yield improved lift-to-drag ratios.

Lacasse, James M.↗

Methodology for simulation of unsteady flows to determine the time-dependent interference between stationary and moving boundaries

A new methodology is developed to simulate unsteady flows about prescribed and aerodynamically determined moving boundary problems. The method couples the fluid dynamics and rigid-body dynamics equations to capture the time-dependent interference between stationary and moving boundaries. The unsteady, compressible, inviscid (Euler) equations are solved on dynamic, unstructured grids by an explicit, finite-volume, upwind method. For efficiency, the grid adaptation is performed within a window around the moving object. The Eulerian equations of the rigid-body dynamics are solved by a Runge-Kutta method in a non-inertial frame of reference. The two-dimensional flow solver is validated by computing the flow past a sinusoidally-pitching airfoil and comparing these results with the experimental data. The overall methodology is used for two two-dimensional examples: the flow past an airfoil which is performing a three-degrees-of-freedom motion in a transonic freestream, and the free-fall of a store after separation from a wing-section. Then the unstructured mesh methodology is extended to three-dimensions to simulate unsteady flow past bodies in relative motion, where the trajectory is determined from the instantaneous aerodynamics. The flow solver and the adaptation scheme in three dimensions are validated by simulating the transonic, unsteady flow around a wing undergoing a forced, periodic, pitching motion, and comparing the results with the experimental data. To validate the trajectory code, the six-degrees-of-freedom motion of a store separating from a wing was computed using the experimentally determined force and moment fields, then comparing with an independently generated trajectory. Finally, the overall methodology was demonstrated by simulating the unsteady flowfield and the trajectory of a store dropped from a wing. The methodology, its computational cost notwithstanding, has proven to be accurate, automated, easy for dynamic gridding, and relatively efficient for the required man-hours.

Singh, Kamakhya Prasad↗

Adaptively-refined overlapping grids for the numerical solution of systems of hyperbolic conservation laws

Adaptive mesh refinement (AMR) in conjunction with higher-order upwind finite-difference methods have been used effectively on a variety of problems in two and three dimensions. In this paper we introduce an approach for resolving problems that involve complex geometries in which resolution of boundary geometry is important. The complex geometry is represented by using the method of overlapping grids, while local resolution is obtained by refining each component grid with the AMR algorithm, appropriately generalized for this situation. The CMPGRD algorithm introduced by Chesshire and Henshaw is used to automatically generate the overlapping grid structure for the underlying mesh.

Brislawn, Kristi D.↗

Adaptive unstructured triangular mesh generation and flow solvers for the Navier-Stokes equations at high Reynolds number

A method for generating high quality unstructured triangular grids for high Reynolds number Navier-Stokes calculations about complex geometries is described. Careful attention is paid in the mesh generation process to resolving efficiently the disparate length scales which arise in these flows. First the surface mesh is constructed in a way which ensures that the geometry is faithfully represented. The volume mesh generation then proceeds in two phases thus allowing the viscous and inviscid regions of the flow to be meshed optimally. A solution-adaptive remeshing procedure which allows the mesh to adapt itself to flow features is also described. The procedure for tracking wakes and refinement criteria appropriate for shock detection are described. Although at present it has only been implemented in two dimensions, the grid generation process has been designed with the extension to three dimensions in mind. An implicit, higher-order, upwind method is also presented for computing compressible turbulent flows on these meshes. Two recently developed one-equation turbulence models have been implemented to simulate the effects of the fluid turbulence. Results for flow about a RAE 2822 airfoil and a Douglas three-element airfoil are presented which clearly show the improved resolution obtainable.

Ashford, Gregory A.↗

Unstructured viscous flow solution using adaptive hybrid grids

A three dimensional finite volume scheme based on hybrid grids containing both tetrahedral and hexahedral cells is presented. The application to hybrid grids offers the possibility to combine the flexibility of tetrahedral meshes with the accuracy of hexahedral grids. An algorithm to compute a dual mesh for the entire computational domain was developed. The dual mesh technique guarantees conservation in the whole flow field even at interfaces between hexahedral and tetrahedral domains and enables the employment of an accurate upwind flow solver. The hybrid mesh can be adapted to the solution by dividing cells in areas of insufficient resolution. The method is tested on different viscous and inviscid cases for hypersonic, transonic and subsonic flows.

Galle, Martin↗

Shock-Induced Separated Structures in Symmetric Corner Flows

Three-dimensional supersonic viscous laminar flows over symmetric corners are considered in this paper. The characteristic features of such configurations are discussed and an historical survey on the past research work is presented. A new contribution based on a numerical technique that solves the parabolized form of the Navier-Stokes equations is presented. Such a method makes it possible to obtain very detailed descriptions of the flowfield with relatively modest CPU time and memory storage requirements. The numerical approach is based on a space-marching technique, uses a finite volume discretization and an upwind flux-difference splitting scheme (developed for the steady flow equations) for the evaluation of the inviscid fluxes. Second order accuracy is reached following the guidelines of the ENO schemes. Different free-stream conditions and geometrical configurations are considered. Primary and secondary streamwise vortical structures embedded in the boundary layer and originated by the interaction of the latter with shock waves are detected and studied. Computed results are compared with experimental data taken from literature.

DAmbrosio, Domenic↗

Polarimetric measurements of sea surface brightness temperatures using an aircraft K-band radiometer

This paper presents the first experimental evidence that the polarimetric brightness temperatures of sea surfaces are sensitive to ocean wind direction in the incidence angle range of 30 to 50 degrees. Our experimental data were collected by a K-band (19.35 GHz) polarimetric wind radiometer (WINDRAD) mounted on the NASA DC-8 aircraft. A set of aircraft radiometer flights was successfully completed in November 1993. We performed circle flights over National Data Buoy Center (NDBC) moored buoys deployed off the northern California coast, which provided ocean wind measurements. The first WINDRAD flight was made on November 4, 1993. There was clear weather with a wind speed of 12 m/s at 330 degrees around the Pt. Arena buoy. We circled the buoy at three incidence angles, and all data when plotted as functions of azimuth angles show clear modulations of several Kelvin. At 40 degrees incidence angle, there is a 5 Kelvin peak-to-peak signal in the second Stokes parameter Q and the third Stokes parameter U. The Q data maximum is in the upwind direction and U has a 45 degrees phase shift in azimuth as predicted by theory. There is also an up/downwind asymmetry of 2 Kelvin in the Q data, and 1 Kelvin in the U data. At 50 degrees incidence angle, the collected data show very similar wind direction signatures to the SSM/I model function. Additional flights were made on other days under cloudy conditions. Data taken at a wind speed of 8 m/s show that at 40 degrees incidence Q and U have a smaller azimuthal modulation of 3 Kelvin, probably due to the lower wind speed. Additionally, the simultaneously recorded video images of sea surfaces suggested that Q and U data were less sensitive to unpolarized geophysical variations, such as clouds and whitecaps, while the T(v) and T(h) increased by a few Kelvin when the radiometer beam crossed over clouds, or there was a sudden increase of whitecaps in the radiometer footprint. The results of our aircraft flights indicate that passive polarimetric radiometry has a strong potential for global ocean wind speed and direction measurements from space.

Yueh, Simon H.↗

Assessment of an Unstructured-Grid Method for Predicting 3-D Turbulent Viscous Flows

A method is presented for solving turbulent flow problems on three-dimensional unstructured grids. Spatial discretization is accomplished by a cell-centered finite-volume formulation using an accurate linear reconstruction scheme and upwind flux differencing. Time is advanced by an implicit backward-Euler time-stepping scheme. Flow turbulence effects are modeled by the Spalart-Allmaras one-equation model, which is coupled with a wall function to reduce the number of cells in the sublayer region of the boundary layer. A systematic assessment of the method is presented to devise guidelines for more strategic application of the technology to complex problems. The assessment includes the accuracy in predictions of skin-friction coefficient, law-of-the-wall behavior, and surface pressure for a flat-plate turbulent boundary layer, and for the ONERA M6 wing under a high Reynolds number, transonic, separated flow condition.

Frink, Neal T.↗

Coupled Rolling and Pitching Oscillation Effects on Transonic Shock-Induced Vortex-Breakdown Flow of a Delta Wing

Unsteady, transonic vortex dominated flow over a 65 deg. sharp edged, cropped-delta wing of zero thickness undergoing forced coupled pitching and rolling oscillations is investigated computationally. The wing mean angle of attack is 20 deg. and the free stream Mach number and Reynolds number are 0.85 and 3.23 x 10(exp 6), respectively. The initial condition of the flow is characterized by a transverse terminating shock and vortex breakdown of the leading edge vortex cores. The computational investigation uses the time-accurate solution of the laminar, unsteady, compressible, full Navier-Stokes equations with the implicit, upwind, Roe flux-difference splitting, finite volume scheme. The main focus is to analyze the effects of coupled motion on the wing response and vortex breakdown flow by varying oscillation frequency and phase angle while the maximum pitch and roll amplitude is kept constant at 4.0 deg. Four cases demonstrate the following: simultaneous motion at a frequency of 1(pi), motion with a 90 deg. phase lead in pitch, motion with a rolling frequency of twice the pitching frequency, and simultaneous motion at a frequency of 2(pi). Comparisons with single mode motion at these frequencies complete this study and illustrate the effects of coupling the oscillations.

Kandil, Osama A.↗

Effects of Coupled Rolling and Pitching Oscillations on Transonic Shock-Induced Vortex-Breakdown Flow of a Delta Wing

Unsteady, transonic vortex-breakdown flow over a 65 deg. sharp edged, cropped-delta wing of zero thickness undergoing forced coupled pitching and rolling oscillations is investigated computationally. The initial condition of the flow is characterized by a transverse terminating shock which induces of the leading edge vortex cores to breakdown. The computational investigation uses the time-accurate solution of the laminar, unsteady, compressible, full Navier-Stokes equations with the implicit, upwind, Roe flux-difference splitting, finite-volume scheme. The main focus is to analyze the effects of coupled motion on the wing response and vortex-breakdown flow by varying oscillation frequency and phase angle while keeping the maximum pitch and roll amplitude equal.

Kandil, Osama A.↗

Forced Rolling Oscillation of a 65 deg-Delta Wing in Transonic Vortex-Breakdown Flow

Unsteady, transonic, vortex dominated flow over a 65 deg. sharp-edged, cropped-delta wing of zero thickness undergoing forced rolling oscillations is investigated computationally. The wing angle of attack is 20 deg. and the free stream Mach number and Reynolds number are 0.85 and 3.23 x 10(exp 6), respectively. The initial condition of the flow is characterized by a transverse terminating shock which induces vortex breakdown of the leading edge vortex cores. The computational investigation uses the time accurate solution of the laminar, unsteady, compressible, full Navier-Stokes equations with the implicit, upwind, Roe flux difference splitting, finite-volume scheme. While the maximum roll amplitude is kept constant at 4.0 deg., both Reynolds number and roll frequency are varied covering three cases of forced sinusoidal rolling. First, the Reynolds number is held at 3.23 x 10(exp 6) and the wing is forced to oscillate in roll around the axis of geometric symmetry at a reduced frequency of 2(pi). Second, the Reynolds number is reduced to 0.5 x 10(exp 6) to observe the effects of added viscosity on the vortex breakdown. Third, with the Reynolds number held at 0.5 x 10(exp 6), the roll frequency is reduced to 1(pi) to complete the study.

Menzies, Margaret A.↗

Pick-up ion pressure gradients modulating the solar wind dynamics

Neutral interstellar atoms penetrate deeply into the inner heliosphere before they become ionized by various processes. As ions they are picked-up by the frozen-in magnetic fields and are convected outwards with the solar wind plasma. Thereby the primary plasma flow is mass, momentum, and energy-loaded. The dynamics of the distant multi-constituent solar wind is, however, not solely determined by these loading processes, but is also affected by the wave-mediated pick-up ion pressure gradients derivable from the pick-up ion distribution function. The action of the radial components of these pressures essentially counter balances the decelerating effect of the solar wind momentum loading, diminishing strongly the deceleration of the distant solar wind. Furthermore the latitudinal components of the pick-up ion pressures induce latitudinal forces acting on the multiconstituent solar plasma outflow and inducing nonradial bulk flow components. The enforced nonradial outflow geometry on the upwind hemisphere may partly be responsible for the magnetic flux deficit which was claimed since several years in the PIONEER-10 magnetic flux data.

Fahr, Hans J.↗

Investigation of supersonic jet plumes using an improved two-equation turbulence model

Supersonic jet plumes were studied using a two-equation turbulence model employing corrections for compressible dissipation and pressure-dilatation. A space-marching procedure based on an upwind numerical scheme was used to solve the governing equations and turbulence transport equations. The computed results indicate that two-equation models employing corrections for compressible dissipation and pressure-dilatation yield improved agreement with the experimental data. In addition, the numerical study demonstrates that the computed results are sensitive to the effect of grid refinement and insensitive to the type of velocity profiles used at the inflow boundary for the cases considered in the present study.

Lakshmanan, B.↗

NASA flow fields analysis

The objectives of the present research are to improve design capabilities for low thrust rocket engines through understanding the detailed mixing and combustion processes in a representative combustor. Of particular interest is a small gaseous hydrogen-oxygen thruster which is considered as a coordinated part of an on-going experimental program at NASA LERC. Detailed computational modeling involves the solution of both the two- and three-dimensional Navier Stokes equations, coupled with chemical reactions and the species diffusion equations. Computations of interest include both steady state and time-accurate flowfields and are obtained by means of LU approximate factorization in time and flux split upwinding differencing in space. The emphasis in the research is focused on using numerical analysis to understand detailed combustor flowfields, including the shear layer dynamics created between fuel film cooling and the core gas in the vicinity on the nearby combustor wall; the integrity and effectiveness of the coolant film; and three-dimensional fuel and oxidizer jet injection/mixing/combustion characteristics in the primary combustor along with their joint impacts on global engine performance.

Merkle, Charles L.↗