Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “local minimum”

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 217 records · Page 12
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

SPiRaL: a multiresolution global tomography model of seismic wave speeds and radial anisotropy variations in the crust and mantle

SUMMARY SPiRaL is a joint global-scale model of wave speeds (P and S) and anisotropy (vertical transverse isotropy, VTI) variations in the crust and mantle. The model is comprised of >2.1 million nodes with five parameters at each node that capture velocity variations for P- and S-waves travelling at arbitrary directions in transversely isotropic media with a vertical symmetry axis (VTI). The crust (including ice, water, sediments and crystalline layers) is directly incorporated into the model. The default node spacing is approximately 2° in the lower mantle and 1° in the crust and upper mantle. The grid is refined with ∼0.25° minimum node spacing in highly sampled regions of the crust and upper mantle throughout North America and Eurasia. The data considered in the construction of SPiRaL includes millions of body wave traveltimes (crustal, regional and teleseismic phases with multiples) and surface wave (Rayleigh and Love) dispersion. A multiresolution inversion approach is employed to capture long-wavelength heterogeneities commonly depicted in global-scale tomography images as well as more localized details that are typically resolved in more focused regional-scale studies. Our previous work has demonstrated that such global-scale models with regional-scale detail can accurately predict both teleseismic and regional body wave traveltimes, which is necessary for more accurate location of small seismic events that may have limited signal at teleseismic distances. SPiRaL was constructed to predict traveltimes for event location and long-period waveform dispersion for seismic source inversion applications in regions without sufficiently tuned models. SPiRaL may also serve as a starting model for full-waveform inversion (FWI) with the goal of fitting waves with periods 10–50 s over multiple broad regions (thousands of kilometres) and potentially the globe. To gain insight to this possibility, we simulated waveforms for a small set of events using SPiRaL and independent waveform-based models for comparison. For the events tested, the performance of the traveltime-based SPiRaL model is shown to be generally on par with regional 3-D waveform-based models in three regions (western United States, Middle East, Korean Peninsula) suggesting SPiRaL may serve as a starting model for FWI over broad regions.

58 GEOSCIENCES↗
DOE OSTI2026

Carbon Capture on Air Liquide United States Gulf Coast Steam Methane Reformer using the Cryocap TM Flue Gases Process

The objective of this DOE-funded project is to execute and complete a Front-End Engineering and Design (FEED) study for a commercial-scale carbon capture project for Air Liquide’s Steam Methane Reformer (SMR) located at La Porte, near Houston, Texas. The Host Site is one of Air Liquide’s largest SMR assets, and supplies hydrogen to the Air Liquide Gulf Coast Hydrogen System, which serves refiners and petrochemical manufacturers through an existing H2 pipeline spanning roughly 200 miles from Port Arthur, TX to Bay City, TX. The Host Site produces up to 116.5 million standard cubic feet per day (MMSCFD) of gaseous hydrogen and emits approximately 950 ktpy CO 2 at the nameplate capacity. The evolving importance of decarbonized H2 for both existing customers and new markets, as well as the availability and intent of prominent third-party partners to transport and sequester captured CO 2 in locally abundant geological sites for 45Q credits, make the AL US Gulf Coast SMR an ideal site for the proposed FEED study. The envisioned carbon capture system is based on Air Liquide’s proprietary Cryocap™ Flue Gas (FG) process. The integration of the Cryocap™ FG technology to the existing SMR would enable the capture of 900 ktpy of CO 2 , with a net carbon capture rate of >95% and with minimum impact on the levelized cost of hydrogen produced at 99.97% purity.

03 NATURAL GAS↗
DOE OSTI2022

Modeling the distribution of the endangered Jemez Mountains salamander (Plethodon neomexicanus) in relation to geology, topography, and climate

The Jemez Mountains salamander (Plethodon neomexicanus; hereafter JMS) is an endangered salamander restricted to the Jemez Mountains in north-central New Mexico, United States. This strictly terrestrial and lungless species requires moist surface conditions for activities such as mating and foraging. Threats to its current habitat include fire suppression and ensuing severe fires, changes in forest composition, habitat fragmentation, and climate change. Forest composition changes resulting from reduced fire frequency and increased tree density suggest that its current aboveground habitat does not mirror its historically successful habitat regime. However, because of its limited habitat area and underground behavior, we hypothesized that geology and topography might play a significant role in the current distribution of the salamander. We modeled the distribution of the JMS using a machine learning algorithm to assess how geology, topography, and climate variables influence its distribution. The best habitat suitability model indicates that geology type and maximum winter temperature (November to March) were most important in predicting the distribution of the salamander (23.5% and 50.3% permutation importance, respectively). Minimum winter temperature was also an important variable (21.4%), suggesting this also plays a role in salamander habitat. Our habitat suitability map reveals low uncertainty in model predictions, and we found slight discrepancies between the designated critical habitat and the most suitable areas for the JMS. Because geological features are important to its distribution, we recommend that geological and topographical data are considered, both during survey design and in the description of localities of JMS records once detected.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Improved representation of laminar and turbulent sheet flow in subglacial drainage models

Subglacial hydrology models struggle to reproduce seasonal drainage patterns that are consistent with observed subglacial water pressures and surface velocities. We modify the standard sheet-flow parameterization within a coupled sheet–channel subglacial drainage model to smoothly transition between laminar and turbulent flow based on the locally computed Reynolds number in a physically consistent way (the ‘transition’ model). We compare the transition model to standard laminar and turbulent models to assess the role of the sheet-flow parameterization in reconciling observed and modelled water pressures under idealized and realistic forcing. Relative to the turbulent model, the laminar and transition models improve seasonal simulations by increasing winter water pressure and producing a more prominent late-summer water pressure minimum. In contrast to the laminar model, the transition model remains consistent with its own internal assumptions across all flow regimes. Based on the internal consistency of the transition model and its improved performance relative to the standard turbulent model, we recommend its use for transient simulations of subglacial drainage.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Commercial Zero Code Plug-In: Zero Energy and Operational Emissions Overlay for Model Energy Codes (Technical Brief)

Model energy codes (MEC) describe requirements needed to demonstrate building energy performance compliance. They can be readily adopted by states and local jurisdictions to support desirable energy efficiency investment in new buildings and major renovations. This report provides commercial building energy code language as an overlay to current MEC, which is recognized by the U.S. DOE as ASHRAE Standard 90.1 2022. The code language provides a performance-based compliance path for achieving net zero energy or net zero operational energy emissions buildings with MEC. The path requires meeting two compliance metric target values: 1) a required minimum level of efficiency and 2) a measure of zero energy or emissions. The supporting documentation illustrates the magnitude of needed efficiency improvement and the additional offsets to be achieved from renewable energy sources.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2026

Insights into Rotational and Translational Dynamics in Mixtures of Ethylene Glycol and Choline Chloride Using Nuclear Magnetic Resonance Techniques

This work examines molecular dynamics and interactions in ethylene glycol–choline chloride (EG–ChCl) mixtures across 0–33 mol % ChCl, spanning the true eutectic region near 17–20 mol % and the commonly used 1:2 formulation. We combine pulsed-field-gradient (PFG) diffusion, fast-field-cycling (FFC) relaxometry, temperature-dependent 13 C T 1 , and nuclear Overhauser effect spectroscopy (NOESY) to disentangle local from macroscopic dynamics. PFG and FFC show that both translational and average rotational motions largely track the strong increase in viscosity with ChCl content, with ethylene glycol consistently diffusing faster than the choline cation and no global dynamical anomaly at the eutectic composition. More subtle, site-specific composition effects nevertheless emerge. The ratio of the diffusion coefficient of the hydroxyl group of choline to the diffusion coefficient of the methyl group of choline displays a shallow minimum in the 17–25 mol % region, indicating a modest change in how the hydroxyl-bearing end of choline samples the underlying translational motion relative to the methyl groups. 13 C T 1 analysis shows that rotational correlation times at 25 °C generally increase with ChCl, reflecting viscosity-coupled slowing, while the CH 2 –N α site exhibits a small but reproducible deviation from this monotonic trend near the eutectic. NOESY spectra at similar compositions reveal enhanced cross-relaxation between EG and choline protons, consistent with increased headgroup–solvent contact density rather than a wholesale structural rearrangement. Overall, our multitechnique study demonstrates that EG–ChCl dynamics are predominantly viscosity-dominated, with the eutectic region acting as a subtle dynamical crossover where specific choline segments become maximally coupled to the hydrogen-bond network. These insights refine the structure–dynamics picture of choline-chloride DESs and provide practical guidance for tuning composition in electrochemical, separation, and catalytic applications.

diffusion↗
DOE OSTI2021

Pericentric passage-driven star formation in satellite galaxies and their hosts: CLUES from local group simulations

Local Group satellite galaxies show a wide diversity of star formation histories (SFHs) whose origin is yet to be fully understood. Using hydrodynamical simulations from the Constrained Local UniversE project, we study the SFHs of satellites of Milky Way-like galaxies in a cosmological context: while in the majority of the cases the accretion on to their host galaxy causes the satellites to lose their gas, with a subsequent suppression in SF, in about 25 per cent of our sample we observe a clear enhancement of SF after infall. Peaks in SF clearly correlate with the satellite pericentric passage around its host and, in one case, with a satellite–satellite interaction. We identify two key ingredients that result in enhanced SF after infall: galaxies must enter the host’s virial radius with a reservoir of cold gas M gas , inf/M vir , inf ≳ 10 -2 and with a minimum pericentric distance ≳10 kpc (mean distance ~50 kpc for the full sample), in order to form new stars due to compression of cold gas at pericentric passage. On the other hand, satellites that infall with little gas or whose pericentric distance is too small, have their gas ram-pressure stripped and subsequent SF quenched. The pericentric passage of satellites likewise correlates with SF peaks in their hosts, suggesting that this mechanism induces bursts of SF in satellites and central galaxies alike, in agreement with recent studies of our Galaxy’s SFH. Our research findings can explain the recently reported multiple stellar populations observed in dwarf galaxies such as Carina and Fornax, and should be taken into account in semi-analytic models of galaxy formation and satellite quenching.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Optimum Injector Parameters for Thermoacoustic Stability in a Multi-Nozzle Can Combustion System

Abstract High-frequency transverse instabilities are an important concern in can combustor configurations. In these configurations which are typically operated with multiple injectors around a central injector, each injector is subjected to different parts of the acoustic mode shape and thus respond differently for the same instability mode. Recent work by the author has modeled the response of premixed flames to excitation by natural high-frequency transverse modes in a can combustor both in the center and outer nozzles. The stability of these acoustically non-compact flames was assessed using the Rayleigh criterion (Rayleigh Integral denoted as RI) and not the overall unsteady heat release as is the case for compact flames. Several key control parameters were studied, namely — flame angle, swirling strength, nozzle location. For non-axisymmetric modes such as the commonly occurring 1-T mode, both radial and azimuthal offsets of the nozzle location affected stability. The framework was applied to an optimization study to identify the optimal combination of parameters that minimizes RI for the different nozzles in the multi-nozzle system. In this study, a N-around-1 configuration was studied, and the results indicated that the different nozzles needed to be operated at different flame angles and swirl numbers to result in an overall minimum RI. However, the specific response of the different injectors was not considered. The helical mode distribution at each injector varies as we azimuthally go around the combustor’s injector distribution and thus the most amplified mode and the resulting flame response would be different. To minimize RI, it is important to determine the injector configurations that result in a hydrodynamic profile that minimizes the individual RI for each nozzle. The resulting relationship between the injector’s flow and local hydrodynamics can then be used in a hydrodynamics study of an individual injector so that the most optimal injector is chosen depending on its location in the combustor dump plane.

Acharya, Vishal↗
DOE OSTI2026

Shuttling Majorana zero modes in disordered and noisy topological superconductors

The braiding of Majorana zero modes (MZMs) forms the fundamental building block for topological quantum computation. Braiding protocols which involve the physical exchange of MZMs are typically envisioned on a network of topological superconducting wires. An important component of these protocols is the transport of MZMs, which can be performed by using electric gates to locally tune sections of the wire between topologically trivial and nontrivial phases. In this work, we numerically simulate this transport by tuning a single section of a superconducting wire which contains either disorder (uncorrelated and correlated) or noise. We focus on the impact of these additional effects on the diabatic error, which describes unwanted transitions between the ground state and excited states. We show that the behavior of the average diabatic error is predominantly controlled by the statistics of the minimum bulk energy gap which is suppressed in the presence of disorder. The increase in diabatic error can be several orders of magnitude and is most deleterious when the disorder correlation length is a finite fraction of the transport distance and negligible when these lengths are far apart. In the presence of noise, the diabatic error is significantly enhanced due to optical transitions which depend on the minimum bulk energy gap as well as the frequency modes present in the noise. The results presented here serve to further characterize the diabatic error in disordered and noisy settings, which are important considerations in practical implementations of physical braiding schemes.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Void-free Copper Electrodeposition in High Aspect Ratio, Full Wafer Thickness Through-Silicon Vias with Endpoint Detection

High density interconnects are required for increased input/output for microelectronics applications, incentivizing the development of Cu electrochemical deposition (ECD) processes for high aspect ratio through-silicon vias (TSVs). This work outlines Cu ECD processes for 62.5 μ m diameter TSVs, etched into a 625 μ m thick silicon substrate, a 10:1 aspect ratio. Cu ECD in high aspect ratio features relies on a delicate balance of electrolyte composition, solution replenishment, and applied voltage. Implementing a CuSO 4 -H 2 SO 4 electrolyte, which contains suppressor and a low chloride concentration, allows for a tunable relationship between applied voltage and localized deposition in the vias. A stepped potential waveform was applied to move the Cu growth front from the bottom of the via to the top. Sample characterization was performed through mechanical cross-sections and X-ray computed tomography (CT) scans. The CT scans revealed small seam voids in the Cu electrodeposit, and process parameters were tuned accordingly to produce void-free Cu features. During the voltage-controlled experiments, measured current data showed a characteristic current minimum, which was identified as an endpoint detection method for Cu deposition in these vias. We believe this is the first report of this novel endpoint detection method for TSV filling.

42 ENGINEERING↗
DOE OSTI2024

The ASSIST trial: Acute effects of manipulating strength exercise volume on insulin sensitivity in obese adults: A protocol for a randomized controlled, crossover, clinical trial

Type 2 diabetes mellitus is a disease in which insulin action is impaired, and an acute bout of strength exercise can improve insulin sensitivity. Current guidelines for strength exercise prescription suggest that 8 to 30 sets could be performed, although it is not known how variations in exercise volume impact insulin sensitivity. Additionally, this means an almost 4-fold difference in time commitment, which might directly impact an individual’s motivation and perceived capacity to exercise. This study will assess the acute effects of high- and low-volume strength exercise sessions on insulin sensitivity. After being thoroughly familiarized, 14 obese individuals of both sexes (>40 year old) will undergo 3 random experimental sessions, with a minimum 4-day washout period between them: a high-volume session (7 exercises, 3 sets per exercise, 21 total sets); a low-volume session (7 exercises, 1 set per exercise, 7 total sets); and a control session, where no exercise will be performed. Psychological assessments (feeling, enjoyment, and self-efficacy) will be performed after the sessions. All sessions will be held at night, and the next morning, an oral glucose tolerance test will be performed in a local laboratory, from which indexes of insulin sensitivity will be derived. We believe this study will aid in strength exercise prescription for individuals who claim not to have time to exercise or who perceive high-volume strength exercise intimidating to adhere to. This trial was prospectively registered.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2025

Methodology to Establish Performance Targets for Building Energy Codes

This report explains the underlying development of the U.S. NZE building sector performance targets reported in Federal Register Notice XX. This document also includes NZE performance targets for U.S. states, which were developed following the same procedures used to develop the national values. These national and state targets are intended to support the establishment and attainment of NZE goals over time. Where residential and commercial building energy code compliance mechanisms apply to new buildings and major renovations and dictate performance requirements, these targets reflect the recommended, minimum levels of performance. Appendix A provides the performance values and targets by U.S. state for residential and commercial buildings. These include each state’s adopted energy code performance, current published MEC performance value, and a nominal NZE code compliance performance target. The current code and NZE values can be compared to the state adopted code value to indicate the level of advancement needed to move towards recommended NZE building performance levels. Additional supplemental resources are available to support states and local governments ready to move from setting NZE policy goals to adopting NZE codes. More information about these resources and the actions states can take to achieve NZE buildings with energy codes are provided in the body of this document.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2025

Propagation of EMIC Waves From Shabansky Orbits in the Dayside Magnetosphere

We explore the characteristics of EMIC waves generated in a non-dipole, compressed magnetic field at the minimum of the magnetic field. We conducted 2D full-wave simulations using the Petra-M code, focusing on a compressed magnetic field in the outer dayside magnetosphere for a range of L values (L = 8 - 10). By comparing the simulation results with MMS observations, we aim to understand how the observed wave characteristics are affected by a shifting source region across different L-shells. Our findings indicate that the direction of the Poynting vector systematically changes depending on the local source location of the wave, which is consistent with the observations. EMIC waves propagate along the magnetic field line and reach both the northern and southern hemispheres; however, there is a notable difference in the power of EMIC waves between the two hemispheres, indicating seasonal asymmetries in their occurrence.

58 GEOSCIENCES↗
DOE OSTI2020

Temperature-driven BCS-BEC crossover and Cooper-paired metallic phase in coupled boson-fermion systems

Motivated by strongly correlated and frustrated systems, in this work we propose an effective model that describes the dynamics of pairs of opposite spin fermions scattering from localized bosons. Integrating out one of the degrees, either the bosons or fermions, generates temperature-dependent long-range effective interactions between the entities that we investigate using Monte Carlo techniques. The behavior of bosons is dominated by vortex-antivortex unbinding, with effective interboson interactions beyond the nearest-neighbor Josephson coupling of phases. Remarkably, in the fermion sector we observe a temperature-driven phase transition from a SC phase with a “BCS” spectral function that shows a gap minimum on the underlying Fermi surface to a conducting phase of pairs but with a “BEC” spectral function with a gap minimum at k = 0. Tunneling and angle-resolved photoemission spectroscopy on Bose-Fermi mixtures in cold atomic systems and superconducting islands on graphene are some of the promising experimental platforms to test our predictions.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Local Thermomechanical Processing for Improving Formability of High Strength Aluminum Sheets

Limited room temperature formability hinders the wide-spread use of high strength aluminum alloys in structural body-in-white parts. Stamping or extrusion at warm temperatures or from softer tempers are the current solutions. In this work, our approach is to start with age-hardened sheets from 7xxx and 6xxx family of alloys and improve their formability using local thermomechanical processing only in the regions demanding highest ductility in the forming processes. We achieved local formability improvements with friction stir processing and introduce another process named roller-bending as a concept and showed its feasibility through finite element simulations. Initial results from FSP indicated significant deformation in the processed zones with minimal sheet distortion. FSP also resulted in dynamically recrystallized, fine grained (d < 5 µm) microstructures in the processed regions with textures significantly different from the base material. This resulted in formability improvements of > 30% in 7085-T76 and > 50% in 6111-T6 alloys as measured by minimum bend radius (r/t) ratios established with V-bending tests performed at 90° and 180° bend angles. Based on the encouraging initial results, we envision applicability of these processes in both pre- and post-stamping stages through integration of the process tools as end-effectors to robotic arms.

Nasim, Wahaz↗
DOE OSTI2025

Creep ductility limiting mechanisms in an additively manufactured Al-Ce-Ni-Mn-Zr alloy

Tensile creep response and cavitation damage evolution in an additively manufactured Al-7.5Ce-4.5Ni-0.4Mn-0.7Zr (wt%) alloy with peak-aging and overaging treatments were investigated in the 300–400 ºC range. Microstructural heterogeneity and its response to heat treatment and subsequent creep deformation were studied to understand the interplay between cavity formation, creep lifetime and ductility. Increasing the applied stress activated the nucleation of more cavities, an experimental observation that is well described using the vacancy accumulation model. Cavities nucleated prematurely due to localized plasticity in the denuded zones that formed at/near melt-pool or grain boundaries. Microstructure/deformation heterogeneity with consequent evolution of stress triaxiality, especially at lower stresses, causes accelerated cavitation, thus producing low creep ductility (∼ 0.2–2.4 %), compared to (∼12–21 %) ductility of the alloy measured by regular tensile tests at equivalent temperatures. A constrained diffusional cavity growth mechanism with continuous cavity nucleation during creep is established as the dominant mechanism, implying that cavitation involves vacancy diffusion, yet its growth rate is dictated by the minimum creep rate. In conclusion, the ductility-limiting creep and cavitation mechanisms discussed here provide new insight into the creep behavior of 3D-printed metallic alloys.

Additive manufacturing↗
DOE OSTI2020

ELM Suppression by Boron Powder Injection and Comparison with Lithium Powder Injection on EAST

Type I edge-localized modes (ELMs) in the Experimental Advanced Superconducting Tokamak (EAST) were completely suppressed via boron powder injection into the X-point region of an upper-single null configuration over a wide range of operating conditions (2.8 < P aux < 7.5 MW, 3.8 × 10 19 < n e < 6 × 10 19 m -3 , RF-only and RF + NBI heating scenarios, both grad-B drift directions, and even He ion majority plasmas) (Sun et al. in Nucl. Fusion, 2020). A window of edge B concentration for stable long pulse operation was identified: too low and ELMs return, too high and the discharge suffers radiative collapse. The injection of boron powder above the minimum for ELM suppression coincided with the occurrence of an edge harmonic oscillation detected in magnetics (both on the high-field side and low-field side), in AXUV diodes near the upper X-point, divertor Dα emission, and in a range of other diagnostics (Diallo et al. in: Proceedings of 2020 IAEA fusion energy conference, 2021). Furthermore, no harmonic oscillation was observed when ELMs were present, and stored energy was slightly increased at constant density during ELM suppression. Core tungsten emission during ELM suppression either increased or decreased relative to ELMy H-mode, but the W emission was maintained at acceptable levels. The threshold B injection rate was measured for several conditions, and found to increase with heating power. Li powder injection into comparable discharges also resulted in a short phase of ELM suppression, but density and stored energy both decreased due to the strong pumping effect of lithium; no edge harmonic oscillation was observed with Li injection, indicating that the ELM suppression mechanisms differ. The new set of B-seeded, ELM-suppressed discharges exhibited certain characteristics of quiescent H-mode (Burrell et al. in Phys Plasmas 8:2153, 2001), but did not require high shear, counter beams, etc. The wide operating window and compatibility with RF-only discharges paves the way for future experiments targeting long pulse H-mode discharges with complete ELM suppression.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Optimizing the physical design and layout of a resilient wind, solar, and storage hybrid power plant

We report as wind and solar technologies improve and their costs decrease, the share of power produced by these sources will increase. As the market penetration increases, these power sources will need to provide grid services, such as dispatchability, in addition to providing energy. One way to reduce variability, provide higher quality power to the grid, and address local grid stability issues is through colocating wind and solar power plants. In addition to operating reliably during normal operating conditions, in scenarios with high penetrations of renewable generation, it is important that these hybrid plants can withstand production disruptions and continue to supply power despite prolonged resource reduction, extreme weather events, or other disruptions. In this paper, we present a methodology to optimize a wind-solar-battery hybrid power plant down to the component level that is resilient against production disruptions and that can continually produce some minimum required power. We introduce the models and assumptions we used to simulate a hybrid power plant as well as the design variable parameterization and specific methods we used to optimize the plant. We demonstrate the performance of our method by comparing a plant optimized for different objectives, generation outage durations, minimum power requirements, and power purchase agreements. Although the plant design is sensitive to model parameters and various other assumptions, our results demonstrate some of the optimal designs that occur in different scenarios and what one should expect when designing a hybrid wind-solar-storage power plant.

14 SOLAR ENERGY↗