Search NASA⌕ Search

SEARCH · Search NASA

Results for “Upwinding”

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

An improved PNS scheme for predicting complex three-dimensional hypersonic flows

Upwinding is incorporated into a numerical technique for predicting hypersonic viscous flows over lifting configurations at moderate angles of attack. A general real-gas flux-vector-splitting technique based on Van Leer's (1982) approach is employed to model upwinding, and three techniques are examined for flux-vector differencing. The three methods are evaluated by applying them to an axisymmetric configuration with a 10-deg afterbody flare. The results indicate that an oscillation-free shock front can be described by using first-order full upwinding across the embedded shock and central-differencing for the other zones. This combined approach is found to be highly convergent for the near-wall region, and its performance is examined for predicting a Mach 15 flow over a finned missile. Attention is given to the effects of gas chemistry which can significantly affect the flows over the missile configurations.

Bhutta, Bilal A.↗
NASA NTRS1992

Numerical simulation of steady and unsteady asymmetric vortical flow

The unsteady, compressible, thin-layer, Navier-Stokes (NS) equations are solved to simulate steady and unsteady, asymmetric, vortical laminar flow around cones at high incidences and supersonic Mach numbers. The equations are solved by using an implicit, upwind, flux-difference splitting (FDS), finite-volume scheme. The locally conical flow assumption is used and the solutions are obtained by forcing the conserved components of the flowfield vector to be equal at two axial stations located at 0.95 and 1.0. Computational examples cover steady and unsteady asymmetric flows around a circular cone and its control using side strakes. The unsteady asymmetric flow solution around the circular cone has also been validated using the upwind, flux-vector splitting (FVS) scheme with the thin-layer NS equations and the upwind FDS with the full NS equations. The results are in excellent agreement with each other. Unsteady asymmetric flows are also presented for elliptic- and diamond-section cones, which model asymmetric vortex shedding around round- and sharp-edged delta winds.

Kandil, Osama A.↗
NASA NTRS1992

Convection treatment for high Rayleigh number, laminar, natural convection calculation

The problem of natural convection in a square cavity of a Boussinesq incompressible fluid is computed using a consistent central-difference approximation, a consistent second-order upwind scheme, and the hybrid scheme for the energy equation. It was found that, for this problem, the second-order upwind scheme yields satisfactory solutions and that an apparently low-order combination of schemes, with the second-order upwind scheme for the momentum equations and the hybrid scheme for the energy equation, can produce satisfactory predictions for this problem.

Shyy, Wei↗
NASA NTRS2010

Updates to Multi-Dimensional Flux Reconstruction for Hypersonic Simulations on Tetrahedral Grids

The quality of simulated hypersonic stagnation region heating with tetrahedral meshes is investigated by using an updated three-dimensional, upwind reconstruction algorithm for the inviscid flux vector. An earlier implementation of this algorithm provided improved symmetry characteristics on tetrahedral grids compared to conventional reconstruction methods. The original formulation however displayed quantitative differences in heating and shear that were as large as 25% compared to a benchmark, structured-grid solution. The primary cause of this discrepancy is found to be an inherent inconsistency in the formulation of the flux limiter. The inconsistency is removed by employing a Green-Gauss formulation of primitive gradients at nodes to replace the previous Gram-Schmidt algorithm. Current results are now in good agreement with benchmark solutions for two challenge problems: (1) hypersonic flow over a three-dimensional cylindrical section with special attention to the uniformity of the solution in the spanwise direction and (2) hypersonic flow over a three-dimensional sphere. The tetrahedral cells used in the simulation are derived from a structured grid where cell faces are bisected across the diagonal resulting in a consistent pattern of diagonals running in a biased direction across the otherwise symmetric domain. This grid is known to accentuate problems in both shock capturing and stagnation region heating encountered with conventional, quasi-one-dimensional inviscid flux reconstruction algorithms. Therefore the test problems provide a sensitive indicator for algorithmic effects on heating. Additional simulations on a sharp, double cone and the shuttle orbiter are then presented to demonstrate the capabilities of the new algorithm on more geometrically complex flows with tetrahedral grids. These results provide the first indication that pure tetrahedral elements utilizing the updated, three-dimensional, upwind reconstruction algorithm may be used for the simulation of heating and shear in hypersonic flows in upwind, finite volume formulations.

Gnoffo, Peter A.↗
NASA NTRS2016

On the Synergy Between Numerics and Subgrid Scale Modeling in LES of Stratified Flows: Grid Convergence of a Stratocumulus-Topped Boundary Layer

The effectiveness of a linear upwinding scalar advection scheme to suppress numerical dispersion errors near sharp inversions in large-eddy simulations of a nocturnal stratocumulus-topped boundary layer is assessed. Linear upwinding is a trade-off between non-dissipative and non-linear positive definite advection schemes. It is shown that linear upwinding does not negatively impact the model's grid convergence properties and a sharp inversion free of numerical artifacts is maintained. Even though mean profiles and turbulence fluxes show good grid convergence characteristics the liquid water amount varies significantly with grid resolution. The entrainment rate is identical for all resolutions and independent of the liquid water amount. For the present stratocumulus case, the impact of cloud-top radiative cooling is negligible and turbulence is largely driven by convection emanating from the surface.

Matheou, Georgios↗
NASA NTRS1975

A study of the efficiency of various Navier-Stokes solvers

A comparative study of the efficiency of some finite difference methods for the solution of the Navier-Stokes equations was conducted. The study was restricted to the two-dimensional steady, uniform property vorticity-stream function equations. The comparisons were drawn by recording the CPU time required to obtain a solution as well as the accuracy of this solution using five numerical methods: central differences, first order upwind differences, second order upwind differences, exponential differences, and an ADI solution of the central difference equations. Solutions were obtained for two test cases: a recirculating eddy inside a square cavity with a moving top, and an impinging jet flow. The results show that whenever the central difference method is stable it generates results with a given accuracy for less CPU time than any other method.

Atias, M.↗
NASA NTRS1976

Time-asymptotic solutions of the Navier-Stokes equation for free shear flows using an alternating-direction implicit method

An uncoupled time asymptotic alternating direction implicit method for solving the Navier-Stokes equations was tested on two laminar parallel mixing flows. A constant total temperature was assumed in order to eliminate the need to solve the full energy equation; consequently, static temperature was evaluated by using algebraic relationship. For the mixing of two supersonic streams at a Reynolds number of 1,000, convergent solutions were obtained for a time step 5 times the maximum allowable size for an explicit method. The solution diverged for a time step 10 times the explicit limit. Improved convergence was obtained when upwind differencing was used for convective terms. Larger time steps were not possible with either upwind differencing or the diagonally dominant scheme. Artificial viscosity was added to the continuity equation in order to eliminate divergence for the mixing of a subsonic stream with a supersonic stream at a Reynolds number of 1,000.

Rudy, D. H.↗
NASA NTRS1977

Simulation model of erosion and deposition on a barchan dune

Erosion and deposition over a barchan dune near the Salton Sea, California, are modeled by bookkeeping the quantity of sand in saltation following streamlines of transport. Field observations of near surface wind velocity and direction plus supplemental measurements of the velocity distribution over a scale model of the dune are combined as input to Bagnold type sand transport formulas corrected for slope effects. A unidirectional wind is assumed. The resulting patterns of erosion and deposition compare closely with those observed in the field and those predicted by the assumption of equilibrium (downwind translation of the dune without change in size or geometry). Discrepancies between the simulated results and the observed or predicted erosional patterns appear to be largely due to natural fluctuations in the wind direction. The shape of barchan dunes is a function of grain size, velocity, degree of saturation of the oncoming flow, and the variability in the direction of the oncoming wind. The size of the barchans may be controlled by natural atmospheric scales, by the age of the dunes, or by the upwind roughness. The upwind roughness can be controlled by fixed elements or by sand in the saltation. In the latter case, dune scale is determined by grain size and wind velocity.

Howard, A. D.↗
NASA NTRS1978

Effects of rotor location, coning, and tilt on critical loads in large wind turbines

Several large (1500 kW) horizontal rotor configurations were analyzed to determine the effects on dynamic loads of upwind downwind rotor locations, coned and radial blade positions, and tilted and horizontal rotor axis positions. Loads were calculated for a range of wind velocities at three locations in the structure: (1) the blade shank; (2) the hub shaft; and (3) the yaw drive. Blade axis coning and rotor axis tilt were found to have minor effects on loads. However, locating the rotor upwind of the tower significantly reduced loads at all locations analyzed.

Spera, D. A.↗
NASA NTRS1977

Effects of rotor location, coning, and tilt on critical loads in large wind turbines

Three large (1500 kW) horizontal rotor configurations were analyzed to determine the effects on dynamic loads of upwind and downwind rotor locations, coned and radial blade positions, and tilted and horizontal rotor axis positions. Loads were calculated for a range of wind velocities at three locations in the structure: the blade shank, the hub shaft, and the yaw drive. Blade axis coning and rotor axis tilt were found to have little effect on loads. However, locating the rotor upwind of the tower significantly reduced loads at all locations analyzed.

Spera, D. A.↗
NASA NTRS1977

An empirical model for ocean radar backscatter and its application in inversion routine to eliminate wind speed and direction effects

Several regression models were tested to explain the wind direction dependence of the 1975 JONSWAP (Joint North Sea Wave Project) scatterometer data. The models consider the radar backscatter as a harmonic function of wind direction. The constant term accounts for the major effect of wind speed and the sinusoidal terms for the effects of direction. The fundamental accounts for the difference in upwind and downwind returns, while the second harmonic explains the upwind-crosswind difference. It is shown that a second harmonic model appears to adequately explain the angular variation. A simple inversion technique, which uses two orthogonal scattering measurements, is also described which eliminates the effect of wind speed and direction. Vertical polarization was shown to be more effective in determining both wind speed and direction than horizontal polarization.

Dome, G. J.↗
NASA NTRS1978

Electron heating of inflowing interstellar gas

The heating of inflowing interstellar gas by solar electrons is calculated. In the present work it is assumed that the thermal energy of the electrons is constant with distance from the Sun and that the inflowing gas follows Keplerian trajectories, and the contours of constant temperature so obtained are given. At a heliocentric distance of 1 AU the estimated temperature is 500 K upwind and 2000 K downwind for hydrogen atoms, and is 240 K downwind and 100 K upwind for helium.

Wu, F.-M.↗
NASA NTRS1979

A fast, conservative algorithm for solving the transonic full-potential equation

A fast, fully implicit approximate factorization (AF) algorithm designed to solve the conservative transonic full-potential equation in either two or three dimensions is described. The algorithm uses an upwind bias of the density coefficient for stability in supersonic regions. This provides an effective upwind difference of the streamwise terms for any orientation of the velocity vector (i.e., 'rotated differencing'), and thereby greatly enhances the reliability of the present algorithm. A numerical transformation is used to establish an arbitrary body-fitted finite-difference mesh. Computed results for both airfoils and simplified wings demonstrate substantial improvement in convergence speed for the new algorithm relative to standard successive-line overrelaxation algorithms.

Holst, T. L.↗
NASA NTRS1978

Steady state spurious errors in shock-capturing numerical schemes

The behavior of the steady state spurious error modes of the MacCormack scheme and the upwind scheme of Warming and Beam was obtained from a linearized difference equation for the steady state error. It was shown that the spurious errors can exist either as an eigensolution of the homogeneous part of this difference equation or because of excitation from large discretization errors near oblique shocks. It was found that the upwind scheme does not permit spurious oscillations on the upstream side of shocks. Examples are given for the inviscid Burgers' equation and for one and two dimensional gasdynamic flows.

Karlsen, L. K.↗
NASA NTRS1980

An implicit algorithm for the transonic full-potential equation in conservative form

A fast, implicit approximate factorization algorithm for the solution of the conservative full-potential equation for transonic flow in two and three dimensions is presented. Stability in supersonic regions is maintained by the use of an upwind evaluation of the density coefficient along all coordinate directions, providing an effective upwind difference of the streamwise terms for any orientation of the velocity vector and thereby enhancing the reliability of the algorithm. The algorithm is shown to provide rapid convergence for the computation of certain difficult two-dimensional test cases, including cases with fishtail shock patterns, demonstrating the reliability and efficiency of the procedure. Surface pressure coefficient distributions obtained by the present method are also found to be in good agreement with those computed by successive-line overrelaxation and a hybrid direct-solver/successive-line overrelaxation scheme, with significant reductions in CPU time required. A three-dimensional solution for a swept wing mounted between parallel walls is also presented which demonstrates the high convergence rate of the algorithm in three dimensions as well as two.

Holst, T.↗
NASA NTRS1981

Measured and calculated characteristics of wind turbine noise

The results of an analytical and experimental investigation of wind turbine noise are presented. Noise calculations indicate that for configurations with the rotor downwind of the support tower, the primary source of noise is the rapid change in rotor loadings which occurs as the rotor passes through the tower wake. Noise measurements are presented for solid and truss type tower models with both upwind and downwind rotors. Upwind rotor configurations are shown to be significantly quieter than downwind configurations. The model data suggest that averaged noise measurements and noise calculations based on averaged tower wake characteristics may not accurately represent the impulsive noise characteristics of downwind rotor configurations.

Greene, G. C.↗
NASA NTRS1982

Conservative supra-characteristics method for splitting the hyperbolic systems of gasdynamics with computed boundaries for real and perfect gases

Implicit methods developed by Beam and Warming (1978 and Briley and McDonald (1977) make it possible to overcome the hyperbolic stiffness of the conservative compressible Navier-Stokes equations in the fine wall region computational mesh for high Reynold's number flow. Certain difficulties related to the use of these methods could be overcome by employing an approach reported by Roe (1981). In the present investigation Roe's conceptual framework has been adopted for constructing globally conservative finite difference methods. A globally conservative upwind finite difference method (CSCM) consisting of both implicit interior point and boundary point equations is constructed from a new characteristics based flux difference splitting. It is found that the employed upwind eigenvector split scheme which combines fully coupled implicit interior point and boundary point approximations has the desired properties of robust stability and accuracy.

Lombard, C. K.↗
NASA NTRS1983

The movement and decay of ice edge bands in the winter Bering Sea

A study of the movement and decay of ice bands in the Bering Sea carried out in March 1981 is described. A pair of radio transponders mounted on ice floes within a band was tracked and the band motion was compared with that of a satellite-tracked buoy deployed in the ice interior. The results show that the bands moved away from the interior pack ice at a speed 30 percent greater than that of the interior ice, that the cause of this band acceleration was very likely the wind-wave radiation stress on the upwind side of the band, and that as the bands moved into warmer water, they decayed both by wind-wave erosion of the upwind edge and by bottom melting. The working of the buoys and the method of deployment are described, as are the band shape, trajectory, and decay. By a calculation of the steady state stress balance on the band, it is shown that the radiation stress can account for the velocity increase of the band relative to the ice interior.

Martin, S.↗