Search NASA⌕ Search

SEARCH · Search NASA

Results for “exit time”

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 577 records · Page 32

Measurements in a Transitional Boundary Layer Under Low-Pressure Turbine Airfoil Conditions

This report presents the results of an experimental study of transition from laminar to turbulent flow in boundary layers or in shear layers over separation zones on a convex-curved surface which simulates the suction surface of a low-pressure turbine airfoil. Flows with various free-stream turbulence intensity (FSTI) values (0.5%, 2.5% and 10%), and various Reynolds numbers (50,000, 100,000 200,000 and 300,000) are investigated. Reynold numbers in the present study are based on suction surface length and passage exit mean velocity. Flow separation followed by transition within the separated flow region is observed for the lower-Re cases at each of the FSTI levels. At the highest Reynolds numbers and at elevated FSn, transition of the attached boundary layer begins before separation, and the separation zone is small. Transition proceeds in the shear layer over the separation bubble. For both the transitional boundary layer and the transitional shear layer, mean velocity, turbulence intensity and intermittency (the fraction of the time the flow is turbulent) distributions are presented. The present data are compared to published distribution models for bypass transition, intermittency distribution through transition, transition start position, and transition length. A model developed for transition of separated flows is shown to adequately predict the location of the beginning of transition, for these cases, and a model developed for transitional boundary layer flows seems to adequately predict the path of intermittency through transition when the transition start and end are known. These results are useful for the design of low-pressure turbine stages which are known to operate under conditions replicated by these tests.

Simon, Terrence W.↗

Unsteady Ejector Performance: an Experimental Investigation Using a Pulsejet Driver

An experimental investigation is described in which thrust augmentation and mass entrainment were measured for a variety of simple cylindrical ejectors driven by a gasoline-fueled pulsejet. The ejectors were of varying length, diameter, and inlet radius. Measurements were also taken to determine the effect on performance of the distance between pulsejet exit and ejector inlet. Limited tests were also conducted to determine the effect of driver cross-sectional shape. Optimal values were found for all three ejector parameters with respect to thrust augmentation. This was not the case with mass entrainment, which increased monotonically with ejector diameter. Thus, it was found that thrust augmentation is not necessarily directly related to mass entrainment, as is often supposed for ejectors. Peak thrust augmentation values of 1.8 were obtained. Peak mass entrainment values of 30 times the driver mass flow were also observed. Details of the experimental setup and results are presented. Preliminary analysis of the results indicates that the enhanced performance obtained with an unsteady jet (primary source) over comparably sized ejectors driven with steady jets is due primarily to the structure of the starting vortex-type flow associated with the former.

Paxson, Daniel E.↗

Disturbance of hibernating bats due to researchers entering caves to conduct hibernacula surveys

Estimating population changes of bats is important for their conservation. Population estimates of hibernating bats are often calculated by researchers entering hibernacula to count bats; however, the disturbance caused by these surveys can cause bats to arouse unnaturally, fly, and lose body mass. We conducted 17 hibernacula surveys in 9 caves from 2013 to 2018 and used acoustic detectors to document cave-exiting bats the night following our surveys. We predicted that cave-exiting flights (i.e., bats flying out and then back into caves) of Townsend’s big-eared bats (Corynorhinus townsendii) and western small-footed myotis (Myotis ciliolabrum) would be higher the night following hibernacula surveys than on nights following no surveys. Those two species, however, did not fly out of caves more than predicted the night following 82% of surveys. Nonetheless, the activity of bats flying out of caves following surveys was related to a disturbance factor (i.e., number of researchers × total time in a cave). We produced a parsimonious model for predicting the probability of Townsend’s big-eared bats flying out of caves as a function of disturbance factor and ambient temperature. That model can be used to help biologists plan for the number of researchers, and the length of time those individuals are in a cave to minimize disturbing bats.

59 BASIC BIOLOGICAL SCIENCES↗

The response of an axial compressor to unsteady disturbances

The disturbances of flow in a uniform duct are examined, taking into account a disturbance which is sinusoidal in time. The equations of motion for a plane disturbance with blade forces represented by body forces are considered along with the stage configuration of the compressor and the temperature rise for two-dimensional flow through the stage. An investigation is conducted for the case in which the disturbances are restricted to small amplitudes about the design operating condition for all points within the compressor. The three independent solutions of the disturbance equations obtained for the system of compressor and associated inlet and exit ducts can be used with three arbitrary constants to satisfy three external conditions.

Rannie, W. D.↗

Operation of a Modular 3D-Pixelated Liquid Argon Time-Projection Chamber in a Neutrino Beam

The 2x2 Demonstrator, a prototype for the Deep Underground Neutrino Experiment (DUNE) liquid argon (LAr) Near Detector, was exposed to the Neutrinos from the Main Injector (NuMI) neutrino beam at Fermi National Accelerator Laboratory (Fermilab). This detector is a prototype of a new modular design for a liquid argon time-projection chamber (LArTPC), comprising a two-by-two array of four modules, each further segmented into two optically isolated LArTPCs. The 2x2 Demonstrator features a number of pioneering technologies, including a low-profile resistive field shell to establish drift fields, native 3D ionization pixelated imaging, and a high-coverage dielectric light readout system. The 2.4-tonne active mass detector is flanked upstream and downstream by supplemental solid-scintillator tracking planes, repurposed from the MINERvA experiment, which track ionizing particles exiting the argon volume. The antineutrino beam data collected by the detector over a 4.5 day period in 2024 include over 30,000 neutrino interactions in the LAr active volume—the first neutrino interactions reported by a DUNE detector prototype. During its physics-quality run, the 2x2 Demonstrator operated at a nominal drift field of 500 V/cm and maintained good LAr purity, with a stable electron lifetime of approximately 1.25 ms. This paper describes the detector and supporting systems, summarizes the installation and commissioning, and presents the initial validation of collected NuMI beam and off-beam self-triggers. In addition, it highlights observed interactions in the detector volume, including candidate muon antineutrino events.

47 OTHER INSTRUMENTATION↗

Mechanism of Supercooled Water Droplet Breakup near the Leading Edge of an Airfoil

This work presents the results of an experimental study on supercooled droplet deformation and breakup near the leading edge of an airfoil. The results are compared to prior room temperature droplet deformation results to explore the effects of droplet supercooling. The experiments were conducted in the Adverse Environment Rotor Test Stand (AERTS) at The Pennsylvania State University. An airfoil model placed at the end of the rotor blades mounted onto the hub in the AERTS chamber was moved at speeds ranging between 50 and 80 m/sec. The temperature of the chamber was set at -20˚C. A monotonic droplet generator was used to produce droplets that fell from above, perpendicular to the path of the airfoil. The supercooled state of the droplets was determined by measurement of the temperature of the drops at various locations below the droplet generator exit. A temperature prediction code was also used to estimate the temperature of the droplets based on vertical velocity and the distance traveled by droplets from the droplet generator to the airfoil stagnation line. High speed imaging was employed to observe the interaction between the droplets and the airfoil. The high speed imaging provided droplet deformation information as the droplet approached the airfoil near the stagnation line. A tracking software program was used to measure the horizontal and vertical displacement of the droplet against time. It was demonstrated that to compare the effects of water supercooling on droplet deformation, the ratio of the slip velocity and the initial droplet velocity must be equal. A case with equal slip velocity to initial velocity ratios was selected for room temperature and supercooled droplet conditions. The airfoil velocity was 60 m/s and the slip velocity for both sets of data was 40 m/s. In these cases, the deformation of the weakly supercooled and warm droplets did not present different trends. The similar behavior for both environmental conditions indicates that water supercooling has no effect on particle deformation for the limited range of the weak supercooling of the droplets tested and the selected impact velocity. The assumption of a constant surface tension value was further supported by the equal trend of the Bond number obtained for supercooled and room temperature droplets.

in-flight icing↗

Performance of a short annular dump diffuser using wall trailing-edge suction

A short annular dump (abrupt flow area change) diffuser was tested with suction through wall trailing-edge slots at inlet Mach numbers of 0.19 and 0.27 and at near ambient inlet temperature and pressure, with suction flow varied from zero to 10 percent of the inlet air mass-flow rate. The overall ratio of diffuser exit area to inlet area was 4.0, and the ratio of length to inlet height was 2.0. By applying suction flow separately on either wall or to both walls simultaneously, the original annular jet profile could be altered to either a hub- or tip-biased profile. Diffuser effectiveness was increased from about 25 percent with no suction to 50 percent at 6 percent outer-wall suction and to 52 percent at a combined suction rate on both walls of 10.25 percent. At the same time, diffuser total pressure loss was reduced by one-fourth.

Juhasz, A. J.↗

Inside out: Measuring the effect of wood anatomy on the efflux and assimilation of xylem–transported CO 2

Carbon dioxide concentrations (including aqueous CO 2 , carbonic acid, bicarbonate, and carbonate) in woody tissues can be up to 750 times higher than atmospheric CO 2 concentrations, ranging from <1% to >26% vs atmospheric ~0.04% (Fig. 1) (Teskey et al. 2008). CO 2 formed through respiration is generally assumed to diffuse to the atmosphere from tissues adjacent to where it is produced. Here, this CO 2 buildups in the stem due to the diffusional barriers in woody and bark tissues. CO 2 in the stem has three fates: 1) it can be refixed for photosynthesis, 2) it can be used for anaplerotic reactions, or 3) it can exit the plant either adjacent to where it was produced (radial diffusion) or in an area remote from its point of origin (xylem-transported CO 2 ). Not accounting for assimilation of xylem-transported CO 2 may result in underestimating total plant photosynthesis. Alternatively, overlooking the transport of CO 2 away from its point of origin complicates the estimation of respiration in stems, branches, or even leaves. For example, CO 2 efflux from the stem may not solely represent stem respiration, but it may represent CO 2 generated through respiration in areas remote from the point of efflux (Stutz et al., 2017). Additionally, wood anatomy and branching architecture likely influence how much and how far xylem-transported CO 2 travels (Fig. 2). Few studies have combined measurements of both the efflux of xylem-transported CO 2 from trees together and the assimilation of xylem-transported CO 2 . Excitingly in this issue of Plant, Cell & Environment, Salomón et al. (pp. ) demonstrate the importance of xylem-transported CO 2 in plants with different wood anatomies by comparing the amount of xylem-transported CO 2 used for photosynthesis to the amount effluxed.

54 ENVIRONMENTAL SCIENCES↗

Tropical and extratropical forcing of the large-scale circulation - A diagnostic study

Seasonal anomalies in the large-scale circulation are studied in relation to seasonal-mean tropical forcing and extratropical transient-eddy forcing using nine years of tropical-wind and outgoing-longwave-radiation (OLR) data. Results show that in the interannual time scale, extreme swings in anomalous convection over the maritime continent of Indonesia/Borneo and the equatorial central Pacific are linked, respectively, by an east-west seesaw or a dipole pattern, identified with major fluctuations of the Pacific Walker circulation. It is found that extratropical circulation anomalies are related to the fluctuation of this dipole heat source. The overall extratropical circulation seems to be more responsive to the dipole heating with enhanced convection over the maritime continent. While the characteristics of the transient eddies, especially near the exit region of the East Asia jet stream, change substantially between maximum or minimum convection episodes, these transient eddies (defined as transients from the seasonal mean) appear to be unimportant in forcing the seasonal mean circulation anomalies.

Lau, K. M.↗

Three-dimensional calculations for underexpanded and overexpanded supersonic jet flows

A three-dimensional code is developed to simulate the characteristics of the jet exhaust plume issuing from nonaxisymmetric nozzles at different flight conditions. In the present paper, the time-iterative Space Marching Scheme (SMS) is used to simulate the shock-cell structure of underexpanded supersonic jet issued from round and square nozzles. Also, the computational results are compared with the experimental data from underexpanded and overexpanded supersonic round jets. The agreement with the experimental data validates the applicability of the new technique for predicting the shock-cell structure of supersonic jet flows. The computational results reveals that the azimuthal structure of the supersonic square jet changed as the jet progressed downstream of the exit nozzle.

Abdol-Hamid, Khaled S.↗

Integrated Display System for Low Visibility Landing and Surface Operations

This report summarizes the software products and system architectures developed by Lockheed Martin in support of the Low Visibility Landing and Surface Operations (LVLASO) program at NASA Langley Research Center. It presents an overview of the technical aspects, capabilities, and system integration issues associated with an integrated display system (IDS) that collects, processes and presents information to an aircraft flight crew during all phases of landing, roll-out, turn-off, inbound taxi, outbound taxi and takeoff. Communications hardware, drivers, and software provide continuous real-time data at varying rates and from many different sources to the display programs for presentation on a head-down display (HDD) and/or a head-up display (HUD). An electronic moving map of the airport surface is implemented on the HDD which includes the taxi route assigned by air traffic control, a text messaging system, and surface traffic and runway status information. Typical HUD symbology for navigation and control of the aircraft is augmented to provide aircraft deceleration guidance after touchdown to a pilot selected exit and taxi guidance along the route assigned by ATC. HUD displays include scene-linked symbolic runways, runway exits and taxiways that are conformal with the actual locations on the airport surface. Display formats, system architectures, and the various IDS programs are discussed.

Beskenis, Sharon Otero↗

Ground impingement of supersonic jets from nozzles with various exit geometries

Aerodynamc and acoustic measurements of single and twin free and impinging jets were made for circular convergent and CD nozzles, rectangular CD nozzles, and D-shaped convergent nozzles. Measured shock cell lengths and screech tone frequencies compare well with existing theory. Noise due to jet impingement correlates on the basis of the nozzle-to-ground distance relative to the free jet shock structure, with the highest sound pressure levels generally due to jet impingement tones. All nozzle types experienced near-field OASPL of 160 dB at some nozzle-to-ground distances, including the CD nozzles operating at their design pressure. The effect of the second jet was to sometimes increase levels due to twin jet resonance and at other times reduce levels due either to suppression of the screech/impingement feedback process or more simply by acoustic interference of the noise produced by the two jets.

Norum, Thomas D.↗

A numerical investigation of a subsonic jet in a crossflow

The flowfield induced by a single, subsonic jet exhausting perpendicularly from a flat plate into a subsonic crossflow has been numerically investigated. The test case was chosen to match available experimental data where the jet Mach number was 0.78, and the freestream Mach number was 0.13. Time-averaged solutions were obtained using the thin-layer Navier-Stokes equations and two overlapping grids. The solutions were sensitive to the radial grid clustering near the edge of the jet and to the far-field boundary conditions. Experimental data comparisons were required to determine the most appropriate jet grid and satisfactory boundary conditions. Globally, the solutions converged in about 6000 iterations. The computational results accurately showed the deflected jet and associated contrarotating vortices. The fine clustered grid in the region upstream of the jet exit allowed the horseshoe vortex in the boundary layer near the jet exit to be captured. Most importantly for aircraft applications, the computed plate pressure distributions compared favorably with the experimental data over most of the surface. However, in the wake region immediately downstream of the jet exit, where there is extensive flow separation, some discrepancies with experimental data were observed. Two turbulence models were used in this study: (1) the zero-equation, two layer Baldwin-Lomax turbulence model; and (2) one-equation Baldwin-Barth turbulence model. The turbulence models gave results which generally compared no better with experimental data than the laminar computation results.

Chiu, Stephen H.↗

Sensitivity of the Tropical Atmosphere Energy Balance to ENSO-Related SST Changes: How Well Can We Quantify Hydrologic and Radiative Responses?

The continuing debate over feedback mechanisms governing tropical sea surface temperatures (SSTs) and tropical climate in general has highlighted the diversity of potential checks and balances within the climate system. Competing feedbacks due to changes in surface evaporation, water vapor, and cloud long- and shortwave radiative properties each may serve critical roles in stabilizing or destabilizing the climate system. It is also intriguing that even those climate variations having origins internal to the climate system-- changes in ocean heat transport for example, apparently require complementary equilibrating effects by changes in atmospheric energy fluxes. Perhaps the best observational evidence of this is the relatively invariant nature of tropically averaged net radiation exiting the top-of-atmosphere (TOA) as measured by broadband satellite sensors over the past two decades. Thus, analyzing how these feedback mechanisms are operating within the context of current interannual variability may offer considerable insight for anticipating future climate change. In this paper we focus on how fresh water and radiative fluxes over the tropical oceans change during ENSO warm and cold events and how these changes affect the tropical energy balance. At present, ENSO remains the most prominent known mode of natural variability at interannual time scales. Although great advances have been made in understanding this phenomenon and realizing prediction skill over the past decade, our ability to document the coupled water and energy changes observationally and to represent them in climate models seems far from settled (Soden, 2000 J Climate). Our analysis makes use a number of data bases, principally those derived from space-based measurements, to explore systematic changes in rainfall, evaporation, and surface and top-of-atmosphere (TOA) radiative fluxes, A reexamination of the Langley 8-Year Surface Radiation Budget data set reveals errors in the surface longwave emission due to use of biased SSTs. Subsequent correction allows subsequent use of this data set along with ERBE TOA fluxes to infer net atmospheric radiative beating. Further analysis of recent rainfall algorithms provides new estimates for precipitation variability in line with interannual evaporation changes inferred from the da Silva, Young, Levitus COADS analysis. The overall results from our analysis suggest an increase (decrease) of the hydrologic cycle during ENSO warm (cold) events at the rate of about 5 Wm-2 per K of SST change. This rate is slightly less than that which would be expected for constant relative humidity over the tropical oceans. Corresponding radiative fluxes seem systematically smaller resulting in a enhanced (suppressed) export of energy from the tropical ocean regions during warm (cold) SST events. Discussion of likely errors due to sampling and measurement strategies are discussed along with their impacts on our conclusions.

Robertson, Franklin R.↗

Data Communications Performance of AOCDL and AUTOMET Over a VDL Mode 2 Link

Airline Operational Control Data Link (AOCDL) is a service that allows for downlink of in-flight data to help an airline operate efficiently. Automated Meteorological Transmission (AUTOMET) is a weather reporting system that provides forecasts and atmospheric conditions to pilots during flight. These are two examples of new or extended services that will be enabled by implementation of the ATN and its digital air/ground communications data links. One of the possible data links identified for these services is VDL Mode 2. This research focuses on the data link performance of AOCDL and AUTOMET services in a 2015 scenario using VDL Mode 2 as the data link. In the near term, it is expected that the Controller-Pilot Data Link Communication (CPDLC) and Decision Support System Data Link (DSSDL) services will also use the VDL Mode 2 link. CPDLC and DSSDL are high-priority services that communicate aircraft control information and require low end-to-end delays for safety. VDL Mode 3, which is the preferred link for these services because it provides priority-based link access, will not be deployed in time for earlier implementations of these services to use. As such, the research also investigates the impact of CPDLC and DSSDL services on the same VDL link as AOCDL and AUTOMET. Performance of these services over the data link will be separately examined for the airport, terminal, and en route domains. OPNET simulations are performed for each domain. The number of aircraft for each domain is representative of the expected peak aircraft in 2015. Results are collected during a steady-state condition, with aircraft neither entering nor exiting the domain. The study explores the services in terms of the traffic load they present to the VDL Mode 2 radio. The amount of load each service contributes to the data link will be shown for each domain. The end-to-end delays for each service will be discussed. Data link metrics such as subnetwork and queuing delays, channel utilization, and protocol overhead will be used to show the impact of these services on VDL Mode 2. All results will be presented both with and without the inclusion of CPDLC and DSSDL services.

Bretmersky, Steven↗

Aerojet - AJ10-137 Apollo Service Module Engine

The general configuration of the SPS engine was 20,000 pounds of thrust, with a chamber pressure of 100 psi and specific impulse (Isp) of 314.5. The very large nozzle had an area ratio of 62.5:1 (exit area to throat area). The propellants were nitrogen tetroxide (also known as N2O4 and nitrous oxide) and A-50. A-50 was a hydrazine family fuel. Aerojet developed it for the Titan Missile Program when they went with Titan II, to store it in the launch silos. They wanted the highest performance they could get. N2H4 was just pure hydrazine, which doesn't take low temperature very well. In fact, it freezes about like water. We started adding unsymmetrical-dimethylhydrazine (UDMH) to the hydrazine until such time as it would meet the environmental specifications the Air Force needed for Titan II. It turned out it s roughly a fifty-fifty mix. We still had to be careful with that fuel because the two fluids didn't mix very well chemically. We had to spray the two fluids through some special nozzles to get them to emulsify with each other into a single fluid. If we ever got it too cold or froze it, the hydrazine separated back out. Then, if we tried to run the engine, things could go boom in the night. The inlet pressure was only 165 pounds per square inch absolute (psia), but we needed at least forty psi pressure drop across the injector just to get some kind of stable flow. It was a whole new game for some of us. We didn't have much supply pressure to work with. It had the aluminum injector to keep the weight down. That was a couple feet in diameter, and we didn't have a lot of propellant to cool it. In fact, we had to use both propellants to keep the injector cool. There were twenty-two ring channels in the injector. Specification required 750 seconds duration, or fifty engine restarts during a flight. There were several first flight things we accomplished with the engine. It was the first ablative thrust chamber of any size to fly. (See Slide 6, Appendix G) There were no liners in it. It was just straight ablative material. It took us a while to figure that out. It was a throat-gimbaled engine, and it was the first engine to fly with columbium (also known as niobium, used as an alloying element in steels and superalloys) in the nozzle.

Boyce, Clay↗

Particle Scale Impacts on Deconstruction Energy of Pine Residues

The goal of this Case Study was to quantify the impacts of variable moisture and ash on hammer mill throughput and energy consumption and on generation of fines that are not able to be fed to conversion, as compared to a status quo Base Case system. Also considered was convertible carbon content (minimum carbon specification) and maximum ash content and the delivered feedstock cost impacts of not being able to feed residue not meeting both specifications to the conversion reactor. Laboratory data on the impacts of input particle size and moisture content on the exit particle size were received from FCIC Subtask 5.2: Preprocessing, High Temperature Conversion Preprocessing from their single particle impact population balance modeling study (Tiasha Bhattacharjee, INL). Additional throughput and energy consumption data were obtained from FCIC Subtask 5.2 (Jordan Klinger, INL) for the same grinder with a 6 mm screen in place. These data were utilized to develop the necessary response surface equations to perform throughput analysis using discrete event simulation. Because the ash contents in the separated fines had not been analyzed in the laboratory at the time of the model runs, we chose to assume that the ash distributed proportionally with total mass into the overs and unders in the disk screen following grinding.

energy consumption↗

A reduced-order modeling of a tubular solar reactor for long duration thermochemical energy storage

The storage of solar energy in a solid form, referred to as a “solar fuel”, can be achieved through a process known as endothermic solar thermochemistry. This process transforms the absorbed solar energy into a stable and retrievable form that can be stored for extended periods of time. This paper presents a low–order heat transfer model of a counter–current tubular falling bed reactor designed to produce thermally reduced magnesium manganese oxide pellets for long duration thermochemical energy storage. The energy required for the endothermic reduction was supplied by concentrated solar energy or renewable electricity via indirect heating of the gas and solid reactants flowing in a ceramic tube. The counter-current gas flow enhances the mixing of the solid particles with the heat recuperation zone, allowing the gas and particles to enter and exit the tubular reactor close to room temperature. Further, the reactor was vertically oriented and was heated circumferentially by an adjustable level heat flux along a finite segment of its length. The temperature distribution of the reactor in response to transient changes along the tube was modeled by considering conduction, convection, and radiation heat transfer. Governing equations for the heat transfer model were solved by discretizing the reactor tube into a finite number of control volumes and using an energy balance for the heat exchange between the reactor wall, gas, and particles within the control volume. The energy absorbed during this endothermic reaction was modeled numerically by fitting the data of the chemical conversion rate with the corresponding temperature of particles in the heating zone. The numerical model has been experimentally validated using a reactor prototype made of a 121.92 cm alumina tube heated by a 7kW electric tube–furnace. The alumina tube receives magnesium manganese oxide pellets of 3.66±0.516 mm in diameter from the top, and a counter–current gas flow from the bottom. The reactor wall temperature was monitored by six thermocouples installed along the reactor tube length. The experimental procedure was numerically simulated, and the temperature variation along the reactor tube was compared with a matrix of experimental runs for a range of particles mass flowrates (0.75–1.25g/s) and corresponding gas flowrates (36–65 SLPM). The reactor system was heated gradually from room temperature to a steady state temperature of 1673K, and then cooled down to room temperature. The heating and cooling processes were simulated, and the numerical and experimental results were compared throughout processes. The numerical model showed similar trends to the experimental results, with an error of 0.69 to 7.9% for the particle inlet and 0.7 to 7.9% for the gas inlet during steady-state operation. The proposed numerical model can be implemented as a simplified physical model to design a feedback control system to regulate reactor temperature.

14 SOLAR ENERGY↗