Search NASA⌕ Search

SEARCH · Search NASA

Results for “Linear filters”

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.

132 records · Page 5
All sourcesNASA NTRSDOE OSTI
NASA NTRSDate not supplied

Softening the Gap between Wöhler and Paris – New Approaches for Fatigue Analysis –

Fatigue analysis tools can vary across industries. For example, automotive engineers often use the Wöhler (S-N) approach to design for safe-life, while aerospace engineers prioritize damage tolerance and inspection intervals, relying instead on crack growth models such as Paris’ law. Although both approaches may deal with the control of cracks in similar materials, their analysis tools and material characterizations are fundamentally distinct. This divide mirrors the classic split between stress-based strength analysis and linear elastic fracture mechanics. However, modern nonlinear models that incorporate material softening, such as cohesive laws, blur this boundary and capture fracture behaviors across scales. This presentation describes the CF23 fatigue model, which uses cohesive softening to link S-N crack initiation with crack propagation rates. CF23 spans the full fatigue spectrum, from initial propagation transients to steady-state growth and threshold conditions, offering a unified framework that bridges Wöhler and Paris-based methodologies. Example applications include fatigue crack propagation transients in adhesive interfaces and skin/stiffener separation.

cohesive elements↗
DOE OSTI2024

The cluster decomposition of the configurational energy of multicomponent alloys

Abstract The cluster expansion method (CEM) is a widely used lattice-based technique in the study of multicomponent alloys. Despite its prevalent use, a clear understanding of expansion terms is lacking. We present a modern mathematical formalism of the CEM and introduce thecluster decomposition—a unique and basis-independent decomposition for functions of the atomic configuration in a crystal. We identify the cluster decomposition as an invariant ANOVA decomposition; and demonstrate how functional analysis of variance and sensitivity analysis can be used to interpret interactions among species. Furthermore, we show how the mathematical structure of the cluster decomposition enables numerical evaluation that scales with the number of clusters and is independent of the number of species. Overall, our work enables rigorous interpretations of interactions among species, provides opportunities to explore parameter estimation beyond linear regression, introduces a numerical efficient implementation, and enables analysis of cluster expansions based on established mathematical and statistical principles.

Chemistry↗
DOE OSTI2025

The role of fast and slow dynamics in nonlinear resonant ultrasound spectroscopy of consolidated granular materials

Abstract Elastic nonlinearity observed in consolidated granular media can be attributed to the combination of slow and fast effects, which give rise to hysteresis and relaxation of both modulus and damping after the sample is perturbed. A consequence is a high level of complexity in the measurements of the sample linear and nonlinear elastic parameters. The results of experiments are dependent on the experimental protocol that is adopted to measure the relevant quantities and it is hard to quantify parameters with accuracy and repeatability. Here we focus on examining Nonlinear Resonant Ultrasound Spectroscopy, showing experimentally the role of slow dynamics in the process and quantifying/discussing its influence on the quantification of nonlinearity. We also propose a model to describe the process, which shows that different contributions to nonlinearity (e.g., classical and hysteretic) could be due to physical features (defects) relaxing with different relaxation times.

Science & Technology - Other Topics↗
NASA NTRS2026

Parabolic Deformable Mirror Project Development and Device Characterization

Existing state-of-the-art wavefront sensing and control setups for the imaging of exoplanets utilize two high-actuator count flat deformable mirrors (DMs) to control both phase and amplitude aberrations. One of the objectives of the Exoplanet Spectroscopy (ExoSpec) project at Goddard Space Flight Center (GSFC) is to improve upon this regime through the development of parabolic deformable mirrors (PDMs), which can address several of the key limitations in existing setups. Specifically, the need for flat DMs to be located in conjugate pupil planes drives unique packaging challenges that can be overcome by allowing the off-axis imaging elements to be controllable as well. In addition, simulations show that the use of deformable powered elements will increase the overall controllable bandwidth, up to 35% over 5-12 lambda/D. To pursue this effort, the ExoSpec group has begun testing a first-generation parabolic DM, which is already informing a second generation design. This paper will provide a comprehensive update on its characterization, including influence function measurements, linearity analysis, closed-loop wavefront control performance, and thermal stability assessment.

spectroscopy↗
NASA NTRSDate not supplied

Stationary Testing of a Pulse Tube Cryocooler Designed to Rotate at 6800 RPM with the Superconducting Rotor of an Electric Machine

Electrified aircraft are being developed to increase the efficiency and reduce the cost of operating subsonic transport aircraft. Achieving a substantial impact necessitates focusing on single- and twin-aisle aircraft which use propulsion systems with >20 MW ratings. For these aircraft, multi-MW superconducting electric machines are being developed due to their high specific power and high efficiency. Most of these electric machines employ superconductors on their rotor, thereby requiring cryogenic cooling of the rotating system. An attractive solution is to conductively cool the rotor using a cryocooler that rotates with the rotor. NASA has been developing a 50 K pulse tube cryocooler for its High Efficiency Megawatt Motor (HEMM) that can operate while rotating at 6,800 rpm. The design of this cryocooler and testing of its linear motor have been discussed in prior work. This presentation presents the first measurements of the complete HEMM cryocooler. The design and results of a stationary (non-rotating) test of the cryocooler are described.

electrified aircraft propulsion↗
DOE OSTI2024

System for controller area network payload decoding

A system for decoding an unknown automotive controller area network (“CAN”) message definitions. CAN data vehicle signal mappings are typically held in secret and varied by automotive model and year. Without knowledge of the mappings, the wealth of real-time vehicle data hidden in the automotive CAN packets is uninterpretable—impeding research, after-market tuning, efficiency and performance monitoring, fault diagnosis, and privacy-related technologies. This system can ascertain the CAN signals' boundaries (start bit and length), endianness (byte ordering), signedness (binary-to-integer encoding) from raw CAN data. This allows conversion of CAN data to time series. Interpreting the translated CAN data's physical meaning and finding a linear mapping to standard units (e.g., knowing the signal is speed and scaling values to represent units of miles per hour) can be achieved for many signals by leveraging diagnostic standards to obtain real-time measurements of in-vehicle systems. The system can be integrated into lightweight hardware enabling an OBD-II plugin for real-time in-vehicle CAN decoding or run on standard computers. The system can output a standard DBC file with the signal definition information.

Verma, Kiren E.↗
NASA NTRS1985

Summary of recent NASA propeller research

Advanced high speed propellers offer large performance improvements for aircraft that cruise in the Mach 0.7 to 0.8 speed regime. At these speeds, studies indicate that there is a 15 to near 40 percent block fuel savings and associated operating cost benefits for advanced turboprops compared to equivalent technology turbofan powered aircraft. Recent wind tunnel results for five eight to ten blade advanced models are compared with analytical predictions. Test results show that blade sweep was important in achieving net efficiencies near 80 percent at Mach 0.8 and reducing nearfield cruise noise about 6 dB. Lifting line and lifting surface aerodynamic analysis codes are under development and some results are compared with propeller force and probe data. Also, analytical predictions are compared with some initial laser velocimeter measurements of the flow field velocities of an eight bladed 45 swept propeller. Experimental aeroelastic results indicate that cascade effects and blade sweep strongly affect propeller aeroelastic characteristics. Comparisons of propeller nearfield noise data with linear acoustic theory indicate that the theory adequately predicts nearfield noise for subsonic tip speeds, but overpredicts the noise for supersonic tip speeds.

Daniel C Mikkelson↗
DOE OSTI2026

Toward Fully Soft and Multifunctional Shape Sensing via Optical Waveguide Arrays

For soft bodies, surface deformation and pressure provide proprioceptive and exteroceptive information, including body configuration, body compliance, and external forces. We develop a sheet sensor with a fully soft sensing surface that provides surface shape reconstruction using optical waveguide arrays. The waveguides are fabricated to achieve a tunable linear response to bi‐directional bending curvature, and the waveguide arrays are configured to differentiate between ambiguous shapes. We characterize the waveguide performance, relating curvature sensitivity to the core's surface roughness. Synergy of waveguide responses reduces the number of sensing elements required and achieves damage resilience. Using waveguide sensitivity to pressure, we also demonstrate feasibility for exteroception. We demonstrate the multifunctional sensing capability by wrapping the sheet around an upper arm, showcasing joint motion and external force sensing. Integrated into or applied onto surfaces of robotic or living systems, this design can be implemented in applications such as virtual reality, teleoperation, physical therapy, and soft robotics.

Yu, Qifan [Department of Mechanical Engineering Ma↗
DOE OSTI2024

Isotropic negative thermal expansion in ZrW 2 O 8 and HfW 2 O 8 from 1100 to 1275°C

Abstract ZrW 2 O 8 (ZrO 2 •2WO 3 ) and HfW 2 O 8 (HfO 2 •2WO 3 ) have been the focus of thermal expansion studies due to their isotropic negative thermal expansion (NTE) measured previously at temperatures below 775°C. This work presents measurements of these materials at their thermodynamically stable temperature ranges of 1105 and 1257°C for ZrW 2 O 8 and 1105–1276°C for HfW 2 O 8 , where they were characterized with in situ, powder X‐ray diffraction. The linear coefficients of thermal expansion were measured to be −5.52 × 10 −6 and −4.87 × 10 −6 °C −1 for ZrW 2 O 8 and HfW 2 O 8 , respectively. The mechanism leading to this NTE is discussed. Powder samples were synthesized by a solution‐based process called the organic–inorganic steric entrapment method. In situ characterization in air was carried out at the National Synchrotron Light Source II using a hexapole lamp, optical furnace and the Advanced Photon Source using a quadrupole lamp, optical furnace to achieve elevated temperatures.

Materials Science↗
NASA NTRSDate not supplied

Preliminary Assessment of the Impact on the V-Band Oxygen Channels From Satellite Communication Uplinks

We calculate the percentage of time that an ATMS-like instrument [1] will be illuminated by the uplink beam of one of the proposed V-band communication system and estimate the damage resulting from such exposure. Using a combination of openly available information and educated guesses about the location and characteristics of the up/down link terminals, we constructed the ground segment of a hypothetical high-speed communication network. The space segment of the network was constructed from the orbital data of the existing Starlink constellation [2] of 6223 communication satellites (comsats) which is used as strawman to represent any other possible constellation of communication satellites. It is shown that without a very delicate balance of frequency allocations (science vs telecommunications), coupled with extremely steep and deep bandpass-defining filters, and strict adherence to the agreed limits (i.e. no out-of-band transmissions) the deployment of the telecommunication network leads to almost-complete loss of some important geophysical data. For the analysis we use the spectral characteristics of the ATMS instrument with the ephemeris for the NOAA-21 satellite [3]. The analysis is conducted for the USA and the simulation covers 8 consecutive days in July 2024. Effective and accurate microwave remote sensing of the atmosphere depends on the availability of interference-free spectrum windows at frequencies which are prescribed by physical processes [e.g. 4]. The family of resonant lines of the oxygen molecule near 60 GHz provides a unique opportunity to sample the vertical distribution of temperature and density from space, and it has been exploited for weather and climate studies from polar-orbiting satellites since 1978 (MSU on TIROS-N [5]). It remains a staple in the payloads operated by Russia, China, USA, Japan, France, India, UK, Ukraine [6] which are built around a common blueprint: a few wide-band (hundreds of MHz) channels around 50 GHz to sample the atmosphere and the surface while several more channels with high spectral resolution (few MHz) sample the individual resonant lines. Accurate retrieval of the environmental parameters depends upon the data provided by both sets of channels, and the their location in frequency space is not arbitrary and cannot be altered at will [7, 8]. The introduction of 5G technology in 2019 has driven telecommunication companies to request more bandwidth to be dedicated to their devices. This additional bandwidth is only available in spectral regions traditionally reserved for environmental and astrophysical research, such as the V-band between 50 and 60 GHz for up/downlink between satellites in low-earth orbits and terminals connected to fiberoptics network for distribution to high-speed local internet services. The power broadcast by the uplink communication leg is many orders of magnitude greater than the natural thermal signal emitted from the Earth scene. If the ground antenna were to perfectly align with the passive instrument’s antenna, the spaceborne receiver would suffer permanent, irreparable damage. While a direct boresight-to-boresight conjunction is extremely unlikely (even with a large constellation of satellites the fraction of the celestial sphere occupied by the satellites remains minuscule) the finite size of the ground station’s antenna beam in the sky suggests that the ATMS will be in the near background (as seen from the ground station) of one of the communication satellites and will be illuminated by either the main lobe or the near sidelobes of the uplink antenna more often than it is desirable. For our analysis we first calculate the position of the ATMS with respect to each of the ground stations at a resolution of 0.2 sec, then calculate the position of each of the comsats which are at least 25 deg above the station’s local horizon; finally we calculate the angle between the line-of-sight of the ATMS and the line-of-sight of the comsat. We assume that the gain pattern of the ground station is circularly symmetric; the angle-off-station-boresight then provides an attenuation of the uplink power which we use to assess the likely effect upon the passive instrument’s operations. We assume that each ground station can communicate with all the comsats in its field of view; this implies that, on average, a ground station can engage with 46 comsats simultaneously. The analysis is repeated for the case when the uplink broadcast within the ATMS passive channels (in-band scenario) and for the case when the uplink is limited to frequencies adjacent to the ATMS channels (out-of-band scenario). The antenna of the ground station is modelled as having a HPBW (Half-Power Beam Width) of 0.16 deg and EIRP (Equivalent Isotropic Radiated Power) of 70 dBW. We account for the geometric dissipation of the signal caused by the satellite orbital altitude, the attenuation induced by atmospheric gasses at 51 GHz and the mismatch between the circular polarization of the ground-based transmitting antenna and the linear polarization of the satellite-borne receiving antenna. The damages on ATMS are estimated from bench-level measurement conducted at the ATMS’ manufacturer facilities [unpublished].

passive microwave↗
DOE OSTI2026

Collision Tracking in OpenMC: Methods and Applications in Neutron Noise, Neutron Imaging, Time-of-Flight, and Multiplicity Counting

We present the development and application of a collision tracking feature within the OpenMC Monte Carlo particle transport code, designed for diverse applications such as neutron spectroscopy, scatter camera system, neutron noise, and multiplicity counting simulations. This feature enables the tracking of individual particle collisions, with potential applications in nuclear nonproliferation, reactor physics, and nuclear security. Additionally, the feature holds potential for the calibration of neutron detectors, specifically in converting light output into energy deposited within the detectors. The implementation consists of a set of filters—such as reaction type, energy, cell, and material—that constrain the set of collisions that are tracked, extensions to the Python API to enable simple input specification, and support for writing either OpenMC’s native HDF5-based format or the Monte Carlo particle list format. This feature was added to the official OpenMC release in version 0.15.3. In this work, the feature will be applied to showcase scenarios such as time-of-flight simulations, scatter-camera imaging for neutron source localization, neutron-noise analysis to extract integral kinetic parameters such as the prompt decay constant α, and multiplicity counting to estimate the mass of special nuclear materials. Ultimately, this feature aims to expand the application scope of open-source Monte Carlo particle transport codes such as OpenMC.

Monte Carlo code↗
DOE OSTI2024

Chiral kinematic theory and converse vortical effects

Response theories in condensed matter typically describe the response of an electron fluid to external electromagnetic fields, while perturbations on neutral particles are often designed to mimic such fields. Here, we study the response of fermions to a space-time-dependent velocity field, thereby sidestepping the issue of gauge charge. First, we use a semiclassical chiral kinematic theory to obtain the local density of current and extract the orbital magnetization. The theory immediately predicts a "converse vortical effect," defined as an orbital magnetization driven by linear velocity. It receives contributions from magnetic moments on the Fermi surface and the Berry curvature of the occupied bands. Then, transcending semiclassics via a complementary Kubo formalism reveals that the uniform limit of a clean system receives only the Berry curvature contribution while other limits sense the Fermi surface magnetic moments too. We propose CoSi as a candidate material and suggest magnetometry of a sample under a thermal gradient to detect the effect. Overall, our study sheds light on the effects of a space-time-dependent velocity field on electron fluids and paves the way for exploring quantum materials using new probes and perturbations.

Chen, Kai [University of Houston]↗
DOE OSTI2023

AutoSourceID-Classifier: Star-galaxy classification using a convolutional neural network with spatial information

Aims.Traditional star-galaxy classification techniques often rely on feature estimation from catalogs, a process susceptible to introducing inaccuracies, thereby potentially jeopardizing the classification’s reliability. Certain galaxies, especially those not manifesting as extended sources, can be misclassified when their shape parameters and flux solely drive the inference. We aim to create a robust and accurate classification network for identifying stars and galaxies directly from astronomical images. Methods.The AutoSourceID-Classifier (ASID-C) algorithm developed for this work uses 32x32 pixel single filter band source cutouts generated by the previously developed AutoSourceID-Light (ASID-L) code. By leveraging convolutional neural networks (CNN) and additional information about the source position within the full-field image, ASID-C aims to accurately classify all stars and galaxies within a survey. Subsequently, we employed a modified Platt scaling calibration for the output of the CNN, ensuring that the derived probabilities were effectively calibrated, delivering precise and reliable results. Results.We show that ASID-C, trained on MeerLICHT telescope images and using the Dark Energy Camera Legacy Survey (DECaLS) morphological classification, is a robust classifier and outperforms similar codes such as SourceExtractor. To facilitate a rigorous comparison, we also trained an eXtreme Gradient Boosting (XGBoost) model on tabular features extracted by SourceExtractor. While this XGBoost model approaches ASID-C in performance metrics, it does not offer the computational efficiency and reduced error propagation inherent in ASID-C’s direct image-based classification approach. ASID-C excels in low signal-to-noise ratio and crowded scenarios, potentially aiding in transient host identification and advancing deep-sky astronomy.

Astronomy & Astrophysics↗
NASA NTRS2025

On the Stability Analysis of Astrophysical Cooling Functions

To model the temperature evolution of optically thin astrophysical environments at MHD scales, radiative and collisional cooling rates are typically either pretabulated or fit into a functional form and then input into MHD codes as a radiative loss function. Thermal balance requires estimates of the analogous heating rates, which are harder to calculate, and due to uncertainties in the underlying dissipative heating processes these rates are often simply parameterized. The resulting net cooling function defines an equilibrium curve that varies with density and temperature. Such cooling functions can make the gas prone to thermal instability (TI), which will cause departures from equilibrium. There has been no systematic study of thermally unstable parameter space for nonequilibrium states. Motivated by our recent finding that there is a related linear instability, catastrophic cooling instability, that can dominate over TI, here we carry out such a study. We show that Balbus instability criteria for TI can be used to define a critical cooling rate, Λc, that permits a nonequilibrium analysis of cooling functions through the mapping of TI zones. We furthermore illustrate how thermal conduction modifies the shape of TI zones. Upon applying a Λc-based stability analysis to coronal loop simulations, we find that loops undergoing periodic episodes of coronal rain formation are linearly unstable to catastrophic cooling instability, while TI is stabilized by thermal conduction.

Amanda Stricklan↗
NASA NTRSDate not supplied

The Aeolus-Earth DWTS Heliophysics Mission

The Aeolus-Earth Mission (DWTS-Helio) is a recently selected NASA Heliophysics investigation that will apply the innovative Doppler Wind Temperature Sounder (DWTS) technique to the Space Atmosphere Interaction Region (SAIR) at ~105-180 km. The technique, invented by Gordley & Marshall (2011), applies the Doppler Scanning with Gas Filter (DSGF) technique using a cryo-cooled MWIR radiometer to view the limb of the atmosphere from LEO (orbital velocity creates the Doppler shifts required by the technique). For this implementation of DWTS, a Nitric Oxide (NO) low pressure gas cell is used. The limb of the atmosphere is observed from a spacecraft at orbital velocity, and as an observed air volume passes through the field of view, this velocity creates a set of Doppler-shifted images. These images are used to create Doppler Integrated Pass (DIP) signals, from which atmospheric temperature and wind profiles are extracted (it may be considered a converse of the successful in-flight calibration of gas cells used for the HALOE instrument). The technique is relative, removing the need for absolute radiances, and the signal-to-noise ratio for this technique is high, resulting in good measurement precision. Different gas cells can be used to create different altitude coverage and may be combined to measure a more continuous altitude range. For example, using a NO cell captures the lower atmosphere from 30-50km as well as the SAIR, while using the 13 isotopologue of CO2 would capture 50-120km. Due to the desired data volume as well as power and thermal considerations, this initial Aeolus-Earth flight will be more suited for a ‘hosted’ payload on an Orbital Maneuvering Vehicle platform instead of the ‘free-flyer’ option originally investigated. Lastly, the technique is sufficiently powerful that a two-channel system could produce full Mars atmosphere data not currently possible.

DWTS↗
DOE OSTI2026

A hybrid Monte Carlo-deterministic second moment method with efficient variance reduction

In this work, we present a hybrid method that combines Monte Carlo with deterministic finite element methods to solve a linear Boltzmann transport equation. Our hybrid method runs orders of magnitude faster than Monte Carlo, without sacrificing accuracy, for a proxy problem from radiative transfer that contains both optically-thick and optically-thin material. We believe that this is the first demonstration of a hybrid Second Moment Method in more than one spatial dimension, the first to consider more than one material, and the first to use variance reduction. Our variance reduction approach arises from an asymptotic analysis in which we show that the magnitude of the scattering source grows without bound. We transform the problem to compute the deviation of the radiation intensity from isotropy. The magnitude of the source in the transformed problem is bounded, and the quality of the hybrid method solution is dramatically improved by a substantial reduction in the variance.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2026

Evolution of the Antarctic Ice Sheet from 2000–2300 and beyond: model sensitivity and uncertainty analysis using MPAS-Albany Land Ice

We present a description of the Antarctic Ice Sheet model configuration submitted to the ISMIP6-Antarctica-2300 experiment using the MPAS-Albany Land Ice model, along with three new sets of simulations: (1) a set of extended simulations to 2500 for three forced experiments and to 2775 for the control experiment; (2) a sensitivity analysis of our model configuration to parameters controlling basal sliding and sub-shelf melt, and to model structural choices including the choice of the energy and stress balances; and (3) a 72-member ensemble run on graphics processing units (GPUs) and analysis of variance to determine the primary sources of uncertainty in our ice-sheet model projections. Our extended simulations predict rapid retreat beginning after 2300 for SSP1-2.6 forcing and after 2500 for present-day (control) forcing, primarily in the Amundsen Sea Embayment. We find that varying the sub-shelf melt parameter between the 5th to 95th percentile values for a mean-Antarctic calibration target results in an up to ∼ ± 40 % change in sea-level contribution relative to our baseline simulations that used the median value. Using a linear basal sliding law reduces sea-level contribution by 51 %–73 % relative to our baseline nonlinear sliding law with an exponent of 1/5. When using basal sliding law exponents of 1/3 and 1/10, the overall difference from our baseline simulations at 2300 is on the order of 10 %. The Amundsen Sea Embayment region displays a strongly non-linear dependence of mass loss on the sliding law exponent, with no discernible relationship between the sliding law exponent and the mass loss by 2300, while the sectors feeding the Ross and Filchner-Ronne ice shelves exhibit more mass loss with a more-plastic sliding law. Our model fidelity sensitivity experiments reveal a 9 %–31 % increase in sea-level contribution when using a depth-integrated stress balance approximation relative to our three-dimensional solver, while using a fixed-in-time temperature field increases sea-level contribution by 14 %–88 % relative to two thermomechanically coupled configurations. Our 72-member ensemble and analysis of variance show that the uncertainty in long-term projections is dominated by the choice of Earth system model forcing and the presence or absence of hydrofracture forcing, rather than uncertainty in sliding and sub-shelf melt parameters.

58 GEOSCIENCES↗
DOE OSTI2026

A Critical Review of Heliostat Design for Concentrating Solar Thermal Technologies

Heliostat designs have undergone a widespread and eclectic development process, with many unique designs demonstrated. However, recent developments do not show the overall cohort converging toward a globally accepted universal design. Here, this study characterizes heliostats by breaking down and evaluating design traits based on emergent patterns from a comprehensive compilation of known heliostat designs spanning several decades. Four main categories for evaluation emerged: heliostat base, heliostat primary axis, heliostat drive, and facet support. Each of these four categories is further defined by four subtypes so that all heliostats fall into a single subtype within each category. The classified heliostats are ranked, yielding several view slices into the heliostat compilation or a breakdown of heliostats by type. An analysis of the breakdown shows several trends: a scale-up of established designs, new approaches at small to medium scales, and a movement toward greater adoption of linear drives. These trends reflect the most meaningful contributing factors to a proposed trajectory for a new era of heliostat designs striving to meet widely considered cost targets of $\$$50/m 2 or $\$$75/m 2 .

14 SOLAR ENERGY↗