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 307 records · Page 17
All sourcesNASA NTRSDOE OSTI
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.↗
NASA NTRS1984

TVD finite difference schemes and artificial viscosity

The total variation diminishing (TVD) finite difference scheme can be interpreted as a Lax-Wendroff scheme plus an upwind weighted artificial dissipation term. If a particular flux limiter is chosen and the requirement for upwind weighting is removed, an artificial dissipation term which is based on the theory of TVD schemes is obtained which does not contain any problem dependent parameters and which can be added to existing MacCormack method codes. Numerical experiments to examine the performance of this new method are discussed.

Davis, S. F.↗
NASA NTRS1985

Time-dependent Navier-Stokes computations of separated flows over airfoils

The application of both a central-difference and an upwind implicit approximate factorization Navier-Stokes algorithm to highly separated flow is described. Using the thin-layer approximation, both algorithms are employed to solve the low Reynolds number laminar flow around a circular cylinder with periodic shedding. Results agree with experimental Strouhal numbers to within about 10 percent. A grid refinement study is conducted using the upwind code for both the thin-layer and complete Navier-Stokes approximations. The circular cylinder solution is shown to be highly dependent upon both grid density outward from the body and grid extent. Finally, the thin-layer central-difference code is used to predict the time-accurate unsteady flow about two stalled airfoils. The NACA 0012 airfoil turbulent solutions are periodic with Strouhal numbers of about 0.1; a strong vortex is periodically shed from the leading edge. The LS(1)-0417 airfoil turbulent solution shows a high frequency oscillatory trailing edge separation. The 20 degree NACA 0012 airfoil solution is massively separated, but shows no evidence of periodicity.

Rumsey, C. L.↗
NASA NTRS1985

Large-scale computations in fluid mechanics; Proceedings of the Fifteenth Summer Seminar on Applied Mathematics, University of California, La Jolla, CA, June 27-July 8, 1983. Parts 1 & 2

Papers are presented on such topics as the use of semi-Lagrangian advective schemes in meteorological modeling; computation with high-resolution upwind schemes for hyperbolic equations; dynamics of flame propagation in a turbulent field; a modified finite element method for solving the incompressible Navier-Stokes equations; computational fusion magnetohydrodynamics; and a nonoscillatory shock capturing scheme using flux-limited dissipation. Consideration is also given to the use of spectral techniques in numerical weather prediction; numerical methods for the incorporation of mountains in atmospheric models; techniques for the numerical simulation of large-scale eddies in geophysical fluid dynamics; high-resolution TVD schemes using flux limiters; upwind-difference methods for aerodynamic problems governed by the Euler equations; and an MHD model of the earth's magnetosphere.

Engquist, B. E.↗
NASA NTRS1986

Application of program LAURA to three-dimensional AOTV flowfields

Program LAURA (Langley Aerothermodynamic Upwind Relaxation Algorithm) is a robust, finite volume, single-level storage, implicit upwind differencing algorithm which has been documented and tested on several three-dimensional blunt-body flows. The algorithm can run at unlimited Courant numbers (relaxing the steady-state equations) but requires the inversion of only a 5 x 5 matrix per computational cell. An alternating directional sweep Gauss-Seidel substitution strategy is used to relax the governing equations. At present, the Euler and thin-layer Navier-Stokes equations using Sutherland's law for viscosity have been modeled for a perfect gas, equilibrium air, and nonequilibrium air chemistry neglecting diffusion. The equilibrium and nonequilibrium air chemistry options have been described in a companion paper. Good comparisons with experimental data and another calculation method for pressure distributions, aerodynamic coefficients, and heat-transfer distributions have been demonstrated for three-dimensional blunt-body flows.

Gnoffo, P. A.↗
NASA NTRS1986

An efficient method for solving the steady Euler equations

An efficient numerical procedure for solving a set of nonlinear partial differential equations is given, specifically for the steady Euler equations. Solutions of the equations were obtained by Newton's linearization procedure, commonly used to solve the roots of nonlinear algebraic equations. In application of the same procedure for solving a set of differential equations we give a theorem showing that a quadratic convergence rate can be achieved. While the domain of quadratic convergence depends on the problems studied and is unknown a priori, we show that firstand second-order derivatives of flux vectors determine whether the condition for quadratic convergence is satisfied. The first derivatives enter as an implicit operator for yielding new iterates and the second derivatives indicates smoothness of the flows considered. Consequently flows involving shocks are expected to require larger number of iterations. First-order upwind discretization in conjunction with the Steger-Warming flux-vector splitting is employed on the implicit operator and a diagonal dominant matrix results. However the explicit operator is represented by first- and seond-order upwind differencings, using both Steger-Warming's and van Leer's splittings. We discuss treatment of boundary conditions and solution procedures for solving the resulting block matrix system. With a set of test problems for one- and two-dimensional flows, we show detailed study as to the efficiency, accuracy, and convergence of the present method.

Liou, M. S.↗
NASA NTRS1987

Numerical simulation by TVD schemes of complex shock reflections from airfoils at high angle of attack

The shock-capturing capability of total variation diminishing (TVD) schemes is demonstrated for a more realistic complex shock-diffraction problem for which the experimental data are available. Second-order explicit upwind and symmetric TVD schemes are used to solve the time-dependent Euler equations of gas dynamics for the interaction of a blast wave with an airfoil at high angle-of-attack. The test cases considered are a time-dependent moving curved-shock wave and a contant moving planar-shock wave impinging at an angle-of-attack 30 deg on a NACA 0018 airfoil. Good agreement is obtained between isopycnic contours computed by the TVD schemes and those from experimental interferograms. No drastic difference in flow-field structure is found between the curved- and planar-shock wave cases, except for a difference in density level near the lower surface of the airfoil. Computation for cases with higher shock Mach numbers is also possible. Numerical experiments show that the symmetric TVD scheme is less sensitive to the boundary conditions treatment than the upwind scheme.

Moon, Young J.↗
NASA NTRS1988

Numerical simulation of axisymmetric turbulent flow in combustors and diffusors

A method for predicting turbulent flow in combustors and diffusers is developed. The Navier-Stokes equations, incorporating a turbulence kappa-epsilon model equation, were solved in a nonorthogonal curvilinear coordinate system. The solution applied the finite volume method to discretize the differential equations and utilized the SIMPLE algorithm iteratively to solve the differenced equations. A zonal grid method, wherein the flow field was divided into several subsections, was developed. This approach permitted different computational schemes to be used in the various zones. In addition, grid generation was made a more simple task. However, treatment of the zonal boundaries required special handling. Boundary overlap and interpolating techniques were used and an adjustment of the flow variables was required to assure conservation of mass, momentum and energy fluxes. The numerical accuracy was assessed using different finite differencing methods, i.e., hybrid, quadratic upwind and skew upwind, to represent the convection terms. Flows in different geometries of combustors and diffusers were simulated and results compared with experimental data and good agreement was obtained.

Yung, Chain Nan↗
NASA NTRS1987

Recent progress in the development and understanding of SUPG methods with special reference to the compressible Euler and Navier-Stokes equations

The current status of streamline-upwind/Petrov-Galerkin (SUPG) methods for the analysis of flow problems is surveyed in an analytical review. Problem areas addressed include classical Galerkin, upwind, artificial-diffusion, SUPG, discontinuous Galerkin, space-time FEM, and discontinuity-capturing approaches to the scalar advection-diffusion equation; incompressible flows; advective-diffusive systems; and the compressible Euler and Navier-Stokes equations. Graphs and diagrams are provided, and the good stability properties of state-of-the-art SUPG methods are pointed out.

Hughes, Thomas J. R.↗
NASA NTRS1988

A point implicit unstructured grid solver for the Euler and Navier-Stokes equations

An upwind finite element technique that uses cell centered quantities and implicit and/or explicit time marching has been developed for computing hypersonic laminar viscous flows using adaptive unstructured triangular grids. A structured grid of quadrilaterals is laid out near the body surface. For inviscid flows the method is stable at Courant numbers of over 100,000. A first order basic scheme and a higher order flux corrected transport (FCT) scheme have been implemented. This technique has been applied to the problem of predicting type III and IV shock wave interactions on a cylinder, with a view of simulating the pressure and heating rate augmentation caused by an impinging shock on the leading edge of a cowl lip of an engine inlet. The predictions of wall pressure and heating rates compare very well with experimental data. The flow features are very distinctly captured with a sequence of adaptively generated grids. The adaptive mesh generator and the upwind Navier-Stokes solver are combined in a set of programs called LARCNESS, an acronym for Langley Adaptive Remeshing Code and Navier-Stokes Solver.

Thareja, Rajiv R.↗
NASA NTRS1988

Third-order multi-dimensional Euler/Navier-Stokes solver

A nonoscillatory multidimensional advection code is presented which is applicable to the study of very high Reynolds number steady flow, and results are obtained for the oblique advection of a step profile. The scheme is based on an exponential upwinding or linear extrapolation refinement of a previous QUICK scheme, allowing sharp monotonic resolution of the step while retaining third-order accuracy. The resulting simple flux-limiting criteria have been applied to a higher order method and a formally third-order supercompressive algorithm. Results are presented for first-order, second-order, and third-order upwinding.

Leonard, B. P.↗