Search NASA⌕ Search

SEARCH · Search NASA

Results for “voltage balancing”

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 361 records · Page 20

Robust Disturbance Rejection Rotor Current Control of Doubly-Fed Induction Generators

This article proposes a robust transient disturbance rejection current controller (TDRCC) to improve the transient performance of doubly-fed induction generators (DFIGs) in wind turbines during grid disturbances, such as short circuits. Here, the proposed robust TDRCC replaces the proportional-integral-derivative (PID) current controllers conventionally used in the vector control scheme of the DFIG rotor side converter (RSC). The TDRCC estimates the transient disturbances caused by grid faults or other external disturbances and provides compensation to the DFIG rotor voltage in the current control loops of the RSC to improve the robustness of the controller to disturbance. A sliding-mode current controller (SMCC) is also designed to highlight the transient performance improvement of the DFIG using the TDRCC over the state-of-the-art control methods during grid short circuit faults. Simulation studies are conducted in PSCAD/EMTDC for a 3.6-MW DFIG wind turbine with the proposed TDRCC, the conventional PI controller, and the SMCC, respectively during the most severe balanced three-phase grid short circuit fault specified by the U.S. grid code as well as an unbalance single-phase grid short circuit fault. Hardware experiments are conducted on a 200-W DFIG wind turbine emulator with the three different controllers for the same three-phase short circuit fault. Simulation and hardware experiment results show that the TDRCC reduces the peak rotor current significantly when compared with the PI and SMCC controlled DFIG system and, therefore, would help prolong the lifespan of the DFIG's power electronics and reduce maintenance costs.

42 ENGINEERING↗

Control of multiple resonant power processors in a multi-source system

Analysis and test results show that phasor-regulated, Mapham-derived resonant inverters can be paralleled to provide standardizing interfaces for multiple sources on a utility-type, aerospace power distribution bus. The basic sources do not require matching in any way, and may have grossly different characteristics. Fully stable system architectures with multiple sources, parallel/redundant distribution buses, and a wide variety of loads can be easily constructed and controlled. The commands and parameters available for system control allow for tight tolerance bus voltage control, and absolute power-sharing control from the various sources over the full range of possible source and load variations. That level of control enables simplified load power processing hardware and the distribution of losses to optimally load the source thermal control system. Positive control of all system performance and allocation of losses are not required by all missions or vehicles, and overall vehicle considerations do not always require the loads on vehicle energy sources and thermal control systems to be balanced. In those cases, power system control can be simplified, and a hierarchical set of defaults can be substituted for computer-generated or supervisory input commands to allow for stable, fully autonomous system operation.

Mildice, James↗

As‐Cast Ternary Organic Solar Cells Based on an Asymmetric Side‐Chains Featured Acceptor with Reduced Voltage Loss and 14.0% Efficiency

Abstract A new small‐molecule nonfullerene acceptor based on the benzo[1,2‐ b :4,5‐ b ′]dithiophene (BDT) fused central core with asymmetrical alkoxy and thienyl side chains, namely TOBDT , is designed and synthesized. The alkoxy unit helps narrow the bandgap, and thienyl side chain helps enhance the intermolecular interaction. As a result, TOBDT is suitable to match the deep‐lying highest occupied molecular orbital (HOMO) of polymer donor PM6 . Then, a strong crystalline acceptor IDIC is introduced as the third component to fabricate as‐cast nonfullerene ternary devices to achieve absorption and morphology control. Addition of IDIC not only mixes well with TOBDT but modulates the morphology of the blend film, which helps to balance the charge transport properties and reduce the photovoltage loss of ternary devices. All these contribute to synergetic improvement of J sc , V oc , and fill factor parameters, leading to a power conversion efficiency of 14.0% for the as‐cast fullerene‐free ternary device.

Chen, Xuebin↗

Thin-Film Air-Mass-Flow Sensor of Improved Design Developed

Researchers at the NASA Glenn Research Center have developed a new air-mass-flow sensor to solve the problems of existing mass flow sensor designs. NASA's design consists of thin-film resistors in a Wheatstone bridge arrangement. The resistors are fabricated on a thin, constant-thickness airfoil to minimize disturbance to the airflow being measured. The following photograph shows one of NASA s prototype sensors. In comparison to other air-mass-flow sensor designs, NASA s thin-film sensor is much more robust than hot wires, causes less airflow disturbance than pitot tubes, is more accurate than vane anemometers, and is much simpler to operate than thermocouple rakes. NASA s thin-film air-mass-flow sensor works by converting the temperature difference seen at each leg of the thin-film Wheatstone bridge into a mass-flow rate. The following figure shows a schematic of this sensor with air flowing around it. The sensor operates as follows: current is applied to the bridge, which increases its temperature. If there is no flow, all the arms are heated equally, the bridge remains in balance, and there is no signal. If there is flow, the air passing over the upstream legs of the bridge reduces the temperature of the upstream legs and that leads to reduced electrical resistance for those legs. After the air has picked up heat from the upstream legs, it continues and passes over the downstream legs of the bridge. The heated air raises the temperature of these legs, increasing their electrical resistance. The resistance difference between the upstream and downstream legs unbalances the bridge, causing a voltage difference that can be amplified and calibrated to the airflow rate. Separate sensors mounted on the airfoil measure the temperature of the airflow, which is used to complete the calculation for the mass of air passing by the sensor. A current application for air-mass-flow sensors is as part of the intake system for an internal combustion engine. A mass-flow sensor is used to provide accurate information about the amount of air entering the engine so that the amount of fuel can be adjusted to give the most efficient combustion. The ideal mass-flow sensor would be a rugged design that minimizes the disturbance to the flow stream and provides an accurate reading of both smooth and turbulent flows; NASA's design satisfies these requirements better than any existing design. Most of the mass-flow sensors used today are the hot wire variety. Hot wires can be fragile and cannot accurately measure a turbulent or reversing flow, which is often encountered in an intake manifold. Other types of mass-flow sensors include pitot tubes, vane anemometers, and thermocouple rakes-all of which suffer from some type of performance problem. Because it solves these performance problems while maintaining a simple design that lends itself to low-cost manufacturing techniques, NASA s thin-film resistance temperature detector air-mass-flow sensor should lead to more widespread use of mass-flow sensors.

Fralick, Gustave C.↗

Development of Fixtures and Methods to Assess the Durability of Balance of Systems Components

The degradation of photovoltaic (PV) balance of systems (BoS) components is not well studied, but the consequences include offline modules, strings, and inverters; system shutdown; arc faults; and fires. A utility provider experienced a ~30% failure rate in their power transfer chain, originally attributed to branch connectors. Field-failed specimen assemblies were, therefore, examined, consisting of cable connector, branch connector, and discrete fuse components. In this study, unused field-vintage specimens are examined using a benchtop prototype fixture to identify the most influential environmental stressors on BoS components as well as the effect of external mechanical perturbation. The prototype fixture was used to develop a perturbation capability for future use in the combined-accelerated stress testing chamber. The benchtop experiments were also used to develop the in-situ data acquisition of specimen current, voltage, and temperature. A significant increase in operating temperature (~100 °C from ~40 °C) and a different failure mode (arcing at the metal pins rather than overheating of the fuse filament) were observed promptly once periodic mechanical perturbation was applied. The current at failure was decreased from 35 A (measured for static specimens, with failure occurring in the fuses) to 15 A (for tests with mechanical perturbation, with failure at the male/female metal pin connection). After initial examination using X-ray computed tomography, the external plastic was machined away from failed specimens to allow for failure analysis, including the extraction of the internal convolute springs for morphological examination (optical and electron microscopy). Chemical composition analysis included energy-dispersive X-ray spectroscopy, differential scanning calorimetry, and Fourier transform infrared spectroscopy.

14 SOLAR ENERGY↗

A Difluoro‐Methoxylated Ending‐Group Asymmetric Small Molecule Acceptor Lead Efficient Binary Organic Photovoltaic Blend

Abstract Developing a new end group for synthesizing asymmetric small molecule acceptors (SMAs) is crucial for achieving high‐performance organic photovoltaics (OPVs). Herein, an asymmetric small molecule acceptor, BTP‐BO‐4FO, featuring a new difluoro‐methoxylated end‐group is reported. Compared to its symmetric counterpart L8‐BO, BTP‐BO‐4FO exhibits an upshifted energy level, larger dipole moment, and more sequential crystallinity. By adopting two representative and widely available solvent additives (1‐chloronaphthalene (CN) and 1,8‐diiodooctane (DIO)), the device based on PM6:BTP‐BO‐4FO (CN) photovoltaic blend demonstrates a power conversion efficiency (PCE) of 18.62% with an excellent open‐circuit voltage (V OC ) of 0.933 V, which surpasses the optimal result of L8‐BO. The PCE of 18.62% realizes the best efficiencies for binary OPVs based on SMAs with asymmetric end groups. A series of investigations reveal that optimized PM6:BTP‐BO‐4FO film demonstrates similar molecular packing motif and fibrillar phase distribution as PM6:L8‐BO (DIO) does, resulting in comparable recombination dynamics, thus, similar fill factor. Besides, it is found PM6:BTP‐BO‐4FO possesses more efficient charge generation, which yields betterV OC –J SC balance. This study provides a new ending group that enables a cutting‐edge efficiency in asymmetric SMA‐based OPVs, enriching the material library and shed light on further design ideas.

Chemistry↗

Scalability and Effectiveness of Smart Charge Management

The rise in electric vehicle (EV) adoption presents growing challenges for power grids, particularly from simultaneous residential charging, which can cause voltage fluctuations and increase feeder peak loads. Baltimore Gas and Electric (BGE), with support from the U.S. Department of Energy, initiated a pilot program to evaluate managed residential EV charging through Smart Charge Management (SCM). This study analyzes real-world charging behavior data from the pilot and feeder-level base loads from BGE to simulate residential charging scenarios through 2035 across the Washington, DC–Baltimore region. Grid impacts under unmanaged charging are compared to three SCM strategies: TOU-immediate, TOU-distributed, and Load Balancing. Results show that the magnitude of peak reduction is highly feeder-dependent. Some feeders achieve reductions of more than 40% at high enrollment levels, while others show improvements closer to 10–15%. This heterogeneity reflects differences in baseline feeder load shapes, EV penetration, and plug-in behavior across customers. Results also highlight trade-offs between shifting load away from peak periods and minimizing secondary demand peaks, offering practical insights for future utility program design.

Electric vehicle↗

Reactive Oxygen Ion Beam-Induced Deposition for Concurrent Purification of Platinum Nanostructures

Oxygen-focused ion beam induced deposition (O-FIBID) enables the direct-write fabrication of Pt nanostructures while simultaneously enhancing purity concurrently through reactive oxygen–deposit interactions. By systematically varying the dwell time, accelerating voltage, and precursor pressure, the Pt content and conductivity can be controlled. Under optimum conditions, the Pt content reached 63 at.%. Across the dwell-time range used for resistivity measurements, the Pt content increased from 20 to 33 at.%, while the resistivity decreased from 2.9 × 10 4 μΩ·cm to 1.2 × 10 3 μΩ·cm, which is consistent with enhanced percolation through Pt grains and the lower intrinsic resistivity of the purer Pt deposit. The simulation results support a purification mechanism driven by the beam-induced activation of implanted oxygen balanced against the preferential sputtering of Pt. These results demonstrate O-FIBID as a viable method for the nanoscale direct write of conductive Pt without post-processing, and some deviations from conventional FIBID wisdom are observed. These results serve as a foundation for exploring nascent, reactive focused ion beam-induced deposition processes.

deposit purification↗

Definition study for photovoltaic residential prototype system

A parametric sensitivity study and definition of the conceptual design is presented. A computer program containing the solar irradiance, solar array, and energy balance models was developed to determine the sensitivities of solar insolation and the corresponding solar array output at five sites selected for this study as well as the performance of several solar array/battery systems. A baseline electrical configuration was chosen, and three design options were recommended. The study indicates that the most sensitive parameters are the solar insolation and the inverter efficiency. The baseline PST selected is comprised of a 133 sg m solar array, 250 ampere hour battery, one to three inverters, and a full shunt regulator to limit the upper solar array voltage. A minicomputer controlled system is recommended to provide the overall control, display, and data acquisition requirements. Architectural renderings of two photovoltaic residential concepts, one above ground and the other underground, are presented. The institutional problems were defined in the areas of legal liabilities during and after installation of the PST, labor practices, building restrictions and architectural guides, and land use.

Imamura, M. S.↗

A Transformer Flux Balancing Scheme Based on Magnetizing Current Harmonic in Dual-Active-Bridge Converters

The Dual-Active-Bridge (DAB) converter transformer flux unbalance, or dc bias, could be a critical issue due to the mismatches of converter circuit parameters and control transients. In this paper, a steady-state flux balancing scheme, based on the second-order harmonic of magnetizing current has been proposed. The relationship between steady-state flux/magnetizing current unbalance and its harmonic contents has been analyzed, with which a new current sensing method and control strategy are introduced. The proposed magnetizing current sensing method can be implemented for both unity and complicated fractional turns ratios, without high sensitivity to dc offset from the analog circuitry. The proposed scheme has been validated with experimental results on 850-V/850-V DAB converters within a system, as well as in an 850-V/6700-V medium voltage DAB unit.

Gao, Zihan↗

Transient Studies of Nonequilibrium Transport in Two-Dimensional Semiconductors

The overarching objective of research under DOE support has been to investigate emergent carrier phenomena in different two-dimensional (2D) systems, under conditions where they are driven far from equilibrium. The description of nonequilibrium transport has long represented one of the most enduring challenges in condensed-matter physics, in large part due to the need to invoke the role of complex many-body phenomena (electron-phonon and electron-electron coupling) in this regime. In this project, we have explored aspects of this problem in a variety of different two-dimensional systems. These include the modulation-doped two-dimensional electrons gas, or 2DEG, in GaAs/AlGaAs heterostructures; the electron and hole systems in graphene and its related hetero-structures, and; the strongly-confined electron layers in transition-metal dichalcogenide (MoS 2 and WS 2 ) and trichalcogenide (TiS 3 ) materials. Significant understanding has been achieved in a number of different areas, as we elaborate upon in subsequent sections of this proposal: Nanoscale spectroscopy of defect-related barriers in MoS2 Negative differential conductance due to inter-valley scattering in WS 2 under high fields Incipient localization due to scattering from frozen-phonon disorder in quantum point contacts High-field transient transport in graphene/hexagonal BN heterostructures Universal voltage scaling of quantum corrections in graphene under nonequilibrium A defining feature of much of this research has been an approach to problem solving in which novel experimental techniques (such as transient electrical pulsing and scanning-probe microscopy) are combined with the development (in collaboration with colleagues from condensed-matter theory) of sophisticated models, a strategy that we propose to continue in our future research. In recent years, we have developed an increasing capacity for handling atomically-thin materials, and their van der Waals systems, a capacity that is reflected in the balance of topics listed above.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Generation of excited species in a streamer discharge

At or near atmospheric pressure, most transient discharges, particularly in molecular gases or gas mixture containing molecular gases result in a space-charge dominated transport called a streamer discharge. The excited species generation in such discharges forms the basis for plasma chemistry in most technological applications. In this paper we simulate the propagation of streamers in atmospheric pressure N 2 to understand the energy partitioning in the formation of various excited species and compare the results to a uniform Townsend discharge. The model is fully two-dimensional with azimuthal symmetry. Results show a significantly larger fraction of the energy goes into vibrational excitation of the N 2 ground state in a streamer-type discharge in comparison to a Townsend discharge. For lower applied voltages, the anode-directed (negative) steamer is slightly more efficient in channeling energy for excited species production in comparison to a cathode-directed (positive) streamer. Near 70% overvoltage both types of streamers show very similar energy partitioning but quite different from a Townsend discharge.

42 ENGINEERING↗

Ensemble voting-based fault classification and location identification for a distribution system with microgrids using smart meter measurements

This study presents an ensemble learning approach for fault classification and location identification in a smart distribution network containing photovoltaics (PV)-based microgrid. Lack of available data points and the unbalanced nature of the distribution system make fault handling a challenging task for utilities. The proposed method uses event-driven voltage data from smart meters to classify and locate faults. The ensemble voting classifier is composed of three base learners; random forest, k-nearest neighbours, and artificial neural network. The fault location (FL) task has been formulated as a classification problem where the fault type is classified in the first step and based on the fault type, the faulty bus is identified. The method is tested on IEEE-123 bus system modified with added PV-based microgrid along with dynamic loading conditions and varying fault resistances from 0 to 20 Ω for both unbalanced and balanced fault types. A further sensitivity analysis has been done to test the robustness of the proposed method under various noise levels and data loss errors in the smart meter measurements. The ensemble method shows improved performance and robustness compared to some previously proposed FL methods. Finally, the proposed method has been experimentally validated on a real-time simulation-based testbed using a state-of-the-art digital real-time simulator, industry standard DNP3 communication protocol and a cpu-based control centre running the FL algorithm.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Thermal Control using Electrochromism

The applicability of a charge balanced electrochromic device to modulate the frequencies in the thermal infrared region is examined in this study. The device consisted of a transparent conductor, WO3 anode, PMMA/LiClO4 electrolyte, V2O5 cathode and transparent conductor. The supporting structure in the device is SnO2 coated glass and the edges are sealed with epoxy to reduce moisture absorption. The performance evaluation comprised of cyclic voltammetric measurements and determination of transmittance at various wavelengths. The device was subjected to anodic and cathodic polarization by sweeping the potential at a rate of 10 mV/sec from -0.8V to 1.8V. The current versus voltage profile indicated no reaction between -0.5 and +0.5 V. The device is colored green at 1.8V with a transmittance of 5% at a wavelength, lambda=900 nm and colorless at -0.8V with a transmittance of 74% at lambda=500 nm. The optical modulation is limited to 400-1500 nm and there is no activity in the thermal infrared. The switching time is 75 seconds for transmittance to decrease from 74% to 50%. The device yielded reproducible values for transmittance when cycled between colored and bleached states by application of 1.8V and -0.8V, respectively.

Vaidyanathan, Hari↗

Thermal Control Using Electrochromism

The applicability of a charge balanced electrochromic device to modulate the frequencies in the thermal infrared region is examined in this study. The device consisted of a transparent conductor, WO3, anode, PMMA/LiClO4 electrolyte, V2O5 cathode and transparent conductor. The supporting structure in the device is SnO2 coated glass and the edges are sealed with epoxy to reduce moisture absorption. The performance evaluation comprised of cyclic voltammetric measurements and determination of transmittance at various wavelengths. The device was subjected to anodic and cathodic polarization by sweeping the potential at a rate of 10 mV/sec from -0.8V to 1.8V. The current versus voltage profile indicated no reaction between -0.5 and +0.5 V. The device is colored green at 1.8 V with a transmittance of 5% at a wavelength, lambda = 900 nm and colorless at -0.8 V with a transmittance of 74% at lambda = 500 nm. The optical modulation is limited to 400-1500 nm and there is no activity in the thermal infrared. The switching time is a function of temperature and time for coloring reaction was slower than the bleaching reaction. The device yielded reproducible values for transmittance when cycled between colored and bleached states by application of 1.8V and -0.8V, respectively.

Vaidyanathan, Hari↗

Robust Medium-Voltage Distribution System State Estimation using Multi-Source Data

Due to the lack of sufficient online measurements for distribution system observability, pseudo-measurements from short-term load or distributed renewable energy resources (DERs) forecasting are used. However, the accuracy of them is low and thus significantly limits the performance of distribution system state estimation (DSSE). In this paper, a robust DSSE that integrates multi-source measurement data is proposed. Specifically, the historical low-voltage (LV) side smart meters are used to forecast load and DERs injections via the support vector machine (SVM) with optimally tuned parameters. By contrast, the online smart meters at LV side are utilized to derive equivalent power injections at the MV/LV transformers, yielding more accurate pseudo-measurements compared to the forecasted injections. Furthermore, to deal with bad data caused by communication loss, instrumental errors and cyber attacks, robust DSSE that relies on generalized maximum-likelihood (GM)-estimation criterion is developed. The projection statistics are developed to adjust the weights of each measurement, leading to better balance between pseudo- and real-time measurements. Numerical results conducted on modified IEEE 33-bus system with DG integration demonstrate the effectiveness and robustness of the proposed method.

distribution system state estimation↗

Generalized Power Flow Model of an Extra High-Power Multi-Terminal HVdc Transmission Grid with Parallel-Connected and Voltage-Stacked Converters

To increase power transfer capacity of high-voltage direct current (HVdc) transmission, a new extra high-power HVdc architecture with multiple standard modular multilevel converters (MMCs) per substation has recently been introduced. This paper proposes a power flow model for a multi-terminal HVdc (MTdc) grid with this innovative substation architecture. The proposed MTdc model can be integrated seamlessly with existing ac-dc power flow algorithms with minimal modifications. The model is applicable to various MTdc grid types and topologies, different numbers of dc buses, dc lines, and MMCs per substation, along with diverse control parameters. In addition, it accurately captures both balanced and unbalanced operations of the MTdc grid. The proposed model is verified using a 5-terminal bipole MTdc grid that spans 4 areas in the Eastern Interconnection system of the USA. The numerical solutions obtained from unified and sequential ac-dc power flow algorithms under different operating conditions closely match the time-domain simulation results in PSCAD, validating the accuracy and versatility of the proposed MTdc power flow model.

Nguyen, Quan H.↗

Direct Air Capture Using Trapped Small Amines in Hierarchical Nanoporous Capsules on Porous Electrospun Fibers (Final Technical Report)

This report summarizes the carbon capture research and development conducted by The State University of New York at Buffalo (UB) and GTI Energy (GTI) for award “DE-FE0031969: Direct Air Capture Using Trapped Small Amines in Hierarchical Nanoporous Capsules on Porous Electrospun Fibers” sponsored by the U.S. Department of Energy (DOE). The objective of this project is to develop an innovative sorbent structure of trapped small amines in HNC embedded in PEF for DAC. This involves tailoring both sorbent and PEF materials to achieve a compact system for DAC with high capacity for CO 2 at concentrations typically available in air and at near ambient conditions. An innovative sorbent structure of trapped small amines in hierarchical nanoporous capsules (HNC) embedded in porous electrospun fibers (PEF) was developed for direct air capture (DAC). This involves tailoring both sorbent and PEF materials to achieve a compact system for DAC with high capacity for CO 2 at concentrations typically available in air and at near ambient conditions. An interfacial polymerization process was developed, which utilized loaded amines inside mesoporous silica and trimesoyl chloride (TMC) dissolved in organic solvents as the precursors, to generate a polyamide (PA) coating layer on mesoporous silica and thus trap amines. Reaction conditions, including TMC concentration, organic solvents, reaction time, etc., for interfacial polymerization were optimized to effectively trap loaded amines, and cyclic heating-cooling operation was conducted to evaluate the coating quality. Larger pore volume mesoporous silica was also synthesized to increase amine loading and thus increase CO 2 capacity. The optimized sorbent material exhibited CO 2 capacity as high as 4.88 mmol/g under humid DAC conditions and negligible loss (<1%) during 10 cyclic heating-cooling operations. The optimized PA-coated sorbent also showed fast adsorption and desorption kinetics, with <20% t1/2 increase compared to uncoated sorbent. PEF fabrication conditions, including organic solvents for dissolving core and shell polymers, voltage, distance from the nozzle to the collection panel, etc. were adjusted to better incorporate HNC. After incorporating the optimized sorbent material into PEF, the structured sorbent had a CO 2 capacity of approximately 4.0 mmol/g under humid DAC conditions, with capacity loss of 0.17% per cycle and t1/2 increase less than 10%. A techno-economic analysis (TEA) for the process design for a DAC system based on our developed sorbent structure of trapped small amines in HNC embedded in PEF was conducted. The process design included process description and major equipment sizing and energy and mass balances in addition to scale-up research results and estimated capture cost. Aspen Adsorption Simulator was used to fit the experimentally measured breakthrough curves and extract equilibrium and kinetic data of the optimized sorbent. Our results indicated that for a DAC plant with CO 2 productivity of 3,000 tonne/year, the levelized cost of CO 2 capture was $\$$612/tonne, with the largest contribution of 44.33% from the fixed operation cost. Increasing CO 2 productivity, while maintaining similar fixed operation cost, is expected to significantly reduce the CO 2 capture cost. A sensitivity study was also conducted to understand the influence of total plant cost, sorbent cost, CO 2 concentration in the feed, sorbent mat lifetime, sorbent regeneration electricity, and adsorption blower pressure drop on the levelized cost of CO 2 capture, revealing a capture cost range of $\$$520-870/tonne.

36 MATERIALS SCIENCE↗