Search NASA⌕ Search

SEARCH · Search NASA

Results for “Engine Simulation”

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 1,333 records · Page 74

SISGR: Chemomechanics of Far-From-Equilibrium Interfaces (COFFEI)

Portable, reliable, and deployable devices for energy storage and conversion require fundamental changes in design of solid-state composites comprising ceramics and metals. These materials comprise the electrodes and electrolytes of next-generation solid-oxide fuel cells and all solid-state batteries, forming solid-state functional composites. The advent of solid-state batteries – which replace liquid electrolytes with solid electrolytes capable of lithium ion transport for reliable energy storage in portable batteries – and the increased demand space for all solid-state fuel cells capable of oxygen reduction at intermediate temperatures remain important challenges for improved material stability and decreased system cost. However, little is understood about three fundamental facets of materials that enable such solid-state energy applications. First, how do such materials deform, fracture, or delaminate under operando conditions? Second, how does such mechanical deformation limit or facilitate electronic and ionic transport within and across such material interfaces? Third, how we can predictably design interface-rich composites to engineer both structural and electrochemical stability? This COFFEI Group comprised expertise from Materials Science & Engineering and Nuclear Science & Engineering to integrate unique in situ experiments, simulations, and fabricated interfaces that address these fundamental questions in solid-state interfaces of nanoscale composites that will guide solid-state electrochemistry, transport kinetics, and mechanical deformation for nonstoichiometric materials that enable such applications. In particular, we built on COFFEI’s understanding of chemomechanical coupling among defect concentrations, ionic transport, electron transport, and stored elastic energy that is particularly acute in the far-from-equilibrium conditions typical of energy device applications. By tailoring our focus to solid-state interfaces, we addressed these important issues by (a) developing and applying advanced in situ and ex situ characterization tools to characterize model materials and interfaces synthesized with molecular-level control, under both laboratory-controlled and extreme environments representative of energy device operation; and (b) employing computational modeling and simulation frameworks to predict transport mechanisms, reactivity and stability of these model materials and interfaces under significant chemical strains typical of energy device operation. Recent progress provided insights to additional materials systems and electrochemomechanical fatigue and fracture that were not fully envisioned when the program was initiated. Specifically, in the final three years of COFFEI we pursued two integrated thrusts, with complementary focus. Thrust I focused on failure-resistant electrochemomechanical composites, while Thrust II focused on strain-modulated conductivity and reactivity across interfaces. In contrast to our initial COFFEI focus, these thrusts concentrated wholly on solid-state material interfacial interactions and included greater integration of multiscale visualization including in situ electron microscopy of strained structures/interactions and mesoscale simulations. Successful development of functionally superior and long-lived battery and fuel cell systems and stress adaptable oxides requires a deeper, fundamental understanding of the coupling among the historically important subfields of solid-state electrochemistry, transport kinetics, and mechanical deformation for nonstoichiometric metal oxide electrodes. In this program, the understanding and the application of chemomechanical coupling of defect concentrations, ionic transport, electro-catalytic activity and stored elastic energy, particularly acute in the far-from-equilibrium conditions typical of energy device applications, are being refined and implications for device operation clarified, including for miniaturized solid-state batteries and fuel cells.

36 MATERIALS SCIENCE↗

Simulating target tests at AP0 facility for Mu2e

The predicted production of muons and energy deposition for the production target of the Mu2e experiment are dependent on the documentation of cross-sections values, the properties of the target’s material, etc. Hence, making real tests on similar targets before its use on Mu2e will be a good indicator of the accuracy of the simulations done for the experiment. The main goal of the project is to simulate a target test at the AP0 Target Facility, previously used to create the muon beam for the Muon g-2 experiment. Results from these simulations will be used as input to engineering models of the target to produce predicted stress and surface temperature distributions, parameters that will later be measured in future beam runs to validate underlying cross section models.

Delgado, Amii Matamoros↗

Technical Assessment of Proposed LAW Simulant for Cold Commissioning

Prior to radioactive operations, the Hanford Tank Waste Treatment and Immobilization Plant (WTP) will implement a cold-commissioning plan or strategy that will include water runs and other melter feed simulants to not only test the various engineering systems, unit operations, and processing strategies with respect to safety and performance, but to also allow facility operations to gain operational experience before introducing radioactive material into the vitrification facilities. Based on direction from Bechtel National, Inc. to eliminate nitrogen oxide (NOx) generation, EnergySolutions has recommended the use of a simple NaOH solution with glass formers to support this stage of testing during cold commissioning.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Waveform Simulation Framework: User Manual with Tutorials

This manuscript describes the Waveform Simulation Framework (WSF), a Python-based framework that provides a unified, programmable interface for generating synthetic seismograms for applications such as seismic array design, method development, and special event analysis. WSF standardizes how users define sources, receivers, and velocity models while abstracting simulator-specific configuration details, enabling workflows that are largely independent of the underlying numerical engine. The document provides installation guidance and tutorial-driven examples for three WSF simulator wrappers—WSF PyFK, WSF SW4, and WSF SPECFEM2D—illustrating end-to-end workflows from forward waveform simulation to common post-processing tasks (e.g., visualization and backprojection) using consistent data products (e.g., ObsPy Stream objects and SAC files).

97 MATHEMATICS AND COMPUTING↗

Fabricated devices for performing bacterial-fungal interaction experiments across scales

Diverse and complex microbiomes are found in virtually every environment on Earth. Bacteria and fungi often co-dominate environmental microbiomes, and there is growing recognition that bacterial-fungal interactions (BFI) have significant impacts on the functioning of their associated microbiomes, environments, and hosts. Investigating BFI in vitro remains a challenge, particularly when attempting to examine interactions at multiple scales of system complexity. Fabricated devices can provide control over both biotic composition and abiotic factors within an experiment to enable the characterization of diverse BFI phenotypes such as modulation of growth rate, production of biomolecules, and alterations to physical movements. Engineered devices ranging from microfluidic chips to simulated rhizosphere systems have been and will continue to be invaluable to BFI research, and it is anticipated that such devices will continue to be developed for diverse applications in the field. This will allow researchers to address specific questions regarding the nature of BFI and how they impact larger microbiome and environmental processes such as biogeochemical cycles, plant productivity, and overall ecosystem resilience. Devices that are currently used for experimental investigations of bacteria, fungi, and BFI are discussed herein along with some of the associated challenges and several recommendations for future device design and applications.

59 BASIC BIOLOGICAL SCIENCES↗

BINARY NEUTRON STAR MERGERS: A JET ENGINE FOR SHORT GAMMA-RAY BURSTS

We perform magnetohydrodynamic simulations in full general relativity (GRMHD) of quasi-circular, equal-mass, binary neutron stars that undergo merger. The initial stars are irrotational, n = 1 polytropes and are magnetized. We explore two types of magnetic-field geometries: one where each star is endowed with a dipole magnetic field extending from the interior into the exterior, as in a pulsar, and the other where the dipole field is initially confined to the interior. In both cases the adopted magnetic fields are initially dynamically unimportant. The merger outcome is a hypermassive neutron star that undergoes delayed collapse to a black hole (spin parameter a / M {sub BH} ∼ 0.74) immersed in a magnetized accretion disk. About 4000 M ∼ 60( M {sub NS}/1.625 M {sub ⊙}) ms following merger, the region above the black hole poles becomes strongly magnetized, and a collimated, mildly relativistic outflow—an incipient jet—is launched. The lifetime of the accretion disk, which likely equals the lifetime of the jet, is Δ t ∼ 0.1 ( M {sub NS}/1.625 M {sub ⊙}) s. In contrast to black hole–neutron star mergers, we find that incipient jets are launched even when the initial magnetic field is confined to the interior of the stars.

79 ASTRONOMY AND ASTROPHYSICS↗

Numerical Modeling of Supersonic Combustion in a Non- Premixed Rotating Detonation Engine

A computational fluid dynamics (CFD) model was developed to simulate a representative non-premixed RDE combustor. Unsteady Reynolds-Averaged Navier-Stokes (RANS) simulations were performed for the full combustor geometry (including the separate fuel and air injection ports), with hydrogen as the fuel and air as the oxidizer. Adaptive mesh refinement (AMR) was incorporated to achieve a good trade-off between model accuracy and computational expense. A finite-rate chemistry model was employed along with a 10-species detailed kinetic mechanism to describe the H2-Air combustion chemistry. Two different operating conditions were simulated, corresponding to the same global equivalence ratio of unity but different fuel and air mass flow rates. Numerical results were compared against available experimental data on detonation wave frequency/height, reactant fill height, oblique shock angle and axial pressure distribution in the channel. The CFD model predicted the sensitivity of these wave characteristics to operating conditions with good accuracy, both qualitatively and quantitatively. The present CFD model offers a potential capability to perform rapid design space exploration and/or performance optimization studies for realistic full-scale RDE configurations.

Pal, Pinaki↗

Auto-ignitive deflagration of CH4 blended DME/Air mixtures with composition stratification

Front propagation speeds from one dimensional simulations of DME/Air and DME/CH4/Air mixtures under engine relevant conditions are presented with reduced kinetics and transport. Different time scales of monochromatic inhomogeneities in DME concentration and varying DME-CH4 blending ratios are explored to understand its subsequent effect on the propagation speed as well as the ensuing transition from deflagration to spontaneous propagation. It is found that for a given timescale, the instantaneous peak propagation speed achieved by a front as well as its overall variation is significantly affected by the level of CH4 concentration in the mixture. Specifically, at relatively smaller time scales, the magnitude of variation in the instantaneous propagation speed is found to subside as the level of CH4 concentration in the mixture is increased. An opposite trend is observed at comparatively larger time scales, wherein, a rapid change in the instantaneous propagation speed is observed with an increase in the level of CH4 concentration. This study provides evidence as to how the propagation speed of a single fuel flame varies compared to a dual fuel flame under the influence of composition stratification.

Desai, Swapnil S.↗

Thermal and Mechanical Energy Performance Analysis of Closed-loop Systems in Hot-Dry-Rock and Hot-Wet-Rock Reservoirs

To understand the potential and limitations for recovering thermal and mechanical energy from closed-loop geothermal systems a collaborative study is underway that will investigate an array of system configurations, working fluids, geothermal reservoir characteristics, operational periods, and heat transfer enhancements. Closed-loop geothermal systems are distinguished from hydrothermal or enhanced geothermal systems (EGS) in that the working fluid only circulates through drilled boreholes. Principal objectives of this study are to determine upper limits for thermal and mechanical energy recovery and optimal operational and configuration parameters for each scenario. Teams of scientists and engineers are applying a suite of numerical simulation and analytical tools to model the heat recovery from closed-loop geothermal systems, and then optimizing operational and configuration parameters to maximize the thermal and mechanical energy recovery. Results from the suite of numerical simulators and analytical tools, such as outlet and inlet states and temperature profiles in the geothermal reservoir over time are intercompared to increase confidence in the analysis. This paper documents the study findings for closed-loop systems in hot-dry-rock and hot-wet-rock reservoirs, where water is the working fluid. The characteristics of the hot-dry-rock reservoir were based on the U.S. Department of Energy's Utah Frontier Observatory for Research in Geothermal Energy (FORGE) site, near Milford Utah. Two objective functions are defined to optimize the operational and configuration parameters of the system, one each for the recovery of mechanical and thermal energy over the period of operation. For both objective functions, a surface plant thermal to mechanical energy conversion factor and an energy drilling cost is required. In keeping with the study objectives the surface plant conversion factor is determined from a second-law of thermodynamics analysis of a generic binary plant, and drilling costs are based on those from the Utah FORGE site and current national electrical costs.

closed-loop geothermal systems↗

Thermal and Mechanical Energy Performance Analysis of Closed-loop Systems in Hot-Dry-Rock and Hot-Wet-Rock Reservoirs

To understand the potential and limitations for recovering thermal and mechanical energy from closed-loop geothermal systems a collaborative study is underway that will investigate an array of system configurations, working fluids, geothermal reservoir characteristics, operational periods, and heat transfer enhancements. Closed-loop geothermal systems are distinguished from hydrothermal or enhanced geothermal systems (EGS) in that the working fluid only circulates through drilled boreholes. Principal objectives of this study are to determine upper limits for thermal and mechanical energy recovery and optimal operational and configuration parameters for each scenario. Teams of scientists and engineers are applying a suite of numerical simulation and analytical tools to model the heat recovery from closed-loop geothermal systems, and then optimizing operational and configuration parameters to maximize the thermal and mechanical energy recovery. Results from the suite of numerical simulators and analytical tools, such as outlet and inlet states and temperature profiles in the geothermal reservoir over time are intercompared to increase confidence in the analysis. This paper documents the study findings for closed-loop systems in hot-dry-rock and hot-wet-rock reservoirs, where water is the working fluid. The characteristics of the hot-dry-rock reservoir were based on the U.S. Department of Energy’s Utah Frontier Observatory for Research in Geothermal Energy (FORGE) site, near Milford Utah. Two objective functions are defined to optimize the operational and configuration parameters of the system, one each for the recovery of mechanical and thermal energy over the period of operation. For both objective functions, a surface plant thermal to mechanical energy conversion factor and an energy drilling cost is required. In keeping with the study objectives the surface plant conversion factor is determined from a second-law of thermodynamics analysis of a generic binary plant, and drilling costs are based on those from the Utah FORGE site and current national electrical costs.

Closed-loop geothermal systems, hot-dry-rock, hot-↗

Novel Solver Algorithms for Nearly Singular Linear Systems Arising in Combustion Modelling

Direct Numerical Simulations of realistic combustion devices are extremely challenging due to the wide separation of scales in the simulation, for example an internal combustion (IC) engine chamber, and the flame thickness of a high-pressure flame. The PeleLMeX solver uses adaptive mesh refinement (AMR) to evolve multi-species reacting flows in the low Mach number limit at the Exascale and relies on an embedded boundary (EB) approach to represent complex geometries. In that framework, the EB geometries often give rise to very small cut-cells along the boundary, which translate into extreme ill-conditioning of the pressure-projection, with eigenvalues that span 15-16 orders of magnitude. In this talk, we focus on the case of a typical IC piston bowl geometry for which we present on a novel approach towards solving these nearly singular linear systems with ILU-based, C-AMG smoothers on massively parallel architectures. In particular, we use scaling and equilibration algorithms to handle the non-normality of the upper triangular factors. This enables us to approximate the highly sequential triangular solve algorithm, embedded in the AMG smoothing-solve phase, with Jacobi iterations. This approximation can be written as a convergent Neumann series whose terms are composed of highly parallel sparse matrix vector multiplications. The result is an algorithm that substantially decreases setup and solve time, compared to state-of-the-art, for these challenging linear systems.

combustion modelling↗

Engineering scale demonstration of a prospective cast Stone process - 15584

A cementitious waste form is one of the alternatives being evaluated for Supplemental Immobilization of Hanford Low Activity Waste (LAW). Washington River Protection Solutions (WRPS) is collecting data on cementitious or pozzolanic waste forms such as Cast Stone. The goal of this project for WRPS is to obtain data on the performance of the Cast Stone waste form for immobilizing LAW. As part of the data package being assembled, an engineering-scale demonstration with non-radioactive simulants was performed at the Savannah River National Laboratory (SRNL) using the Scaled Continuous Processing Facility (SCPF) to fill a container with simulated Cast Stone grout. Legacy salt solution from previous Hanford salt waste testing was trimmed to resemble the average composition generated from the Hanford Tank Waste Operation Simulator (HTWOS) used in the screening tests. The dry blend materials, ordinary portland cement (OPC), Class F fly ash, and ground granulated blast furnace slag (GGBFS or BFS) were obtained from Lafarge North America in Pasco, WA. Over three days, the SCPF was used to fill a 6,056 liter (1600 gallon) container with simulated Cast Stone grout. The tank was instrumented with x-, y-, and z-axis thermocouples to monitor curing temperature and two formed core sampling vials. The target production rate was 5.7 L/min (1.5 gpm). This required a salt solution flow rate of approximately 3.78 L/min (1 gpm) and a premix feed rate of approximately 263 kh/h (580 lb/h). The final surface slope at a fill height of 100 cm was 2.5-3.8 cm across the 260 cm (8.5 ft) diameter tank. During processing, grout was collected from both the mixer discharge and the discharge into the tank. These samples were stored in a humid environment either in a closed box proximal to the tank or inside the laboratory. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Numerical Simulation of Recirculating Flow and Physical Model of Slag–Metal Behavior in an RH Reactor: Application to Desulfurization

Computational fluid dynamics (CFD) techniques and a 1:7.5 physical model of an RH degasser have been used to evaluate the flow of gas and metal inside an RH reactor for vacuum degassing of liquid steel. The effect of gas injection on the gas spatial distribution, steel circulation rate and flow field inside the ladle, snorkels and vacuum chamber have been assessed. N-pentane oil was employed to evaluate the average residence time as well as the slag droplet size distribution. The predicted radial gas distribution and liquid circulation rate have been validated against experimental data from a physical model. The results with incorporation of the virtual mass force coefficient of 0.25 and the turbulence dispersion force showed better predictions of gas distribution in the up-snorkel as well as circulation rate. Full-scale simulations were performed, and the predicted circulation rate was significantly affected by the argon bubble expansion. Data from these simulations were used to analyze the degree of desulfurization performed by the addition of desulfurizing agents inside the vacuum chamber. A model of the kinetics of desulfurization based on the results from the physical model and CFD simulation and on slag dispersion inside liquid steel yields degrees of desulfurization similar to the industrial trials reported in the literature.

36 MATERIALS SCIENCE↗

Investigation of Cycle-to-Cycle Variations in Internal Combustion Engine Using Proper Orthogonal Decomposition

The understanding, modeling and control of the cycle-to-cycle variation (CCV) in the modern internal combustion engine (ICE) is a key scientific challenge to achieve stable engine operation. High CCV in the engine combustion chamber may contribute to partial burn, misfire and knock, which adversely affects the engine performance and may potentially damage the engine. The objective of the current study is to leverage high-fidelity numerical simulations to improve the understanding of the causes of CCV. Using the massively parallel code, Nek5000, multi-cycle, wall-resolved large-eddy simulations (LES) were performed for the General Motors (GM), Transparent Combustion Chamber (TCC-III) optical engine under motored operating conditions. Further, the large-scale structures of the in-cylinder flow were investigated using a triple proper orthogonal decomposition (POD) technique to explore the characteristics of different parts of the flow and their contributions to CCV. The kinetic energy of the subset of flow structures were determined and correlated between the intake and compression strokes. The insights from the analysis of the large-scale flow structures will be used to assist the development of improved engine designs with reduced CCV and enhance the engine performance.

33 ADVANCED PROPULSION SYSTEMS↗

High Performance Computing Facility Operational Assessment 2015: Oak Ridge Leadership Computing Facility

Oak Ridge National Laboratory’s (ORNL’s) Leadership Computing Facility (OLCF) continues to surpass its operational target goals: supporting users; delivering fast, reliable systems; creating innovative solutions for high-performance computing (HPC) needs; and managing risks, safety, and security aspects associated with operating one of the most powerful computers in the world. The results can be seen in the cutting-edge science delivered by users and the praise from the research community. Calendar year (CY) 2015 was filled with outstanding operational results and accomplishments: a very high rating from users on overall satisfaction that ties the highest-ever mark set in CY 2014; the greatest number of core-hours delivered to research projects; the largest percentage of capability usage since the OLCF began tracking the metric in 2009; and success in delivering on the allocation of 60, 30, and 10% of core hours offered for the INCITE (Innovative and Novel Computational Impact on Theory and Experiment), ALCC (Advanced Scientific Computing Research Leadership Computing Challenge), and Director’s Discretionary programs, respectively. These accomplishments, coupled with the extremely high utilization rate, represent the fulfillment of the promise of Titan: maximum use by maximum-size simulations. The impact of all of these successes and more is reflected in the accomplishments of OLCF users, with publications this year in notable journals Nature, Nature Materials, Nature Chemistry, Nature Physics, Nature Climate Change, ACS Nano, Journal of the American Chemical Society, and Physical Review Letters, as well as many others. The achievements included in the 2015 OLCF Operational Assessment Report reflect first-ever or largest simulations in their communities; for example Titan enabled engineers in Los Angeles and the surrounding region to design and begin building improved critical infrastructure by enabling the highest-resolution Cybershake map for Southern California to date. The Titan system provides the largest extant heterogeneous architecture for computing and computational science. Usage is high, delivering on the promise of a system well-suited for capability simulations for science. This success is due in part to innovations in tracking and reporting the activity on the compute nodes, and using this information to further enable and optimize applications, extending and balancing workload across the entire node. The OLCF continues to invest in innovative processes, tools, and resources necessary to meet continuing user demand. The facility’s leadership in data analysis and workflows was featured at the Department of Energy (DOE) booth at SC15, for the second year in a row, highlighting work with researchers from the National Library of Medicine coupled with unique computational and data resources serving experimental and observational data across facilities. Effective operations of the OLCF play a key role in the scientific missions and accomplishments of its users. Building on the exemplary year of 2014, as shown by the 2014 Operational Assessment Report (OAR) review committee response in Appendix A, this OAR delineates the policies, procedures, and innovations implemented by the OLCF to continue delivering a multi-petaflop resource for cutting-edge research. This report covers CY 2015, which, unless otherwise specified, denotes January 1, 2015, through December 31, 2015.

97 MATHEMATICS AND COMPUTING↗

Model Reduction by Generalized Falk Method for Efficient Field-Circuit Simulations

The Generalized Falk Method (GFM) for coordinate transformation, together with two model-reduction strategies based on this method, are presented for efficient coupled field-circuit simulations. Each model-reduction strategy is based on a decision to retain specific linearly-independent vectors, called trial vectors, to construct a vector basis for coordinate transformation. The reduced-order models are guaranteed to be stable and passive since the GFM is a congruence transformation of originally symmetric positive definite systems. We also show that, unlike the Pad´e-via-Lanczos (PVL) method, the GFM does not generate unstable positive poles while reducing the order of circuit problems. Further, the proposed GFM is also faster when compared to methods of the type Lanczos (or Krylov) that are already widely used in circuit simulations for electrothermal and electromagnetic problems. The concept of response participation factors is introduced for the selection of the trial vectors in the proposed model-reduction methods. Further, we present methods to develop simple equivalent circuit networks for the field component of the overall field-circuit system. The implementation of these equivalent circuit networks in circuit simulators is discussed. With the proposed model-reduction strategies, significant improvement on the efficiency of the generalized Falk method is illustrated for coupled field-circuit problems.

42 ENGINEERING↗

Predicting the 3D fatigue crack growth rate of small cracks using multimodal data via Bayesian networks: In-situ experiments and crystal plasticity simulations

Small crack propagation accounts for most of the fatigue life of engineering structures subject to high cycle fatigue loading conditions. Determining the fatigue crack growth rate of small cracks propagating into polycrystalline engineering alloys is critical to improving fatigue life predictions, thus lowering cost and increasing safety. In this work, cycle-by-cycle data of a small crack propagating in a beta metastable titanium alloy is available via phase and diffraction contrast tomography. Crystal plasticity simulations are used to supplement experimental data regarding the micromechanical fields ahead of the crack tip. Experimental and numerical results are combined into a multimodal dataset and sampled utilizing a non-local data mining procedure. Furthermore, to capture the propensity of body-centered cubic metals to deform according to the pencil-glide model, a non-local driving force is postulated. The proposed driving force serves as the basis to construct a data-driven probabilistic crack propagation framework using Bayesian networks as building blocks. The spatial correlation between the postulated driving force and experimental observations is obtained by analyzing the results of the proposed framework. Results show that the above correlation increases proportionally to the distance from the crack front until the edge of the plastic zone. Moreover, the predictions of the propagation framework show good agreement with experimental observations. Finally, we studied the interaction of a small crack with grain boundaries (GBs) utilizing various slip transmission criteria, revealing the tendency of a crack to cross a GB by propagating along the slip directions minimizing the residual Burgers vector within the GB.

36 MATERIALS SCIENCE↗

Full-System Linearization for Floating Offshore Wind Turbines in OpenFAST

The wind engineering community relies on multiphysics engineering software to run nonlinear time-domain simulations (e.g., for design-standards-based loads analysis). Although most physics involved in wind energy are nonlinear, linearization of the underlying nonlinear system equations is often advantageous to understand the system properties and exploit well-established methods and tools for analyzing linear systems. This paper presents the development of the new linearization functionality of the open-source engineering tool OpenFAST for floating offshore wind turbines, as well as the concepts and mathematical background needed to understand and apply it.

17 WIND ENERGY↗