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 271 records · Page 15

Numerical experiments with a symmetric high-resolution shock-capturing scheme

Characteristic-based explicit and implicit total variation diminishing (TVD) schemes for the two-dimensional compressible Euler equations have recently been developed. This is a generalization of recent work of Roe and Davis to a wider class of symmetric (non-upwind) TVD schemes other than Lax-Wendroff. The Roe and Davis schemes can be viewed as a subset of the class of explicit methods. The main properties of the present class of schemes are that they can be implicit, and, when steady-state calculations are sought, the numerical solution is independent of the time step. In a recent paper, a comparison of a linearized form of the present implicit symmetric TVD scheme with an implicit upwind TVD scheme originally developed by Harten and modified by Yee was given. Results favored the symmetric method. It was found that the latter is just as accurate as the upwind method while requiring less computational effort. Currently, more numerical experiments are being conducted on time-accurate calculations and on the effect of grid topology, numerical boundary condition procedures, and different flow conditions on the behavior of the method for steady-state applications. The purpose here is to report experiences with this type of scheme and give guidelines for its use.

Yee, H. C.↗

Application of ray theory to propagation of low frequency noise from wind turbines

Ray theory is used to explain data from two experiments (1985 and 1985) on the propagation of low frequency sound generated by the WTS-4 wind turbine. Emphasis is on downwind data, but some upwind measurements taken during the 1985 experiment are also considered. General ray theory for a moving medium is reviewed and ray equations obtained. Restrictions are introduced simplifying the equations and permitting the use of a ray theory program MEDUSA, the computed propagation loss curve of which is compared to the measurements. Good qualitative agreement is obtained with 1984 downwind data. The results indicate that the downwind sound field is that of a near-ground sound channel. Although more scatter is seen in the 1985 data, agreement between theory and data is also good. In particular, the position and magnitude of the jump in the sound levels associated with the beginning of the sound channel is correctly predicted. The theoretical explanation of the upwind data is less successful. Ray theory calculations indicate the formation of a shadow zone that, in fact, does not occur. While no sharp shadow zone is apparent in the data, the general expectation (based on ray theory) that sound levels should be much reduced upwind is confirmed by the data.

Hawkins, James A.↗

Distorted-wave born approximation calculations for turbulence scattering in an upward-refracting atmosphere

Weiner and Keast observed that in an upward-refracting atmosphere, the relative sound pressure level versus range follows a characteristic 'step' function. The observed step function has recently been predicted qualitatively and quantitatively by including the effects of small-scale turbulence in a parabolic equation (PE) calculation. (Gilbert et al., J. Acoust. Soc. Am. 87, 2428-2437 (1990)). The PE results to single-scattering calculations based on the distorted-wave Born approximation (DWBA) are compared. The purpose is to obtain a better understanding of the physical mechanisms that produce the step-function. The PE calculations and DWBA calculations are compared to each other and to the data of Weiner and Keast for upwind propagation (strong upward refraction) and crosswind propagation (weak upward refraction) at frequencies of 424 Hz and 848 Hz. The DWBA calculations, which include only single scattering from turbulence, agree with the PE calculations and with the data in all cases except for upwind propagation at 848 Hz. Consequently, it appears that in all cases except one, the observed step function can be understood in terms of single scattering from an upward-refracted 'skywave' into the refractive shadow zone. For upwind propagation at 848 Hz, the DWBA calculation gives levels in the shadow zone that are much below both the PE and the data.

Gilbert, Kenneth E.↗

Development of new flux splitting schemes

Maximizing both accuracy and efficiency has been the primary objective in designing a numerical algorithm for computational fluid dynamics (CFD). This is especially important for solutions of complex three dimensional systems of Navier-Stokes equations which often include turbulence modeling and chemistry effects. Recently, upwind schemes have been well received for their capability in resolving discontinuities. With this in mind, presented are two new flux splitting techniques for upwind differencing. The first method is based on High-Order Polynomial Expansions (HOPE) of the mass flux vector. The second new flux splitting is based on the Advection Upwind Splitting Method (AUSM). The calculation of the hypersonic conical flow demonstrates the accuracy of the splitting in resolving the flow in the presence of strong gradients. A second series of tests involving the two dimensional inviscid flow over a NACA 0012 airfoil demonstrates the ability of the AUSM to resolve the shock discontinuity at transonic speed. A third case calculates a series of supersonic flows over a circular cylinder. Finally, the fourth case deals with tests of a two dimensional shock wave/boundary layer interaction.

Liou, Meng-Sing↗

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.↗

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.↗

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↗

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.↗

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↗

Floating platform effects on power generation in spar and semisubmersible wind turbines

The design and financing of commercial-scale floating offshore wind projects require a better understanding of how power generation differs between newer floating turbines and well-established fixed-bottom turbines. In floating turbines, platform mobility causes additional rotor motion that can change the time-averaged power generation. In this work, OpenFAST simulations examine the power generated by the National Renewable Energy Laboratory's 5-MW reference turbine mounted on the OC3-UMaine spar and OC4-DeepCWind semisubmersible floating platforms, subjected to extreme irregular waves and below-rated turbulent inflow wind from large-eddy simulations of a neutral atmospheric boundary layer. For these below-rated conditions, average power generation in floating turbines is most affected by two types of turbine displacements: an average rotor pitch angle that reduces power, caused by platform pitch; and rotor motion upwind-downwind that increases power, caused by platform surge and pitch. The relative balance between these two effects determines whether a floating platform causes power gains or losses compared to a fixed-bottom turbine; for example, the spar creates modest (3.1%–4.5%) power gains, whereas the semisubmersible creates insignificant (0.1%–0.2%) power gains for the simulated conditions. Furthermore, platform surge and pitch motions must be analyzed concurrently to fully capture power generation in floating turbines, which is not yet universal practice. Finally, a simple analytical model for predicting average power in floating turbines under below-rated wind speeds is proposed, incorporating effects from both the time-averaged pitch displacement and the dynamic upwind-downwind displacements.

17 WIND ENERGY↗

Control co-design of 13 MW downwind two-bladed rotors to achieve 25% reduction in levelized cost of wind energy

Wind energy is recognized worldwide as cost-effective and environmentally friendly and is among the fastest-growing sources of electrical energy. To further decrease the cost of wind energy, wind turbines are being designed at ever larger scales, which is challenging due to greater structural loads and deflections. Large-scale systems such as modern wind turbines increasingly require a control co-design approach, whereby the system design and control design are performed in a more integrated fashion. We overview a two-bladed downwind morphing rotor concept that is expected to lower the cost of energy at wind turbine sizes beyond 13 megawatts (MW) compared with continued upscaling of traditional three-bladed upwind rotor designs. We describe an aero-structural-control co-design process that we have used in designing such extreme-scale wind turbines, and we discuss how we were able to achieve a 25% reduction in levelized cost of energy for our final turbine design compared to a conventional upwind three-bladed rotor design.

17 WIND ENERGY↗

Estimation of wintertime cloud radiative effects in the Western Arctic, a function of cloud-moisture-coupling and sea ice conditions

Cloud radiative effects are estimated using up- and downward ground-based observations of broadband thermal infrared sensors in cloudy conditions which are provided by the ARM data repository. The ARM Surface Radiation measurement is analyzed. To study the cloud radiative effects, the observed thermal infrared measurements under cloudy conditions are analyzed together with microphysical cloud properties like liquid and ice water path obtained from the CloudNet algorithms applied to the remote-sensing instrumentation suite at the ARM site in Utqiagvik. Cloud observations have been sorted depending on the coupling of the cloud with upwind water vapor transport events that can interact with sea ice upwind. Statistics of surface net thermal fluxes and cloud water/ice properties show a distinguished distribution for coupled and decoupled clouds. It was found that coupled clouds contain more liquid water and are more efficient reduce the cooling of surface LW radiation.

54 ENVIRONMENTAL SCIENCES↗

Virtual tower measurements during the American WAKE ExperimeNt (AWAKEN)

Dual-Doppler lidar measurements were made during the American WAKE ExperimeNt to provide height-resolved measurements of wind speed and direction at multiple locations immediately south of the leading row turbines in the King Plains wind farm in Oklahoma. These so-called virtual tower measurements were performed to characterize the inflow into the wind farm and to assess possible upwind blockage effects due to the collective action of the wind farm. The campaign was conducted from 12 November 2022 to 17 October 2023, during which time 14 unique virtual tower locations were sampled with heights ranging from 240 to 490 m AGL. The wind retrieval algorithm provided estimates of the horizontal winds and their uncertainties with a vertical resolution of about 10 m, while also accounting for the tilt of the lidar platform. The virtual tower results are compared to collocated lidar wind profiling data at the A1 site, which was located roughly 2.4 rotor diameters south of the nearest turbine. The wind speed difference between the wind profiler and the virtual tower was found to be quite sensitive to atmospheric stability and wind direction below 250 m AGL. The largest differences were observed for inflow under stable conditions, where the profiler wind speeds were observed to be about 22% lower than the virtual tower near hub height. These results suggest that there are persistent horizontal gradients in the flow upwind of the wind farm which result in biased estimates using standard ground-based lidar wind profiling methods.

17 WIND ENERGY↗

Diffusion Limit–Preserving Lumped DFEMs on AMR Meshes

Here, we present sweep-compatible, novel upwinding recipes for the bilinear discontinuous (BLD) finite element method (FEM) that allows lumped BLD to be used on adaptive mesh refinement (AMR) meshes for thick transport applications without adding additional degrees of freedom at hanging nodes that exist on refinement boundaries. We analyze the properties of the upwinding and lumping that are needed for BLD to get the thick diffusion limit on such meshes, present results demonstrating locking with the wrong recipe, and present results showing error convergence and robustness properties for two diffusive problems on a variety of AMR meshes.

42 ENGINEERING↗

A Fast Algebraic Multigrid Solver and Accurate Discretization for Highly Anisotropic Heat Flux I: Open Field Lines

We present a novel solver technique for the anisotropic heat flux equation, aimed at the high level of anisotropy seen in magnetic confinement fusion plasmas. Such problems pose two major challenges: (i) discretization accuracy and (ii) efficient implicit linear solvers. We simultaneously address each of these challenges by constructing a new finite element discretization with excellent accuracy properties, tailored to a novel solver approach based on algebraic multigrid (AMG) methods designed for advective operators. We pose the problem in a mixed formulation, introducing the directional temperature gradient as an auxiliary variable. The temperature and auxiliary fields are discretized in a scalar discontinuous Galerkin space with upwinding principles used for discretizations of advection. We demonstrate the proposed discretization’s superior accuracy over other discretizations of anisotropic heat flux, achieving error 1000x smaller for anisotropy ratio of 10 9 , for closed field lines. The block matrix system is reordered and solved in an approach where the two advection operators are inverted using AMG solvers based on approximate ideal restriction, which is particularly efficient for upwind discontinuous Galerkin discretizations of advection. To ensure that the advection operators are nonsingular, in this paper we restrict ourselves to considering open (acyclic) magnetic field lines for the linear solvers. We demonstrate fast convergence of the proposed iterative solver in highly anisotropic regimes where other diffusion-based AMG methods fail.

97 MATHEMATICS AND COMPUTING↗

Performance Evaluation of a Two-Dimensional Flood Model on Heterogeneous High-Performance Computing Architectures

This paper describes the implementation of a two-dimensional hydrodynamic flood model with two different numerical schemes on heterogeneous high-performance computing architectures. Both schemes were able to solve the nonlinear hyperbolic shallow water equations using an explicit upwind first-order approach on finite differences and finite volumes, respectively, and were conducted using MPI and CUDA. Four different test cases were simulated on the Summit supercomputer at Oak Ridge National Laboratory. Both numerical schemes scaled up to 128 nodes (768 GPUs) with a maximum 98.2x speedup of over 1 GPU. The lowest run time for the 10 day Hurricane Harvey event simulation at 5 meter resolution (272 million grid cells) was 50 minutes. GPUDirect communication proved to be more convenient than the standard communication strategy. Both strong and weak scaling are shown.

Sharif, Md Bulbul↗

Loads Response That is Due to Wake Steering on a Pair of Utility-Scale Wind Turbines

The goal of this report is to present the resulting turbine loads of a full-scale wind turbine wake steering field campaign. In this experiment, wake steering controls were applied to an upwind turbine (T2) and the downwind turbine (T3) was instrumented to measure mechanical loads. The test subjects (T2 & T3) were GE 1.5SLE CWE turbines located on a wind farm with a prevailing wind direction coming from the northwest. The data collection strategy involved toggling the wake steering controls of T2 between on and off. This resulted in two databases (baseline and wake steered) with similar turbulence intensities. Valid data was extracted and processed. The analysis involved scaling data to engineering units, applying coordinate transformations where applicable, calculating 10-minutes statistics, and calculating damage equivalent loads (DELs) to assess fatigue. Figures are shown as statistics and DELs binned by wind speed with supplemental scatter plots provided in the appendices. Results between the baseline and wake steered databases were compared. Overall, the binned statistics showed minimal differences in loading between the two cases. However, the DELs of the wake steered case were observed to be consistently smaller than the baseline for most of the turbine components. A potential for future work was recognized based on the findings of this report. Some areas of further study include an analysis for more granular loads sensitivity to different yaw offset angles and rotor wake overlap, an in-depth comparison of loads between the upwind and downwind turbines, and further fatigue analysis to quantify differences observed in this experiment.

17 WIND ENERGY↗

Isolated Building Wake Experiment At Texas Tech University’s Wind Engineering Research Field Lab

Texas Tech University’s (TTU’s) Wind Engineering Research Field Lab (WERFL) building is an idealized rectangular building on the scale of a typical suburban house. It is surrounded almost entirely by unobstructed planes and is built on a turn table. WERFL is ideal for providing measurements of an idealized building under real-world conditions including atmospheric stability and variable relative wind angles. The intent of this data set is for verification and validation of computational fluid dynamics model simulations of flow around buildings. The primary focus of this dataset features 24 3D sonic anemometers (20 downwind and 4 upwind) and 10 2D sonic anemometers (9 downwind and 1 upwind) located in the immediate vicinity of the WERFL building. It also includes a variety of instruments intended to measure the ambient meteorological parameters and the undisturbed wind profile.

17 WIND ENERGY↗