Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “continuous reconstruction”

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 73 records · Page 4

Continuous scanning mode for ptychography

We outline how ptychographic imaging can be performed without the need for discrete scan positions. Through an idealized experiment, we demonstrate how a discrete-position scan regime can be replaced with a continuously scanned one with suitable modification of the reconstruction scheme based on coherent modes. The impact of this is that acquisition times can be reduced, significantly aiding ptychographic imaging with x rays, electrons, or visible light.

Clark, Jesse N.↗

Improving reproducibility in synchrotron tomography using implementation-adapted filters

For reconstructing large tomographic datasets fast, filtered backprojection-type or Fourier-based algorithms are still the method of choice, as they have been for decades. These robust and computationally efficient algorithms have been integrated in a broad range of software packages. The continuous mathematical formulas used for image reconstruction in such algorithms are unambiguous. However, variations in discretization and interpolation result in quantitative differences between reconstructed images, and corresponding segmentations, obtained from different software. This hinders reproducibility of experimental results, making it difficult to ensure that results and conclusions from experiments can be reproduced at different facilities or using different software. In this paper, a way to reduce such differences by optimizing the filter used in analytical algorithms is proposed. These filters can be computed using a wrapper routine around a black-box implementation of a reconstruction algorithm, and lead to quantitatively similar reconstructions. Use cases for this approach are demonstrated by computing implementation-adapted filters for several open-source implementations and applying them to simulated phantoms and real-world data acquired at the synchrotron. Our contribution to a reproducible reconstruction step forms a building block towards a fully reproducible synchrotron tomography data processing pipeline.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Attitude and translation control of a low-altitude Gravsat

Dynamic simulation of the Gravsat's attitude and translation control system is used to provide an upper bound for the fuel supply and test the feasibility of the preliminary design. A preliminary design is made for the disturbance compensation system (DISCOS) sensor, the thruster control laws, reaction wheel control laws, and the onboard state estimators. The sensor analysis and noise measurements show no problems in scaling the Triad navigation satellite sensor design up to meet the Gravsat requirements, except for proof mass center-of-mass offset. A promising technique is proposed to measure and eliminate this error. The covariance analysis confirms that a sophisticated post-flight data fit will be necessary to reconstruct a scientifically useful proof mass state. The DISCOS sensor will have to be continuously calibrated from the inflight data to achieve this reconstruction.

Ray, J. C.↗

ED-cPSD: Fast Phase-Size Distribution via Sequential Erosion-Dilation

The Erosion-Dilation continuous Phase-Size Distribution, ED-cPSD, is an application for calculating continuous pore and particle-size distribution from digital reconstructions and/or image-based structural data. It is based on the erosion-dilation continuous phase-size distribution method. A continuous size distribution is a measure of the probability density of finding a particle or pore of a certain size. These distributions are of interest in any field of study involving porous media, including but not limited to electrochemistry, petroleum engineering, geology, and food science. The algorithm behind the software provides a computationally efficient way to calculate phase-size distributions for large domains. For a 3D battery electrode reconstruction with 1.3 x 10 8 voxels, the particle size distribution is derived in under 2 min on a desktop, while also retaining flexibility and computational efficiency for HPC-scale multi-threading. The software can handle structures with over 10 9 voxels. The algorithm is roughly 280 times faster than a previous version on the same task.

Characterization↗

The Design of Mechanically Compatible Fasteners for Human Mandible Reconstruction

Mechanically compatible fasteners for use with thin or weakened bone sections in the human mandible are being developed to help reduce large strain discontinuities across the bone/implant interface. Materials being considered for these fasteners are a polyetherertherketone (PEEK) resin with continuous quartz or carbon fiber for the screw. The screws were designed to have a shear strength equivalent to that of compact/trabecular bone and to be used with a conventional nut, nut plate, or an expandable shank/blind nut made of a ceramic filled polymer. Physical and finite element models of the mandible were developed in order to help select the best material fastener design. The models replicate the softer inner core of trabecular bone and the hard outer shell of compact bone. The inner core of the physical model consisted of an expanding foam and the hard outer shell consisted of ceramic particles in an epoxy matrix. This model has some of the cutting and drilling attributes of bone and may be appropriate as an educational tool for surgeons and medical students. The finite element model was exercised to establish boundary conditions consistent with the stress profiles associated with mandible bite forces and muscle loads. Work is continuing to compare stress/strain profiles of a reconstructed mandible with the results from the finite element model. When optimized, these design and fastening techniques may be applicable, not only to other skeletal structures, but to any composite structure.

Roberts, Jack C.↗

The scaling of complex craters

The empirical relation between the transient crater diameter (Dg) and final crater diameter (Dr) of complex craters and basins is estimated using cumulative terrace widths, central uplift diameters, continuous ejecta radii, and transient crater reconstructions determined from lunar and terrestrial impact structures. The ratio Dg/Dr is a power law function of Dr, decreasing uniformly from unity at the diameter of the simple-complex crater morphology transition to about 0.5 for large multiring basins like Imbrium on the moon. The empirical constants in the Dg/Dr relation are interpreted physically to mean that the position of the final rim relative to the transient crater, and hence the extent of collapse, is controlled or greatly influenced by the properties of the zone of dissociated material produced by the impact shock. The continuity of the Dg/Dr relation over the entire spectrum of morphologic types from complex craters to multiring basins implies that the rims of all these structures form in the same tectonic environment despite morphologic differences.

Croft, S. K.↗

Black hole and galaxy coevolution from continuity equation and abundance matching

We investigate the coevolution of galaxies and hosted supermassive black holes (BHs) throughout the history of the universe by a statistical approach based on the continuity equation and the abundance matching technique. Specifically, we present analytical solutions of the continuity equation without source terms to reconstruct the supermassive BH mass function from the active galactic nucleus (AGN) luminosity functions. Such an approach includes physically motivated AGN light curves tested on independent data sets, which describe the evolution of the Eddington ratio and radiative efficiency from slim- to thin-disk conditions. We nicely reproduce the local estimates of the BH mass function, the AGN duty cycle as a function of mass and redshift, along with the Eddington ratio function and the fraction of galaxies with given stellar mass hosting an AGN with given Eddington ratio. We exploit the same approach to reconstruct the observed stellar mass function at different redshift from the ultraviolet and far-IR luminosity functions associated with star formation in galaxies. These results imply that the build-up of stars and BHs in galaxies occurs via in situ processes, with dry mergers playing a marginal role at least for stellar masses ≲3×10{sup 11} M{sub ⊙} and BH masses ≲10{sup 9} M{sub ⊙}, where the statistical data are more secure and less biased by systematic errors. In addition, we develop an improved abundance matching technique to link the stellar and BH content of galaxies to the gravitationally dominant dark matter (DM) component. The resulting relationships constitute a testbed for galaxy evolution models, highlighting the complementary role of stellar and AGN feedback in the star formation process. In addition, they may be operationally implemented in numerical simulations to populate DM halos or to gauge subgrid physics. Moreover, they may be exploited to investigate the galaxy/AGN clustering as a function of redshift, mass, and/or luminosity. In fact, the clustering properties of BHs and galaxies are found to be in full agreement with current observations, thus further validating our results from the continuity equation. Finally, our analysis highlights that (i) the fraction of AGNs observed in the slim-disk regime, where most of the BH mass is accreted, increases with redshift; and (ii) already at z≳6 a substantial amount of dust must have formed over timescales ≲10{sup 8} yr in strongly star-forming galaxies, making these sources well within the reach of ALMA surveys in (sub)millimeter bands.

79 ASTRONOMY AND ASTROPHYSICS↗

SU-G-IeP4-03: Cone Beam X-Ray Luminescence Computed Tomography Based On Generalized Gaussian Markov Random Field

Purpose: Cone beam X-ray luminescence computed tomography (CB-XLCT), which aims to achieve molecular and functional imaging by X-rays, has recently been proposed as a new imaging modality. However, the inverse problem of CB-XLCT is seriously ill-conditioned, hindering us to achieve good image quality. In this work, a novel reconstruction method based on Bayesian theory is proposed to tackle this problem Methods: Bayesian theory provides a natural framework for utilizing various kinds of available prior information to improve the reconstruction image quality. A generalized Gaussian Markov random field (GGMRF) model is proposed here to construct the prior model of the Bayesian theory. The most important feature of GGMRF model is the adjustable shape parameter p, which can be continuously adjusted from 1 to 2. The reconstruction image tends to have more edge-preserving property when p is slide to 1, while having more noise tolerance property when p is slide to 2, just like the behavior of L1 and L2 regularization methods, respectively. The proposed method provides a flexible regularization framework to adapt to a wide range of applications. Results: Numerical simulations were implemented to test the performance of the proposed method. The Digimouse atlas were employed to construct a three-dimensional mouse model, and two small cylinders were placed inside to serve as the targets. Reconstruction results show that the proposed method tends to obtain better spatial resolution with a smaller shape parameter, while better signal-to-noise image with a larger shape parameter. Quantitative indexes, contrast-to-noise ratio (CNR) and full-width at half-maximum (FWHM), were used to assess the performance of the proposed method, and confirmed its effectiveness in CB-XLCT reconstruction. Conclusion: A novel reconstruction method for CB-XLCT is proposed based on GGMRF model, which enables an adjustable performance tradeoff between L1 and L2 regularization methods. Numerical simulations were conducted to demonstrate its performance.

60 APPLIED LIFE SCIENCES↗

Effects of Catalyst Processing on the Activity and Stability of Pt–Ni Nanoframe Electrocatalysts

Pt-based alloys have shown great promise as cathodic catalysts for cost-effective proton exchange membrane fuel cells. Post-synthesis treatment has been recognized as a critical step to improve the catalytic performance of Pt-based alloys. Here we present the effects of catalyst processing on the catalytic behavior of Pt-Ni nanoframe electrocatalysts in oxygen reduction reaction. The Pt-Ni nanoframes were made by corroding Ni-rich phase from solid rhombic dodecahedral particles. Three different corrosion procedures were compared. Among them, electrochemical corrosion led to the highest initial specific activity (1.35 mA cm -2 at 0.95 VRHE) by retaining more Ni in the nanoframes. However, the high activity gradually went down in a subsequent stability test due to continuous Ni loss and concomitant surface reconstruction. On the other hand, the best stability was achieved by a more aggressive corrosion using oxidative nitric acid. Although the initial activity was compromised, this procedure imparted a less-defective surface and thus the specific activity dropped by only 7% over 30,000 potential cycles. These results indicate a delicate trade-off between the activity and stability of Pt-Ni nanoframe electrocatalysts. The obtained understandings on how to balance the activity-stability trade-off via catalyst processing can be generalized to other Pt-based alloys.

36 MATERIALS SCIENCE↗

Bayesian inference of electron density and ion temperature profiles from neutral beam and halo Balmer- α emission at Wendelstein 7-X

By employing Bayesian inference techniques, the full electron density profile from the plasma core to the edge of Wendelstein 7-X (W7-X) is inferred solely from neutral hydrogen beam and halo Balmer- α (H α ) emission data. The halo is a cloud of neutrals forming in the vicinity of the injected neutral beam due to multiple charge exchange reactions. W7-X is equipped with several neutral hydrogen beam heating sources and an H α spectroscopy system that views these sources from different angles and penetration depths in the plasma. As the beam and halo emission form complex spectra for each spatial point that are non-linearly dependent on the plasma density profile and other parameters, a complete model from the neutral beam injection and halo formation through to the spectroscopic measurements is required. The model is used here to infer electron density profiles for a range of common W7-X plasma scenarios. The inferred profiles show good agreement with profiles determined by the Thomson scattering and interferometry diagnostics across a broad range of absolute densities without any changes to the input or fitting parameters. The time evolution of the density profile in a discharge with continuous core density peaking is successfully reconstructed, demonstrating sufficient spatial resolution to infer strongly shaped profiles. Furthermore, it is shown as a proof of concept that the model is also able to infer the main ion temperature profile using the same data set.

forward modeling↗

EJFAT Scientific Perspective

Presented new computing model to the test by deploying the EJFAT system alongside a data-stream processing framework running the production-level CLAS12 event reconstruction application. In this experiment, a continuous stream of CLAS12 Level-1 identified events was processed in real-time using the EJFAT load balancer, distributing the workload across 90 computing nodes located across the U.S. This marks the first-ever large-scale, real-time distributed data stream processing experiment, demonstrating that scientific data-streaming pipelines can efficiently scale across four dimensions, thanks to EJFAT’s advanced hardware and software capabilities.

Gyurjyan, Vardan [Thomas Jefferson National Accele↗

Forty-Year Calibrated Record of Earth-Surface Reflected Radiance from Landsat: A Review

Sensors on Landsat satellites have been collecting images of the Earth's surface for nearly 40 years. These images have been invaluable for characterizing and detecting changes in the land cover and land use of the world. Although initially conceived as primarily picture generating sensors, even the early sensors were radiometrically calibrated and spectrally characterized prior to launch and incorporated some capabilities to monitor their radiometric calibration once on orbit. Recently, as the focus of studies has shifted to monitoring Earth surface parameters over significant periods of time, serious attention has been focused toward bringing the data from all these sensors onto a common radiometric scale over this 40-year period. This effort started with the most recent systems and then was extended back in time. Landsat-7 ETM+, the best-characterized sensor of the series prior to launch and once on orbit, and the most stable system to date, was chosen to serve as the reference. The Landsat-7 project was the first of the series to build an image assessment system into its ground system, allowing systematic characterization of its sensors and data. Second, the Landsat-5 TM (still operating at the time of the Landsat-7 launch and continues to operate) calibration history was reconstructed based on its internal calibrator, vicarious calibrations, pseudo-invariant sites and a tie to Landsat-7 ETM+ at the time of the commissioning of Landsat-7. This process was performed in two iterations: the earlier one relied primarily on the TM internal calibrator. When this was found to have some deficiencies, a revised calibration was based more on pseudo-invariant sites, though the internal calibrator was still used to establish the short-term variations in response due to icing build up on the cold focal plane. As time progressed, a capability to monitor the Landsat-5 TM was added to the image assessment system. The Landsat-4 TM, which operated from 1982-1992, was the third system to which the radiometric scale was extended. The limited and broken use of the Landsat-4 TM made this analysis more difficult. Eight-day separated image pairs from Landsat-5 combined with analysis of pseudo invariant sites established this history. The fourth and most challenging effort was making the Landsat-1 to -5 MSS sensors' data internally radiometrically consistent. This effort was particularly complicated by the age of the MSS data, varying formats and processing levels in the archive, limited datasets, and limited documentation available. Ultimately, pseudo-invariant sites were identified in North America and used for this effort. Note that most of the Landsat-MSS archived data had already been calibrated using the MSS internal calibrators, so this processing was imbedded in the result. The final effort was developing an absolute scale for Landsat MSS similar to what was already established for the "TM" sensors. This was accomplished by using simultaneous data from Landsat-5 MSS and Landsat-5 TM, accounting for spectral differences between the sensors using EO-1 Hyperion data. The recalibrated history of the Landsat data and implications to users are discussed. The key result from this work is a consistently calibrated Landsat data archive that spans nearly 40 years with total uncertainties on the order of 10% or less for most sensors and bands.

Markham, Brian↗

High-resolution wave-theory-based ultrasound reflection imaging using the split-step fourier and globally optimized fourier finite-difference methods

Methods for enhancing ultrasonic reflection imaging are taught utilizing a split-step Fourier propagator in which the reconstruction is based on recursive inward continuation of ultrasonic wavefields in the frequency-space and frequency-wave number domains. The inward continuation within each extrapolation interval consists of two steps. In the first step, a phase-shift term is applied to the data in the frequency-wave number domain for propagation in a reference medium. The second step consists of applying another phase-shift term to data in the frequency-space domain to approximately compensate for ultrasonic scattering effects of heterogeneities within the tissue being imaged (e.g., breast tissue). Results from various data input to the method indicate significant improvements are provided in both image quality and resolution.

47 OTHER INSTRUMENTATION↗

Kinetic equilibrium reconstruction and the impact on stability analysis of KSTAR plasmas

High fidelity kinetic plasma equilibrium reconstructions are an essential requirement for accurate stability and disruption prediction analyses to support continuous operation of high beta tokamak plasmas. The present kinetic equilibrium reconstructions of plasmas in the KSTAR device include plasma density and temperature profiles from Thomson scattering and ion temperature from charge exchange spectroscopy diagnostics, and allowance for fast particle pressure. In addition, up to 25 channels of motional Stark effect diagnostic data are used to constrain the magnetic field pitch angle profile in the plasma to produce a reliable computation of the safety factor, q, profile. In this work, H-mode plasmas exhibit clear pedestal characteristics in the reconstructed pressure profile compared to internal transport barrier or L-mode plasmas. The plasma configuration and vertical position of inner strike points are validated by CCD and infrared camera images. Ideal and resistive magnetohydrodynamic (MHD) stability analyses using the DCON and resistive DCON codes utilize these kinetic equilibrium reconstructions to compare to experimental plasma stability. Equilibria with sufficiently low convergence error can provide reliable computation of ideal and resistive MHD stability analysis.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Conditional deep generative models for simultaneous simulation and reconstruction of entire events

We extend the particle-flow neural assisted simulations (arnassus) framework of fast simulation and reconstruction to entire collider events. In particular, we use two generative artificial intelligence tools, continuous normalizing flows and diffusion models, to create a set of reconstructed particle-flow objects conditioned on truth-level particles from CMS Open Simulations. While previous work focused on jets, our updated methods now can accommodate all particle-flow objects in an event along with particle-level attributes like particle type and production vertex coordinates. This approach is fully automated, entirely written in Python, and GPU-compatible. Using a variety of physics processes at the LHC, we show that the extended arnassus is able to generalize beyond the training dataset and outperforms the standard, public tool elphes.

Dreyer, Etienne [Weizmann Institute of Science, Re↗

On the Reconstruction of Palaeo-Ice Sheets: Recent Advances and Future Challenges

Reconstructing the growth and decay of palaeo-ice sheets is critical to understanding mechanisms of global climate change and associated sea-level fluctuations in the past, present and future. The significance of palaeo-ice sheets is further underlined by the broad range of disciplines concerned with reconstructing their behaviour, many of which have undergone a rapid expansion since the 1980s. In particular, there has been a major increase in the size and qualitative diversity of empirical data used to reconstruct and date ice sheets, and major improvements in our ability to simulate their dynamics in numerical ice sheet models. These developments have made it increasingly necessary to forge interdisciplinary links between sub-disciplines and to link numerical modelling with observations and dating of proxy records. The aim of this paper is to evaluate recent developments in the methods used to reconstruct ice sheets and outline some key challenges that remain, with an emphasis on how future work might integrate terrestrial and marine evidence together with numerical modelling. Our focus is on pan-ice sheet reconstructions of the last deglaciation, but regional case studies are used to illustrate methodological achievements, challenges and opportunities. Whilst various disciplines have made important progress in our understanding of ice-sheet dynamics, it is clear that data-model integration remains under-used, and that uncertainties remain poorly quantified in both empirically-based and numerical ice-sheet reconstructions. The representation of past climate will continue to be the largest source of uncertainty for numerical modelling. As such, palaeo-observations are critical to constrain and validate modelling. State-of-the-art numerical models will continue to improve both in model resolution and in the breadth of inclusion of relevant processes, thereby enabling more accurate and more direct comparison with the increasing range of palaeo-observations. Thus, the capability is developing to use all relevant palaeo-records to more strongly constrain deglacial (and to a lesser extent pre-LGM) ice sheet evolution. In working towards that goal, the accurate representation of uncertainties is required for both constraint data and model outputs. Close cooperation between modelling and data-gathering communities is essential to ensure this capability is realised and continues to progress.

Ice sheet reconstruction↗

The potential of positron emission tomography for intratreatment dynamic lung tumor tracking: A phantom study

Purpose: This study aims to evaluate the potential and feasibility of positron emission tomography for dynamic lung tumor tracking during radiation treatment. The authors propose a center of mass (CoM) tumor tracking algorithm using gated-PET images combined with a respiratory monitor and investigate the geometric accuracy of the proposed algorithm. Methods: The proposed PET dynamic lung tumor tracking algorithm estimated the target position information through the CoM of the segmented target volume on gated PET images reconstructed from accumulated coincidence events. The information was continuously updated throughout a scan based on the assumption that real-time processing was supported (actual processing time at each frame ≈10 s). External respiratory motion and list-mode PET data were acquired from a phantom programmed to move with measured respiratory traces (external respiratory motion and internal target motion) from human subjects, for which the ground truth target position was known as a function of time. The phantom was cylindrical with six hollow sphere targets (10, 13, 17, 22, 28, and 37 mm in diameter). The measured respiratory traces consisted of two sets: (1) 1D-measured motion from ten healthy volunteers and (2) 3D-measured motion from four lung cancer patients. The authors evaluated the geometric accuracy of the proposed algorithm by quantifying estimation errors (Euclidean distance) between the actual motion of targets (1D-motion and 3D-motion traces) and CoM trajectories estimated by the proposed algorithm as a function of time. Results: The time-averaged error of 1D-motion traces over all trajectories of all targets was 1.6 mm. The error trajectories decreased with time as coincidence events were accumulated. The overall error trajectory of 1D-motion traces converged to within 2 mm in approximately 90 s. As expected, more accurate results were obtained for larger targets. For example, for the 37 mm target, the average error over all 1D-motion traces was 1.1 mm; and for the 10 mm target, the average error over all 1D-motion traces was 2.8 mm. The overall time-averaged error of 3D-motion traces was 1.6 mm, which was comparable to that of the 1D-motion traces. There were small variations in the errors between the 3D-motion traces, although the motion trajectories were very different. The accuracy of the estimates was consistent for all targets except for the smallest. Conclusions: The authors developed an algorithm for dynamic lung tumor tracking using list-mode PET data and a respiratory motion signal, and demonstrated proof-of-principle for PET-guided lung tumor tracking. The overall tracking error in phantom studies is less than 2 mm.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Safer Foundation Solar Demands Skill Collaborative

The milestones and accomplishments achieved in Safer’s Solar Demands Skill Training program required a collective effort of stakeholders. The ability to achieve the desired results and impact lives of the population we serve works best when collaborating with community stakeholders such as community-based organizations, faith-based organizations, community activists, local law enforcement, community residents, employers, and elected and appointed officials. We are working harder than ever to place clients in-demand careers and high-growth industry sectors. Our employer engagement with high growth sectors continues to increase, we have been able to deepen our relationships in this space by providing in demand stackable credential training. Launched our initiative in the green job space by partnering with organizations like ComEd, the Department of Energy, community-based organizations, faith-based organizations, and manufacturers that helped to train clients in the green jobs and renewable energy space. By the end of the year, we will have completed our tenth cohort of photovoltaic solar installation training. The training has been used as a foundation for solar panel installers to acquire new skills in in the green jobs career pathways, such as sales and customer service. It has also served as a pipeline to union-level trade positions through the skills participants obtained through partnership engagement. Safer Foundation's policy and advocacy team plays a major role in expanding opportunities for people with records to over one hundred occupations, including the trades. This reform has allowed many to secure living wage employment, reducing the high recidivism rate within Illinois. Safer continues to build our social enterprise with Reconstructive Technology Partners (RTP). RTP is introducing people with records to the construction trades such as solar, carpentry and electrical through residential remodeling of homes on smaller construction projects. Understanding the need to have a greater community presence, outreach was extended to a boots on the ground model, including, but not limited to door-to-door engagement, DE-EE0008571 Safer Foundation Page 4 of 29 local radio broadcasting and print ads in local newspapers. We continue to seek innovative ways to increase Safer Foundations' presence in the renewable energy space. Outreach efforts gave us an opportunity to share critical reentry information and opportunities about our PV installer program. Giving us reach with local and national audiences alike. We launched a systematic approach to improving client data and client tracking through an evidence-based practices initiative. That implemented an agency wide cross-functional data management system. The system supported our goal of reviewing, evaluating, and making recommendations to improve operational processes, program delivery, information sharing, and more. Played a significant part in the success of the solar program. Finally, economists around the nation agree that there is a significant labor shortage. The shortage directly threatens our ability to sustain our economic growth; if employers do not have access to the workers they need, it can lead to a shutdown in the economic recovery. The demand for Safer Foundation services is more important and impactful than ever. We will continue to build upon our 50 years of experience to address the challenges ahead and serve more people in a better way. We are confident we will accomplish beautiful things that benefit everyone involved because together, we are powerful.

14 SOLAR ENERGY↗