Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “validation”

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,207 records · Page 67

Validating Greater Sage-Grouse Individual-based Model (IBM) Tool (Final Report)

The project focused on validating the previously developed Greater Sage-Grouse Individual-based Model (GrSG IBM; LaGory et al. 2012, 2021). The objective was to transform this predictive, spatially and temporally explicit model into a portable resource to assist siting/resource managers in proactively assessing the cumulative impacts of wind energy development on the greater sage-grouse. Utilizing a bottom-up, individual-based approach, the GrSG IBM accounts for landscape context and species behavior, aiming to reduce uncertainty in estimating development impacts and support ecologically mindful land-based wind energy development. The validation effort covered approximately 6,540 km 2 near the Seven Mile Hill Wind Project in Wyoming. The GrSG IBM tool, built on the NetLogo platform (Tisue and Wilensky 2004), was executed over a 50-year period, with the analysis focusing on years following a 10-year initialization phase. Key results demonstrated the tool’s biological soundness across five key biological metrics: non-chick age class distribution (older than 10 weeks), adult sex ratio, life expectancy, population size, and overall population growth. For instance, the tool estimated that 58.6% of the non-chick population was reproductively immature, while the reference ranges from 51.4% to 57.8% (Patterson 1952, Rogers 1964). Experts confirmed the tool’s estimate was within a reasonable range for the species. The tool estimated average life expectancy of 1.43 years, while the reference ranges from 0.9 years to 1.1 years (Ammann 1957, Hamerstrom 1949). Experts also supported the model’s life-expectancy estimate as ecologically sound for the species in the study area. In terms of population change, the model estimated an annual shift between a 0.6% decline and a 1.0% increase over 50 years. While the reference suggests 2.9% annual decline in range-wide populations (Cortes et al. 2023), that includes many at-risk populations in South Dakota and Washington, for example. Our study area—in the northeastern part of Carbon County and western-edge of Albany County, Wyoming—is one of the remaining greater sage-grouse habitats supporting some of the most stable populations. Experts confirmed that the range of the annual population change spanning from a 0.6% decline to a 1.0% increase estimated by the tool was reasonable for our study area for this reason and confirmed that aligned with population estimates from existing studies on the greater sage-grouse and wind energy development in the study area (LeBeau et al. 2017a, Smith et al. 2024). Furthermore, the project showed that temporally explicit biological metrics generated by the GrSG IBM tool can complement the USGS’ Prioritizing Restoration of Sagebrush Ecosystems Tool (PReSET; Duchardt et al. 2021) by incorporating habitat restoration strategies into seasonal habitat suitability models to visualize population responses over time.

17 WIND ENERGY↗

Microreactor Core Transportation Cask Model Description for Criticality Safety Validation Basis Assessment (Rev. 1)

Criticality safety analyses are completed on transportation casks used for microreactor whole core shipment to provide examples of models and analyses to industry, regulators, and nuclear community at large to be used in verification and validation analyses of similar applications. The microreactors considered are based on a Gas-Cooled Microreactor (GCMR) and a Heat-Pipe Microreactor (HPMR), both utilize HALEU fuel in the form of TRISO particles and various other design options considered in industry microreactor designs. Variant design options of GCMR and HPMR were also investigated to provide a wider application range for each technology. Criticality safety analyses for the GCMR and HPMR packages were performed using the CSAS6 sequence of SCALE 6.3.2 with the ENDF/B-VII.1-based continuous energy neutron libraries. Different scenarios were investigated, including normal operation and water flooded conditions to represent nominal and hypothetical accident scenarios. Sensitivity and similarity analyses are also performed using the TSUNAMI sequence of SCALE 6.3.2, and the similarity analysis uses all the experiments from the ICSBEP Handbook with Intermediate and mixed Enriched Uranium (IEU) and Low Enriched Uranium (LEU) systems, together with additional experiments that are sponsored by the DNCSH program. Many experiments were found marginally similar to the GCMR and HPMR models, with similarity index (ck) values greater than 0.8 but less than 0.9. Among all the experiments analyzed, no case has a ck value greater than 0.9, indicating that additional critical experiments might be needed to further validate the criticality safety models for microreactors transport packages.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Data Compilation and Analysis from the Sirius-1 Experiment at TREAT for Transient Simulation Validation

This report assembles comprehensive data from the Sirius-1 experiment conducted by Idaho National Laboratory in collaboration with the National Aeronautics and Space Administration. The primary goal is to provide a robust data set that external users can utilize for the validation of computational methods for transient multiphysics simulations. By compiling all relevant data, including experiment design calculations, detailed engineering drawings for the experiment and data from reactor and fuel specimen measurement, this report is intended to serves as a reference for researchers and engineers working on the development and validation of computational models for transient nuclear behavior.

21 - SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLAN↗

Microreactor Core Transportation Cask Model Description for Criticality Safety Validation Basis Assessment (Rev. 2)

Criticality safety analyses are completed on transportation casks used for microreactor whole core shipment to provide examples of models and analyses to industry, regulators, and nuclear community at large to be used in verification and validation analyses of similar applications. The microreactors considered are based on a Gas-Cooled Microreactor (GCMR) and a Heat-Pipe Microreactor (HPMR), both utilize HALEU fuel in the form of TRISO particles and various other design options considered in industry microreactor designs. Variant design options of GCMR and HPMR were also investigated to provide a wider application range for each technology. Criticality safety analyses for the GCMR and HPMR packages were performed using the CSAS6 sequence of SCALE 6.3.2 with the ENDF/B-VII.1-based continuous energy neutron libraries. Different scenarios were investigated, including normal operation and water flooded conditions to represent nominal and hypothetical accident scenarios. Sensitivity and similarity analyses are also performed using the TSUNAMI sequence of SCALE 6.3.2, and the similarity analysis uses all the experiments from the ICSBEP Handbook with High Enriched Uranium (HEU), Intermediate and mixed Enriched Uranium (IEU) and Low Enriched Uranium (LEU) systems. Many ICSBEP experiments were found marginally similar, with similarity index (ck) values greater than 0.8, to the GCMR and HPMR cask models, especially in flooded conditions that may be more constraining due to reduced reactivity margins. Some experiments using TRISO fuel particles and graphite moderator, including the recent THETA and Deimos experiments from LANL, together with some IRPhEP experiments, were also used for the similarity analyses. These TRISO-fueled and graphite moderated experiments behaved similarly to the GCMR models, with maximum ck value greater than 0.9, but they are less similar to the HPMR models. This observation indicates that additional critical experiments might be needed to further validate the criticality safety models, especially for heat-pipe based microreactor transport packages and for non-flooded transportation configurations.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Draft Impacts: Modeling and Lab Validation (CRADA Final Report)

Chimney draft is a significant source of variability in real-world wood heater performance, yet laboratory tests conducted for certification do not capture this variability. Because draft varies with climate, chimney height, and home conditions, a heater that performs well in the lab can perform quite differently when installed in a home. Until now, wood heater manufacturers and installers have lacked the tools needed to anticipate how draft will vary across real installations or to advise on corrective measures when draft is too low or too high. The goal of this project was to develop and validate an open-source draft prediction tool for cordwood heater chimney systems. Over 18 months, Lawrence Berkeley National Laboratory (LBNL) and the Hearth, Patio & Barbecue Association (HPBA) collaborated to build a physics-based, Python tool that predicts how chimney draft evolves during operation, from ignition through steady burning. HPBA convened a stakeholder group of manufacturers and venting experts who provided feedback throughout development. LBNL validated the tool against laboratory measurements from a catalytic and a non-catalytic cordwood heater operated across a range of chimney heights and room-pressure conditions.

42 ENGINEERING↗

CMS Endcap Timing Layer: System Validation and Assembly

The High-Luminosity Large Hadron Collider (HL-LHC) will operate at unprecedented luminosities, resulting in up to 200 simultaneous interactions per bunch crossing. To mitigate the resulting pileup challenges, the CMS experiment is implementing the MIP Timing Detector (MTD), designed to provide precise timing information with a resolution of 30 – 40 picoseconds per track. The MTD consists of the Barrel Timing Layer (BTL) and the Endcap Timing Layer (ETL), each optimized for different regions of the detector. The ETL, comprising two double-sided disks, utilizes Low-Gain Avalanche Diode (LGAD) sensors coupled with the Endcap Timing Readout Chip (ETROC) to achieve high-precision timing measurements in the forward region. Significant progress has been made toward the realization of the ETL through extensive system-level validation of the ETROC readout chain and the development of scalable module assembly procedures. This presentation will provide a comprehensive overview of the ETL and highlight recent advances in system validation and module assembly, emphasizing their roles in ensuring the detector readiness for HL-LHC operation.

Lee, Dongyub [Kyungpook Natl. U.] (ORCID:000000034↗

Simulation Validation for the DUNE Near Detector

DUNE will study neutrino oscillations as they are beamed from Fermilab in Batavia, Illinois, 1300 km away to Sanford Underground Research Facility in South Dakota. However, when preparing for an experiment of this scale, it is crucial that we run simulations to validate the geometry and overall design of the project. This poster is comprised of my work during a 2026 SULI internship at Fermilab, involving reading the simulated data and assessing its validity.

Vershaw, Andre [Unlisted, US, IL] (ORCID:000900084↗

Additional Data Focused on Phase 1 Geared Toward Computational Fluid Dynamics (CFD) validation

A new validation campaign was developed within the Offshore Code Comparison Collaboration, Continued, with Correlation and unCertainty (OC6) to better understand the complex interactions between components of a floating wind system (e.g., columns, pontoons, etc.) in a more generic sense, and to validate modeling approaches for a variety of floating designs. A set of cylinders based on the Offshore Code Comparison Collaboration Continuation (OC4)-DeepCwind semidesign were tested under wave loading, both individually and in different combinations (including the attachment of heave plates of different sizes). The cylinders were held rigid, and pressure sensors were used to measure the distributed loading on the structures, while load cells measured the total integrated hydrodynamic loads.

17 WIND ENERGY↗

Validation of the Community Land Model Version 5 over the Contiguous United States (CONUS) using in-situ and remote sensing datasets

The Community Land Model (CLM) is an effective tool to simulate the biophysical and biogeochemical processes and their interactions with the atmosphere. Although CLM Version 5 (CLM5) constitutes various updates in these processes, its performance in simulating energy, water and carbon cycles over the Contiguous United States (CONUS) at scales which land surface changes and hydrometeorological and hydroclimatological applications are more locally relevant is yet to be assessed. In this study, we conducted three simulations at 0.125? during 1979-2018 over the CONUS using different configurations of CLM, namely CLM5-biogeochemistry (CLM5BGC), CLM4.5BGC, and CLM5-satellite phenology (CLM5SP). We validated and compared their simulations against multiple remote-sensed and in-situ datasets. Overall, the parametric and structural updates (e.g., carbon cost for nitrogen uptake, variable soil thickness, dry surface layer) in CLM5 improve its ability in capturing terrestrial biogeochemical dynamics. The low evapotranspiration in CLM5BGC is associated with biases in simulating vegetation phenological characteristics rather than soil water limitations. The mismatch between CLM5BGC-simulated peak leaf area index and reference data can be attributed to CLM5BGC's inability in simulating phenology of trees and grasses. The differences between CLM-simulated irrigation and reference estimates can be attributed to differences between processes represented in models and in reality, and uncertainties in input and validation datasets. Evaluation against observations at small catchments suggest that hydrologic parameters needed to be calibrated to improve simulations of runoff, especially subsurface runoff. Additional efforts are needed to incorporate spatially-distributed plant phenology and physiology parameters and regional-specific agricultural management practices (e.g., planting, harvest).

54 ENVIRONMENTAL SCIENCES↗

Development and Validation of a Process Model and Open-Source Process Simulator for Microalgae-Based Tertiary Phosphorus Recovery

Microalgae-based tertiary wastewater treatment has the potential to meet stringent effluent phosphorus limits, with the added benefit of producing a marketable feedstock. However, the lack of validated mechanistic models and their implementation in process simulators have limited the adoption of this technology. In this study, an updated lumped pathway metabolic model (Phototrophic-Mixotrophic Process Model, PM 2 ), including both photoautotrophic and heterotrophic metabolisms of microalgae, was developed to predict effluent phosphorus concentration and biomass yield in response to dynamic influent and varying environmental conditions. The model was implemented in QSDsan – an open-source, Python-based design and simulation platform – for robust simulation under uncertainty. A global sensitivity analysis was performed to prioritize model parameters for calibration. The model was then calibrated and validated using batch experimental data and 45 days of continuous online monitoring data from a full-scale (568 m 3 ·d -1 ) microalgae-based tertiary wastewater treatment plant (EcoRecover process). In particular, along with dynamic influent composition, temperature and light intensity data with diel variation were provided as model inputs to reflect the microalgal behavior under day-night cycling. Overall, the QSDsan-based microalgae process simulator was able to predict effluent phosphorus within 0.02–0.04 mg-P·L -1 , while also capturing the general trends of state variables according to nutrient availability.

Lumped pathway metabolic model↗

Simulation-Based Validation of An Open-Source, Scalable Framework for Building Energy Management in Small and Medium-Sized Commercial Buildings

Abstract: Small and medium-sized commercial buildings (SMCBs) represent 94% of U.S. commercial buildings but encounter substantial obstacles in adopting Building Energy Management (BEM) systems. Current approaches exhibit fundamental limitations: vendor-specific API platforms restrict interoperability through proprietary ecosystems; commercial automation software demands extensive technical expertise and licensing costs; open-source IoT solutions lack native support for building automation protocols and semantic models. This paper introduces a configuration-driven web interface framework addressing the gap between smart device advancements and accessible BEM software infrastructure for SMCBs. The framework leverages VOLTTRON middleware integrated with an automated converter that processes unified YAML configurations into heterogeneous system files, reducing required configuration artifacts from six separate files to a single unified specification. The system architecture enables vendor-agnostic operation through BACnet and Modbus protocols while supporting semantic building model integration via automated Brick Schema parsing. Configuration-driven interfaces automatically adapt to diverse HVAC types without custom development. Simulation-based validation using BOPTEST demonstrates automatic interface generation between fan coil and hydronic systems, with the automated converter successfully generating all platform-specific outputs from the single YAML input. The result demonstrates the framework's capability to streamline BEM system deployment through reduced configuration complexity. This work bridges simulation capabilities with operational deployment, demonstrating how virtual testbeds validate generalizable software frameworks for real-world building automation.

Chung, Jihoon [ORNL] (ORCID:0000000184880815)↗

Validation of a Proteomic Signature of Lung Cancer Risk from Bronchial Specimens of Risk-Stratified Individuals

A major challenge in lung cancer prevention and cure hinges on identifying the at-risk population that ultimately develops lung cancer. Previously, we reported proteomic alterations in the cytologically normal bronchial epithelial cells collected from the bronchial brushings of individuals at risk for lung cancer. The purpose of this study is to validate, in an independent cohort, a selected list of 55 candidate proteins associated with risk for lung cancer with sensitive targeted proteomics using selected reaction monitoring (SRM). Bronchial brushings collected from individuals at low and high risk for developing lung cancer as well as patients with lung cancer, from both a subset of the original cohort (batch 1: n = 10 per group) and an independent cohort of 149 individuals (batch 2: low risk (n = 32), high risk (n = 34), and lung cancer (n = 83)), were analyzed using multiplexed SRM assays. ALDH3A1 and AKR1B10 were found to be consistently overexpressed in the high-risk group in both batch 1 and batch 2 brushing specimens as well as in the biopsies of batch 1. Validation of highly discriminatory proteins and metabolic enzymes by SRM in a larger independent cohort supported their use to identify patients at high risk for developing lung cancer.

60 APPLIED LIFE SCIENCES↗

Analyzing the Validity of Brazilian Testing Using Digital Image Correlation and Numerical Simulation Techniques

Characterizing the mechanical behavior of rocks plays a crucial role to optimize the fracturing process in unconventional reservoirs. However, due to the intrinsic anisotropy and heterogeneity in unconventional resources, fracture process prediction remains the most significant challenge for sustainable and economic hydrocarbon production. During the deformation tracking under compression, deploying conventional methods (strain gauge, extensometer, etc.) is insufficient to measure the deformation since the physical attachment of the device is restricted to the size of the sample, monitoring limited point-wise deformation, producing difficulties in data retrieval, and a tendency to lose track in failure points, etc. Where conventional methods are limited, the application of digital image correlation (DIC) provides detailed and additional information of strain evolution and fracture patterns under loading. DIC is an image-based optical method that records an object with a camera and monitors the random contrast speckle pattern painted on the facing surface of the specimen. To overcome the existing limitations, this paper presents numerical modeling of Brazilian disc tests under quasi-static conditions to understand the full-field deformation behaviors and finally, it is validated by DIC. As the direct tensile test has limitations in sample preparation and test execution, the Brazilian testing principle is commonly used to evaluate indirectly the tensile strength of rocks. The two-dimensional numerical model was built to predict the stress distribution and full-field deformation on Brazilian disc under compression based on the assumptions of a homogenous, isotropic and linear elastic material. The uniaxial compression test was conducted using the DIC technique to determine the elastic properties of Spider Berea sandstone, which were used as inputs for the simulation model. The model was verified by the analytical solution and compared with the digital image correlation. The numerical simulation results showed that the solutions matched reasonably with the analytical solutions where the maximum deviation of stress distribution was obtained as 14.59%. The strain evolution (normal and shear strains) and displacements along the central horizontal and vertical planes were investigated in three distinguishable percentages of peak loads (20%, 40%, and 90%) to understand the deformation behaviors in rock. The simulation results demonstrated that the strain evolution contours consistently matched with DIC generated contours with a reasonable agreement. The changes in displacement along the central horizontal and vertical planes showed that numerical simulation and DIC generated experimental results were repeatable and matched closely. In terms of validation, Brazilian testing to measure the indirect tensile strength of rocks is still an issue of debate. The numerical model of fracture propagation supported by digital image correlation from this study can be used to explain the fracturing process in the homogeneous material and can be extended to non-homogeneous cases by incorporating heterogeneity, which is essential for rock mechanics field applications.

58 GEOSCIENCES↗

Highly Accurate Experimental Heave Decay Tests with a Floating Sphere: A Public Benchmark Dataset for Model Validation of Fluid–Structure Interaction

Highly accurate and precise heave decay tests on a sphere with a diameter of 300 mm were completed in a meticulously designed test setup in the wave basin in the Ocean and Coastal Engineering Laboratory at Aalborg University, Denmark. The tests were dedicated to providing a rigorous benchmark dataset for numerical model validation. The sphere was ballasted to half submergence, thereby floating with the waterline at the equator when at rest in calm water. Heave decay tests were conducted, wherein the sphere was held stationary and dropped from three drop heights: a small drop height, which can be considered a linear case, a moderately nonlinear case, and a highly nonlinear case with a drop height from a position where the whole sphere was initially above the water. The precision of the heave decay time series was calculated from random and systematic standard uncertainties. At a 95% confidence level, uncertainties were found to be very low—on average only about 0.3% of the respective drop heights. Physical parameters of the test setup and associated uncertainties were quantified. A test case was formulated that closely represents the physical tests, enabling the reader to do his/her own numerical tests. The paper includes a comparison of the physical test results to the results from several independent numerical models based on linear potential flow, fully nonlinear potential flow, and the Reynolds-averaged Navier–Stokes (RANS) equations. A high correlation between physical and numerical test results is shown. The physical test results are very suitable for numerical model validation and are public as a benchmark dataset.

42 ENGINEERING↗

Distributed Energy Management for Networked Microgrids with Hardware-in-the-Loop Validation

For the cooperative operation of networked microgrids, a distributed energy management considering network operational objectives and constraints is proposed in this work. Considering various ownership and privacy requirements of microgrids, utility directly interfaced distributed energy resources (DERs) and demand response, a distributed optimization is proposed for obtaining optimal network operational objectives with constraints satisfied through iteratively updated price signals. The alternating direction method of multipliers (ADMM) algorithm is utilized to solve the formulated distributed optimization. The proposed distributed energy management provides microgrids, utility-directly interfaced DERs and responsive demands the opportunity of contributing to better network operational objectives while preserving their privacy and autonomy. Results of numerical simulation using a networked microgrids system consisting of several microgrids, utility directly interfaced DERs and responsive demands validate the soundness and accuracy of the proposed distributed energy management. The proposed method is further tested on a practical two-microgrid system located in Adjuntas, Puerto Rico, and the applicability of the proposed strategy is validated through hardware-in-the-loop (HIL) testing.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Modeling Sepsis: Establishment and Validation of a 72-Hour Swine Model of Penetrating Abdominal Trauma

Background/Objectives: Fecal peritonitis following penetrating abdominal trauma is a serious condition that often results in sepsis and organ failure. The aim of our study was to develop a novel conscious porcine model of sepsis and organ dysfunction caused by multiple penetrating injuries to the small and large intestines. Methods: Twelve female Yorkshire pigs (average weight 50.6 ± 6.5 kg) were divided into two groups: Penetrating Abdominal Trauma (PAT) (n = 8) and Control (n = 4). All surgical procedures were performed under anesthesia with adequate analgesia. In the PAT group, the small and large intestines were punctured, and feces mixed with saline were introduced into the abdominal cavity to induce peritonitis. The Control group received sham surgery with only saline solution. The animals were observed in a conscious state over a period of 72 h, vital parameters were recorded, and blood samples were taken regularly. We adapted a pig-specific SOFA score and developed pig-specific SIRS criteria and NEWS2 score to assess organ function. The model was validated by independent investigators. Results: The survival rate in the PAT group was 75%, with an average survival time of 58.5 h, while all animals in the Control group survived to euthanasia. Monitoring showed pathophysiological changes, such as tachycardia, leucopenia, and thrombocytopenia, indicative of sepsis and organ dysfunction. Blinded investigators independently confirmed the model’s validity. Conclusions: A new swine model of penetrating abdominal trauma and sepsis has been successfully developed that demonstrates significant physiological and immunologic changes comparable to human sepsis. This new model provides a realistic platform for future research into sepsis, its diagnostics, and the evaluation of therapeutic strategies.

General & Internal Medicine↗

Atmospheric Transmittance Model Validation for CSP Tower Plants

In yield analysis and plant design of concentrated solar power (CSP) tower plants, increased uncertainties are caused by the mostly unknown solar attenuation between the concentrating heliostat field and the receiver on top of the tower. This attenuation is caused mainly by aerosol particles and water vapor. Various on-site measurement methods of atmospheric extinction in solar tower plants have been developed during recent years, but during resource assessment for distinct tower plant projects in-situ measurement data sets are typically not available. To overcome this lack of information, a transmittance model (TM) has been previously developed and enhanced by the authors to derive the atmospheric transmittance between a heliostat and receiver on the basis of common direct normal irradiance (DNI), temperature, relative humidity and barometric pressure measurements. Previously the model was only tested at one site. In this manuscript, the enhanced TM is validated for three sites (CIEMAT’s Plataforma Solar de Almería (PSA), Spain, Missour, Morocco (MIS) and Zagora, Morocco (ZAG)). As the strongest assumption in the TM is the vertical aerosol particle profile, three different approaches to describe the vertical profile are tested in the TM. One approach assumes a homogeneous aerosol profile up to 1 kilometer above ground, the second approach is based on LIVAS profiles obtained from Lidar measurements and the third approach uses boundary layer height (BLH) data of the European Centre for Medium-Range Weather Forecasts (ECMWF). The derived broadband transmittance for a slant range of 1 km (T 1$km$ ) time series is compared with a reference data set of on-site absorption- and broadband corrected T 1$km$ derived from meteorological optical range (MOR) measurements for the temporal period between January 2015 and November 2017. The absolute mean bias error (MBE) for the TM’s T 1$km$ using the three different aerosol profiles lies below 5% except for ZAG and one profile assumption. The MBE is close to 0 for PSA and MIS assuming a homogeneous extinction coefficient up to 1 km above ground. The root mean square error (RMSE) is around 5–6% for PSA and ZAG and around 7–8% for MIS. The TM performs better during summer months, during which more data points have been evaluated. This validation proves the applicability of the transmittance model for resource assessment at various sites. It enables the identification of a clear site with high T 1$km$ with a high accuracy and provides an estimation of the T 1$km$ for hazy sites. Thus it facilitates the decision if on-site extinction measurements are necessary. The model can be used to improve the accuracy of yield analysis of tower plants and allows the site adapted design.

14 SOLAR ENERGY↗