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 865 records · Page 48
All sourcesNASA NTRSDOE OSTI
NASA NTRS1992

Navier-Stokes simulations of unsteady transonic flow phenomena

Numerical simulations of two classes of unsteady flows are obtained via the Navier-Stokes equations: a blast-wave/target interaction problem class and a transonic cavity flow problem class. The method developed for the viscous blast-wave/target interaction problem assumes a laminar, perfect gas implemented in a structured finite-volume framework. The approximately factored implicit scheme uses Newton subiterations to obtain the spatially and temporally second-order accurate time history of the blast-waves with stationary targets. The inviscid flux is evaluated using either of two upwind techniques, while the full viscous terms are computed by central differencing. Comparisons of unsteady numerical, analytical, and experimental results are made in two- and three-dimensions for Couette flows, a starting shock-tunnel, and a shock-tube blockage study. The results show accurate wave speed resolution and nonoscillatory discontinuity capturing of the predominantly inviscid flows. Viscous effects were increasingly significant at large post-interaction times. While the blast-wave/target interaction problem benefits from high-resolution methods applied to the Euler terms, the transonic cavity flow problem requires the use of an efficient scheme implemented in a geometrically flexible overset mesh environment. Hence, the Reynolds averaged Navier-Stokes equations implemented in a diagonal form are applied to the cavity flow class of problems. Comparisons between numerical and experimental results are made in two-dimensions for free shear layers and both rectangular and quieted cavities, and in three-dimensions for Stratospheric Observatory For Infrared Astronomy (SOFIA) geometries. The acoustic behavior of the rectangular and three-dimensional cavity flows compare well with experiment in terms of frequency, magnitude, and quieting trends. However, there is a more rapid decrease in computed acoustic energy with frequency than observed experimentally owing to numerical dissipation. In addition, optical phase distortion due to the time-varying density field is modelled using geometrical constructs. The computed optical distortion trends compare with the experimentally inferred result, but underpredicts the fluctuating phase difference magnitude.

Atwood, C. A.↗
NASA NTRS1992

Three-dimensional computational study of asymmetric flows using Navier-Stokes equations

The unsteady, compressible, thin-layer Navier-Stokes equations are used to obtain three-dimensional, asymmetric, vortex-flow solutions around cones and cone-cylinder configurations. The equations are solved using an implicit, upwind, flux-difference splitting, finite-volume scheme. The computational applications cover asymmetric flows around a 5 semi-apex angle cone of unit length at various Reynolds number. Next, a cylindrical afterbody of various length is added to the conical forebody to study the effect of the length of cylindrical afterbody on the flow asymmetry. All the asymmetric flow solutions are obtained by using a short-duration side-slip disturbance.

Cheung, Y. K.↗
NASA NTRS1992

Passive control of supersonic asymmetric vortical flows around cones

The unsteady, compressible, thin-layer Navier-Stokes equations are used to numerically study the passive control of steady and unsteady supersonic asymmetric flows around circular and noncircular cones. The main computational scheme of the present study is an implicit upwind, flux-difference splitting, finite-volume scheme. Passive control of flow asymmetry is studied by using a vertical fin in the leeward plane of geometric symmetry and side strakes with and without thickness at different orientations. The study focuses on circular-section cones since they are the most likely section-shapes for strong flow asymmetry. Side-strake passive control is shown to be more efficient and practical than vertical-fin passive control.

Liu, C. H.↗
NASA NTRS1993

Preconditioned conjugate-gradient methods for low-speed flow calculations

An investigation is conducted into the viability of using a generalized Conjugate Gradient-like method as an iterative solver to obtain steady-state solutions of very low-speed fluid flow problems. Low-speed flow at Mach 0.1 over a backward-facing step is chosen as a representative test problem. The unsteady form of the two dimensional, compressible Navier-Stokes equations is integrated in time using discrete time-steps. The Navier-Stokes equations are cast in an implicit, upwind finite-volume, flux split formulation. The new iterative solver is used to solve a linear system of equations at each step of the time-integration. Preconditioning techniques are used with the new solver to enhance the stability and convergence rate of the solver and are found to be critical to the overall success of the solver. A study of various preconditioners reveals that a preconditioner based on the Lower-Upper Successive Symmetric Over-Relaxation iterative scheme is more efficient than a preconditioner based on Incomplete L-U factorizations of the iteration matrix. The performance of the new preconditioned solver is compared with a conventional Line Gauss-Seidel Relaxation (LGSR) solver. Overall speed-up factors of 28 (in terms of global time-steps required to converge to a steady-state solution) and 20 (in terms of total CPU time on one processor of a CRAY-YMP) are found in favor of the new preconditioned solver, when compared with the LGSR solver.

Ajmani, Kumud↗
NASA NTRS1993

Basic studies of baroclinic flows

A fully nonlinear 3-dimensional numerical model (GEOSIM), previously developed and validated for several cases of geophysical fluid flow, has been used to investigate the dynamical behavior of laboratory experiments of fluid flows similar to those of the Earth's atmosphere. The phenomena investigated are amplitude vacillation, and the response of the fluid system to uneven heating and cooling. The previous year's work included hysteresis in the transition between axisymmetric and wave flow. Investigation is also continuing of the flows in the Geophysical Fluid Flow Cell (GFFC), a low-gravity Spacelab experiment. Much of the effort in the past year has been spent in validation of the model under a wide range of external parameters including nonlinear flow regimes. With the implementation of a 3-dimensional upwind differencing scheme, higher spectral resolution, and a shorter time step, the model has been found capable of predicting the majority of flow regimes observed in one complete series of baroclinic annulus experiments of Pfeffer and co-workers. Detailed analysis of amplitude vacillation has revealed that the phase splitting described in the laboratory experiments occurs in some but not all cases. Through the use of animation of the models output, a vivid 3-dimensional view of the phase splitting was shown to the audience of the Southeastern Geophysical Fluid Dynamics Conference in March of this year. A study on interannual variability was made using GEOSIM with periodic variations in the thermal forcing. Thus far, the model has not predicted a chaotic behavior as observed in the experiments, although there is a sensitivity in the wavenumber selection to the initial conditions. Work on this subject, and on annulus experiments with non-axisymmetric thermal heating, will continue. The comparison of GEOSIM's predictions will result from the Spacelab 3 GFFC experiments continued over the past year, on a 'back-burner' basis. At this point, the study (in the form of a draft of a journal article) is nearly completed. The results from GEOSIM compared very well with the experiments, and the use of the model allows the demonstration of flow mechanics that were not possible with the experimental data. For example, animation of the model output shows that the forking of the spiral bands is a transient phenomenon, due to the differential east-west propagation of convection bands from different latitudes.

Miller, Timothy L.↗
NASA NTRS1993

Methodology for sensitivity analysis, approximate analysis, and design optimization in CFD for multidisciplinary applications

In this study involving advanced fluid flow codes, an incremental iterative formulation (also known as the delta or correction form) together with the well-known spatially-split approximate factorization algorithm, is presented for solving the very large sparse systems of linear equations which are associated with aerodynamic sensitivity analysis. For smaller 2D problems, a direct method can be applied to solve these linear equations in either the standard or the incremental form, in which case the two are equivalent. Iterative methods are needed for larger 2D and future 3D applications, however, because direct methods require much more computer memory than is currently available. Iterative methods for solving these equations in the standard form are generally unsatisfactory due to an ill-conditioning of the coefficient matrix; this problem can be overcome when these equations are cast in the incremental form. These and other benefits are discussed. The methodology is successfully implemented and tested in 2D using an upwind, cell-centered, finite volume formulation applied to the thin-layer Navier-Stokes equations. Results are presented for two sample airfoil problems: (1) subsonic low Reynolds number laminar flow; and (2) transonic high Reynolds number turbulent flow.

Taylor, Arthur C., III↗
NASA NTRS1993

Computation of incompressible viscous flows through turbopump components

Flow through pump components, such as an inducer and an impeller, is efficiently simulated by solving the incompressible Navier-Stokes equations. The solution method is based on the pseudocompressibility approach and uses an implicit-upwind differencing scheme together with the Gauss-Seidel line relaxation method. the equations are solved in steadily rotating reference frames and the centrifugal force and the Coriolis force are added to the equation of motion. Current computations use a one-equation Baldwin-Barth turbulence model which is derived from a simplified form of the standard k-epsilon model equations. The resulting computer code is applied to the flow analysis inside a generic rocket engine pump inducer, a fuel pump impeller, and SSME high pressure fuel turbopump impeller. Numerical results of inducer flow are compared with experimental measurements. In the fuel pump impeller, the effect of downstream boundary conditions is investigated. Flow analyses at 80 percent, 100 percent, and 120 percent of design conditions are presented.

Kiris, Cetin↗
NASA NTRS1993

Wing flutter boundary prediction using an unsteady Euler aerodynamic method

Modifications to an existing three-dimensional, implicit, upwind Euler/Navier-Stokes code (CFL3D Version 2.1) for the aeroelastic analysis of wings are described. These modifications, which were previously added to CFL3D Version 1.0, include the incorporation of a deforming mesh algorithm and the addition of the structural equations of motion for their simultaneous time-integration with the government flow equations. The paper gives a brief description of these modifications and presents unsteady calculations which check the modifications to the code. Euler flutter results for an isolated 45 degree swept-back wing are compared with experimental data for seven freestream Mach numbers which define the flutter boundary over a range of Mach number from 0.499 to 1.14. These comparisons show good agreement in flutter characteristics for freestream Mach numbers below unity. For freestream Mach numbers above unity, the computed aeroelastic results predict a premature rise in the flutter boundary as compared with the experimental boundary. Steady and unsteady contours of surface Mach number and pressure are included to illustrate the basic flow characteristics of the time-marching flutter calculations and to aid in identifying possible causes for the premature rise in the computational flutter boundary.

Lee-Rausch, Elizabeth M.↗
NASA NTRS1993

Spatial adaptation procedures on tetrahedral meshes for unsteady aerodynamic flow calculations

Spatial adaptation procedures for the accurate and efficient solution of steady and unsteady inviscid flow problems are described. The adaptation procedures were developed and implemented within a three-dimensional, unstructured-grid, upwind-type Euler code. These procedures involve mesh enrichment and mesh coarsening to either add points in high gradient regions of the flow or remove points where they are not needed, respectively, to produce solutions of high spatial accuracy at minimal computational cost. A detailed description of the enrichment and coarsening procedures are presented and comparisons with experimental data for an ONERA M6 wing and an exact solution for a shock-tube problem are presented to provide an assessment of the accuracy and efficiency of the capability. Steady and unsteady results, obtained using spatial adaptation procedures, are shown to be of high spatial accuracy, primarily in that discontinuities such as shock waves are captured very sharply.

Rausch, Russ D.↗
NASA NTRS1992

ERS-1 scatterometer calibration and validation activities at ECMWF. B: From radar backscatter characteristics to wind vector solutions

Calibration and validation activities for the ERS-1 scatterometer were carried out at ECMWF (European Center for Medium range Weather Forecast) complementary to the 'Haltenbanken' field campaign off the coast of Norway. At a Numerical Weather Prediction (NWP) center a wealth of verifying data is available both in time and space. This data is used to redefine the wind retrieval procedure given the instrumental characteristics. It was found that a maximum likelihood estimation procedure to obtain the coefficients of a reformulated sigma deg to wind relationship should use radar measurements in logarithmic rather than physical space, and use winds as the wind components rather than wind speed and direction. Doing this, a much more accurate transfer function than the one currently operated by ESA was derived. Sigma deg measurement space shows no signature of a separation in an upwind solution cone and a downwind solution cone. As such signature was anticipated in ESA's wind direction ambiguity removal algorithm, reconsideration of the procedure is necessary. Despite the fact that revisions have to be made in the process of wind retrieval; a grid potential is shown for scatterometry in meteorology and climatology.

Stoffelen, AD↗
NASA NTRS1993

The addition of algebraic turbulence modeling to program LAURA

The Langley Aerothermodynamic Upwind Relaxation Algorithm (LAURA) is modified to allow the calculation of turbulent flows. This is accomplished using the Cebeci-Smith and Baldwin-Lomax eddy-viscosity models in conjunction with the thin-layer Navier-Stokes options of the program. Turbulent calculations can be performed for both perfect-gas and equilibrium flows. However, a requirement of the models is that the flow be attached. It is seen that for slender bodies, adequate resolution of the boundary-layer gradients may require more cells in the normal direction than a laminar solution, even when grid stretching is employed. Results for axisymmetric and three-dimensional flows are presented. Comparison with experimental data and other numerical results reveal generally good agreement, except in the regions of detached flow.

Cheatwood, F. Mcneil↗
NASA NTRS1992

CSU FIRE 2 cirrus field experiment: Description of field deployment phase

The Colorado State University (CSU) surface observing systems are described. These systems were deployed at the Parsons, Kansas site during the FIRE 2 Cirrus Special Observing Period (SOP) from 13 Nov. - 7 Dec. 1991. The geographical coordinates of the site containing most of the CSU instrumentation are 37 deg. 18 min N. latitude and 96 deg. 30 min. W. longitude; site elevation was 269 meters. In addition, one surface meteorological and broadband flux observing site was maintained at the Tri City Airport which is approximately 18 miles due west of Parsons (37 deg. 20 min. N. latitude, 95 deg. 30 min. 30 sec. W. longitude). A map of the locations of the CSU deployment sites is presented. At the main Parsons site, the instrumentation was located directly adjacent to and north of a lake. Under most cirrus observing conditions, when the wing had a significant southernly component, the lake was upwind of the observing site. The measurements and observations collected during the experiment are listed. These measurements may be grouped into five categories: surface meteorology; infrared spectral and broadband measurements; solar spectral and broadband measurements; upper air measurements; and cloud measurements. A summary of observations collected at the Parsons site during the SOP are presented. The wind profiler, laser ceilometer, surface meteorology and surface broadband radiation instrumentation were operated on a continuous basis. All other systems were operated on an 'on demand' basis when cloud conditions merited the collection of data.

Cox, S.↗
NASA NTRS1992

Efficient simulation of incompressible viscous flow over multi-element airfoils

The incompressible, viscous, turbulent flow over single and multi-element airfoils is numerically simulated in an efficient manner by solving the incompressible Navier-Stokes equations. The computer code uses the method of pseudo-compressibility with an upwind-differencing scheme for the convective fluxes and an implicit line-relaxation solution algorithm. The motivation for this work includes interest in studying the high-lift take-off and landing configurations of various aircraft. In particular, accurate computation of lift and drag at various angles of attack, up to stall, is desired. Two different turbulence models are tested in computing the flow over an NACA 4412 airfoil; an accurate prediction of stall is obtained. The approach used for multi-element airfoils involves the use of multiple zones of structured grids fitted to each element. Two different approaches are compared: a patched system of grids, and an overlaid Chimera system of grids. Computational results are presented for two-element, three-element, and four-element airfoil configurations. Excellent agreement with experimental surface pressure coefficients is seen. The code converges in less than 200 iterations, requiring on the order of one minute of CPU time (on a CRAY YMP) per element in the airfoil configuration.

Rogers, Stuart E.↗
NASA NTRS1992

Investigation of forced unsteady separated flows using velocity-vorticity form of Navier-Stokes equations

The phenomenon of forced unsteady separation and eruption of boundary-layer vorticity is a highly-complex, high-Reynolds number flow phenomenon, which abruptly leads to the formation of a dynamic stall vortex as demonstrated earlier by the authors for a NACA 0015 airfoil undergoing constant rate pitch-up motion. This, as well as the results of other researchers, have convincingly demonstrated a complex vortical structure within the state of unsteady separation prior to the evolution of dynamic stall. This phenomenon of vortex eruption, although observed in studying dynamic stall phenomena, is also associated with transition from laminar to turbulence flow and its generic nature has been stressed by many researchers including the present investigators. An unsteady Navier-Stokes (NS) analysis is developed for arbitrarily maneuvering bodies using velocity-vorticity variables; this formulation is nearly form-invariant under a generalized non-inertial coordinate transformation. A fully-implicit uniformly second-order accurate method is used, with the nonlinear convective terms approximated using a biased third-order upwind differencing scheme to be able to simulate higher-Re flows. No explicit artificial dissipation is added. The numerical method is fully vectorized and currently achieves a computational index of 7 micro-seconds per time step per mesh point, using a single processor on a CRAY Y-MP. The simulation results show that the energetic free shear from the leading edge is responsible for the wall viscous layer to abruptly erupt near the center of the counterclockwise rotating eddy in the unsteady boundary layer. Primary, secondary, tertiary and quaternary vortices have been observed before the dynamic stall vortex evolves and gathers its maximum strength. This study will discuss the simulation results of Reynolds number up to Re = 45,000 and will also discuss the efforts of initial acceleration in a specific maneuver, on the evolution of the stall vortex.

Ghia, K. N.↗
NASA NTRS1993

Navier-Stokes analysis of three-dimensional flow and heat transfer inside turbine blade rows

A numerical method for solving the three-dimensional, Navier-Stokes equations for unsteady, viscous flow and heat transfer through multiple turbomachinery blade rows is presented. The method solves the fully three-dimensional Navier-Stokes equations with an implicit scheme which is based on a control volume approach. A two-equation turbulence model with a low Reynolds number modification is employed. A third-order accurate upwinding scheme is used to approximate convection terms while a second order accurate central difference scheme is used for the discretization of viscous terms. A second-order accurate scheme is employed for the temporal discretization. The numerical method is applied to study the unsteady flow and heat transfer field of the High Pressure Fuel side Turbo-Pump (HPFTP) of the Space Shuttle Main Engine (SSME). The stage calculation is performed by coupling the stator and the rotor flow fields at each time step through an over-laid grid. Numerical results for the complete geometry with the vane trailing edge cutback are presented and compared with the available experimental data.

Hah, C.↗
NASA NTRS1993

Numerical simulations of three-dimensional laminar flow over a backward facing step; flow near side walls

This paper reports the results of numerical simulations of steady, laminar flow over a backward-facing step. The governing equations used in the simulations are the full 'compressible' Navier-Stokes equations, solutions to which were computed by using a cell-centered, finite volume discretization. The convection terms of the governing equations were discretized by using the Advection Upwind Splitting Method (AUSM), whereas the diffusion terms were discretized using central differencing formulas. The validity and accuracy of the numerical solutions were verified by comparing the results to existing experimental data for flow at identical Reynolds numbers in the same back step geometry. The paper focuses attention on the details of the flow field near the side wall of the geometry.

Steinthorsson, Erlendur↗
NASA NTRS1992

An implicit multigrid scheme for hypersonic strong-interaction flowfields

An implicit multigrid method for the compressible Navier-Stokes equations is applied to a series of two-dimensional, steady, laminar hypersonic flows over a compression ramp, including both attached and separated flow conditions. The algorithm uses upwind-biased differencing for the convective and pressure terms and central differencing for the shear-stress and heat-transfer terms. An implicit spatially factored scheme is used to advance the solution in terms on a given mesh, and a V-cycle multigrid strategy is used to accelerate convergence. The multigrid algorithm provided an order of magnitude decrease in the computational time required to obtain a converged position of the separation location and enabled a substantially grid-independent result to be obtained for each of the flows considered.

Thomas, James L.↗
NASA NTRS1992

Viscous equilibrium computations using program LAURA

Modifications have been made to the Langley Aerothermodynamic Upwind Relaxation Algorithm (LAURA) that enable it to compute viscous airflows under the assumption of thermal and chemical equilibrium. Equilibrium thermodynamic and transport property information are input to the code via curve fits. The periodic updating of this information enables the equilibrium algorithm to perform at a computational rate that is only a small percentage larger than the rate associated with the perfect-gas algorithm. Presented in this article are the results of the initial validation of the modified code. Solutions for surface pressure and heating are presented for the flow over slender and blunt cones at realistic reentry conditions. LAURA solutions are compared with those produced by a viscous shock-layer method, and, for one case considered, with heat transfer data from a flight experiment. For both pressure and heating, the agreement is good. In general, differences in pressure of a few percent were noted, while differences in heating rates were in the 5-10 percent range.

Greene, Francis A.↗