Search NASA⌕ Search

SEARCH · Search NASA

Results for “Initialization method”

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
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Determining the N -Representability of a Reduced Density Matrix via Unitary Evolution and Stochastic Sampling

The N-representability problem consists in determining whether, for a given p-body matrix, there exists at least one N-body density matrix from which the p-body matrix can be obtained by contraction, that is, if the given matrix is a p-body reduced density matrix (p-RDM). The knowledge of all necessary and sufficient conditions for a p-body matrix to be N-representable allows the constrained minimization of a many-body Hamiltonian expectation value with respect to the p-body density matrix and, thus, the determination of its exact ground state. However, the number of constraints that complete the N-representability conditions grows exponentially with system size, and hence, the procedure quickly becomes intractable for practical applications. This work introduces a hybrid quantum-stochastic algorithm to effectively replace the N-representability conditions. The algorithm consists of applying to an initial N-body density matrix a sequence of unitary evolution operators constructed from a stochastic process that successively approaches the reduced state of the density matrix on a p-body subsystem, represented by a p-RDM, to a target p-body matrix, potentially a p-RDM. The generators of the evolution operators follow the well-known adaptive derivative-assembled pseudo-Trotter method (ADAPT), while the stochastic component is implemented by using a simulated annealing process. The resulting algorithm is independent of any underlying Hamiltonian, and it can be used to decide whether a given p-body matrix is N-representable, establishing a criterion to determine its quality and correcting it. We apply the proposed hybrid ADAPT algorithm to alleged reduced density matrices from a quantum chemistry electronic Hamiltonian, from the reduced Bardeen–Cooper–Schrieffer model with constant pairing, and from the Heisenberg XXZ spin model. In all cases, the proposed method behaves as expected for 1-RDMs and 2-RDMs, evolving the initial matrices toward different targets.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Elastic Changepoint Detection for Globally-indexed Functional Time Series Data with Climate Applications

Changepoint detection is a vital tool in the application of climate data analysis. Numerous types of climate observation data are most properly represented by functional time series, implying a need for accurate changepoint detection methods applicable to functional time series data. Such data taken at a global scale often contain both spatial heterogeneity and dependence as well as phase (time) misalignment. In this report, we present methods which can detect spatially-dependent changepoints while allowing different estimates of change time and change strength depending on location. Additionally, we provide extensions to this spatially-predicted model which controls for phase variability among observations. Our methods provide the ability to detect a single change, or control for epidemic changes (where a “return-to-normal” change is more likely to be detected than the initial change). We showcase results analyzing the June 1991 eruption of Mt. Pinatubo, where our methods demonstrate the ability to accurately detect both single and epidemic changepoints even in the presence of strong seasonal variability. We find that our spatially-predicted model improves the detection of relevant changepoints versus methods which do not take spatial information into account, and we find that controlling for phase variability helps to control the false discovery rate during the detection process.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Application of physics-informed neural networks (PINNs) solution to coupled thermal and hydraulic processes in silty sands

Abstract The accurate modeling of water and heat transport in soils is crucial for both geo-environmental and geothermal engineering. Traditional modeling methods are problematic because they require well-defined boundaries and initial conditions. Recently, physics-informed neural networks (PINNs), which incorporate partial differential equations (PDEs) to solve forward and inverse problems, have attracted increasing attention in machine learning research. In this study, we applied PINNs to tackle hydraulic and thermal transport coupling forward problems in silty sands. A fully connected deep neural network was utilized for training. This neural network model leverages automatic differentiation to apply the governing equations as constraints, based on the mathematical approximations established by the neural network itself. We conducted forward problems and compared the solutions derived from PINNs with those from Finite Element Method (FEM) simulations. The forward problem results demonstrate the PINNs model’s capability in predicting hydraulic transport, heat transport, and thermal–hydraulic coupling in silty sands under various boundary conditions. The PINNs exhibited great performance in simulating the thermal–hydraulic coupling problem. The accuracy of the PINNs solutions shows its potential for simulation in geotechnical engineering.

Feng, Yuan↗
DOE OSTI2026

A time-parallel method for scalable heat transfer simulations of additive manufacturing

Here, a major challenge in simulating the thermal behavior in additive manufacturing processes is the disparate length and time scales between transport phenomena occurring in the melt pool and the component. A common simulation approach relies on spatial decomposition for parallel computing, but due to the nature of heat transfer in AM, where most of the computational expenditure is localized near the melt pool, the computational speedup from spatial parallelization saturates quickly. Therefore, additional parallelism by means of time-domain decomposition is needed to fully take advantage of high-performance computing (HPC) resources. This work introduces a time-parallel method to improve the computational scalability of additive manufacturing simulations on HPC systems, while maintaining high temporal resolution of heat transfer near the melt pool. The method, inspired by the nonlinear paraexp formalism, performs an iterative superposition of nonlinear solutions to the initial value problem, integrating the heat equation across overlapping time-parallel intervals. For a single layer of the NIST AMB2018–01 L7 benchmark problem, the method achieves a 38.51x speedup in wall-clock time with a maximum error in the global temperature solution of 0.99%. This reduces the total solution time from 196.72 min to 5.11 min on 128 nodes of the ORNL Frontier supercomputer. The tradeoff between accuracy and total wall-clock time is investigated and recommendations for time-parallel deployment for AM problems are made.

Additive manufacturing↗
DOE OSTI2024

Fast permeability measurement for tight reservoir cores using only initial data of the one chamber pressure pulse decay test

Here, in this study, a mathematical model for fast determination of the permeabilities of tight rocks using measurements taken from the initial period of the One Chamber Pressure Pulse Decay (OC-PPD) test is presented. The model applies to measurements taken both before and after the pressure pulse front has reached the downstream end of the specimen. The analytical solutions for the pressure decay in the upstream chamber are derived based on a parabolic arc approximation of pore pressure distribution along the test specimen. This approximation allows converting the initial–boundary value problem of fluid diffusion in the specimen, governed by partial differential equations, to a system of ordinary differential equations that can be easily solved by explicit formulae. Thus, an explicit formula for the pressure decay rate is obtained, which enables inverse analysis of the initial experimental data to estimate the rock permeability. The proposed method expedites the pulse decay test as it does not require the system to reach equilibrium. The method is validated with three sets of experimental data of the OC-PPD test using helium as the diffusing fluid, for which the relative error of the permeability is found to be less than 6%. This method is particularly useful if the equilibrium time of the pulse decay test for rock specimens with permeabilities in the range of nano-Darcy takes hours or days.

early-time solution↗
DOE OSTI2025

A theoretical study of the upper bound of surface elevation variance in the Phillips initial stage during wind-wave generation

The resonance mechanism in the initial of wind-wave generation proposed by Phillips is a foundation of wind-wave generation theory, but a precise theoretical quantification of wave energy growth in this initial stage has not been obtained yet after more than six decades of research. In this study, we aim to address this knowledge gap by developing an analytical approach based on a novel complex analysis method to theoretically investigate the temporal evolution of the wave energy in the Phillips initial stage. We quantitatively derive and analyse the growth behaviour of the surface wave energy and obtain an analytical solution for its upper bound. Our result highlights the crucial effects of surface tension. Because the phase velocity of gravity–capillary waves has a minimal value at a critical wavenumber, gravity–capillary waves and gravity waves (which neglect surface tension) exhibit distinct resonance curve properties and wave energy growth behaviours. For gravity waves, the resonance curve extends indefinitely; for gravity–capillary waves, it either forms a finite-length curve or does not exist, depending on the wind speed. The leading-order term of the upper-bound solution of the energy of gravity waves increases linearly over time, while for gravity–capillary waves, the term increases linearly over time under strong wind conditions but remains finite under weak wind conditions. This theoretical study provides an analytical framework for the generation of wind-waves in the Phillips initial stage, which may inspire further theoretical, numerical and experimental research.

58 GEOSCIENCES↗
DOE OSTI2025

Supramolecular Support of Cuprophilic Network Bonding in 2-D Copper n -Alkanethiolates

Here, the development of heterogeneous materials, catalysts, and semiconductors is often reliant on precise control of self-assembly and crystal packing. Many new materials are initially synthesized as microcrystalline powders, making them incompatible with typical methods of structure determination, such as single-crystal X-ray diffraction. This resultant lack of structural information has made thorough investigation into the effect of metal substitution on crystal structure in metal-organic chalcogenolates (MOChas) challenging. Here, we use small molecule serial femtosecond crystallography (smSFX) to present the structures of four copper n-alkanethiolates: CuSC4, CuSC5, CuSC6, and CuSC7. Divergent patterns of alkyl chain packing are identified from microcrystalline powders via smSFX. An odd-even effect in crystal packing has been identified and attributed to different orientations of symmetry elements in the even- and odd-numbered chains. This results in minute changes in the azimuthal organization of the even-numbered chains and the network of cuprophilic interactions. Additionally, we present a synthesis of crystalline gold n-alkanethiolates to provide the first comparison between three d 10 coinage metals (Cu, Ag, and Au) and their resultant n-alkanethiolates.

Willson, Maggie C. [Univ. of Connecticut, Storrs, ↗
DOE OSTI2024

Imaging shapes of atomic nuclei in high-energy nuclear collisions

Atomic nuclei are self-organized, many-body quantum systems bound by strong nuclear forces within femtometre-scale space. These complex systems manifest a variety of shapes, traditionally explored using non-invasive spectroscopic techniques at low energies. However, at these energies, their instantaneous shapes are obscured by long-timescale quantum fluctuations, making direct observation challenging. Here we introduce the collective-flow-assisted nuclear shape-imaging method, which images the nuclear global shape by colliding them at ultrarelativistic speeds and analysing the collective response of outgoing debris. This technique captures a collision-specific snapshot of the spatial matter distribution within the nuclei, which, through the hydrodynamic expansion, imprints patterns on the particle momentum distribution observed in detectors. We benchmark this method in collisions of ground-state uranium-238 nuclei, known for their elongated, axial-symmetric shape. Our findings show a large deformation with a slight deviation from axial symmetry in the nuclear ground state, aligning broadly with previous low-energy experiments. This approach offers a new method for imaging nuclear shapes, enhances our understanding of the initial conditions in high-energy collisions and addresses the important issue of nuclear structure evolution across energy scales.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2026

Accelerating catalytic advancements through the precision of high-throughput experiments & calculations

The growing demand for energy-efficient processes to support a sustainable future drives the need for research to rapidly explore chemical and material space through accelerated catalyst discovery initiatives. Recent breakthroughs in high-throughput experimental and computational methods are transforming the catalysis field, surpassing traditional approaches to manipulating variables in catalytic processes. Key advancements in innovation include the integration of machine learning for efficient catalyst screening, high-throughput experimentation, data-driven methodologies employing comprehensive databases, and in situ and in operando techniques for realistic observations. This progress has undoubtedly been intertwined with a collaborative framework across disciplines, reshaping catalyst discovery methods in both industry and academia. This Opinion article presents a multifaceted perspective from coauthors with expertise spanning various stages of the Technology Readiness Level spectrum, highlighting both opportunities and persistent challenges in integrating computational and experimental approaches in catalysis. These challenges span from obtaining high-quality experimental data, scaling simulations to industrially relevant materials and process conditions to navigating the complexity and predictive accuracy of computational models.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Stochastic Optimization to Find Optimum Beginning-of-Life Core Configuration of Stable Salt Reactor with Online Refueling

A stochastic optimization method has been developed to find an optimum equilibrium cycle core configuration of the waste-burning stable salt reactor, which is a fast-spectrum molten salt reactor with frequent online refueling. An optimum core configuration was determined with the goal of minimizing radial power peaking. Because of the vast number of potential candidate core configurations, stochastic optimization was applied based on simulated annealing and an additional acceleration method, which screened out unpromising core configurations. It has been demonstrated that the developed stochastic optimization method successfully finds the optimal core configuration regardless of the initial guess and outperforms the gradient descent approach. In addition, it has been observed that the use of a so-called out-in core configuration as the initial guess speeds up convergence of the iterative solution more than five times. Based on the searched optimum equilibrium cycle core configuration, new beginning-of-life (BOL) core configurations have been developed. In conclusion, the new BOL core configurations will be used in developing optimum refueling strategies.

Moltex static salt reactor↗
DOE OSTI2025

Measurement of photonuclear jet production in ultraperipheral Pb + Pb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV with the ATLAS detector

In ultrarelativistic heavy ion collisions at the LHC, each nucleus acts a sources of high-energy real photons that can scatter off the opposing nucleus in ultraperipheral photonuclear (𝛾 + 𝐴) collisions. Hard scattering processes initiated by the photons in such collisions provide a novel method for probing nuclear parton distributions in a kinematic region not easily accessible to other measurements. ATLAS has measured production of dijet and multijet final states in ultraperipheral Pb + Pb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV using a dataset recorded in 2018 with an integrated luminosity of 1.72 nb −1 . Photonuclear final states are selected by requiring a rapidity gap in the photon direction; this selects events where one of the outgoing nuclei remains intact. Jets are reconstructed using the anti-𝑘 t algorithm with radius parameter, 𝑅 = 0.4. Triple-differential cross sections, unfolded for detector response, are measured and presented using two sets of kinematic variables. The first set consists of the total transverse momentum (𝐻 T ), rapidity, and mass of the jet system. The second set uses 𝐻 T and particle-level nuclear and photon parton momentum fractions, 𝑥 A and 𝑧 𝛾 , respectively. The results are compared with leading-order perturbative QCD calculations of photonuclear jet production cross sections, where all leading order predictions using existing fits fall below the data in the shadowing region. More detailed theoretical comparisons will allow these results to strongly constrain nuclear parton distributions, and these data provide results from the LHC directly comparable to early physics results at the planned Electron-Ion Collider.

Parton distribution functions↗
DOE OSTI2026

An Adaptive Multiparameter Penalty Selection Method for Multiconstraint and Multiblock ADMM

This work presents a new method for online selection of multiple penalty parameters for the alternating direction method of multipliers (ADMM) algorithm applied to optimization problems with multiple constraints or functions with block matrix components. ADMM is widely used for solving constrained optimization problems in a variety of fields, including signal and image processing. Implementations of ADMM often utilize a single hyperparameter, referred to as the penalty parameter, which needs to be tuned to control the rate of convergence. However, in problems with multiple constraints, ADMM may demonstrate slow convergence regardless of penalty parameter selection due to scale differences between constraints. Accounting for scale differences between constraints to improve convergence in these cases requires introducing a penalty parameter for each constraint. The proposed method is able to adaptively account for differences in scale between constraints, providing robustness with respect to problem transformations and initial selection of penalty parameters. It is also simple to understand and implement. Our numerical experiments demonstrate that the proposed method performs favorably compared to a variety of existing penalty parameter selection methods.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

SARS-CoV-2 wastewater variant surveillance: pandemic response leveraging FDA’s GenomeTrakr network

ABSTRACT Wastewater surveillance has emerged as a crucial public health tool for population-level pathogen surveillance. Supported by funding from the American Rescue Plan Act of 2021, the FDA‘s genomic epidemiology program, GenomeTrakr, was leveraged to sequence SARS-CoV-2 from wastewater sites across the United States. This initiative required the evaluation, optimization, development, and publication of new methods and analytical tools spanning sample collection through variant analyses. Version-controlled protocols for each step of the process were developed and published on protocols.io. A custom data analysis tool and a publicly accessible dashboard were built to facilitate real-time visualization of the collected data, focusing on the relative abundance of SARS-CoV-2 variants and sub-lineages across different samples and sites throughout the project. From September 2021 through June 2023, a total of 3,389 wastewater samples were collected, with 2,517 undergoing sequencing and submission to NCBI under the umbrella BioProject,PRJNA757291. Sequence data were released with explicit quality control (QC) tags on all sequence records, communicating our confidence in the quality of data. Variant analysis revealed wide circulation of Delta in the fall of 2021 and captured the sweep of Omicron and subsequent diversification of this lineage through the end of the sampling period. This project successfully achieved two important goals for the FDA’s GenomeTrakr program: first, contributing timely genomic data for the SARS-CoV-2 pandemic response, and second, establishing both capacity and best practices for culture-independent, population-level environmental surveillance for other pathogens of interest to the FDA. IMPORTANCE This paper serves two primary objectives. First, it summarizes the genomic and contextual data collected during a Covid-19 pandemic response project, which utilized the FDA’s laboratory network, traditionally employed for sequencing foodborne pathogens, for sequencing SARS-CoV-2 from wastewater samples. Second, it outlines best practices for gathering and organizing population-level next generation sequencing (NGS) data collected for culture-free, surveillance of pathogens sourced from environmental samples.

Microbiology↗
DOE OSTI2026

Ion Trapping Studies and Mitigation Strategies for the EIC ERL-Based Strong Hadron Cooler

An Energy Recovery Linac based strong hadron cooler was previously considered for the Electron-Ion Collider. The required electron beam parameters for variable-energy strong hadron cooling place significant constraints on ion trapping and collective effects. This paper presents initial studies of these constraints through a combination of analytical modelling and numerical simulations of ion production, trapping behaviour, and mitigation strategies. A multi-bunch tracking framework based on ELEGANT with the ionEffects module is used to simulate machine operation over millisecond time scales, corresponding to more than 3 × 10^5 electron bunches. The simulations include modelling of ionisation processes together with transverse electron–ion dynamics, allowing the evolution and accumulation of ions to be investigated. Analytical expressions based on Gaussian beam distributions are used to estimate ion trapping conditions and benchmark the simulation results. A bi-periodic bunch spacing scheme is also investigated as a possible mitigation method by detuning the ion oscillation frequency. These studies provide an initial assessment of ion trapping in the strong hadron cooler and demonstrate possible approaches for reducing beam–ion effects.

Bi, R. [Lancaster University, Cockcroft Institute]↗
DOE OSTI2025

A New Vehicle-to-Vehicle Communication System: Visual-Enhanced Cooperative Traffic Operations

The advent of Connected and Autonomous Vehicles (CAVs) has highlighted the necessity for robust communication systems between vehicles and their environment. This study introduces a novel vehicle-to-vehicle (V2V) communication system, termed the Visual-Enhanced Cooperative Traffic Operations (VECTOR) system. The VECTOR system addresses the need for robust communication by converting dynamic data (including velocity and yaw angle data) into binary code, which is displayed on an LED panel mounted on the top of the vehicle. Following vehicles detect this panel and decode the information using a camera, implementing a visual-based communication method. VECTOR system employs a comprehensive five-module process. Initially, polynomial fitting techniques are applied to velocity data over fixed time intervals using third-degree polynomials, with validation via R² and MSE metrics. The second module converts velocity and yaw angle data into binary form, thereby enhancing detection and processing efficiency. The third module focuses on improving detection stability across various environmental conditions to enhance traffic safety. The fourth module decodes the binary data back into trajectory information, ensuring the fidelity of velocity and yaw angles. The final module integrates eco-control through the VECTOR system, employing advanced control algorithms to minimize energy consumption in CAVs. Experimental evaluations conducted using a modified CAV test platform based on the Lincoln MKZ demonstrate the feasibility and efficiency of the VECTOR system, achieving a 75% R-squared accuracy rate in replicating original velocity data. This methodology not only highlights potential applications but also underscores significant implications for advancing CAV technology.

Ma, Ke↗
DOE OSTI2025

Standard Library Grid-Forming Hybrid Control Inverter-Based Resource Model Specification (REGFM_C1)

This document describes a standard library grid-forming (GFM) hybrid control inverter-based resource (IBR) model. The GFM hybrid control approach implements both a typical GFM control and a typical grid-following (GFL) control inside one single inverter simultaneously, so that it can take advantage of both methods without comprising the benefits of a typical GFM. The initial version of model specification was jointly developed by Pacific Northwest National Laboratory (PNNL), Tesla Energy, and EPRI, and it was revised multiple times later to incorporate suggestions from WECC MVS members. Tesla Energy provided main control blocks to support the development of this model specification. This standard library model is developed to help the utility industry better understand the GFM technology. The model could be used to represent equipment for long-term planning studies where vendor-specific models are not available. As equipment matures and improves, standard library models will be updated to capture the new functionalities of GFMs. It is not intended that these models will always remain representative of all future GFM technologies.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2024

Ferrocene Bis(Sulfonate) Salt as Redoxmer for Fast and Steady Redox Flow Desalination

Desalination is considered a promising solution to alleviate water shortages, yet current methods are often restricted, due to challenges like high energy consumption, significant cost, or limited desalination capacity. In this study, we present a novel approach of redox flow desalination (RFD) utilizing the highly aqueous-soluble and reversible redox-active compound, potassium 1,1'-bis(sulfonate) ferrocene (1,1'-FcDS). This water-soluble organic compound yielded stable and rapid desalination, sustaining extended operation without notable decay and achieving an impressive desalination rate of up to 457.5 mmol·h -1 ·m -2 and energy consumption as low as 40.2 kJ·mol NaCl -1 . Specifically, the RFD device effectively desalinated a 50 mM NaCl solution to potable standards within 6000 s using 1,1'-FcDS. It maintained an average energy consumption of 178.16 kJ·mol NaCl -1 and exhibited negligible deterioration in desalination rate, energy efficiency, and charge efficiency throughout a rigorous 12,000 s cycling test. Furthermore, the versatility of this method was demonstrated by effectively treating saline water with varying initial concentrations from 10 mM to 50 mM, showcasing its potential across a broad spectrum of applications.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Advanced Perovskite Solar Cells and Modules

The “Advanced perovskite Cells and Modules” research project was the final agreement focused on enhancing perovskite solar cell (PSC) technologies funded by the US Department of Energy's Solar Energy Technologies Office. The project was designed to address three crucial areas in PSC development: stability, manufacturability, and efficiency. The project was then structured around three main tasks, each targeting one of these strategic goals. The team of experienced researchers in these materials worked collaboratively to address the targets outlined in the technical work plan. building on existing PSC research while also exploring promising new concepts arising in the field. An overview of each primary task is summarized below: Task 1 Stability: This first task, aims to identify material characteristics and metrics that can help predict the primary degradation mechanisms impacting PSC stability. This involved developing specific device tests based on hypotheses regarding mechanisms impacting stability, including fast failure procedures to speed up PSC development and improvement. Various strategies to enhance stability, like incorporating additives, post-treatments, novel contact materials etc. were developed using this fast feedback approach. The relationships between indoor and outdoor stresses were also validated. Task 2 Manufacturability: This second task, focused on creating a scalable production process for PSCs. Initially the objective is to establish a best-known method for a 182 cm2 minimodule. However, given resource limitations, these metrics were modified to focus on the other goal of outlined in the TWP. Specifically, this task worked to demonstrate the transferability of this best-known method to another research institution. Work scope in this area was expanded to material purity and understanding of reagent/process relationships. Examination of other difficulties in PSC production and potential solutions for large-scale production were also evaluated. Given challenges observed in process transfer, work to develop data infrastructure and recording tools for processing of material and devices was then also prioritized in this task. Task 3 Efficiency: This task was focused on improvements to PCE, while still considering Task 1 and Task 2 goal. The efforts targeted a PCE greater than 22% with a T95 exceeding 1000 hours at 25°C in a nitrogen environment for lab-scale devices (approximately 0.1 cm2 devices) across a range of solar-relevant perovskite compositions, including wide-gap (around 1.7 eV) and low-gap (around 1.3 eV) materials, using standard metal contacts. This work then provides a foundation for MHP-based tandem efforts undertaken in other projects and the All-MHP tandem efforts outlined in this projects TWP. Work in this project emphasized disseminating its findings through peer-reviewed publications (PRP), conference presentations, and industrial collaborations. Significant products were produced in all these areas, over 53 peer reviewed publications, 32 conference presentations and industrial investment based on NLR assistance on precompetitive challenges. The team also developed significant intellectual property and awards for their technical excellence, innovations and leadership. The team also leveraged traditional and social media platforms to engage with stakeholders and the public.

14 SOLAR ENERGY↗