Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “space operations”

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 181 records · Page 10

Dualities among massive, partially massless and shift symmetric fields on (A)dS

We catalog all the electromagnetic-like dualities that exist between free dynamical bosonic fields of arbitrary symmetry type and mass on (anti-) de Sitter space in all dimensions, including dualities among the partially massless and shift symmetric fields. This generalizes to all these field types the well known fact that a massless p-form is dual to a massless (D − p − 2)-form in D spacetime dimensions. In the process, we describe the structure of the Weyl modules (the spaces of local operators linear in the fields and their derivative relations) for all the massive, partially massless and shift symmetric fields.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Effects of 18 O Exchange on Neutron Emission Rates of Aging 238 PuO 2 LWRHUs

The Lightweight Radioisotope Heater Unit (LWRHU) provides about 1 Watt of thermal power to critical electrical and mechanical systems in the extreme temperatures of space, maintaining the operating temperature of the spacecraft’s scientific and general operating equipment. Los Alamos National Laboratory (LANL) designed and tested the LWRHU in the 1970s through early 1980s. The first production campaign at LANL was in the early 1980s in support of the Galileo mission in which over 100 LWRHUs were manufactured. In the mid-1990s, almost 200 LWRHUs were manufactured in support of the Cassini mission. After the primary Cassini manufacturing campaign, the Department of Energy (DOE) approved the fabrication of more than 10 additional LWRHUs as spares in 1998. The LWRHU assembly contains the following components: a fuel pellet, a vented capsule, a pyrolytic graphite insulator, and a fine-weave pierced fabric (FWPF) graphite aeroshell. The fuel pellet is a hot-pressed 238 PuO 2 cylindrical pellet that is sintered to create a ceramic pellet. The pellet is encapsulated in a platinum-rhodium encapsulation. The welded fueled clad is placed inside a pyrolytic graphite thermal insulation to protect the fuel from the heat produced in reentry events. The exterior layer of the LWRHU assembly is the FWPF aeroshell. Before sintering, PuO 2 granules are heated in a furnace at temperatures greater than or equal to 700°C under an enriched 16 O atmosphere in order to reduce trace quantities of 18 O. Because 18 O undergoes an (α,n) reaction, 18 O needs to be removed as much as possible to meet radiation emission requirements set by the Jet Propulsion Laboratory (JPL) and the National Aeronautics and Space Administration (NASA), as well as other agencies. The relative abundances of the oxygen isotopes are indicated by neutron emission rate (NER) measurements. LWRHUs are vented with a sintered platinum frit vent (George 1986). During extraterrestrial operation, the vent in the cladding allows helium from the decay of plutonium to escape and relieves pressure from the system while preventing the solids from escaping (Rinehart 1996, Tate 1982 & 1985). This vent is protected during manufacture by a platinum-30 rhodium cover. The capsule vent is activated by milling a 0.025-inch diameter hole through the protective cover to a nominal depth of 0.015 inches, generally immediately before the LWRHU is loaded into the aeroshell (George 1986). Prior to long-term storage, LWRHUs are vented and thus potentially allow infiltration of 18 O from the ambient atmosphere. Although NER decreases with age due to the reduction in emitted alpha particles, it might be expected that NER reduction is mitigated by the exchange of 18 O. However, this is not necessarily the case and no trend was observed relating NER to age other than the natural decay of 238 Pu (Mulford 2021). The most probable and dominant factor for determining specific NER was manufacturing discrepancies between individual units. It is important to understand the lack of observed 16 O- 18 O exchange and to compare LWRHU pellets as much as possible.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Development and field demonstration of residential air source integrated heat pump using a three-stage compressor

To promote decarbonization and all electrification at residential sectors, it is necessary to use air source heat pumps (ASHPs) to replace natural gas for space heating and water heating. ASHPs are widely utilized for residential space cooling, heating, and water heating due to their simplicity and cost-effectiveness. However, their performance can be compromised in cold climates, where they may experience reduced heating capacity. A multi-functional heat pump, using a single compressor, to meet all home space conditioning and water heating demands, is an emerging technology. To address this limitation, we have developed and demonstrated an air source integrated heat pump to fulfill comprehensive home comfort requirements. This system employs a three-stage compressor and a single set of heat exchangers and valves, optimizing functionality while minimizing costs. The performance of the developed system was rigorously evaluated in both laboratory and field settings. In laboratory conditions, the system achieved a Seasonal Energy Efficiency Ratio of 17.0 (average COP of 4.98) and a Heating Seasonal Performance Factor of 11.0 (3.22). Additionally, in its most efficient operational mode—combining space cooling and water heating—the unit attained a total energy efficiency exceeding 7.0 seasonal COP in the field and could heat a 189-liter tank of water in just 25 min. The field study corroborated the laboratory findings, validating the system’s performance in real-world conditions. Here, this integrated heat pump represents an ideal solution for decarbonizing homes in northern climates by providing efficient space heating and water heating, thereby replacing the need for natural gas.

Integrated heat pump↗

Geothermal Heat Pump System Showcase: Short-Term Validation of Borehole Heat Exchanger Performance from Field Data to Numerical Modeling: Preprint

Since 2011, a geothermal heat pump (GHP) system has been operating to provide space heating and cooling for the Solar Radiation and Research Laboratory building at the National Laboratory of the Rockies (NLR) in Golden, Colorado. The system consists of 23 vertical boreholes, each extending to a depth of 300 ft (91 m), connected to 11 water-to-air heat pump units and four circulation pumps. Between fiscal years 2023 and 2025, additional power meters and temperature sensors were retrofitted to support detailed system performance assessment and model development. This study presents preliminary monitoring results and the development of an initial numerical model of the borehole heat exchanger field. The model incorporated site-specific geometry, ground thermal properties derived from thermal response tests, and ambient temperatures, and simulated system behavior over a representative operating day in September. Model predictions of outlet temperatures were compared against corresponding field measurements. Results showed that modeling initialized with a simplified linear subsurface temperature gradient presents systematic discrepancies in outlet temperature, whereas incorporating depth-resolved borehole temperature measurements for initialization yields substantially improved agreement with observations. The findings highlight the sensitivity of short-term predictive modeling to the representation of initial subsurface thermal conditions and underscore the value of high-resolution field measurements for model calibration and validation. These preliminary results inform ongoing efforts to extend the modeling framework to longer time horizons and to refine monitoring and modeling strategies that support the design guidance and operational optimization of GHP systems in research and commercial buildings.

15 GEOTHERMAL ENERGY↗

Tracer gas test in Net-Zero Energy Residential Test Facility

The Net-Zero Energy Residential Test Facility (NZERTF) at NIST was built in 2012 to support the development and adoption of cost-effective net-zero energy designs, technologies, and construction methods. It is two-story, and has four bedrooms, three bathrooms, and an open living room, dining room, and kitchen space. The NZERTF also has a basement and attic, both located within the conditioned space because the thermal and air-moisture barriers encompass the basement walls and attic roof. Transfer grilles link the living spaces to these two zones. The central heating and cooling system includes an air-toair heat pump, which delivers air to the basement, first and second floors. The heat pump has a cooling capacity of 7.6 kW and a heating capacity of 7.8 kW. The indoor unit is in the basement and ductwork runs along the basement ceiling. A balanced heat recovery ventilator (HRV) is installed in the basement and has its own dedicated ductwork. The HRV supplies 47 L/s of outdoor air to the house, with supplies on the first floor (in the kitchen/dining area) and in each of the three second-floor bedrooms. Air from the first-floor bathroom and both second-floor bathrooms is returned to the HRV before being exhausted, though the HRV was turned off during the tests. The house also has a range hood exhaust and a clothes dryer exhaust, both of which were turned off for the tests. The temperature of the house was measured and the operation of the space conditioning system, including its recirculating air distribution fan, were controlled by a commercially available thermostat.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Mini Micropulse Lidar measurements in vertical and scanning modes

The Atmospheric Radiation Measurement (ARM) program operates ten Sigma Space (currently Droplet Measurement Technologies/DMT) micropulse lidar systems (MPL-4B-IDS/POL-FS series), which detect the backscattered light in two orthogonal linear polarization channels allowing to distinguish between spherical (e.g. water) and non-spherical (e.g. ice) scattering particles. The ARM MPLs have been providing valuable measurements to the aerosol and cloud research community for well over a decade. DMT’s Mini-Micropulse lidar (MiniMPL) is a compact version of the standard MPL. It reduces the power-aperture product to minimize cost, size, weight, power requirements and has all the electronics integrated with optics in a single box. Although as a result, it has a lower detection range than the MPL, limiting the measurements to the troposphere, its performance is promised to match that of the standard MPL within the troposphere. The MiniMPL weighs 13 kg, and the complete system, consisting of an optical transceiver and a laptop running the data acquisition and post-processing software, fits in a storm case, allowing its transport on domestic or international flights. This can be especially valuable for short field campaigns or Intensive Operational Periods, where time is of the essence when replacing a faulty instrument. In addition, the scanning polarization lidar capability (MiniMPL-SCAN option) could be of great interest for satisfying ARM measurement needs for certain research topics related to clouds and precipitation. To evaluate its full potential in comparison with co-located remote sensing instruments including the ARM MPL-4B-IDS/POL-FS at SGP, a miniMPL was hosted at the ARM Southern Great Plains (SGP) mega-site in 2020 for approximately 3 months. During this time the miniMPL was tested in both vertical sampling and scanning modes.

54 ENVIRONMENTAL SCIENCES↗

Computational Fluid Dynamics Simulations to Support Efficiency Improvements in Aluminum Smelting Process

Smelting is broadly described as the extraction of a metal from its ore. In the United States, aluminum is commonly produced by smelting alumina in bauxite using the Hall-Héroult process. Optimization of equipment and processes in conventional smelting is crucial to enhancing process efficiency and productivity, is necessary for improving the techno-economic feasibility, which directly manifests as the growth of the American economy. To achieve optima, insightful data on the multiphysics phenomena that are inherent to the process must be obtained through physical investigation or high-fidelity numerical simulations. The resolution of relevant scales in time and space for smelting operations requires intensive, high-performance computing (HPC) simulations. Hostile operating conditions limit physical data acquisition to specific techniques; therefore, these data do not describe the multiscale interaction of simultaneous effects. Fortunately, in recent decades, significant advancements in computing hardware and computational methods have made the numerical resolution of such a complex process possible. In this study, a high-fidelity simulation of aluminum smelting was performed using an open-source tool, OpenFOAM, which analyzed many parameters characteristic to underlying phenomena. A multiphysics model based on the Eulerian-Eulerian multifluid approach was adopted. This model can resolve critical issues in the electrolytic smelting of aluminum, such as bubbling of carbon dioxide from the anode(s), magnetohydrodynamics from electromagnetic effects, ionic dissolution of the alumina in the electrolyte, and the evolution of thermal profiles. This study provides valuable connectivity for characteristic data that can direct the future designs of efficient smelters. A basic framework to model and simulate the smelting process using OpenFOAM is presented for user modification in keeping with process development. Of relevance to the flow field, a detailed investigation of vortices produced by bubble motion and electromagnetics is discussed, along with their impact on the evolution of thermal profiles. The predictions show small-scale vortices in the clearance between the anode and cathode caused by magnetic forces. Predictions also indicate relatively large-scale vortices in the inter-anode space resulting from carbon dioxide rising through the electrolytic flow field. The formation of vortices at the edges of anodes was shown to direct alumina charged by the feeder to the bottom of the anodes, thus preventing the entrapment of gas bubbles in the periphery of the bottom of the anode. Symmetry was observed in the location of cold spots in the electrolytic mixture in the vicinity of the feeder. Cold spots were also observed in the clearance between the anode and cathode due to the flow’s transmission of unconverted alumina to this region.

36 MATERIALS SCIENCE↗

Recursion relation for general 3d blocks

We derive closed-form expressions for all ingredients of the Zamolodchikov-like recursion relation for general spinning conformal blocks in 3-dimensional conformal field theory. This result opens a path to efficient automatic generation of conformal block tables, which has immediate applications in numerical conformal bootstrap program. Our derivation is based on an understanding of null states and conformally-invariant differential operators in momentum space, combined with a careful choice of the relevant tensor structures bases. This derivation generalizes straightforwardly to higher spacetime dimensions d, provided the relevant Clebsch-Gordan coefficients of Spin (d) are known.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Finite-rate entrainment effects on nitrogen oxide (NO x ) emissions in staged combustors

Axially staged combustors offer the potential for achieving low NO x emissions from gas turbines at the elevated temperatures required to achieve at least 65% combined cycle efficiencies. While there is evidence that premixing between the main burner products and secondary stream reduces NO x , it is important to characterize the specific requirements needed to achieve low emissions at gas turbine conditions. This study examines the sensitivity of NO x to finite-rate, large-scale entrainment of the main and secondary streams under the simplifying assumption of infinitely-fast small-scale mixing. Essentially, this allows us to isolate the effects of limited large-scale entrainment rates by using a homogeneous/uniform condition for the entrained and reacting gases. We use a reduced-order reactor network model to examine a generic staged-combustor whose inputs are physical time scales that embody the finite-rate entrainment characteristics. With simulations conducted over a large parameter space at typical operating conditions (25 atm, 650 K air), the results show that, even when small-scale mixing is infinitely fast and the reaction zone is uniform, entrainment significantly affects NO x emissions due to its influence on the equivalence ratio — and thus the temperature and time — at which the entrained mixture ignites. Fuel-air staging is shown to be vital to NO x reduction in most practical cases where entrainment times exceed roughly 1 ms and the secondary stream finishes entraining before the main burner products. A constrained NO x minimization indicates that using the leanest possible main burner and re-routing as much air as possible to the secondary stage is key to low NO x with finite-rate entrainment; for example, less than 10 ppm NO x can be achieved if the secondary stage contains 20% of the total mass flow.

03 NATURAL GAS↗

A spatially-resolved model of neutron-irradiated tungsten coupling stochastic cluster dynamics and finite deformation plasticity

Structural materials used in nuclear reactors face severe degradation in mechanical properties, such as hardening and embrittlement. At the microscopic scale, this occurs due to creation and accumulation of irradiation-induced defects and their interaction with system dislocations. Although techniques exist which can model evolution of irradiation defects, for instance kinetic transport theory-based models, their interaction with mechanical deformation of the bulk material has not been investigated extensively. In this work, we demonstrate a novel spatially-resolved multiscale coupling between microscopic irradiation defect evolution, modeled using Stochastic Cluster Dynamics (SCD) and macroscopic mechanical deformation modeled using a finite-deformation plasticity model. SCD is used to determine the statistically averaged defect cluster spacing, dependent on operating conditions such as irradiation dose and temperature. This acts as an initial condition that governs the critical resolved shear stress of dislocation glide in the macroscopic plasticity model. This framework is used to predict mechanical behavior in post-mortem test of irradiated Tungsten samples, which has found its importance as structural material used in nuclear reactors. The results obtained using the coupled approach are in good agreement with experimental data of uniaxial tension tests. The model is able to capture the effect of temperature and irradiation dose on the material hardening. Two methods are proposed to estimate hardness – using Tabor's Law relating uniaxial yield stress to hardness and from flat-punch simulations. The results are in reasonable agreement with hardness data from micro-indentation experiments of irradiated Tungsten samples. Finally, the model is also able to reveal microstructural details such as spatial variation in defect density and local stress.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Kinetic closures for unmagnetized and magnetized plasmas

Parallel and perpendicular closures with cyclotron resonance effects retained for the five-moment (density, temperature, and flow velocity) fluid equations are derived by solving the kinetic equation with the Bhatnagar–Gross–Krook operator in Fourier space. For parallel propagation, the parallel closures are reduced to those of Ji et al. [Phys. Plasmas 20, 082121 (2013)]. The closures when combined to the fluid equations reproduce the fully kinetic dispersion relation that can be directly derived from the kinetic equation. The closures for the five-moment fluid system can be utilized to derive closures for the extended fluid system, which is demonstrated by deriving closures for the ten-moment system consisting of density, flow velocity, temperature, and viscosity tensor equations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Progress of the laser ion source upgrade (LION2) for RHIC and NSRL program at Brookhaven National Laboratory

At Brookhaven National Laboratory (BNL), the LION2 ion source is being constructed to replace an existing laser ion ablation ion source (LIS) at the EBIS facility, which provides heavy ion beams of multiple ion species for the operation of NASA Space Radiation Laboratory (NSRL) and Relativistic Heavy Ion Collider (RHIC). The LION1 ion source currently provides singly charged ions of Li, B, C, O, Al, Si, Ca, Ti, Fe, Cu, Zr, Nb, Ag, Tb, Ta, W, Au, Bi, and Th with a rapid-species-change capability. An electron beam ion source, Extended-EBIS captures, confines, and ionizes the ions to high charge state, suitable for injection and acceleration by an RFQ accelerator. Typically, single pulses of the LIS ion species for NSRL are changed sequentially during Galactic Cosmic Ray experiments, while multiple pulses of a given ion beam are provided quasi-simultaneously for RHIC. LION2 will have the same capability of the rapid-species-change with improved beam performance and reliability. LION2 is being constructed in a remote assembly location and is expected to finish in December 2023. The removal of LION1 and installation of LION2 is planned during the December 2023 or summer 2024 shutdown.

43 PARTICLE ACCELERATORS↗

Latent Twins

Over the past decade, scientific machine learning has transformed the development of mathematical and computational frameworks for analyzing, modeling, and predicting complex systems. From inverse problems to numerical partial differential equations (PDEs), dynamical systems, and model reduction, these advances have pushed the boundaries of what can be simulated. Yet they have often progressed in parallel, with representation learning and algorithmic solution methods evolving largely as separate pipelines. With Latent Twins, we propose a unifying mathematical framework that creates a hidden surrogate in latent space for the underlying equations. Whereas digital twins mirror physical systems in the digital world, Latent Twins mirror mathematical systems in a learned latent space governed by operators. Through this lens, classical modeling, inversion, model reduction, and operator approximation all emerge as special cases of a single principle. We establish the fundamental approximation properties of Latent Twins for both ordinary differential equations (ODEs) and PDEs and demonstrate the framework across three representative settings: (i) canonical ODEs, capturing diverse dynamical regimes; (ii) a PDE benchmark using the shallow-water equations, contrasting Latent Twin simulations with deep operator network and forecasts with a four-dimensional variational method baseline; and (iii) a challenging real-data geopotential reanalysis dataset, reconstructing and forecasting from sparse, noisy observations. Latent Twins provide a compact, interpretable surrogate for solution operators that evaluate across arbitrary time gaps in a single-shot, while remaining compatible with scientific pipelines such as assimilation, control, and uncertainty quantification. Looking forward, this framework offers scalable, theory-grounded surrogates that bridge data-driven representation learning and classical scientific modeling across disciplines.

Latent Twins↗

Phonons in complex twisted crystals: Angular momenta, interactions, and topology

We report phonons and their interactions are critically important for a variety of energy-relevant applications ranging from low thermal resistance substrates to thermal barrier coatings. Fundamental insights into the nature of phonons and allowed interactions are governed by the underlying symmetries of the crystal lattice. In this paper, we provide a comprehensive and detailed description of lattice dynamics derived from twist symmetries of chiral and achiral crystals—twist dynamics. Phonon bands naturally carry quantized crystal angular momentum derived from rotational phases, which give insights into the nature of band crossings and avoided crossings, selection rules for phonon interactions, and topological band crossing behaviors. Twist dynamics is demonstrated for a variety of materials covering a range of space groups, symmetry operations, twist axis orientations, and constituent elements. Furthermore, twist symmetry offers insights into peculiar features observed in scattering measurements. In this context, we present inelastic neutron scattering measurements for rutile TiO 2 and explain them using twist dynamics.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Converting three-space matrices to equivalent six-space matrices for Delone scalars in S 6

The transformations from the primitive cells of the centered Bravais lattices to the corresponding centered cells have conventionally been listed as three-by-three matrices that transform three-space lattice vectors. Using those three-by-three matrices when working in the six-dimensional space of lattices represented as Selling scalars as used in Delone (Delaunay) reduction, one could transform to the three-space representation, apply the three-by-three matrices and then back-transform to the six-space representation, but it is much simpler to have the equivalent six-by-six matrices and apply them directly. The general form of the transformation from the three-space matrix to the corresponding matrix operating on Selling scalars (expressed in space S 6 ) is derived, and the particular S 6 matrices for the centered Delone types are listed. (Note: in his later publications, Boris Delaunay used the Russian version of his surname, Delone.)

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Emergent tetragonality in a fundamentally orthorhombic material

Symmetry plays a key role in determining the physical properties of materials. By Neumann’s principle, the properties of a material remain invariant under the symmetry operations of the space group to which the material belongs. Continuous phase transitions are associated with a spontaneous reduction in symmetry. Less common are examples where proximity to a continuous phase transition leads to an increase in symmetry. We find signatures of an emergent tetragonal symmetry close to a charge density wave (CDW) bicritical point in a fundamentally orthorhombic material, ErTe 3 , for which the two distinct CDW phase transitions are tuned via anisotropic strain. We first establish that tension along the a axis favors an abrupt rotation of the CDW wave vector from the c to a axis and infer the presence of a bicritical point where the two continuous phase transitions meet. We then observe a divergence of the nematic elastoresistivity approaching this putative bicritical point, indicating an emergent tetragonality in the critical behavior.

36 MATERIALS SCIENCE↗

Bilevel Parameter Optimization for Nonlocal Image Denoising Models

We propose a bilevel optimization approach for the determination of parameters in nonlocal image denoising. We consider both spatial weights in front of the fidelity term, as well as weights within the kernel of the nonlocal operator. In both cases we investigate the differentiability of the solution operator in function spaces and derive a first order optimality system that characterizes local minima. For the numerical solution of the problems, we propose a second-order optimization algorithm in combination with a finite element discretization of the nonlocal denoising models and a computational strategy for the solution of the resulting dense linear systems. Several experiments are run in order to show the suitability of our approach.

97 MATHEMATICS AND COMPUTING↗

An Integrated Feasibility Study of Reservoir Thermal Energy Storage (RTES) in Portland, OR, USA

In regions with long cold overcast winters and sunny summers, Deep Direct-Use (DDU) can be coupled with Reservoir Thermal Energy Storage (RTES) technology to take advantage of pre-existing subsurface permeability to save summer heat for later use during cold seasons. Many aquifers worldwide are underlain by permeable regions (reservoirs) containing brackish or saline groundwater that has limited beneficial use due to poor water quality. We investigate the utility of these relatively deep, slow flowing reservoirs for RTES by conducting an integrated feasibility study in the Portland Basin, Oregon, USA, developing methods and obtaining results that can be widely applied to RTES systems elsewhere. As a case study, we have conducted an economic and social cost-benefit analysis for the Oregon Health and Science University (OHSU), a teaching hospital that is recognized as critical infrastructure in the Portland Metropolitan Area. Our investigation covers key factors that influence feasibility including 1) the geologic framework, 2) heat and fluid flow modeling, 3) capital and maintenance costs, 4) the regulatory framework, and 5) operational risks. By pairing a model of building seasonal heat demand with an integrated model of RTES resource supply, we determine that the most important factors that influence RTES efficacy in the study area are operational schedule, well spacing, the amount of summer heat stored (in our model, a function of solar array size), and longevity of the system. Generally, heat recovery efficiency increases as the reservoir and surrounding rocks warm, making RTES more economical with time. Selecting a base-case scenario, we estimate a levelized cost of heat (LCOH) to compare with other sources of heating available to OHSU and find that it is comparable to unsubsidized solar and nuclear, but more expensive than natural gas. Additional benefits of RTES include energy resiliency in the event that conventional energy supplies are disrupted (e.g., natural disaster) and a reduction in fossil fuel consumption resulting in a smaller carbon footprint. Key risks include reservoir heterogeneity and a possible reduction in permeability through time due to scaling (mineral precipitation). Lastly, a map of thermal energy storage capacity for the Portland Basin yields a total of 43,400 GWh, suggesting tremendous potential for RTES in the Portland Metropolitan Area.

14 SOLAR ENERGY↗