Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “deep well”

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 469 records · Page 26

Seismic monitoring of a small CO 2 injection using a multi-well DAS array: Operations and initial results of Stage 3 of the CO 2 CRC Otway project

Active time-lapse seismic is widely employed for monitoring CO 2 geosequestration due to its ability to track the distribution of fluids in space and time. However, standard 4D seismic monitoring suffers from several challenges, including high cost, disruption to other land uses, and, consequently, relatively large intervals between monitor surveys. Some of these challenges can be mitigated using permanently installed sources and receivers. Such an approach was tested at the CO 2 CRC Otway site by continuous offset VSP monitoring of 15,000 t of supercritical CO 2 injected into an aquifer 1,500 m deep with nine permanent seismic sources (surface orbital vibrators or SOVs) and five downhole fibre-optic receivers. This continuous monitoring is complemented by multi-well 4D VSP using a mobile vibroseis source and the same DAS receivers, which included one baseline and two monitor surveys after injection of 4,000 and 12,000 t of CO 2 . The continuous DAS-SOV monitoring detected an abrupt increase of travel times below the injection interval on the second day of injection (after injection of 300 t of CO 2 ) and tracked the growth of the areal CO 2 plume by mapping changes of reflection amplitudes. The plume is also detected by time-lapse changes of reflection amplitudes in multi-well 4D VSPs. The plume images obtained from continuous offset VSP and 4D VSP are broadly consistent with each other but with some differences due to differences in illumination, lateral variations of velocities and seismic anisotropy. Furthermore, these differences also serve as a measure of uncertainty of 4D VSP images.

4D VSP↗

GETaLM: A generator for electron tagger and luminosity monitor for electron - proton and ion collisions

The study of elastic bremsstrahlung and electron tagging in electron-proton or ion collisions is gaining importance with the planned construction of several experimental facilities focused on deep-inelastic scattering (DIS) measurements. This paper describes a program which generates bremsstrahlung photons in electron-proton and electron-ion interactions as well as scattered electrons in bremsstrahlung processes and in a quasi-real photon approximation to the general DIS process. Here, the effects of electron beam divergence and the spread of the interaction vertex are implemented. The program can be used as an input to simulations of instrumentation for bremsstrahlung photon detection, luminosity measurements, electron tagging, and the determination of the cross sections of corresponding processes.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Tuned Reactivity at the Lithium Metal–Argyrodite Solid State Electrolyte Interphase

Thin intermetallic Li 2 Te–LiTe 3 bilayer (0.75 µm) derived from 2D tellurene stabilizes the solid electrolyte interphase (SEI) of lithium metal and argyrodite (LPSCl, Li 6 PS 5 Cl) solid-state electrolyte (SSE). Tellurene is loaded onto a standard battery separator and reacted with lithium through single-pass mechanical rolling or transferred directly to SSE surface by pressing. State-of-the-art electrochemical performance is achieved, e.g., symmetric cell stable for 300 cycles (1800 h) at 1 mA cm -2 and 3 mAh cm -2 (25% DOD, 60 µm foil). Cryo-stage focused ion beam (Cryo-FIB) sectioning and Raman mapping demonstrate that the Li 2 Te–LiTe 3 bilayer impedes SSE decomposition. The unmodified Li–LPSCl interphase is electrochemically unstable with a geometrically heterogeneous reduction decomposition reaction front that extends deep into the SSE. Decomposition drives voiding in Li metal due to its high flux to the reaction front, as well as voiding in the SSE due to the associated volume changes. Analysis of cycled SSE found no evidence for pristine (unreacted) lithium metal filaments/dendrites, implying failure driven by decomposition phases with sufficient electrical conductivity that span electrolyte thickness. In conclusion, DFT calculations clarify thermodynamic stability, interfacial adhesion, and electronic transport properties of interphases, while mesoscale modeling examines interrelations between reaction front heterogeneity (SEI heterogeneity), current distribution, and localized chemo-mechanical stresses.

25 ENERGY STORAGE↗

Resolution dependence of the turbulent atmospheric boundary layer in global storm-resolving climate simulations

The current generation of state-of-the-art global climate models are being run at increasingly higher horizontal resolutions, with the goal of resolving organised deep convection explicitly. How a kilometre-scale resolution impacts the representation of the atmospheric boundary layer is, however, not well known. Using statistical analysis on global fields as well as high-frequency data at selected locations, produced with the Integrated Forecasting System (IFS) model for the Next Generation Earth-system Models (nextGEMS) project, we investigate the horizontal resolution dependence of some boundary-layer processes. We find that a change in resolution from 9 to 2.8 km causes no substantial changes to boundary-layer properties and processes at most of the locations studied, although some global changes are detected that indicate circulation changes. Small changes to the boundary-layer depth and structure are found in the Tropics. The short simulation length and lack of data for optimal boundary-layer analysis limits the conclusions, especially in relation to the connection between the boundary layer and the atmosphere general circulation.

54 ENVIRONMENTAL SCIENCES↗

Tracking natural CO 2 migration through a sandstone aquifer using Sr, U and C isotopes: Chimayó, New Mexico, USA

The geochemical and isotopic characteristics of groundwaters in Chimayó, New Mexico, reflect processes that affect water quality in the Tesuque Aquifer, which overlies a leaking natural CO 2 source in a structurally complex region. In this study, select isotopes (δ 13 C, 87 Sr/ 86 Sr, 234 U/ 238 U) are applied to groundwaters to better understand CO 2 transport mechanisms and related water-rock interactions that impact Chimayó groundwater. Carbon stable isotope ratios of dissolved inorganic carbon (DIC; δ 13 CDIC = -15.10‰ to 4.50‰) identify a distinct source of upward-migrating CO 2 that interacts with the groundwater. Additionally, groundwater 87 Sr/ 86 Sr compositions (0.7098 to 0.7154) reflect intrusion of varying amounts of saline water associated with the high CO 2 source, while 234 U/ 238 U ratios corroborate presence of deep groundwater and suggest impacts from distinctive natural uranium sources are affecting groundwater. Previous work proposed two CO 2 transport mechanisms at the site: (1) dissolved in deep brine that underlies the aquifer and (2) in the gas phase; this study uses isotope mixing models to identify wells that are affected by these CO 2 transport mechanisms and demonstrates that both transport mechanisms are associated with impaired groundwaters. Overall, this study demonstrates that applying multiple isotope systems (δ 13 CDIC, 87 Sr/ 86 Sr, 234 U/ 238 U) is a dynamic tool for identifying and measuring the impact of CO 2 leakage from a sequestration site.

58 GEOSCIENCES↗

Evidence of a toroidal magnetic field in the core of 3C 84

The spatial scales of relativistic radio jets, probed by relativistic magneto-hydrodynamic (RMHD) jet launching simulations and by most very long baseline interferometry (VLBI) observations differ by an order of magnitude. Bridging the gap between these RMHD simulations and VLBI observations requires selecting nearby active galactic nuclei (AGN), the parsec-scale region of which can be resolved. The radio source 3C 84 is a nearby bright AGN fulfilling the necessary requirements: it is launching a powerful, relativistic jet powered by a central supermassive black hole, while also being very bright. Using 22 GHz globe-spanning VLBI measurements of 3C 84 we studied its sub-parsec region in both total intensity and linear polarisation to explore the properties of this jet, with a linear resolution of ~0.1 parsec. We tested different simulation set-ups by altering the bulk Lorentz factor Γ of the jet, as well as the magnetic field configuration (toroidal, poloidal, helical). We confirm the persistence of a limb brightened structure, which reaches deep into the sub-parsec region. The corresponding electric vector position angles (EVPAs) follow the bulk jet flow inside but tend to be orthogonal to it near the edges. Our state-of-the-art RMHD simulations show that this geometry is consistent with a spine-sheath model, associated with a mildly relativistic flow and a toroidal magnetic field configuration.

3C 84 (NGC 1275)↗

A search for cosmogenic neutrinos with the ARIANNA test bed using 4.5 years of data

The primary mission of the ARIANNA ultra-high energy neutrino telescope is to uncover astrophysical sources of neutrinos with energies greater than 10 16 eV. A pilot array, consisting of seven ARIANNA stations located on the surface of the Ross Ice Shelf in Antarctica, was commissioned in November 2014. We report on the search for astrophysical neutrinos using data collected between November 2014 and February 2019. A straight-forward template matching analysis yielded no neutrino candidates, with a signal efficiency of 79%. We find a 90% confidence upper limit on the diffuse neutrino flux of E 2 Φ=1.7× 10 -6 GeV cm -2 s -1 sr -1 for a decade wide logarithmic bin centered at a neutrino energy of 1018,eV, which is an order of magnitude improvement compared to the previous limit reported by the ARIANNA collaboration. The ARIANNA stations, including purpose built cosmic-ray stations at the Moore's Bay site and demonstrator stations at the South Pole, have operated reliably. Sustained operation at two distinct sites confirms that the flexible and adaptable architecture can be deployed in any deep ice, radio quiet environment. We show that the scientific capabilities, technical innovations, and logistical requirements of ARIANNA are sufficiently well understood to serve as the basis for large area radio-based neutrino telescope with a wide field-of-view.

79 ASTRONOMY AND ASTROPHYSICS↗

Search for narrow resonances in the b -tagged dijet mass spectrum in proton-proton collisions at s = 13 TeV

A search is performed for narrow resonances decaying to final states of two jets, with at least one jet originating from a b quark, in proton-proton collisions at s = 13 TeV . The data set corresponds to an integrated luminosity of 138 fb - 1 collected with the CMS detector at the LHC. Jets originating from energetic b hadrons are identified through a b -tagging algorithm that utilizes a deep neural network or the presence of a muon inside a jet. The invariant mass spectrum of jet pairs is well described by a smooth parametrization and no evidence for the production of new particles is observed. Upper limits on the production cross section are set for excited b quarks and other resonances decaying to dijet final states containing b quarks. These limits exclude at 95% confidence level models of Z ' bosons with masses from 1.8 TeV to 2.4 TeV and of excited b quarks with masses from 1.8 TeV to 4.0 TeV. This is the most stringent exclusion of excited b quarks to date.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Electric Grid Supply Chain Review: Large Power Transformers and High Voltage Direct Current Systems

The report “America’s Strategy to Secure the Supply Chain for a Robust Clean Energy Transition” lays out the challenges and opportunities faced by the United States in the energy supply chain as well as the federal government plans to address these challenges and opportunities. It is accompanied by several issue-specific deep dive assessments, including this one, in response to Executive Order 14017 “America’s Supply Chains,” which directs the Secretary of Energy to submit a report on supply chains for the energy sector industrial base. The Executive Order is helping the federal government to build more secure and diverse U.S. supply chains, including energy supply chains. The need to modernize and increase the capacity of the U.S. power grid is increasing due to growing population, aging infrastructure, grid resilience requirements, operational flexibility needs, and a growing portfolio of renewable energy. As part of the Paris Climate Agreement the United States has committed to reducing its net greenhouse gas emissions by 50 to 52 percent below 2005 levels in 2030 and has set a goal of 100 percent carbon pollution free electricity by 2035 (The United States of America Nationally Determined Contribution, 2021). Meeting these targets will require a significant expansion of the power system to integrate a large amount of new renewable resources (United States Department of State & United States Executive Office of the President, 2021). However, many critical components supporting the power grid have limited to no domestic manufacturing capacity and face complex challenges in supporting a rapid expansion of the grid to meet multiple objectives, including decarbonization goals. This report focuses on two key grid components: large power transformers (LPTs) and high-voltage direct current (HVDC) transmission.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Deep learning-based spatio-temporal estimate of greenhouse gas emissions using satellite data

Accurate estimation of greenhouse gases (GHGs) emissions is very important for developing mitigation strategies to climate change by controlling and reducing GHG emissions. This project aims to develop multiple deep learning approaches to estimate anthropogenic greenhouse gas emissions using multiple types of satellite data. NO2 concentration is chosen as an example of GHGs to evaluate the proposed approach. Two sentinel satellites (sentinel-2 and sentinel-5P) provide multiscale observations of GHGs from 10-60m resolution (sentinel-2) to ~kilometer scale resolution (sentinel-5P). Among multiple deep learning (DL) architectures evaluated, two best DL models demonstrate that key features of spatio-temporal satellite data and additional information (e.g., observation times and/or coordinates of ground stations) can be extracted using convolutional neural networks and feed forward neural networks, respectively. In particular, irregular time series data from different NO 2 observation stations limit the flexibility of long short-term memory architecture, requiring zero-padding to fill in missing data. However, deep neural operator (DNO) architecture can stack time-series data as input, providing the flexibility of input structure without zero-padding. As a result, the DNO outperformed other deep learning architectures to account for time-varying features. Overall, temporal patterns with smooth seasonal variations were predicted very well, while frequent fluctuation patterns were not predicted well. In addition, uncertainty quantification using conformal inference method is performed to account for prediction ranges. Overall, this research will lead to a new groundwork for estimating greenhouse gas concentrations using multiple satellite data to enhance our capability of tracking the cause of climate change and developing mitigation strategies.

54 ENVIRONMENTAL SCIENCES↗

Refining Principal Stress Measurements in Reservoir Underburden in Regions of Induced Seismicity through Seismological Tools, Laboratory Experiments - Final Technical Report

This project developed methodologies to measure the in-situ principal stress in the deep subsurface through use of multiple independent, but complementary, seismic methods, laboratory verification, and development of theoretical frameworks. By leveraging existing regional and local datasets we developed, tested, and refined a set of diagnostic tools for determining the in-situ stress state with reduced uncertainty at and below reservoir depths (1.5-6 km). A set of novel tools was produced that are scale independent, such that their utility is equivalent on regional, field scale, and near borehole monitoring of principal stresses in reservoir underburden for carbon storage projects. During a 4-year Department of Energy (DOE) and Southern Company funded project, carried out by the Electric Power Research Institute (EPRI), Lawrence Livermore National Laboratory (LLNL), the University of Oklahoma (OU), and the U.S. Geological Survey (USGS), the project team developed methodologies to measure the far-field in-situ principal stress in the deep subsurface, leveraging induced seismicity data from waste-water disposal projects. These methodologies consisted in the use of well-established and technically advanced seismic processing methods, such as virtual seismometer method-moment tensor (VSM-MT) and shear wave splitting (SWS), that are adept at recovering the stress orientation and certain components of the stress tensor. These methods were applied to robust seismicity catalogs created with matched filter techniques near sites of active fluid disposal—a proxy for carbon storage sites where such datasets are more limited. Estimates of the stress orientation made with seismic processing tools were considered along with laboratory acoustic emission experiments conducted on rock samples from the region of interest. Stress orientations in the studied region do not vary significantly across distances of ~100 km, nor are they found to rotate through time as a consequence of local wastewater disposal, as previously speculated. Finally, the project team investigated the trade-offs among the different seismic methods and evaluated the range of uncertainty that is generated with these methodologies, which led to a practical use and refinement of the VSM-MT technique when it is applied to field datasets. Understanding the trade-offs between these different methods highlighted the potential benefits of improved quantification of uncertainties on stress field estimations.

58 GEOSCIENCES↗

Techno-Economic Performance of Eavor Loop 2.0

This project evaluated techno-economic performance for a sample Eavor-Loop 2.0 design for electricity production and direct-use heating. The Eavor-Loop 2.0 design investigated is a 7.5-km deep closed-loop geothermal system consisting of 12 laterals for a total of more than 90 km of downhole well and lateral length. Both a high geothermal gradient scenario of 60 degrees C/km and a low geothermal gradient scenario of 30 degrees C/km were considered. With pure water injected at 60 degrees C and 80 kg/s, reservoir simulations with the Slender-Body Theory simulator indicate average production temperatures over a 30-year lifetime of ~125 degrees C and ~210 degrees C for the low and high geothermal gradient scenario, respectively. These correspond to heat production of ~22 M Wth and ~51 M Wth, respectively. Using IPSEpro simulations, we find average power production of ~2.2 M We and ~8.6 M We, respectively, for a subcritical organic Rankine cycle power plant with air-cooled condensers. Cost estimates indicate the overall capital and levelized costs are dominated by the lateral drilling cost. Obtaining a levelized cost of electricity below $70/M Wh requires a geothermal gradient of 60 degrees C/km, a discount rate below 9%, and lateral drilling cost below $400/m. A well cost model indicates that ~$400/m for the Eavor-Loop 2.0 design investigated can be obtained for a drilling rate of penetration about 40 ft/hr (with bit life of 50 hours), and omitting casing and cement. Traditional (geothermal) well drilling has achieved these drilling rate conditions, including the Utah FORGE project where the rate of penetration has exceeded 50 ft/hr in granite. However, it is unclear if these conditions are still valid for drilling the Eavor-Loop 2.0 laterals (i.e., ~82 km of laterals at 4 to 7.5-km vertical depth with rock temperatures up to 460 degrees C), as such downhole completion has never been developed before. Competitive levelized cost of heat values ($1.2-$8.2/GJ) are calculated, even for the low geothermal gradient scenario (30 degrees C/km) and lateral drilling cost of $600/m.

advanced geothermal system↗

TopTemp: Parsing Precipitate Structure from Temper Topology

Technological advances are in part enabled by the development of novel manufacturing processes that give rise to new materials or material property improvements. Development and evaluation of new manufacturing methodologies is labor-, time-, and resource-intensive expensive due to complex, poorly defined relationships between advanced manufacturing process parameters and the resulting microstructures. In this work, we present a topological representation of temper (heat-treatment) dependent material micro-structure, as captured by scanning electron microscopy, called TopTemp. We show that this topological representation is able to support temper classification of microstructures in a data limited setting, generalizes well to previously unseen samples, is robust to image perturbations, and captures domain interpretable features. The presented work outperforms conventional deep learning baselines and is a first step towards improving understanding of process parameters and resulting material properties.

Kassab, Lara↗

Fe Oxidation and Species Distribution at the Rock–Fluid Interface of Marcellus Shale Reacted with Hydraulic Fracturing Fluid

Hydraulic fracturing of shale reservoirs resulted in significant opportunity for increased oil and gas production in the United States. Rock-fluid interactions can cause mineral dissolution and precipitation reactions that lead to permeability changes in the shale matrix, which ultimately may affect transport pathways and hydrocarbon production. Understanding the distribution of secondary precipitates, such as barite and Fe(III) (hydro)oxides, and cation leaching at the rock-fluid interface is an important step to further investigate how these geochemical processes can change permeability and transport pathways. In this study, thin sections of the fracture-matrix interface were made from reacted Marcellus shale cores. The thin sections were characterized using synchrotron X-ray fluorescence imaging and synchrotron X-ray absorption spectroscopy. Fe species with different oxidation states were identified in the maps, together with barite and Ca distribution. The results show that ferrihydrite, as newly formed Fe(III)-bearing precipitates, aligned well with the border of the Ca (e.g., calcite) leaching region in the reaction front. Some Fe-containing clay also dissolved, but the dissolution region for the clay was not as deep as the calcite. Further, the reaction front is about three times deeper in the direction parallel to the shale bedding than that perpendicular to the bedding. The Ca leaching region can be an index for reaction front detection for Marcellus shale. Reactive transport modeling was conducted and the predicted Ca leaching boarder align well with ferrihydrite precipitation, consistent with the experimental observation. The carbonate mineral dissolution can be crucial to promote fluid access into the shale matrix. Together with our previous study on the shale reactive surface, this follow-up study showed similar Ca leaching region and Fe(III) precipitates distribution in the reaction front regardless of barite precipitation on the surface, indicating that the barite coatings on the surface may not pose a significant impact on reactive transport at the shale-fluid interface.

04 OIL SHALES AND TAR SANDS↗

Thermal and Well Flow Performance of Closed-Loop Geothermal in the Wattenberg Area

Closed-loop geothermal systems provide an alternative to resource-constrained hydrothermal systems and stimulation-intensive enhanced geothermal systems. In this work, we apply the slender-body theory (SBT) model, to simulate the well flow and heat transfer performance of U-loop well designs in the Wattenberg area of the Denver-Julesburg Basin. Three U-loop well patterns are investigated, including a single-, double-, and multi-lateral design. The subsurface within the area of interest is characterized by deep, hot (> 200 degrees C) igneous/metamorphic basement rock underlying multiple sedimentary formations. The lateral section(s) of the U-loop lie(s) within a target depth of 6 km, where temperatures are estimated to approach 300 degrees C. As a base case, conduction-only heat transfer is investigated through simulations with the SBT model within U-loops with open-hole laterals that exchange heat directly with the hot, dry rock using water as a working fluid. The utilization of supercritical CO2 as a heat transfer fluid is also considered. For each scenario, the system performance in terms of annual heat production and temperature profile over a 20-year project lifetime are assessed. Also, the levelized costs of heat and electricity (LCOH and LCOE) are determined using a top-down techno-economic analysis model. The results show that the performance- and cost-optimized U-loop design is one having an injection-production well spacing of 1,000 meters with ten 50-meter-spaced laterals that traverse a subsurface system with a temperature gradient of 60 degrees C/km. By injecting 20 degrees C-water at a rate of 60 kg/s through this loop, an average heat production of 19 MWth (i.e., 2.2 MWe net plant output) can be achieved, resulting in an LCOE and LCOH of $136/MWhe and $1.53/GJ, respectively, over a 20-year project life.

closed-loop geothermal↗

Unified Entrainment and Detrainment Closures for Extended Eddy-Diffusivity Mass-Flux Schemes

The extended Eddy Diffusivity Mass Flux (EDMF) scheme is used as a unied parameterization that simulates subgrid-scale dynamics across a range of regimes.This work presents unied entrainment and detrainment closures that allow this scheme to simulate the dry convective boundary layer, shallow convection and deep convection , while adjacent works present other necessary closures for this unied scheme and test its performance for stable, stratocumulus topped and convective boundary layers. Entrainment and detrainment are modeled as a combination of turbulent and dynamical processes. The turbulent component is taken as downgradient diffusion between subdomains. The dynamical component is based on a ratio of buoyancy and vertical velocity that ensures consistent treatment of the ratio of mass-flux at vanishing area fraction. Dynamical detrainment is controlled by the difference in relative humidity between subdomains when liquid water is present. These closures are rst evaluated with entrainment and detrainment diagnosed from Large Eddy Simulations where decaying tracers are used to identify the different subdomains and mass exchanges between them. When implemented in a single column model (SCM), these closures are tested in simulations of several observationally motivated cases of dry convective boundary layer, shallow convection and deep convection. Results show good agreement between the SCM and Large Eddy Simulations in key properties for climate models and profiles of first, second and third moments. The development of diurnal cycles and the onset of and precipitation (with simplified microphysics) are well predicted by the SCM.

54 ENVIRONMENTAL SCIENCES↗

Disentangling the EMC effect: from free to bound nucleon structure

Deep inelastic scattering experiments using electron beams on fixed targets allow for the study of the structure of nucleons. While the structure of the free proton has been well measured, the free neutron structure and the modification of nucleon structure in the nuclear environment, the EMC effect, remain open questions. Studies suggest that the EMC effect may arise due to strongly-interacting nucleon-nucleon systems, referred to as short-range correlations. The BAND experiment at the Thomas Jefferson National Accelerator Facility seeks to probe the bound proton structure when in a strongly-interacting state, using spectator-tagged deep inelastic scattering. This thesis presents the results of the BAND experiment, demonstrating the proton is modified when strongly-interacting with the spectator neutron in deuterium. This thesis also explores efforts to describe inclusive and tagged deep inelastic scattering data within a uniform theoretical framework to explore nucleon structure modification effects and extract the free neutron structure.

Segarra, Efrain↗

Utah FORGE - Development of a Reservoir Seismic Velocity Model and Seismic Resolution Study

This is data from and a final report on the development of a 3D velocity model for the larger FORGE area and on the seismic resolution in the stimulated fracture volume at the bottom of well 16A-32. The velocity model was developed using RMS velocities of the seismic reflection survey and seismic velocity logs from borehole measurements as an input model. To improve the accuracy of the model in the shallow subsurface, travel times phase arrivals of the direct propagating P-waves were determined from the seismic reflection data, using PhaseNet, a deep-neural-network-based seismic arrival time picking method. The travel times were subsequently inverted using the input velocity model. The seismic resolution study used borehole and surface seismic sensors as well as the seismicity observed during the April 2022 stimulation experiment to estimate the seismic resolution in the activated fracture reservoir. The data contain a 3D P- and S-wave velocity model for the larger FORGE area.

15 GEOTHERMAL ENERGY↗