Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “SST”

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,405 records · Page 78

Heated Supersonic Axisymmetric Jet Cases for the NASA Turbulence Modeling Resource

This paper will describe three new supersonic jet cases for the NASA Turbulence Modeling Resource (TMR). These cases were taken from the 6th AIAA Propulsion Aerodynamics Workshop (PAW) nozzle test problem, which utilized data taken at NASA GRC consisting of non-intrusive measurements of jet plume velocities and temperatures, including both mean values and turbulence statistics. The three jet cases all used a Mach 1.63 axisymmetric nozzle: A heated jet for each at on-design conditions, a temperature-matched on-design condition at Mach 1.63, and an off-design heated condition for the Mach 1.63 nozzle. Computational fluid dynamics (CFD) solutions are obtained using three established Reynolds-averaged Navier-Stokes (RANS) codes: Wind-US, FUN3D, and Vulcan, all using the Menter Shear Stress Transport (SST) k-ω turbulence model. The NASA TMR requires at least two CFD codes yielding essentially identical results to certify that the approaches are verified, that is - providing solutions that may not necessarily provide close comparisons with experimental data, but solve the posed CFD equations as intended.

nozzle↗

Smallsat 2024 - Starling Cubesat Swarm Technology Demonstration Flight Results

The Starling swarm of four 6U CubeSats launched in July 2023 to test four key technologies to enable future swarm missions: 1) Mobile Ad-Hoc Networking (MANET) over a crosslink radio network 2) Autonomous onboard decision-making for operations 3) Optical-based absolute and relative navigation 4) Autonomous maneuver planning and execution The Starling team implemented the Better Approach to Mobile Ad-hoc Networking (B.A.T.M.A.N.) protocol to automatically manage the crosslink network of four satellites. The B.A.T.M.A.N. protocol uses a decentralized approach to managing a multi-hop mesh network of devices, in this case, a satellite swarm. The four satellites were able to successfully establish a network at multiple data rates and demonstrate file transfer and command issuance between spacecraft over the network. Starling incorporated Distributed Spacecraft Autonomy's (DSA) software to demonstrate onboard decision-making. The DSA software takes L1/L2 band GPS measurements and uses them to estimate the relative Total Electron Count (TEC) in the ionosphere. The onboard software then determines if there are any features of interest and provides that information to the other satellites over the crosslink network. The swarm of satellites then reaches a consensus on the optimal TEC observation strategy and adjusts its measurement collection tactics autonomously. The Starling Formation-Flying Optical Experiment (StarFOX), produced by Stanford's Space Rendezvous Laboratory, uses the onboard star trackers to collect images of the other swarm spacecraft and produce angles-only navigation estimates. This system is envisioned to be valuable in applications in which Global Navigation Satellite Systems (GNSS) are not available, such as in cis-lunar or deep space. StarFOX successfully applied its algorithms to multiple simultaneous spacecraft targets using the star tracker imagery. Finally, Starling used Emergent Space's Cluster Flight Application (CFA) software suite for the Reconfiguration and Orbit Maintenance Experiments Onboard (ROMEO) demonstration of autonomously planning and executing propulsive maneuvers. Large swarms will need to be able to maintain formation requirements with minimal operator involvement, especially as the size of the swarm scales up. Results from the ROMEO experiment are presented. Starling is funded by the Small Spacecraft Technology (SST) program out of NASA's Space Technology Mission Directorate (STMD).

distributed systems↗

Comparison of OVERFLOW Computational and Experimental Results of the CobraMRV Mars Entry Vehicle Concept during Supersonic Retropropulsion

The CFD solver Overflow was used to simulate the CobraMRV undergoing supersonic retropropulsion (SRP) in the Langley Unitary Plan Wind Tunnel as part of a CFD and wind tunnel integration effort. Selection of numerical methods and resulting comparisons to selected tunnel conditions of interest are detailed. Comparisons of mean pressures with static taps and pressure-sensitive paint are provided for multiple SST-based turbulence models. URANS predicitons of heat shield pressures at high thrust conditions are poor compared to TMRC and DES predictions. Comparisons to experimental measurements of pressure fluctuations indicate complex relationships between heat shield unsteadiness, flight Mach number, and engine thrust coefficient. Bow shock topology and shock cell structure between experimental and OVERFLOW data show good agreement in Schlieren imaging. Recommendations for future investigation into flow field unsteadiness and turbulence models are suggested.

CFD↗

Novel Turbulence Model Tailored for Complex Rocket Motor Exhaust Jets

This paper describes a new RANS turbulence model intended to give improved accuracy in turbulence that is affected by compressibility, rotation and curvature, and variable density, making it well-suited to complex rocket motor exhaust jets. Four key elements of the model are: the use of simple-averaged (not mass-averaged) turbulence decomposition; a second-order Taylor series for the equilibrated turbulent fluxes; empirical modeling that accounts for the combined effects of aforementioned turbulence modifiers; and model tuning on large-eddy simulation data for mixing layers for a wide range of these turbulence modifiers. The model is implemented as an extension of the k-w SST model. Results include comparisons to detached eddy simulation results for hot and cold co-flowing exhaust jets, and to experimental data for supersonic retropropulsion on the CobraMRV Mars-lander concept vehicle. The model is named TMRC for its constituents elements, the Taylor series with Mach number, rotation and curvature effects.

turbulence modeling↗

Aeroheating Predictions for a Hypersonic, Turbulent Near-Wake

The accuracy of heating predictions using various turbulence models is examined for an axisymmetric near-wake at Mach 6. The CFD predictions are compared with experimental data collected under AGARD Working Group 18 on the wake of a 70-degree sphere-cone. The impact of grid resolution and discretization error is estimated, which allows a comparison of stacked-block and conventional structured meshes. The accuracy of steady Reynolds-averaged Navier-Stokes (RANS) models is contrasted with that of a hybrid RANS/Large-Eddy Simulation model. The predictions are made with three different CFD codes (LAURA, FUN3D, and HyperSolve), to demonstrate the code-to-code variation in the results. Steady SST models capture the qualitative nature of the heating in the wake, including the increasing peak heating with increasing Reynolds number. Spalart-Allmaras models, including SA-Catris, under-predicted the peak heating in the wake. Hybrid RANS/LES models improve upon the SA results but have their own modeling difficulties near the shear layer impingement. These results are generally consistent across solvers and grid topologies.

RANS↗

An Assessment of Turbulence Models for MSR Sample Retrieval Lander Aerothermal Design

Several cases from past entries in the T9 and T5 hypersonic test facilities are examined using various turbulence models and shock-layer thermochemistries with the objectives of verifying/validating turbulence models and determining appropriate multiplicative factors on predicted turbulent heating for use in aerothermal design of the heatshield of the Sample Retrieval Lander (SRL) of the Mars Sample Return (MSR) Flagship mission. The zero-equation turbulence models of Baldwin-Lomax and Cebeci-Smith are shown to replicate measurements to within experimental uncertainty for the T9 cases. For T5 cases, the zero-equation model of Baldwin-Lomax can also replicate measurements provided the shock-layer thermochemistry is modified. The modified thermochemistry resolves a dichotomy from past CFD that required different wall catalycity assumptions in the laminar and turbulent flow regimes. The present work shows that such assumptions are not necessary. Finally, an initial estimate for a multiplicative factor on turbulent heating is estimated for the SST turbulence model – the chosen model for SRL aerothermal design simulations.

Mars entry↗

Aeroheating Predictions for a Hypersonic, Turbulent Near-Wake

The accuracy of heating predictions using various turbulence models is examined for an axisymmetric near-wake at Mach 6. The CFD predictions are compared with experimental data collected under AGARD Working Group 18 on the wake of a 70-degree sphere-cone. The impact of grid resolution and discretization error is estimated, which allows a comparison of stacked-block and conventional structured meshes. The accuracy of steady Reynolds-averaged Navier-Stokes (RANS) models is contrasted with that of a hybrid RANS/Large-Eddy Simulation model. The predictions are made with three different CFD codes (LAURA, FUN3D, and HyperSolve), to demonstrate the code-to-code variation in the results. Steady SST models capture the qualitative nature of the heating in the wake, including the increasing peak heating with increasing Reynolds number. Spalart-Allmaras models, including SA-Catris, under-predicted the peak heating in the wake. Hybrid RANS/LES models improve upon the SA results but have their own modeling difficulties near the shear layer impingement. These results are generally consistent across solvers and grid topologies.

RANS↗

Augmenting RANS Turbulence Models Guided by Field Inversion and Machine Learning

This report investigates the use of a data-driven approach, viz., Field Inversion and Machine Learning (FIML), to improve conventional RANS turbulence models like the Spalart-Allmaras model and the Menter SST k-ω model. One of the crucial aspects of using an ML-based approach with limited training data to produce corrections that are generalizable to a large range of flow configurations is to design appropriate “features” (inputs to the ML model). A model, based on guidance from the FIML methodology, is presented in analytical form. An additional list of potential features is provided. Although these were not used in the present correction, they were considered in the course of its development, and are included to fully document the complete process employed in the present work.

turbulence modeling↗

Numerical Simulations of an Axisymmetric Shock-Wave/Boundary-Layer Interaction

A series of Reynolds-averaged Navier-Stokes (RANS) simulations were performed of two Mach 2.5 axisymmetric shock-wave/boundary-layer interactions (SWBLIs) based off of the NASA Turbulent Computational Fluid Dynamics Validation Experiment (TCFDVE). The goals of the simulations were to (1) determine how well RANS can predict the flow physics of an axisymmetric SWBLI and (2) generate a new validation case for the NASA Turbulence Modeling Resource (TMR). The simulations were run with two different turbulence models and three different computational fluid dynamics (CFD) codes while looking at two different conical shock generator angles, one of which was shown to cause flow separation in the experiment. The results in this paper show that, as expected, the CFD simulations agreed well with each other and the experiment when SWBLI-induced separation was not present. When separation was present, differences in the predicted separation characteristics were noted. These included a stronger interaction and subsequently longer separation extent for the Menter Shear-Stress Transport k-ω (SST) turbulence model simulations compared to both the Spalart-Allmaras (SA) turbulence model simulations and the experiment.

Computational Fluid Dynamics↗

State-of-the-Art Small Spacecraft Technology

The objective of this report is to assess and provide an overview of the state of the art in small spacecraft technologies for use by mission designers, project managers, technologists, and students, connecting current small spacecraft missions to available technologies. This report focuses on the spacecraft system in its entirety, provides current best practices for integration and key considerations for the reader where possible, and presents devices from publicly available sources for each specific spacecraft subsystem. This report is a survey of small spacecraft technologies sourced from open literature; it is not an original source. In addition, this report only considers literature in the public domain. Information presented in this report is limited to SmallSat technology that is publicly available as of September 30, 2024. It does not include information on instrumentation, science payloads, or advances or developments that have not been publicly disclosed. Commonly used sources for data include manufacturer data sheets, press releases, conference papers, journal papers, public filings with government agencies, news articles, presentations, the compendium of databases accessed via NASA’s Small Spacecraft Systems Virtual Institute (S3VI) Information Search, and engagement with companies. Data not appropriate for public dissemination, such as proprietary, export controlled, or otherwise restricted data, are not considered. We encourage the SmallSat community to publish mission outcomes and technology development milestones in publicly available conference papers, press releases, or company websites so they can be reflected in future editions of this report. This report is funded by NASA’s Space Technology Mission Directorate (STMD). It was first commissioned by the Small Spacecraft Technology (SST) program within NASA’s STMD in mid2012 in response to the rapid growth in interest in using small spacecraft for low-Earth orbit, lowcost missions. The report was subsequently updated in 2015, 2018, 2020, 2021, 2022, and 2023 to capture SmallSat technology growth and maturation. In addition to reporting currently available state-of-the-art technologies that have achieved Technology Readiness Level (TRL) 5 or above, a prognosis is provided describing technologies as "on the horizon" if they are being considered for future application.

Sasha V Weston↗

A Smart Silicon Carbide Power Module With Pulse Width Modulation Over Wi-Fi and Wireless Power Transfer-Enabled Gate Driver, Featuring Onboard State of Health Estimator and High-Voltage Scaling Capabilities: Preprint

A wide range of utility applications require controllable switches with features such as high-voltage (HV) blocking and high-current carrying capacity especially at high pulse width modulation (PWM) frequency. Low and medium voltage utility applications such as motor drives and flexible AC transmission systems (FACTS) as well as solid state transformers (SST) could also benefit from a low-cost HV switching module. Wide band gap (WBG) semiconductors such as SiC and GaN MOSFETs are considered to be the present and next generation device choices, although they have their own limitations. For relatively HV applications with demanding thermal management, SiC is still the only choice, and GaN dominates the low voltage regime. This manuscript proposes a new half bridge power MOSFET module that is suitable for conventional H-bridge of multilevel configurations used in HV applications. Constructed from bare SiC dies, this half bridge module takes advantage of (1) optimized MOSFET placement inside the module, (2) customized heat exchanger, manifold and cooling, (3) integrated gate driver module with PWM over wi-fi to eliminate the need for low voltage signals, (4) wireless power transfer (WPT) enabled gate driver and other ancillary circuits, (5) and the option to incorporate an onboard state of health (SOH) estimator module onboard. The entire architecture has been designed and built at National Renewable Energy Laboratory (NREL), Golden, CO.

baseplate design↗

Accelerating Hanford Site Cleanup through Operations Research Modeling - 20238

The Hanford Site cleanup effort will require the integration of dozens of unique facilities and processes, many of which will be first-of-a-kind in implementation and design. Each facility will be governed by its own set of operating logic, configured with a unique array of unit operations, and subject to a set of constraints that will affect its behavior. The collection of facilities have multiple points of interface, making the operations of any one facility potentially significant to the operations of other up- or downstream processes. It is therefore highly desirable to accurately predict these operations, as it allows for Site officials to identify and preempt bottlenecks and vulnerabilities before they unexpectedly inhibit the cleanup mission. With the quantity and complexity of the processes that will be on Site, building a pen-and-paper or even a spreadsheet-assisted model of the cleanup mission quickly becomes overwhelming in scope and inaccurate in execution. The Engineering organization for the Site's Tank Operations Contract (TOC) has therefore implemented the use of operations research (OR) modeling to simulate and predict future operations of Site facilities. These models are created using a discrete event simulation tool that allows for the development of detailed, versatile, and robust models. Not only can these models account for complex logical behaviors, but they can also simulate process details down to the level of vessel sizing, labor utilization, equipment reliability, and resource availability. To date, the TOC has developed OR models for several facilities on Site, including for single-shell tank (SST) farms, double-shell tank (DST) farms, the Effluent Treatment Facility (ETF), and the waste transfer system. These models have focused on identifying bottlenecks and operational constraints, and have been used to quantify the effects of implementing process changes. This latter point is particularly valuable, as it allows for several alternatives to be studied in a virtual setting before committing resources to making a change in the field. The decision to develop OR models has gained tremendous support from the Site's stakeholders and the U.S. Department of Energy (DOE) management, and has prompted the use of the tool to support additional internal and external initiatives. Recently, an initiative was proposed to use the models to help identify and provide quantitative backing for risks and opportunities for the TOC. This application of OR could not only help inform how the TOC manages its risks (e.g. quantities and types of spare parts), but could also help drive process improvements whose benefits might otherwise be hard to quantify. The models have also been used to drive the TOC's cloud computing, artificial intelligence (AI), and machine learning (ML) initiatives. These initiatives will not only improve the ability of the TOC to more rapidly respond to the needs of its customers, but it will also aid in the ability of the TOC to analyze and improve the processes it studies. Partnership with two external software development and consulting companies (Lanner and Ynformed) has furthered not only the application of AI and ML within the TOC, but has also spurred the development of new/improved software tools and platforms used by the companies. These partnerships have proven to be mutually beneficial and productive, and have set a precedent for the types of gains that can be made by exploring such options. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Modeling Efforts to Gain Insight into Historical Leak Events from the Single-Shell Tank A-105 in the Hanford 241-A Tank Farm - 20113

The 241-A Tank Farm is a single-shell tank (SST) farm constructed to store process waste from Hanford nuclear operations. Millions of gallons of nuclear waste were stored in the 241-A Tank Farm SSTs and some of the SSTs leaked in the past. In addition, spills and pipeline leaks during transfers and storage and intentional discharges to cribs and trenches resulted in releasing waste to the ground. Tank liner leaks are referred to as 'leaks' and all other discharges to the soil are referred to as 'releases.' Liquid waste that could be removed by pumping has been removed from all of the SSTs to reduce the potential for future leaks. The most significant historical leak event in 241-A Tank Farm occurred in 1965 when tank A-105 experienced a pressurized steam event. This event resulted in damage to the inner steel liner flooring of A-105, which was separated from the sidewalls over part of its circumference and buckled up to 2.5 m, vertically. During this event up to 7,570 L of contaminated liquid may have leaked to the adjacent soil. This leak occurred under extreme temperature and pressure conditions. Modeling efforts to better understand the nature and extent of historical leaks from subsurface SSTs that were used to store highly radioactive, self-boiling liquid wastes at 241-A Tank Farm have recently been undertaken. These efforts have compiled the available historical information to support a preliminary non-isothermal, multiphase flow and chemical transport modeling effort to re-create the conditions under which the leaks occurred and to formulate a conceptual model as to the extent and distribution of leaked radioactive contaminants in the adjacent soil material. The focus of this work was on developing an understanding of the key features, processes and bounding conditions related to tank A-105 leak events that occurred in the 1960's. The activities include estimation of leak composition, non-isothermal multiphase flow and transport modeling, and geochemical modeling. The STOMP{sup C} Water-Air-Energy modeling code was utilized to implement the three-dimensional representation of the subsurface tank and the surrounding flow field. Once constructed, the model was used to evaluate the conceptual understanding of leaks originating from different parts of the tank under the elevated temperature boundary conditions that were imposed by surrounding tanks at the time of the historical leaks. The chemical evolution of the liquid waste was also evaluated as it leaked from the tank under transient pressure and temperature gradients. The results of the preliminary evaluation suggested the presence of a heat-pipe effect beneath tank A-105 in which water vapor at an elevated temperature is driven away from the base of the tank to a position where water vapor cools, condenses and is then drawn back toward the tank by the strong capillary attraction of the dry soil. A non-sorbing contaminant (Tc-99) was introduced into this flow field to better understand the potential distribution patterns of leaked contaminants that may have occurred during the historical leak events. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Recent Improvements in Pronghorn for Advanced Reactor Modeling

Pronghorn is a thermal-hydraulics computational tool developed using the Idaho National Laboratory's Multiphysics Object-Oriented Simulation Environment (MOOSE). It is designed to support Computational Fluid Dynamics (CFD) modeling, ranging from subchannel and porous media analysis to Reynolds Averaged Navier-Stokes (RANS) turbulence modeling. As an integral part of the MOOSE-based suite of tools, Pronghorn seamlessly couples with other MOOSE-based applications to simulate a variety of physical phenomena. This article highlights recent significant enhancements to Pronghorn's CFD modeling capabilities and demonstrates their application to advanced nuclear reactor designs. The recent improvements in Pronghorn primarily focus on modifications to its turbulence modeling capabilities, near-wall corrections and numerical schemes. In terms of turbulence modeling, the two-equation $k-\epsilon$ and $k-\omega$ SST models have been implemented and validated with both equilibrium and non-equilibrium wall treatments. Additionally, corrections for wall roughness, and curvature, and wall-channeling in pebble beds have been introduced in the near-wall modeling. These developments enable more accurate simulations of advanced nuclear reactors. Two case studies are presented in this work: a pool-type Molten Chloride Reactor and a salt-cooled Pebble-Bed High Temperature Reactor. In both cases, the previous models in Pronghorn are compared with the new implementations, demonstrating the improved accuracy achieved with the updated models.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗

Submodule-Level Decentralized Control for Cascaded Quadruple Active Bridge-Based Medium-Voltage Converters: Preprint

High-power, low-voltage DC (LVDC) loads, such as AI data centers, electric vehicle fast-charging hubs, and hydrogen electrolyzers, are accelerating the deployment of solid-state transformers (SSTs) that directly interface medium-voltage AC (MVAC) to LVDC, streamlining power delivery. This paper proposes a submodule-level decentralized control strategy for an SST comprising a cascaded H-bridge (CHB) active rectifier on the MVAC side and per-module quadruple-active-bridge (QAB) converters on the LVDC side with local controllers. The proposed strategy eliminates the need for high-fidelity, real-time cross-isolation communication by synchronizing all active bridges to a common heartbeat signal and enabling the primary- and secondary-side controllers to generate independent local control signals (phase shifts). The resultant phase shifts naturally govern average and double-line frequency power transfer between active bridges, enabling each side to meet its control objectives using only local measurements. Detailed controller design considerations are presented based on a stability analysis. Comprehensive electromagnetic-transient simulations demonstrate tight DC-link voltage regulation, high-quality grid currents, and robust operation under load transients and grid disturbances, without real-time communication among the submodule controllers. These results indicate that the proposed scheme can significantly enhance the reliability, scalability, and flexibility of cascaded bridge-based medium-voltage converters.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Revisiting ENSO and IOD Contributions to Australian Precipitation

Tropical modes of variability, such as El Niño–Southern Oscillation (ENSO) and the Indian Ocean Dipole (IOD), exert a strong influence on the interannual variability of Australian precipitation. Nevertheless, commonly used indices of ENSO and IOD variability display significant co-variability that prevents a robust quantification of the independent contribution of each mode to precipitation anomalies. This co-variability issue is often addressed by statistically removing ENSO or IOD variability from the precipitation field before calculating teleconnection patterns. However, by performing a suite of coupled and uncoupled modeling experiments in which either ENSO or IOD variability is physically removed, we show that ENSO-only-driven precipitation patterns computed by statistically removing the IOD influence significantly underestimate the impact of ENSO on Australian precipitation variability. Here, inspired by this, we propose a conceptual model that allows one to effectively separate the contribution of each mode to Australian precipitation variability.

54 ENVIRONMENTAL SCIENCES↗

Multiscale System Modeling of Single-Event-Induced Faults in Advanced Node Processors

Integration-technology feature shrink increases computing-system susceptibility to single-event effects (SEE). While modeling SEE faults will be critical, an integrated processor’s scope makes physically correct modeling computationally intractable. Without useful models, presilicon evaluation of fault-tolerance approaches becomes impossible. To incorporate accurate transistor-level effects at a system scope, we present a multiscale simulation framework. Charge collection at the 1) device level determines 2) circuit-level transient duration and state-upset likelihood. Circuit effects, in turn, impact 3) register-transfer-level architecture-state corruption visible at 4) the system level. Furthermore, the physically accurate effects of SEEs in large-scale systems, executed on a high-performance computing (HPC) simulator, could be used to drive cross-layer radiation hardening by design. We demonstrate the capabilities of this model with two case studies. First, we determine a D flip-flop’s sensitivity at the transistor level on 14-nm FinFet technology, validating the model against published cross sections. Second, we track and estimate faults in a microprocessor without interlocked pipelined stages (MIPS) processor for Adams 90% worst case environment in an isotropic space environment.

42 ENGINEERING↗

DC-Link Current Minimization Control for Current Source Converter-Based Solid-State Transformer

This article proposes a fast predictive control method and a small DC-link inductor to minimize the DC-link current in current-source converter (CSC)-based solid-state transformer. The DC-link current minimization can significantly reduce power loss and improve efficiency. The challenge of this problem is on improving both steady-state and dynamic performance. PI control methods and large DC-link inductors are conventionally used in the CSC but have limited dynamic performance. A model predictive control (MPC) method is proposed to achieve switching-cycle-level settling time, and the DC-link inductor is sized for 40% ripple to enable fast current change. Importantly, this article also proposes to minimize the DC-link current by varying the current even within a line cycle under single-phase load to improve the steady-state performance, in contrast with the reduction to a constant value in the literature. The proposed MPC features a constant switching frequency without weighting factors. The MPC does not have a high computational burden and is implemented in a regular digital controller for a prototype of soft-switching solid-state transformer (S4T) with reduced conduction loss. The effectiveness of the proposed method has been experimentally verified on the SiC S4T prototype during steady-state and dynamics under different multiport power flow conditions up to 2 kV peak. Here, the DC-link current in the experiments is close to the minimum current with a short zero-vector duration, which further verifies the performance of the proposed method.

14 SOLAR ENERGY↗