Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “time minimization”

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 973 records · Page 54

PDPTW-DB: MILP-Based Offline Route Planning for PDPTW with Driver Breaks

The Pickup and Delivery Problem with Time Windows (PDPTW) involves optimizing routes for vehicles to meet pickup and delivery requests within specific time constraints, a challenge commonly faced in logistics and transportation. Microtransit, a flexible and demand-responsive service using smaller vehicles within defined zones, can be effectively modeled as a PDPTW. Yet, the need for driver breaks—a key human constraint—is frequently overlooked in PDPTW solutions, despite being necessary for regulatory compliance. This study presents a novel mixed-integer linear programming formulation for the Pickup and Delivery Problem with Time Windows and Driver Breaks (PDPTW-DB). To the best of our knowledge this formulation is the first to consider mandatory periodic driver breaks within optimized Microtransit routes. The proposed model incorporates regulatory compliant break scheduling directly within the vehicle routing optimization framework. By considering driver break requirements as an integral component of the optimization process, rather than as a post-processing step, the model enables the generation of routes that respect hours of service regulations while minimizing operational costs. This integrated approach facilitates the generation of schedules that are operationally efficient and prioritize driver welfare through driver breaks. We work with a public transit agency from the southern USA, and highlight the specific nuances of driver break optimization, and present a Pickup and Delivery Problem with Time Windows formulation for optimizing Microtransit operations and scheduling driver breaks. We validate our approach using real-world data from the transit agency. Our results validate our formulation in producing cost-effective, and regulation-compliant solutions.

Applied Computing, Transportation↗

Probe the Localized Electrochemical Environment Effects and Electrode Reaction Dynamics for Metal Batteries using In Situ 3D Microscopy

Uncontrollable dendrite growth is closely related to non-uniform reaction environments. However, there is a lack of understanding and analysis methods to probe the localized electrochemical environment (LEE). Here the effects of the LEE are investigated, including localized ion concentrations, current density, and electric potential, on metal plating/stripping dynamics and dendrite minimization. A novel in situ 3D microscopy technique is developed to image the morphology dynamics and deposition rate of Zn plating/stripping processes on 3D Zn–Mn anodes. Using the in situ 3D microscope, the electrode morphology changes during the reactions are directly imaged and Zn deposition rate maps at different time points are obtained. It is found that reaction kinetics are highly correlated to LEE and electrode morphology. To further quantify the LEE effects, the digital twin technique is employed that allows the accurate calculation of the electrochemical environments, such as localized ion concentrations, current density, and electric potential, which cannot be directly measured from experiments. We found that the curvature of the 3D electrode surface determines the LEE and significantly influences reaction kinetics. This provides a new strategy to minimize the dendrite formation by designing and optimizing the 3D geometry of the electrode to control the LEE.

36 MATERIALS SCIENCE↗

Models and Strategies for Optimal Demand Side Management in the Chemical Industries

Deregulation and the increase of renewable electricity generation from wind and solar photovoltaics have transformed the U.S. electricity market. Economic and environmental benefits notwithstanding, the presence of renewables has increased variability and uncertainty on the supply side of the grid. Managing demand, rather than generation – a strategy referred to as “demand response (DR)” – is an attractive approach for mitigating this imbalance. DR efforts aim to reduce electricity usage during peak demand times, lessening stress on the grid. Industrial users are particularly attractive entities for DR participation since they present large, localized loads that can provide significant relief on grid demand and –unlike other large loads, such as buildings – are minimally dependent on human needs and preferences. In this project, we accomplished three main objectives. (1) We developed data-driven low-order DR scheduling-relevant dynamic models of chemical processes. Concurrently, we studied the formulation and solution of the associated optimal DR production scheduling problems. (a) A prototype air separation unit (ASU) model was used to generate simulated operating data for initial modeling efforts, which enabled the later use of industrial data for data-driven modeling. (b) We utilized Hammerstein-Wiener (HW) and Finite Step Response (FSR) models to represent nonlinear plant dynamics. (c) The HW models were linearized using exact linearization so they could potentially be embedded in power system models, which are formulated as mixed integer linear programs (MILPs). (d) We solved DR optimization problems under uncertainty and found that even naïve predictions of electricity price and product demand led to significant cost savings benefits. (2) Our DR scheduling optimization problem formulations are amenable to real-time solution. (a) We utilized Lagrangian Relaxation (LR) to efficiently solve the optimization problem by decoupling subproblems linked by complicating constraints. (b) We have achieved computation times for the 3-day DR scheduling problem of an ASU as low as 1.88 minutes. (3) Our representations of the DR behavior of chemical process as grid-level batteries were embedded in power system models. (a) For a small-scale grid, we found that incorporating the dynamics of the chemical plant in the optimal power flow calculations resulted in better resource management leading to up to 15% and 46% cost reduction for the grid and chemical plant operations, respectively, during periods of power line congestion. We have published several works dedicated to modeling and solving DR optimization problems from the user side. These were published in top peer-reviewed journals and are summarized in this report. The most recent work (and papers in preparation) considers DR scheduling from the grid side. Future efforts will consider networked plants (e.g., air separation units operating on a common pipeline) for DR participation, which is expected to amplify the capabilities of industrial DR participants to perform load-shifting. Our consideration of uncertainty in DR has inspired future directions in this area as well: we plan to develop multistage methods to fully account for the effects of uncertainty in DR scheduling.

24 POWER TRANSMISSION AND DISTRIBUTION↗

An adaptive space-time phase field formulation for dynamic fracture of brittle shells based on LR NURBS

Here, we present an adaptive space-time phase field formulation for dynamic fracture of brittle shells. Their deformation is characterized by the Kirchhoff–Love thin shell theory using a curvilinear surface description. All kinematical objects are defined on the shell’s mid-plane. The evolution equation for the phase field is determined by the minimization of an energy functional based on Griffith’s theory of brittle fracture. Membrane and bending contributions to the fracture process are modeled separately and a thickness integration is established for the latter. The coupled system consists of two nonlinear fourth-order PDEs and all quantities are defined on an evolving two-dimensional manifold. Since the weak form requires C 1 -continuity, isogeometric shape functions are used. The mesh is adaptively refined based on the phase field using Locally Refinable (LR) NURBS. Time is discretized based on a generalized-α method using adaptive time-stepping, and the discretized coupled system is solved with a monolithic Newton–Raphson scheme. Finally, the interaction between surface deformation and crack evolution is demonstrated by several numerical examples showing dynamic crack propagation and branching.

42 ENGINEERING↗

Supermassive primordial black holes from inflation

There is controversy surrounding the origin and evolution of our universe's largest supermassive black holes (SMBHs). In this study, we consider the possibility that some of these black holes formed from the direct collapse of primordial density perturbations. Since the mass of a primordial black hole is limited by the size of the cosmological horizon at the time of collapse, these SMBHs must form rather late, and are naively in conflict with constraints from CMB spectral distortions. These limits can be avoided, however, if the distribution of primordial curvature perturbations is highly non-Gaussian. After quantifying the departure from Gaussianity needed to evade these bounds, we explore a model of multi-field inflation — a non-minimal, self-interacting curvaton model — which has all the necessary ingredients to yield such dramatic non-Gaussianities. We leave the detailed model building and numerics to a future study, however, as our goal is to highlight the challenges associated with forming SMBHs from direct collapse and to identify features that a successful model would need to have. This study is particularly timely in light of recent observations of high-redshift massive galaxy candidates by the James Webb Space Telescope as well as evidence from the NANOGrav experiment for a stochastic gravitational wave background consistent with SMBH mergers.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Development of Advanced, Radiation Resistant, Optical-based Detector Technology for Future Experiments.

The primary objective of this project has been to advance the design of high-performance electromagnetic (EM) calorimeters for future particle physics experiments, to identify and measure the timing, position and energy of electrons, positrons and gamma rays, particularly in high-luminosity environments with intense radiation and pileup conditions. To meet such challenges, the proposed research has focused on the development of ultra-compact, radiation-hard calorimeter modules, to provide excellent timing, spatial, and energy resolution. The work aligns with the DOE’s Basic Research Needs (BRN) for High Energy Physics (HEP) Instrumentation and the research team contributes actively to the Coordinating Panel on Advanced Detectors (CPAD) RDC9 calorimetry collaboration in the USA and the European Committee on Future Accelerators (ECFA) DRD-CALO calorimetry collaboration at CERN, the European Laboratory for Particle Physics located in Geneva, Switzerland. The research builds on the RADiCAL (radiation-hard, ultra-compact) modular sampling calorimeter approach, developed by the research team, which employs dense and very bright optical materials such as LYSO:Ce scintillator plates that are interleaved with very dense tungsten plates to minimize detector size while optimizing performance. The modules are comparable in size to a human index finger, dimensionally 14 mm x 14 mm in cross section and 135 mm in length. And despite the small size, the structure is capable of providing excellent timing and energy resolution. This is facilitated through the use of specialized quartz capillaries filled with wavelength-shifting filaments, positioned at various depths along the length of a module, to collect and guide light signals to silicon photomultipliers (SiPMs) which detect and convert the optical signals to electronic signals for analysis. The primary goals of this project have been: (1) Achieve a timing resolution to σ t ≤ 30 ps for high-energy electrons and photons, important for their association with specific events produced in colliding-beam experiments and for the detection of decays-in-flight of long-lived particles. The project has achieved this goal in beam tests of a single RADiCAL module at CERN, during which a timing resolution of σ t = 27 ps was measured for electrons of energy E = 150 GeV. Based upon a mathematical fit to the data measured over a broad energy range from low energy to high energy, a resolution of σ t ≤ 18 ps has been estimated for electrons of very high (TeV) energy. From these measurements and with further expected technical improvements, the timing resolution should reach σ t ≤ 10 ps, important for searches for discovery physics in upcoming and future experiments. (2) Achieve an energy resolution of σ E / E ≤ 10% / $\sqrt{E}$. The project has yet to achieve this goal, but is close to it, having measured a value of σ E / E ≤ 15.9% / $\sqrt{E}$ using a modular array. Ultimately, the resolution goal is expected to be reached by adjustments to material thicknesses within the modules, which will improve the sampling fraction to measure more precisely the shower energy for lower energy particles. The versatility of the modular RADiCAL approach enables the testing of advanced materials, photosensors and electronics, developed in collaboration with CPAD RDC and ECFA DRD-CALO groups. The structure can distinguish electrons, positrons and gamma rays from hadrons and muons and beam-induced backgrounds, making it a valuable tool in a variety of detector environments, including future circular colliders (FCC-ee, FCC-hh) proposed for the European Laboratory for Particle Physics (CERN), the muon-collider proposed for Fermi National Accelerator Laboratory (Fermilab), and searches for new physics in beam-dump, fixed target and forward-physics experiments. And, while designed with particle physics applications in mind, the technologies developed in this project have the potential for application more broadly in particle and nuclear physics, materials science, and medical physics, underscoring the far-reaching potential of this line of instrumentation research and development.

47 OTHER INSTRUMENTATION↗

Parity solution to the strong CP problem and a unified framework for inflation, baryogenesis, and dark matter

It has been known for some time that asymptotic parity invariance of weak interactions can provide a solution to the strong CP problem without the need for the axion. Left-right symmetric theories which employ a minimal Higgs sector consisting of a left-handed and a right-handed doublet is an example of such a theory wherein all fermion masses arise through a generalized seesaw mechanism. In this paper we present a way to understand the origin of matter-antimatter asymmetry as well as the dark matter content of the universe in these theories using the Affleck-Dine (AD) leptogenesis mechanism and inflaton decay, respectively. Three gauge singlet fermions are needed for this purpose, two of which help to implement the Dirac seesaw for neutrino masses while the third one becomes the non-thermal dark matter candidate. A soft lepton number breaking term involving the AD scalar field is used to generate lepton asymmetry which suffers no wash-out effects and maintains the Dirac nature of neutrinos. This framework thus provides a unified description of many of the unresolved puzzles of the standard model that require new physics.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Modulating thermal load through lightweight residential building walls using thermal energy storage and controlled precooling strategy

Precooling is a recognized technique for reducing cooling energy in buildings during peak hours by shifting load to off-peak hours. This technique is particularly effective in buildings with high thermal mass, because of their large thermal energy storage capacity, and in commercial buildings due to their variable electricity pricing based on time-of-use rates. Precooling in residential buildings has been a matter of limited interest in the past because of their low thermal mass and typically uniform electricity pricing rate. While previous studies on precooling primarily focused on cost savings, an important aspect of precooling is the thermal load modulation, which could be very effective in managing peak demand in lightweight residential buildings integrated with thermal energy storage systems. In this study, we examine different precooling strategies to manage the heat gains in lightweight building walls integrated with phase-change materials. We create nine different precooling profiles by controlling the interior temperature, and then evaluate the influence of the precooling profiles on four key building energy performance parameters: total heat gain, peak heat gain, maximum heat gain during peak hours, and time at which peak occurs. To thoroughly understand the fundamental physics, we first consider hypothetical climates and obtain the optimal precooling strategy required to achieve maximum peak shedding and shifting while minimizing the total heat gains. We then extend the model to Baltimore, Maryland, and estimate the benefits of the optimized precooling strategy under real conditions. The optimal precooling strategy proposed in this study can shift the peak heat gain by up to 14 hours, thereby reducing the heat gain during peak period by up to 95%, at the expense of a 23% increase in the total heat gains.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

GMFOLD: Subgraph matching for high-throughput DNA-aptamer secondary structure classification and machine learning interpretability

Aptamers are oligonucleotide receptors that bind to their targets with high affinity. Here, we consider aptamers comprised of single-stranded DNA that undergo target-binding-induced conformational changes, giving rise to unique secondary and tertiary structures. Given a specific aptamer primary sequence, there are well-established computational tools (notably mfold) to predict the secondary structure via free energy minimization algorithms. While mfold generates secondary structures for individual sequences, there is a need for a high-throughput process whereby thousands of DNA structures can be predicted in real-time for use in an interactive setting, when combined with aptamer selections that generate candidate pools that are too large to be experimentally interrogated. We developed a new Python code for high-throughput aptamer secondary structure determination (GMfold). GMfold uses subgraph matching methods to group aptamer candidates by secondary structure similarities. We also improve an open-source code, SeqFold, to incorporate subgraph matching concepts. We represent each secondary structure as a lowest-energy bipartite subgraph matching of the DNA graph to itself. These new tools enable thousands of DNA sequences to be compared based on their secondary structures, using machine-learning algorithms. This process is advantageous when analyzing sequences that arise from aptamer selections via systematic evolution of ligands by exponential enrichment (SELEX). This work is a building block for future machine-learning-informed DNA-aptamer selection processes to identify aptamers with improved target affinity and selectivity and advance aptamer biosensors and therapeutics.

Aptamer↗

Ultra-sensitive radon assay using an electrostatic chamber in a recirculating system

Rare event searches such as neutrinoless double beta decay and Weakly Interacting Massive Particle detection require ultra-low background detectors. Radon contamination is a significant challenge for these experiments, which employ highly sensitive radon assay techniques to identify and select low-emission materials. This work presents the development of ultra-sensitive electrostatic chamber (ESC) instruments designed to measure radon emanation in a recirculating gas loop, for future lower background experiments. Unlike traditional methods that separate emanation and detection steps, this system allows continuous radon transport and detection. This is made possible with a custom-built recirculation pump. A Python-based analysis framework, PyDAn, was developed to process and fit time-dependent radon decay data. Radon emanation rates are given for various materials measured with this instrument. A radon source of known activity provides an absolute calibration, enabling statistically-limited minimal detectable activities of 20 µBq. These devices are powerful tools for screening materials in the development of low-background particle physics experiments.

47 OTHER INSTRUMENTATION↗

Prompt fission neutron multiplicity in the 239 Pu(n, f) reaction and its energy dependence

Measuring prompt fission neutrons to high precision is an experimental challenge, especially for radioactive fissioning nuclides. However, accurate average multiplicities, $\overline{v}_{p}$, and kinetic energy distributions of prompt fission neutrons are essential for fundamental and applied nuclear physics. We present here a recent measurement of the 239 Pu (n,f) $\overline{v}_{p}$ as a function of the incident-neutron energy, over the range 1-700 MeV. The measurement was performed with a cutting-edge setup and an innovative technique, which allowed to minimize and account for the main sources of bias. An unprecedented precision was therefore achieved. Our data are compared to GEF predictions as well as to evaluated libraries. For the first time, at low energies, the ENDF/B-VIII.0 nuclear data evaluation is validated with an independent measurement and the evaluated uncertainty reduced by up to 60%. This work paves the way to precisely measure prompt fission neutron multiplicities on highly radioactive nuclei.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Symbolic pregression: Discovering physical laws from distorted video

In this work, we present a method for unsupervised learning of equations of motion for objects in raw and optionally distorted unlabeled synthetic video (or, more generally, for discovering and modeling predictable features in time-series data). We first train an autoencoder that maps each video frame into a low-dimensional latent space where the laws of motion are as simple as possible, by minimizing a combination of nonlinearity, acceleration, and prediction error. Differential equations describing the motion are then discovered using Pareto-optimal symbolic regression. We find that our pre-regression (“pregression”) step is able to rediscover Cartesian coordinates of unlabeled moving objects even when the video is distorted by a generalized lens. Using intuition from multidimensional knot theory, we find that the pregression step is facilitated by first adding extra latent space dimensions to avoid topological problems during training and then removing these extra dimensions via principal component analysis. An inertial frame is autodiscovered by minimizing the combined equation complexity for multiple experiments.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Streaming Generalized Canonical Polyadic Tensor Decompositions

In this paper, we develop a method which we call OnlineGCP for computing the Generalized Canonical Polyadic (GCP) tensor decomposition of streaming data. GCP differs from traditional canonical polyadic (CP) tensor decompositions as it allows for arbitrary objective functions which the CP model attempts to minimize. This approach can provide better fits and more interpretable models when the observed tensor data is strongly non-Gaussian. In the streaming case, tensor data is gradually observed over time and the algorithm must incrementally update a GCP factorization with limited access to prior data. In this work, we extend the GCP formalism to the streaming context by deriving a GCP optimization problem to be solved as new tensor data is observed, formulate a tunable history term to balance reconstruction of recently observed data with data observed in the past, develop a scalable solution strategy based on segregated solves using stochastic gradient descent methods, describe a software implementation that provides performance and portability to contemporary CPU and GPU architectures and integrates with Matlab for enhanced usability, and demonstrate the utility and performance of the approach and software on several synthetic and real tensor data sets.

97 MATHEMATICS AND COMPUTING↗

AI-Assisted Conceptual Development of a Pre-Geometric Cosmological Model - An Exercise in AI-Assisted Conceptual Framework Generation, Paper I: Foundations and Replication Dynamics

We develop a pre-geometric cosmological framework in which existence is identified with a finite amount of unstructured energy possessing vibration as its only intrinsic property. This vibrational substrate occupies an open, bounded spectral interval $(\omega_{\min},\omega_{\max})$, ensuring finiteness of total energy and excluding infinitely stable configurations. The substrate evolves under two fundamental and competing tendencies---excitation, which amplifies coherence, and randomization, which scrambles it. Their balance produces a metastable unstructured regime in which rare fluctuations may form long-lived self-consistent spectral configurations. Because the substrate is finite and subject to competing order--disorder dynamics, no coherent configuration can be perpetually stable. We show that the only mechanism capable of sustaining long-lived organization is a replication instability: a coherent unit may reproduce into multiple offspring according to a general $1\!\to n$ rule. Replication consumes energy from the finite substrate, breaks the metastable symmetry, and induces a discrete notion of event time through the replication tick $\Delta\tau$. Temporal succession is defined through correlation ordering of spectral microstates, producing an intrinsic pre-causal structure. The compactness of the spectral domain imposes minimal and maximal timescales, bounds the internal coherence of emergent units, and limits their proliferation. These spectral constraints serve as precursors for the emergence of geometry, adjacency, and a limiting propagation speed, developed in subsequent papers of this series. Paper I provides the foundational axioms (PG1--PG9) governing the spectral substrate, its metastable dynamics, the formation of coherent units, and the necessity of replication, establishing a fully pre-geometric stage from which causal and geometric structure naturally emerge.

79 ASTRONOMY AND ASTROPHYSICS↗

Boundary Layer Exploration of Aerosols and Clouds ON Ships (BEACONS)

BEACONS aims to demonstrate autonomous deployment of a shipborne system that routinely observes atmospheric, cloud, and aerosol properties. These observations enable scientists to study aerosol-cloud interactions, advancing energy resilience and the Department of Energy's mission to improve model prediction. Key Project Outcomes Autonomous deployment of aerosol and cloud measurement systems aboard a Pasha Hawaii Marjorie C commercial ship. Six-month campaign for continuous, high-resolution aerosol and cloud data collection with minimal crew intervention. Development and testing of instruments in two phases, advancing from essential to complex systems. AI-ready processed data on the BEACONS Data Platform for near-real-time data access. Validation of aerosol-cloud interaction hypotheses for improved modeling. Establishment of best practices for instrument deployment in marine environments. Scientific collaboration to support the U.S. Department of Energy (DOE) Biological and Environmental Research (BER) mission.

54 ENVIRONMENTAL SCIENCES↗

Dilute Combustion Control Using Spiking Neural Networks

Dilute combustion with exhaust gas recirculation (EGR) in spark-ignition engines presents a cost-effective method for achieving higher levels of engine efficiency. At high levels of EGR, however, cycle-to-cycle variability (CCV) of the combustion process is exacerbated by sporadic occurrences of misfires and partial burns. Previous studies have shown that temporal deterministic patterns emerge at such conditions and certain combustion cycles have a significant influence over future events. Due to the complexity of the combustion process and the nature of CCV, harnessing all the deterministic information for control purposes has remained challenging even with physics based 0-D, 1-D, and high-fidelity computational fluid dynamics (CFD) models. In this study, we present a data-driven approach to optimize the combustion process by controlling CCV adjusting the cycle-to-cycle fuel injection quantity. Readily available data from in-cylinder pressure was used to train a spiking neural network (SNN) which learns the optimal way to manage fuel injection in order to reduce CCV while maintaining acceptable levels of fuel consumption. SNNs are particularly well suited for powertrain control applications due to their ability to be deployed on FPGA-based neuromorphic hardware which are small, inexpensive, and have a low power demand. The high-performance computing (HPC) resources of Oak Ridge National Laboratory were used to run an evolutionary-based training approach for choosing the best SNN configuration that minimizes the size of the network while achieving the desired goal. The neuromorphic hardware with the optimized SNN deployed was connected to the rapid prototyping engine control system for real-time control implementation and tested on a single cylinder version of a GM LNF 4-cylinder engine. The results show a significant reduction of CCV with a small percentage of additional fuel used to stabilize the charge.

33 ADVANCED PROPULSION SYSTEMS↗

Ultra-sensitive radon assay using an electrostatic chamber in a recirculating system

Rare event searches such as neutrinoless double beta decay and Weakly Interacting Massive Particle detection require ultra-low background detectors. Radon contamination is a significant challenge for these experiments, which employ highly sensitive radon assay techniques to identify and select low-emission materials. This work presents the development of ultra-sensitive electrostatic chamber (ESC) instruments designed to measure radon emanation in a recirculating gas loop, for future lower background experiments. Unlike traditional methods that separate emanation and detection steps, this system allows continuous radon transport and detection. This is made possible with a custom-built recirculation pump. A Python-based analysis framework, PyDAn, was developed to process and fit time-dependent radon decay data. Radon emanation rates are given for various materials measured with this instrument. A radon source of known activity provides an absolute calibration, enabling statistically-limited minimal detectable activities of 20 uBq. These devices are powerful tools for screening materials in the development of low-background particle physics experiments.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Feedstock Harvesting and Storage: Post-Harvest Management for Quality Preservation

Research on harvesting, collection, and storage has matured, and our understanding of feedstock supply challenges has led us to the point of our current goals. The goals of this project include using storage time to collect impactful data of critical feedstock properties in storage and to characterize the feedstock before delivery. This data will inform delivery scheduling of biomass and relevant compositional information to minimize variations in the material quality, which are critical to the preprocessing and conversion performance.

09 BIOMASS FUELS↗