Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Operator approximation”

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 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Operation of the H- Linac at FNAL

The Fermi National Accelerator Laboratory (FNAL) Linac has been in operation for 52 years. In approximately four years, it will be replaced by a new 800 MeV superconducting machine, the PIP-II SRF Linac. In the current configuration, the Linac delivers H- ions at 400 MeV and injects protons by charge exchange into the Booster synchrotron. Despite its age, the Linac is the most stable accelerator in the FNAL complex, reliably sending 22 mA in daily operations. We will discuss the status of the operation, beam studies, and plans.

43 PARTICLE ACCELERATORS↗
DOE OSTI2021

Quantum Alternating Operator Ansatz (QAOA) Phase Diagrams and Applications for Quantum Chemistry

Determining Hamiltonian ground states and energies is a challenging task with many possible approaches on quantum computers. While variational quantum eigensolvers are popular approaches for near term hardware, adiabatic state preparation is an alternative that does not require noisy optimization of parameters. Beyond adiabatic schedules, QAOA is an important method for optimization problems. In this work we modify QAOA to apply to finding ground states of molecules and empirically evaluate the modified algorithm on several molecules. This modification applies physical insights used in classical approximations to construct suitable QAOA operators and initial state. We find robust qualitative behavior for QAOA as a function of the number of steps and size of the parameters, and demonstrate this behavior also occurs in standard QAOA applied to combinatorial search. To this end we introduce QAOA phase diagrams that capture its performance and properties in various limits. In particular we show a region in which non-adiabatic schedules perform better than the adiabatic limit while employing lower quantum circuit depth. We further provide evidence our results and insights also apply to QAOA applications beyond chemistry.

Kremenetski, Vladimir↗
DOE OSTI2023

Residual-based error correction for neural operator accelerated infinite-dimensional Bayesian inverse problems

We explore using neural operators, or neural network representations of nonlinear maps between function spaces, to accelerate infinite-dimensional Bayesian inverse problems (BIPs) with models governed by nonlinear parametric partial differential equations (PDEs). Neural operators have gained significant attention in recent years for their ability to approximate the parameter-to-solution maps defined by PDEs using as training data solutions of PDEs at a limited number of parameter samples. The computational cost of BIPs can be drastically reduced if the large number of PDE solves required for posterior characterization are replaced with evaluations of trained neural operators. However, reducing error in the resulting BIP solutions via reducing the approximation error of the neural operators in training can be challenging and unreliable. We provide an a priori error bound result that implies certain BIPs can be ill-conditioned to the approximation error of neural operators, thus leading to inaccessible accuracy requirements in training. To reliably deploy neural operators in BIPs, we consider a strategy for enhancing the performance of neural operators: correcting the prediction of a trained neural operator by solving a linear variational problem based on the PDE residual. We show that a trained neural operator with error correction can achieve a quadratic reduction of its approximation error, all while retaining substantial computational speedups of posterior sampling when models are governed by highly nonlinear PDEs. The strategy is applied to two numerical examples of BIPs based on a nonlinear reaction–diffusion problem and deformation of hyperelastic materials. We demonstrate that posterior representations of the two BIPs produced using trained neural operators are greatly and consistently enhanced by error correction.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Pellet cladding mechanical interaction as a potential failure mechanism during a control rod drop accident in a boiling water reactor

Boiling water reactors (BWRs) represent approximately one-third of the operating fleet in the United States, contributing significantly towards the global effort in reducing carbon emissions. Even though most of the operating fleet has been in operation for quite some time, continued advancements in new nuclear fuel (such as accident tolerant fuel) or operating regimes (such as power up-rates and higher burnup operation) necessitates similar advancements in modeling and simulation capabilities. Bison, a component of the Virtual Environment for Reactor Applications (VERA), is a high-fidelity fuel performance code able to explore the fuel performance of a wide variety of fuel types in one-, two-, and three-dimensions. Until recently, the code had not been used for analyses of BWRs. Modeling capabilities have been added for Gd-bearing UO{sub 2} and pure zirconium liners. New models have been added based upon the U.S. Nuclear Regulatory Commission (NRC) guide-lines for hydrogen pickup in Zircaloy-2 claddings and pellet-clad mechanical interaction (PCMI) failure during a reactivity insertion accident (RIA), known as a rod drop accident (CRDA) in BWRs. Implementing and/or improving Bison modeling capabilities extended its analytical reach to areas beyond its original intended purpose. This paper demonstrates one of these capabilities as a proof of concept. Recently developed models enable Bison to provide an alternate approach for cladding integrity determination in CRDA evaluations, which currently use bounding conservative estimates. As part of this demonstration, NRC guidance on hydrogen-pickup and PCMI failure models were utilized in this research. Even though more research is needed in establishing right inputs and process in this area, this paper demonstrates Bison's ability to determine cladding integrity in a CRDA evaluation. This first of a kind demonstration is a proof of concept in this area, which could potentially be extended to a number of other areas where a more accurate cladding integrity determination would be needed. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2021

Data-Driven Learning for the Mori--Zwanzig Formalism: A Generalization of the Koopman Learning Framework

A theoretical framework which unifies the conventional Mori--Zwanzig formalism and the approximate Koopman learning of deterministic dynamical systems from noiseless observation is presented. In this framework, the Mori--Zwanzig formalism, developed in statistical mechanics to tackle the hard problem of construction of reduced-order dynamics for high-dimensional dynamical systems, can be considered as a natural generalization of the Koopman description of the dynamical system. We next show that, similar to the approximate Koopman learning methods, data-driven methods can be developed for the Mori--Zwanzig formalism with Mori's linear projection operator. We have developed two algorithms to extract the key operators, the Markov and the memory kernel, using time series of a reduced set of observables in a dynamical system. We have adopted the Lorenz `96 system as a test problem and solved for the above operators. These operators exhibit complex behaviors, which are unlikely to be captured by traditional modeling approaches in Mori--Zwanzig analysis. The nontrivial generalized fluctuation-dissipation relationship, which relates the memory kernel with the two-time correlation statistics of the orthogonal dynamics, was numerically verified as a validation of the solved operators. Here we present numerical evidence that the generalized Langevin equation, a key construct in the Mori--Zwanzig formalism, is more advantageous in predicting the evolution of the reduced set of observables than the conventional approximate Koopman operators.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Approximating Trajectory Constraints With Machine Learning – Microgrid Islanding With Frequency Constraints

Here, we introduce deep earning aided constraint encoding to tackle the frequency-constraint microgrid scheduling problem. The nonlinear function between system operating condition and frequency nadir is approximated by using a neural network, which admits an exact mixed-integer formulation (MIP). This formulation is then integrated with the scheduling problem to encode the frequency constraint. With the stronger representation power of the neural network, the resulting commands can ensure adequate frequency response in a realistic setting in addition to islanding success. The proposed method is validated on a modified 33-node system. Successful islanding with a secure response is simulated under the scheduled commands using a detailed three-phase model in Simulink. The advantages of our model are particularly remarkable when the inertia emulation functions from wind turbine generators are considered.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2025

Leveraging operator learning to accelerate convergence of the preconditioned conjugate gradient method

We propose a new deflation strategy to accelerate the convergence of the preconditioned conjugate gradient (PCG) method for solving parametric large-scale linear systems of equations. Unlike traditional deflation techniques that rely on eigenvector approximations or recycled Krylov subspaces, we generate the deflation subspaces using operator learning, specifically the Deep Operator Network (DeepONet). To this aim, we introduce two complementary approaches for assembling the deflation operators. The first approach approximates near-null space vectors of the discrete PDE operator using the basis functions learned by the DeepONet. The second approach directly leverages solutions predicted by the DeepONet. To further enhance convergence, we also propose several strategies for prescribing the sparsity pattern of the deflation operator. Here, a comprehensive set of numerical experiments encompassing steady-state, time-dependent, scalar, and vector-valued problems posed on both structured and unstructured geometries is presented and demonstrates the effectiveness of the proposed DeepONet-based deflated PCG method, as well as its generalization across a wide range of model parameters and problem resolutions.

Deflation↗
DOE OSTI2023

A structural model of the long-term degradation of the concrete biological shield

The concrete biological shield (CBS) of light water reactors is exposed to high neutron radiation dose in the long term, which may lead to the degradation of the concrete’s mechanical properties. Given the important shielding role of the CBS, it is necessary to investigate the irradiation effects at the structural scale and provide estimates of the damage extent from the wall’s inner surface to study potential license renewals. For this purpose, we developed a mechanical model accounting for radiation-induced expansion, creep, and damage in concrete using the Grizzly finite element code, informed by ex-core neutron flux calculations using the VERA tool. The model was applied to a 3D CBS structure represented by the CBS wall, a steel liner, reinforcement bars, and a concrete base mat and evaluated damage at 40, 60, and 80 years of operation. The VERA model predicted a maximum fluence of approximately 2 x 10 19 ncm -2 at 80 years of operation. The results showed that damage is highest at the inner surface of the CBS wall and gradually decreases with depth. It extends beyond the rebar after 60 years and reaches a depth of approximately 12 cm at 80 years.

42 ENGINEERING↗
DOE OSTI2025

JHTDB-wind: a web-accessible large-eddy simulation database of a wind farm with virtual sensor querying

This paper introduces JHTDB-wind (https://turbulence.idies.jhu.edu/datasets/windfarms, last access: 11 November 2025), a publicly accessible database containing large-eddy simulation (LES) data from wind farms. Building on the framework of the Johns Hopkins Turbulence Database (JHTDB), which hosts direct numerical simulation (DNS) and some LES datasets of canonical turbulent flows, JHTDB-wind stores the 4D space–time history of the flow and provides users the ability to access and query the data via a web-based virtual sensor interface. The initial dataset comprises LES results from a large wind farm with 10×6 turbines, modeled using a filtered actuator line method, under conventionally neutral atmospheric conditions. These data comprise 1 h (hour) of flow field data (velocity, pressure, potential temperature deviation, subgrid-scale (SGS) eddy viscosity, and turbine forces, approximately 15 TB (terabytes) and wind turbine data – including both turbine-level operational quantities and blade-level aerodynamic quantities (approximately 1.3 TB) – stored in Zarr and Parquet formats, respectively. Data retrieval is facilitated by the giverny Python package, allowing remote users to query the database in Python or MATLAB (C and Fortran support are available for flow field data). This paper details the simulation setup and demonstrates data access through examples that analyze wind farm flow structures and turbine performance. The framework is extensible to future datasets, including the JHTDB-wind diurnal cycle simulation analyzed in Xiao et al. (2025).

17 WIND ENERGY↗
DOE OSTI2024

Performance Validation of a Thermally Integrated 50 kW High Temperature Electrolyzer System

In the proposed project, INL and OxEon seek to improve the value proposition of hydrogen production by integrating reversible fuel cell operations at relatively small scale for distributed energy applications. This goal will be accomplished by converting a 50 kW solid oxide electrolysis cell or SOEC system into a reversible system that operates at 30 kW in electrolysis mode and approximately 10 kW in fuel cell mode. The reversible SOC system will be operated for over 3,000 hours using an improved catalyst in the fuel electrode. Steam for the electrolysis will come from an electric boiler Thermal Energy Distribution System that will be configured to mimic an industrial source of low-grade heat. Thermodynamic analysis will be performed to demonstrate the potential of the system to achieve >85% system efficiency in electrolysis mode. Finally, a technoeconomic analysis will be completed to show potential to produce hydrogen at a cost of $2/kg. The figure at the right indicates a target cost breakdown to achieve that goal.

08 HYDROGEN↗
DOE OSTI2026

Dynamics of the water-plasma interface in various discharge modes of atmospheric-pressure plasmas

This paper presents an experimental study of the dynamics of the water-plasma interface during the interaction high-voltage (HV) atmospheric pressure discharge with the water surface. Under the influence of the voltage applied to the pin type electrode, discharges are formed at the air layer between the sharp-tip HV electrode and deionized water. Three discharge regimes were identified from synchronized electrical waveforms and imaging: (i) a weak linear regime at 3–10.6 kV, (ii) a branching streamer regime at the maximum applied voltage of 12.6 kV, and (iii) a continuous arc-like regime occurring during the transition to lower voltage and higher current (e.g., U ≈ 3.2 kV, I ≈ 4.07 mA). High-speed shadowgraph images showed a symmetric interfacial cavity whose depth increased nonlinearly with voltage from h 0 ≈ 0.1 mm at 3 kV to h 0 ≈ 2.7 mm at 10.6 kV, reaching a maximum depth of h 0 ≈ 5.9 mm at 12.6 kV, while the cavity disappeared in the continuous-channel regime and was replaced by outward-propagating wave-like motion. The deformation of the water surface during the discharge-water interaction is governed by the balance between electric field forces, surface tension, gravitational forces, and electrohydrodynamic forces, whose relative contributions vary with the applied voltage. This effect was explained by quantitatively estimating the magnitudes of the acting forces and establishing the force balance. Optical emission spectroscopy in the continuous-channel regime indicated air plasma signatures and yielded a gas temperature of approximately 400 K, while prolonged operation (≈ 20 min) increased the water temperature to ~ 70 °C, reduced the water-layer thickness from 6 mm to 3 mm, and decreased pH from 7 to 4. These regime-resolved, quantitative results clarify how the dominant interfacial forcing shifts with discharge mode and provide a mechanistic basis for controlling plasma-water interactions in atmospheric-pressure applications.

atmospheric pressure plasma↗
DOE OSTI2025

Techno-economic and life cycle analyses of the synthesis of a platinum–strontium titanate catalyst

The heterogeneous platinum/strontium titanate (Pt/SrTiO 3 or Pt/STO) catalyst has garnered significant attention as a promising candidate for the hydrogenolysis of polyolefins to hydrocarbon oils. This study evaluates the cost and environmental impacts of a newly developed scalable Pt/STO catalyst production, which includes the synthesis of the STO support and the deposition of Pt onto the support. This two-step synthesis plays a significant role in assessing the commercial feasibility of the catalyst, while energy consumption during the process plays an important role in its environmental impacts. The CatCost and the Research and Development Greenhouse Gases, Regulated Emissions, and Energy use in Technologies (R&D GREET) models were used, respectively, to perform techno-economic analysis (TEA) and life cycle analysis (LCA) of the newly developed catalyst. The TEA showed that the raw materials, accounting for approximately 76% of the total operation cost, has a profound effect on the estimated catalyst cost, mainly due to the platinum precursor. The LCA findings indicated that the catalyst production generates greenhouse gas (GHG) emissions of 66 kg CO 2 e per kg, primarily due to the use of solvents and electricity in the process. The sensitivity analysis indicated that the total operating costs (OpEX), platinum precursor cost, and spent catalyst value (SCV) significantly impact the cost of the synthesized catalyst. Additionally, adopting solvent recovery strategies and using renewable electricity can reduce the GHG emissions of catalyst production to 29 kg CO 2 e per kg.

Heterogeneous Platinum/strontium titanate Catalyst↗
DOE OSTI2021

Bilinear dynamic mode decomposition for quantum control

Abstract Data-driven methods for establishing quantum optimal control (QOC) using time-dependent control pulses tailored to specific quantum dynamical systems and desired control objectives are critical for many emerging quantum technologies. We develop a data-driven regression procedure, bilinear dynamic mode decomposition (biDMD), that leverages time-series measurements to establish quantum system identification for QOC. The biDMD optimization framework is a physics-informed regression that makes use of the known underlying Hamiltonian structure. Further, the biDMD can be modified to model both fast and slow sampling of control signals, the latter by way of stroboscopic sampling strategies. The biDMD method provides a flexible, interpretable, and adaptive regression framework for real-time, online implementation in quantum systems. Further, the method has strong theoretical connections to Koopman theory, which approximates nonlinear dynamics with linear operators. In comparison with many machine learning paradigms minimal data is needed to construct a biDMD model, and the model is easily updated as new data is collected. We demonstrate the efficacy and performance of the approach on a number of representative quantum systems, showing that it also matches experimental results.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

ReachNow EV Driving Data From Seattle, WA, Portland, OR, and New York, NY

ReachNow provided Idaho National Laboratory (INL) with a dataset describing approximately 49,000 trips taken by customers and employees in approximately 100 BMW i3 EVs operating in ReachNow's free-floating car-sharing fleets in Seattle, WA, Portland, OR, and New York, NY between May 2016 and February 2017. Data fields include vehicle rental period start and end timestamps, the location where vehicles were parked at the start and end of rental periods, and distance driven during rental periods. A field categorizing the user during each rental period is also included. This field makes it possible to identify when vehicles were rented by customers and when vehicles were driven by fleet management team employees to reposition, charge, or service the vehicles.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2026

ReachNow EV Driving Data From Seattle, WA, Portland, OR, and New York, NY

ReachNow provided Idaho National Laboratory (INL) with a dataset describing approximately 49,000 trips taken by customers and employees in approximately 100 BMW i3 EVs operating in ReachNow's free-floating car-sharing fleets in Seattle, WA, Portland, OR, and New York, NY between May 2016 and February 2017. Data fields include vehicle rental period start and end timestamps, the location where vehicles were parked at the start and end of rental periods, and distance driven during rental periods. A field categorizing the user during each rental period is also included. This field makes it possible to identify when vehicles were rented by customers and when vehicles were driven by fleet management team employees to reposition, charge, or service the vehicles.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2021

Improving the Cyber and Physical Security Posture of the Electric Sector (Final Report)

Cooperative electric utilities represent an integral part of the larger electric grid and are part of the nation’s critical infrastructure. The National Rural Electric Cooperative Association (NRECA) has a unique relationship with approximately 900 cooperatively owned and operated electric utilities and engaged in a program with the Department of Energy (DOE) to promote a culture of cybersecurity and resiliency within the electric cooperative community. The Rural Cooperative Cybersecurity Capabilities (RC3) Program, funded under a Cooperative Agreement with DOE (Project DE-OE-0000807), focused on improving the cybersecurity and resiliency capabilities of small and mid-sized electric distribution cooperatives. This segment of electric utilities faces many challenges, but also embraces a culture of cooperation that presents opportunities. A customized approach is needed to reach these utilities – one that emphasizes collaboration, more focused and personalized training, use of trusted and familiar experts that can be deployed as needed, software security services that require limited in-house cybersecurity expertise, and shared resource models that enable access to more expensive cybersecurity options.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2023

Theory and Simulation of Ultrafast Multidimensional Nonlinear X-ray Spectroscopy of Molecules (Final Report)

Emerging X-ray free electron laser (XFEL) beam sources offer new types of probes of matter with unprecedented spatial and temporal resolutions. These experimental advances must be met by robust theoretical and computational tools that provide predictive modeling capacity of the underlining electronic and structural dynamics. The latter will be essential for the design of sophisticated multi-pulse experiments and for their interpretation. The proposed research effort will focus on developing cutting- edge simulation tools for nonlinear multidimensional X-ray/optical spectroscopies and aims to address key questions in Priority Research Opportunities 1 (Probing and controlling electron motion within a molecule) and 3 (Capturing rare events and intermediate states in the transformation of matter) as mentioned in the BES Roundtable Report “Opportunities for Basic Research at the Frontiers of XFEL Ultrafast Science”. XFEL multidimensional nonlinear techniques, which combine sequences of X-ray and possibly optical pulses, provide a unique experimental toolbox for probing the dynamics of core and valence electronic excitations, as well as material structure. Predictive modeling of these dynamical processes requires the combination of analytical theory for nonlinear interactions of light and matter, robust quantum-chemical methodologies for the accurate description of electronic structure of various materials, and multiscale ab initio electron and nuclear dynamics techniques operating beyond Born-Oppenheimer approximation. These challenges will be addressed with three research thrusts (i) Develop and implement theoretical apparatus for modeling a broad range of multidimensional spectroscopic techniques enabled by present and upcoming XFEL facilities. This thrust also includes the incorporation of a computational module in the DOE supported open-source NWChem computational chemistry package as well as the development of other open-source codes ready for dissemination across a broad user base; (ii) Propose and design new multi-pulse experiments that make use of the capabilities of the incoming LCLS-II facility; (iii) Perform selected applications to specific molecular systems that can be carried out at LCLS-II and demonstrate how these X-ray sources may be used to study nonadiabatic dynamics through conical interactions, electronic correlations in multi-core excitons, and charge transfer/energy transfer processes. The proposed research will be carried out by a multi-disciplinary four-institution research team which combines academia and national laboratories and spans the broad and necessary expertise in theoretical spectroscopy, nonlinear optics, quantum chemistry, molecular non-adiabatic dynamics and code development. The work will be performed in a highly interactive team environment with junior researchers shared between institutions thus cementing cross- disciplinary interactions. The developed simulation tools will be immediately deployed for XFEL facility users, both experimentalists and theorists, via freely distributed codes and databases. Altogether, this project will facilitate establishing XFEL-based multidimensional spectroscopies as a novel diagnostic tool for monitoring electronic and structural dynamics in molecular materials.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2023

SPUTTERED THIN FILMS FOR VERY HIGH POWER, EFFICIENT, AND LOW-COST COMMERCIAL SOFCS

The aim of this project was to leverage the low area specific resistance (ASR) of Redox’s GDC electrolyte-based cell architecture and increase solid oxide fuel cell (SOFC) efficiency (i.e., open circuit voltage, or OCV) without significantly increasing cell resistance. The key to achieving the increased SOFC efficiency is the introduction of a thin sputtered yttria stabilized zirconia (YSZ) electron-blocking layer and a thin sputtered gadolinia doped ceria (GDC) barrier layer on top of the half-cell substrate (i.e., anode and GDC electrolyte). A key initial effort of the project was to improve the quality of the half-cell substrate surface so that any remaining defects were significantly smaller than the desired film thickness. During lab-scale trials, we had to overcome challenges with film cracking during post-sputtering treatments (e.g., thermal anneals) as well as damage to the sputtering targets. After the initial lab-scale trials, the focus shifted toward the use of commercial scale sputtering equipment with a sputtering equipment manufacturer in the microelectronics industry. Using the commercial sputtering equipment, we deposited films with different variations of power, PO2, sputtering time, platen speed, etc. We then determined the combinations of sputtering conditions and post-sputtering treatments that yielded high-quality films without cracks or other significant defects. While thin films were deposited on cells as large as 10 cm by 10 cm, the processing was optimized using 4 cm by 4 cm cells. Cell performance was evaluated in stainless-steel test fixtures between 500 °C and 700 °C with hydrogen fuel fed to the anode and air fed to the cathode. Extensive studies allowed us to determine that modified cathode and cathode contact firing processes were required to achieve theoretical OCV. Moreover, use of the new firing processes resulted in the need for a modified cathode contact to achieve a low ASR. In summary, the project demonstrated the performance of high OCV (1.13 V at 650 °C) from sputtered layers and a low ASR (~0.25 Ohms-cm2 at 650 °C) resulting from a modification of cathode/contact processing and the introduction of alternative contact layers that are sufficient to yield a Gen-1 cell with a maximum power density of approximately 1.2 W/cm2. At an operating voltage of 0.74 V, this would yield a cell power density of ~1.1 W/cm2. While not utilized in this project, a Redox Gen-2 cell has a catalyst-infiltrated porous anode that reduces the ASR by more than 50% from that of the Gen-1 cells used in this project. Therefore, if the sputtered YSZ electron-blocking layer and GDC barrier layer are added to a Gen-2 half cell with a similar increase in OCV to the theoretical value of ~1.13 V at 650 °C, and if the same improvement in ASR (from that demonstrated in this project) is achieved when using the Gen-2 half-cell architecture as a sputtered cell substrate, then the power density at 0.74 Vop could be as high as ~2.6 W/cm2. The impact of such power density gains, while still maintaining high cell efficiency, and thus high system efficiency, is a dramatic decrease in system cost because the stack represents ~30-40% of the SOFC system cost.

01 COAL, LIGNITE, AND PEAT↗