Search NASA⌕ Search

SEARCH · Search NASA

Results for “wake function”

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

Progress on Optimizing Wind Farms and Rotor Designs Using Adjoints

Modern wind plants are increasingly tasked with multiple performance objectives. In addition to designing plants that maximize power output and minimize the levelized cost of energy (LCOE), the design and operation of wind plants is increasingly influenced by challenges regarding grid integration of variable generation renewables. This places a growing emphasis on making wind plants more controllable and predictable. WindSE is a Reynolds-averaged Navier-Stokes (RANS) model designed around analytical gradient and adjoint methods, with the ability to capture terrain-induced effects, as shown in Figure 1. The recent addition of an unsteady solver with an actuator line method (ALM) and ongoing work to enable massively parallel optimizations gives it a unique niche to explore coupled plant-level controls and design problems. This code is an open source python package built on the FEniCS framework that utilizes fast, parallel PETSc solvers to model fluid flow throughout wind-farm scale domains. Two recent studies performed using WindSE demonstrate the capability to optimize under a wide variety of flow conditions and objective functions. In the first, we present an optimization focused on modifying the layout of a wind farm with a fixed number of turbines for maximum total power output [1]. This study highlights the ability to quickly perform simulations using the steady Navier-Stokes solver combined with rotors represented as actuator disks while also stressing the importance of capturing terrain-induced effects. Gradient-based optimization using the RANS equations is viable due to the inclusion of efficiently computed adjoint derivatives. We interpret the physical results of the optimal layout and also discuss the computational cost of scaling to larger problems. In the second study, we present the capabilities of the unsteady Navier-Stokes solver, where rotor-blade profiles represented by actuator lines are optimized to enhance wake steering effects and overall power production [2]. We quantify the wind plant performance gains obtained from this type of simultaneous control co-design optimization as compared to optimizing the blade design and yaw independently. Figure 2 shows the differences between a baseline two-turbine system and an optimized system where we fine-tune the blade chord profile. Results and challenges from each study are quickly summarized and used to motivate the current development efforts within WindSE. Current and future work is focused on enabling higher-resolution studies with more degrees of freedom through parallelization of both the simulation and optimization algorithms. We present benchmarking results to show that WindSE performs well in both weak- and strong-scaling tests and further demonstrate that the optimizer obtains the same convergence rates in both shared- and distributed-memory environments. Using larger wind farms, we can study deep-array effects within an optimization context, allowing the use of objective functions that have been previously unstudied. As an example, we present ongoing work on a blockage metric which characterizes the loss of available kinetic energy due to wake effects from multiple upstream turbines.

adjoint optimization↗

A generalized framework for reduced-order modeling of a wind turbine wake

A reduced-order model for a wind turbine wake is sought from large eddy simulation data. Fluctuating velocity fields are combined in the correlation tensor to form the kernel of the proper orthogonal decomposition (POD). Proper orthogonal decomposition modes resulting from the decomposition represent the spatially coherent turbulence structures in the wind turbine wake; eigenvalues delineate the relative amount of turbulent kinetic energy associated with each mode. Back-projecting the POD modes onto the velocity snapshots produces dynamic coefficients that express the amplitude of each mode in time. A reduced-order model of the wind turbine wake (wakeROM) is defined through a series of polynomial parameters that quantify mode interaction and the evolution of each POD mode coefficients. The resulting system of ordinary differential equations models the wind turbine wake composed only of the large-scale turbulent dynamics identified by the POD. Tikhonov regularization is used to recalibrate the dynamical system by adding additional constraints to the minimization seeking polynomial parameters, reducing error in the modeled mode coefficients. The wakeROM is periodically reinitialized with new initial conditions found by relating the incoming turbulent velocity to the POD mode coefficients through a series of open-loop transfer functions. The wakeROM reproduces mode coefficients to within 25.2%, quantified through the normalized root-mean-square error. Furthermore, a high-level view of the modeling approach is provided as a platform to discuss promising research directions, alternate processes that could benefit stability and efficiency, and desired extensions of the wakeROM.

17 WIND ENERGY↗

An experimental study of entrainment and transport in the turbulent near wake of a circular cylinder

Attention is given to an experimental investigation of transport processes in the near wake of a circular cylinder, for a Reynolds number of 140,000, in which an X-array of hot wire probes mounted on a pair of whirling arms was used for flow measurement. Rotation of the arms in a uniform flow applies a wide range of relative flow angles to these X-arrays, making them inherently self-calibrating in pitch. A phase signal synchronized with the vortex-shedding process allowed a sorting of the velocity data into 16 populations, each having essentially constant phase. An ensemble average for each population yielded a sequence of pictures of the instantaneous mean flow field in which the vortices are frozen, as they would be on a photograph. The measurements also yield nonsteady mean data for velocity, intermittency, vorticity, stress, and turbulent energy production, as a function of phase. Emphasis is given in a discussion of study results to the nonsteady mean flow, which emerges as a pattern of centers and saddles in a frame of reference that moves with the eddies. The kinematics of the vortex formation process are described in terms of the formation and evolution of saddle points between vortices in the first few diameters of the near wake.

Cantwell, B.↗

Forcing function effects on unsteady aerodynamic gust response. I - Forcing functions

The paper investigates the fundamental gust modeling assumption on the basis of a series of experiments performed in the Purdue Annular Cascade Research Facility. The measured unsteady flow fields are compared to linear-theory gust requirements. The perforated plate forcing functions closely resemble linear-theory forcing functions, with the static pressure fluctuations small and the periodic velocity vectors parallel to the downstream mean-relative flow angle over the entire periodic cycle. The airfoil forcing functions exhibit characteristics far from linear-theory gusts, with the alignment of the velocity vectors and the static pressure fluctuation amplitudes dependent on the rotor-loading condition, rotor solidity, and the inlet mean-relative flow angle. It is shown that airfoil wakes, both compressor and turbine, cannot be modeled with the boundary conditions of current state-of-the-art linear unsteady aerodynamic theory.

Henderson, Gregory H.↗

Aeroacoustic computation of cylinder wake flow

Sound generation by uniform flow over a cylinder at the Reynolds number of 200 is calculated. The incompressible, time dependent flow field is first computed by standard computational fluid dynamics techniques using a stream function/vorticity formulation. The sound radiation is then obtained by integrating over the flow field with the Coriolis acceleration as the source term and a low frequency Green's function technique. Time histories, spectra, and directivity of the radiated sound are obtained. Comparisons with experimental data, where possible, indicate the feasibility of this computational aeroacoustics approach.

Hardin, J. C.↗

Pyroelectric and electrocaloric effects in ferroelectric silicon-doped hafnium oxide thin films

The emergent ferroelectricity in HfO 2 -based systems has attracted significant attention as this simple binary high- k dielectric now offers the possibility of nonvolatile function. Here we employ zero-field and field-dependent pyroelectricity to show that thin films of silicon-doped HfO 2 do exhibit a broken-inversion symmetry and are indeed ferroelectric. In addition, the pyroelectric response is found to exhibit a wake-up behavior akin to the wake-up phenomenon observed in the ferroelectric polarization with electric-field cycling. Using polarization-electric field hysteresis measurements, this wake-up phenomenon is attributed to the presence of defect dipoles which explains the measured pyroelectric response. Finally, direct electrocaloric measurements are performed on these silicon-doped HfO 2 thin films, revealing an electrocaloric coefficient ~ 4 times larger in magnitude than that expected for the measured pyroelectric coefficient. This enhancement is explained using the plausible role played by defect dipoles that contribute to additional configurational or dipolar entropy.

36 MATERIALS SCIENCE↗

Large Eddy Simulation of Aircraft Wake Vortices in a Homogeneous Atmospheric Turbulence: Vortex Decay and Descent

The effects of ambient turbulence on decay and descent of aircraft wake vortices are studied using a validated, three-dimensional: large-eddy simulation model. Numerical simulations are performed in order to isolate the effect of ambient turbulence on the wake vortex decay rate within a neutrally-stratified atmosphere. Simulations are conducted for a range of turbulence intensities, by injecting wake vortex pairs into an approximately homogeneous and isotropic turbulence field. The decay rate of the vortex circulation increases clearly with increasing ambient turbulence level, which is consistent with field observations. Based on the results from the numerical simulations, simple decay models are proposed as functions of dimensionless ambient turbulence intensity (eta) and dimensionless time (T) for the circulation averaged over a range of radial distances. With good agreement with the numerical results, a Gaussian type of vortex decay model is proposed for weak turbulence: while an exponential type of Tortex decay model can be applied for strong turbulence. A relationship for the vortex descent based on above vortex decay model is also proposed. Although the proposed models are based on simulations assuming neutral stratification, the model predictions are compared to Lidar vortex measurements observed during stable, neutral, and unstable atmospheric conditions. In the neutral and unstable atmosphere, the model predictions appear to be in reasonable agreement with the observational data, while in the stably-stratified atmosphere, they largely underestimate the observed circulation decay with consistent overestimation of the observed vortex descent. The underestimation of vortex decay during stably-stratified conditions suggests that stratification has an important influence on vortex decay when ambient levels of turbulence are weak.

Han, Jongil↗

Improving the FLORIS wind plant model for compatibility with gradient-based optimization

The FLORIS (FLOw Redirection and Induction in Steady-state) model, a parametric wind turbine wake model that predicts steady-state wake characteristics based on wind turbine position and yaw angle, was developed for optimization of control settings and turbine locations. This article provides details on changes made to the FLORIS model to make the model more suitable for gradient-based optimization. Changes to the FLORIS model were made to remove discontinuities and add curvature to regions of non-physical zero gradient. Exact gradients for the FLORIS model were obtained using algorithmic differentiation. A set of three case studies demonstrate that using exact gradients with gradient-based optimization reduces the number of function calls by several orders of magnitude. The case studies also show that adding curvature improves convergence behavior, allowing gradient-based optimization algorithms used with the FLORIS model to more reliably find better solutions to wind farm optimization problems.

17 WIND ENERGY↗

Assessment of Scaled Rotors for Wind Tunnel Experiments.

Rotor design and analysis work has been performed to support the conceptualization of a wind tunnel test focused on studying wake dynamics. This wind tunnel test would serve as part of a larger model validation campaign that is part of the Department of Energy Wind and Water Power Program’s Atmosphere to electrons (A2e) initiative. The first phase of this effort was directed towards designing a functionally scaled rotor based on the same design process and target full-scale turbine used for new rotors for the DOE/SNL SWiFT site. The second phase focused on assessing the capabilities of an already available rotor, the G1, designed and built by researchers at the Technical University of München.

42 ENGINEERING↗

Reference floating wind array designs for three representative regions

This work presents the systematic development of three open-source reference floating wind array designs. The designs are tailored to representative site conditions for three regions of the United States: Humboldt Bay off the coast of California, the Gulf of Maine, and the Gulf of America. We adopted existing reference designs for the individual 15 MW turbines, semisubmersible floating platforms, substations, mooring systems, and power cables – integrating and adapting them as needed for each location. We adapted existing dynamic cable designs to use larger conductor sizes to meet the arrays' power transmission requirements, and we set up redundant mooring systems for each substation. The layout of each array is a uniform-grid design optimized to approximately minimize the levelized cost of energy (LCOE) within a square lease area while satisfying spatial constraints. These constraints ensure adequate clearances between adjacent turbines and between underwater components during the layout optimization to prevent clashing and ensure that all components reside within the lease boundaries. Substations are included to allow accounting for intra-array cable costs. They are placed within the uniform grid to maintain the navigability of the arrays. For each feasible layout considered, annual energy production and cable routing costs are calculated and updated in the LCOE objective function. After the optimization, we refined the cable routing with a mix of algorithmic and manual methods to ensure that the cables avoid mooring system components and approach the substation with adequate clearances. We confirmed the suitability of each reference array's layout by comparing the wake losses at each wind heading angle to the wind rose, observing that the optimized layouts largely avoid wake losses in the predominant wind directions. These reference arrays provide open-source baseline designs to enable future research and innovation of floating wind technology at the array scale.

17 WIND ENERGY↗

The Spectral and Statistical Properties of Turbulence Generated by a Vortex/Blade-Tip Interaction

The perpendicular interaction of a streamwise vortex with the tip of a lifting blade was studied in incompressible flow to provide information useful to the accurate prediction of helicopter rotor noise and the understanding of vortex dominated turbulent flows. The vortex passed 0.3 chord lengths to the suction side of the blade tip, providing a weak interaction. Single and two-point turbulence measurements were made using sub-miniature four sensor hot-wire probes 15 chord lengths downstream of the blade trailing edge; revealing the mean velocity and Reynolds stress tensor distributions of the turbulence, as well as its spanwise length scales as a function of frequency. The single point measurements show the flow downstream of the blade to be dominated by the interaction of the original tip vortex and the vortex shed by the blade. These vortices rotate about each other under their mutual induction, winding up the turbulent wakes of the blades. This interaction between the vortices appears to be the source of new turbulence in their cores and in the region between them. This turbulence appears to be responsible for some decay in the core of the original vortex, not seen when the blade is removed. The region between the vortices is not only a region of comparatively large stresses, but also one of intense turbulence production. Velocity autospectra measured near its center suggests the presence quasi-periodic large eddies with axes roughly parallel to a line joining the vortex cores. Detailed two-point measurements were made on a series of spanwise cuts through the flow so as to reveal the turbulence scales as they would be seen along the span of an intersecting airfoil. The measurements were made over a range of probe separations that enabled them to be analyzed not only in terms of coherence and phase spectra but also in terms of wave-number frequency (kappa-omega) spectra, computed by transforming the measured cross-spectra with respect to the spanwise separation of the probes. These data clearly show the influence of the coherent eddies in the spiral wake and the turbulent region between the cores. These eddies produce distinct peaks in the upwash velocity kappa-omega spectra, and strong anisotropy manifested both in the decay of the kappa-omega spectrum at larger wave-numbers and in differences between the kappa-omega spectra of different components. None of these features are represented in the von Karman spectrum for isotropic turbulence that is often used in broadband noise computations. Wave-number frequency spectra measured in the cores appear to show some evidence that the turbulence outside sets tip core waves, as has previously been hypothesized. These spectra also provide for the first time a truly objective method for distinguishing velocity fluctuations produced by core wandering from other motions.

Devenport, William J.↗

The Effects of Gravity on the Circadian Timing System

All vertebrates have a physiological control system that regulates the timing of the rhythms of their daily life. Dysfunction of this system, the circadian timing system (CTS), adversely affects an organism's ability to respond to environmental challenges and has been linked to physiological and psychological disorders. Exposure to altered gravitational environments (the microgravity of space and hyperdynamic environments produced via centrifugation) produces changes in both the functioning of the CTS and the rhythmic variables it controls. The earliest record of primate rhythms in a spaceflight environment come from Biosatellite III. The subject, a pig-tailed macaque, showed a loss of synchronization of the body temperature rhythm and a fragmented sleep-wake cycle. Alterations in the rhythm of body temperature were also seen in rhesus macaques flown on COSMOS 1514. Squirrel monkeys exposed to chronic centrifugation showed an initial decrease in the amplitude and mean of their body temperature and activity rhythms. In a microgravity environment, Squirrel monkeys on Spacelab-3 showed a reduction in the mean and amplitude of their feeding rhythms. Since 1992 we have had the opportunity to participate on three US/Russian sponsored biosatellite missions on which a total of six juvenile male rhesus macaques were flown. These animals uniformly exhibited delays in the phasing of their temperature rhythms, but not their heart rate or activity rhythms during spaceflight. There was also a tendency for changes in waveform mean and amplitude. These data suggest that the spaceflight environment may have a differential effect on the different oscillators controlling different rhythmic variables. Ongoing studies are examining the effects of +G on the CTS. The long-term presence of humans in space highlights the need for effective countermeasures to gravitational effects on the CTS.

Fuller, Charles A.↗

Unsteady Surface Pressure Characterization on Launch Vehicle in Liftoff Configuration

Acquisition of unsteady pressure using remote scanners can provide practical advantages over making the same measurements directly using flush-mounted pressure transducers. However, pneumatic distortion of the signal is introduced by the tubing between the sensing module and surface port. This typically limits the usefulness of such measurements to averaged, steady data only, restricting the application of pressure scanners for unsteady pressure measurements in the wind tunnel and flight-testing applications. In this work, unsteady pressure on the surface of a 1.75%-scale model of the NASA Space Launch System Block 1B Cargo vehicle was remotely measured at the NASA Langley Research Center 14- by 22-Foot Subsonic Tunnel using an ESP pressure scanner. Subsequent reconstruction using an empirically-tuned, Wiener-filtered inverse transfer function produced excellent agreement compared to nearby surface-mounted transducers within a usable bandwidth of approximately 500 Hz. Interactions between the launch vehicle and the support tower produced strong oscillations that were absent in the vehicle-alone configuration, and it is likely that the wake interaction leads to a resonant flow response when the launch tower is immediately downstream of the vehicle.

Unsteady pressure↗

The hot companion and circumbinary disk of the be star HR 2142

We present a spectroscopic investigation of the Be+sdO binary system HR 2142 that is based upon large sets of ultraviolet observations from the International Ultraviolet Explorer (IUE) and ground-based Hα observations. We measured radial velocities for the Be star component from these spectra and computed a revised orbit. In order to search for the spectral signature of the hot subdwarf, we cross-correlated the short wavelength end of each IUE spectrum with a model hot star spectrum, and then we used the predicted Doppler shifts of the subdwarf to shift-and-add all the cross-correlation functions to the frame of the subdwarf. This merged function shows the weak signal from the spectral lines of the hot star, and a best fit is obtained with a mass ratio M{sub 2}/M{sub 1}=0.07±0.02, companion temperature T{sub eff}⩾43±5 kK, projected rotational velocity Vsin i<30 km s{sup −1}, and a monochromatic flux ratio near 1170 Å of f{sub 2}/f{sub 1}>0.009±0.001. This hot subdwarf creates a one-armed spiral, tidal wake in the disk of the Be star, and we present a circumbinary disk model that can explain the occurrence of shell absorption lines by gas enhancements that occur where gas crossing the gap created by the subdwarf strikes the disk boundaries. The faint companion of HR 2142 may be representative of a significant fraction of Be stars with undetected former mass donor companion stars.

79 ASTRONOMY AND ASTROPHYSICS↗

Switching kinetics in ferroelectric zinc magnesium oxide thin films

The switching kinetics of RF-magnetron reactively sputtered ~200 nm thick Zn 1-x Mg x O (ZMO) ferroelectric thin films with x = 0.41 and x = 0.27 prepared on Pt/Ti/SiO 2 /Si substrates were studied at applied fields near the coercive fields, ranging from 3.4 to 5.1 MV cm -1 , and at temperatures ranging from room temperature to 100 °C. Polarization reversal in ZMO followss the Kolmogorov-Avrami-Ishibashi kinetics model for nucleation and growth at applied fields near the coercive field, and obeys an individual column switching (ICS) model at higher applied fields required to switch most of the spontaneous polarization. Switching in the high-applied field regime can be described using either the Gaussian or inverse gamma distribution functions depending on the mole fraction of magnesium in the film. The switching current transients of the high-Mg content ZMO film (x = 0.41) are always bi-modal, whereas the low Mg composition film (x = 0.27) is described by the Gaussian distribution initially following wake-up, and becomes bimodal with continued cycling. Rayleigh-like behavior of the dielectric constant revealed a simultaneous increase of the irreversible and decrease of the reversible contributions to the dielectric constant, which was ascribed to an increase in the density of mobile domain walls with cycling and resulted in faster switching.

30 DIRECT ENERGY CONVERSION↗

STS-78 Flight Day 15

On this fifteenth day of the STS-78 mission, the 4th of July, Cmdr. Terence T. Henricks, Pilot Kevin R. Kregel, Payload Cmdr. Susan J. Helms, Mission Specialists Richard M. Linnehan, Charles E. Brady, Jr., and Payload Specialists Jean-Jacques Favier, Ph.D. and Robert B. Thirsk, M.D., are awakened with Bruce Springsteen's 'Born in the USA' and Lee Greenwood's 'I'm Proud to be an American' to begin another a day on orbit. Mission Commander Tom Henricks responded to Mission Control's wake up call by saying that the five US-born crew members were very proud to be Americans, particularly on the day America celebrates its 220th anniversary. Work in the Spacelab module will continue with investigations into the effects of microgravity on muscle strength and endurance, lung function, and adaptation of the neurovestibular system to a microgravity environment. Henricks and Pilot Kevin Kregel will complete work with a laptop computer designed to test the crew's critical thinking skills and reaction time. They also will test a voice control system that allows them to reposition Columbia's closed-circuit television cameras with verbal cues, keeping their hands free to perform other tasks.

Source record↗

Evaluation of Wind Speed Retrieval from Continuous-Wave Lidar Measurements of a Wind Turbine Wake Using Virtual Lidar Techniques

Wind lidar technology is being extensively used in utility-scale wind turbine wake measurements because of its capabilities to capture the dominant structures of interest in the near-wake region. However, wind lidar devices provide only line-of-sight (LOS) measurements, and retrieval of horizontal wind speed from the LOS measurements of a single lidar in the wind turbine wake could introduce errors due to the assumptions used in the retrieval process. With that in mind, the goals of this paper are to estimate the errors associated with the retrieval processes and provide guidelines for best practices. To achieve these goals, virtual lidar samples are obtained with large-eddy simulations (LES) by sampling the flow field with a model of the Technical University of Denmark (DTU) SpinnerLidar. LES is chosen to simulate the neutral and unstable atmospheric boundary layer (ABL), a benchmark case of the International Energy Agency Wind Technology Collaboration Programme (IEA Wind TCP) Task 31. The precursor ABL simulation is validated with the reference data provided by the IEA Wind Task 31 benchmark case. The data from the precursor simulation is used to drive a simulation of a wind turbine submerged in the ABL. A weighting function is applied to the LES data to mimic the real DTU SpinnerLidar measurements. In the end, the virtual lidar data are processed and compared with the original LES data to estimate the errors. Results show that lidar captures mean ABL profile without any noticeable error, whereas turbulence intensity is underpredicted ~15% at hub height for the neutral case. Volume averaging by the lidar technology has a significant effect in the shear layer of a wind turbine wake.

17 WIND ENERGY↗

Wind tunnel evaluation of air-foil performance using simulated ice shapes

A two-phase wind tunnel test was conducted in the 6 by 9 foot Icing Research Tunnel (IRT) at NASA Lewis Research Center to evaluate the effect of ice on the performance of a full scale general aviation wing. In the first IRT tests, rime and glaze shapes were carefully documented as functions of angle of attack and free stream conditions. Next, simulated ice shapes were constructed for two rime and two glaze shapes and used in the second IRT tunnel entry. The ice shapes and the clean airfoil were tapped to obtain surface pressures and a probe used to measure the wake characteristics. These data were recorded and processed, on-line, with a minicomputer/digital data acquisition system. The effect of both rime and glaze ice on the pressure distribution, Cl, Cd, and Cm are presented.

Bragg, M. B.↗