Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Fluid Properties”

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 199 records · Page 11

A constitutive model for sheared dense fiber suspensions

Here we propose a constitutive model to predict the viscosity of fiber suspensions, which undergoes shear thinning, at various volume fractions, aspect ratios, and shear stresses/rates. We calibrate the model using the data from direct numerical simulation and prove the accuracy by predicting experimental measurements from the literature. We use a friction coefficient decreasing with the normal load between the fibers to quantitatively reproduce the experimentally observed shear thinning in fiber suspensions. In this model, the effective normal contact force, which is directly proportional to the bulk shear stress, determines the effective friction coefficient. A rise in the shear stress reduces the effective friction coefficient in the suspension. As a result, the jamming volume fraction increases with the shear stress, resulting in a shear thinning in the suspension viscosity. Moreover, we extend the model to quantify the effects of fiber volume fraction and aspect ratio in the suspension. We calibrate this model using the data from numerical simulations for the rate-controlled shear flow. Once calibrated, we show that the model can be used to predict the relative viscosity for different volume fractions, shear stresses, and aspect ratios. The model predictions are in excellent agreement with the available experimental measurements from the literature. The findings of this study can potentially be used to tune the fiber size and volume fraction for designing the suspension rheology in various applications.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Thermophysical Property Measurements on Salt Mixture

Thermal conductivity, specific heat, viscosity, and density are critical properties for understanding thermal hydraulics of heat transfer fluids (HTFs) like molten salts. These properties are important for design and safe operation of any HTF loop. These characterizations are required to benchmark computational chemical models that are used to predict properties, but existent data are limited. These measurements are extremely difficult to make due to the operating temperature range and corrosive properties of the salt.Density measurements are needed to predict volumetric behavior. Thermal expansion mismatches can cause mechanical failure of structural materials and important components that are in contact with the molten salt. All of these properties are related to heat transfer in HTFs. The heat generated in a reactor must be high enough to derive power; however it should not be so hot that it cannot be safely removed via the HTF at an appropriate rate. It is important to understand the thermal hydraulics required to operate at the desired power generation level and to avoid situations in which thermal shifts can damage components,cause failure of structural materials, and/or lead to accidents. Understanding these thermal properties can aid in component and process design and can contribute to optimized operational and safety parameters.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Electrically controlled two-dimensional electron-hole fluids

We study the electronic properties of dual-gated electron-hole bilayers in which the two layers are separated by a perfectly opaque tunnel barrier. Combining an electrostatic and thermodynamic analysis with mean-field theory estimates of interacting system chemical potentials, we explain the dependence of the electron and hole densities on the two gate voltages. Because chemical potential jumps occur for both electrons and holes at neutrality, there is a finite area in gate voltage parameter space over which electron and hole densities are equal. In that regime the electron-hole pair density depends only on the sum of the two gate voltages. As a result, we are able to explain a recent experimental study of electrically controlled bilayers by allowing for interlayer tunneling and using a nonequilibrium steady-state transport picture.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

The Effect of Hydrothermal Alteration and Microcracks on Hydraulic Properties and Poroelastic Deformation: A Case Study of the Blue Mountain Geothermal Field

Abstract Geothermal energy plays a vital role in decarbonizing electricity and heat supply. Effective utilization of geothermal resources hinges on identifying or generating permeable reservoir zones and understanding how effective pressure variations affect fluid circulation and reservoir properties by poroelastic deformation. Hydrothermal alteration can modify the petrophysical properties of geothermal reservoir rocks, which may increase or decrease its productivity. Understanding these alteration effects is essential to predict and optimize long‐term sustainable geothermal operations. Here, we investigate the impact of hydrothermal alteration on poroelastic and hydraulic properties of diverse lithologies in a series of deformation tests performed at several confining (0–80 MPa) and pore pressure (10–30 MPa) levels. Experimental results of hydrothermally altered dikes and phyllites obtained from the Blue Mountain geothermal field (Nevada, USA) are compared to thermally cracked La Peyratte granite (France) and correlated with petrophysical properties, mineral composition, and microstructures. Argillic alteration of dikes increases porosity and storage capacity but lowers thermal conductivity and increases pore compressibility. Conversely, silicate precipitation in phyllites increases stiffness and thermal conductivity but also reduces porosity and permeability. Experimentally determined effective pressure coefficients range from 0.1 to 0.9, differ for permeability and volumetric strain and decrease with increasing effective pressure. The presence of compliant microcracks and crack‐like pores significantly increases the stress sensitivity of La Peyratte granite and silicified phyllites. This study demonstrates how thermal and chemical alteration impacts poromechanical and petrophysical characteristics of geothermal targets, which ultimately govern reservoir stability and subsidence, induced seismicity as well as fluid and heat extraction efficiency during geothermal operations.

Schuster, Valerian [Helmholtz Centre Potsdam GFZ G↗

Pool boiling heat transfer evaluation of next-generation dielectric fluid: Opteon™ 2P50

The growing use of artificial intelligence has led to heavy thermal loads and high heat dissipation rates in data centers. Conventional air-cooled technologies are not able to fulfill these requirements. To overcome these challenges, two-phase immersion cooling (2PIC) has emerged as one of the leading technologies for high power-density chips. 2PIC increases the heat dissipation rate and efficiency of the system while reducing the footprint of the cooling equipment. A fluid with adequate dielectric properties, a suitable normal boiling temperature to maintain chip temperatures, and good material compatibility, is desired for 2PIC system. In this study, the pool boiling heat transfer of a new developmental dielectric fluid, Opteon™ 2P50, was experimentally investigated. The heat transfer coefficients at various heat fluxes (20–150 kW/m 2 ) and the critical heat flux were measured using a smooth aluminum surface. Compared with HFE-7100, Opteon™ 2P50 shows higher heat transfer coefficient (up to 59% higher) and a slightly lower value of critical heat flux (around 5.9% lower). The modified Cooper correlation with the optimized leading constant resulted in reliable prediction accuracy with a 5.3% mean absolute error percentage. Overall, these results indicate that the new dielectric fluid provides similar thermal performance to some legacy fluids.

2P50↗

A modified many-body dissipative particle dynamics model for mesoscopic fluid simulation: methodology, calibration, and application for hydrocarbon and water

The many-body dissipative particle dynamics (mDPD) is a prominent mesoscopic multiphase model for fluid transport in mesoconfinement. However, it has been a long-standing challenge for mDPD (and other multiphase-enabled DPD models) to accurately predict real-fluid static and dynamic properties simultaneously. We have developed a modified mDPD model to overcome the issue and a rigorous calibration approach that uses reference data, including experimental and/or molecular dynamics (MD) simulations to parameterise the modified mDPD for real fluids. We choose heptane as a representative example of hydrocarbon in source rocks to demonstrate the model's capability to accurately predict the equation of state (EOS), free surface tension, diffusivity, and viscosity. Our timing test shows that the modified mDPD is 400–500 times faster than its MD counterpart for simulating bulk heptane in equivalent volumes. To further demonstrate the robustness of the model, we revisited the benchmark problem of mesoscopic modelling of water, in which all the previous works on DPD reported only a limited portion of the water properties. Here, we show that the modified mDPD can provide accurate modelling of water static and dynamic properties and an EOS that matches the experimental data to a large range of confinement pressure.

42 ENGINEERING↗

Design and performance of a variable gap system for thermal conductivity measurements of high temperature, corrosive, and reactive fluids

We report high-temperature fluids such as molten salts, liquid metals, and gasses are being proposed for many advanced energy systems including thermal energy storage devices, concentrating solar plants, and advanced nuclear reactor designs. However, the chemical behavior and thermophysical properties of many of these fluids have not been well characterized, which hinders the design, modeling, safety analysis, and deployment of these systems. Thermal conductivity is a property that is especially limited by existing measurement capabilities, which are subject to errors caused by convection, material interaction, radiative heat transfer, and instrument degradation. Therefore, there is a lack of standard, systematic measurement techniques for high-temperature, reactive, and corrosive fluids. In this work, the development of a variable gap thermal conductivity measurement system is detailed. The system is designed to measure the thermal conductivity of highly corrosive and reactive fluids, and survive operation between 100 °C and 800 °C. The effects of convection are minimized by limiting the thickness of the specimen to thin sizes (<0.3 mm). Corrections for radiative heat transfer were included in the working equations to consider specimens with varying optical properties. The design, construction, instrumentation, operating principles, and data analysis techniques are discussed in detail. The system was tested up to 500 °C using helium gas and molten KNO 3 -NaNO 3 to verify the measurement technique and determine the sources of error. At 300 and 400 °C KNO 3 -NaNO 3 , results showed maximum relative error of 6% when compared to results in the literature. The helium results were within 13% of those in the literature at 300 and 400 °C. Higher errors were observed at 500 °C for both fluids, and the sources of these errors are discussed.

42 ENGINEERING↗

Exploration of Heat-Driven Ejector High-Temperature Heat Pumps

Heat-driven ejector heat pumps use a supersonic ejector as a thermo-compressor to replace the mechanical compressor. Supersonic ejectors have many advantages in high-temperature heat pump (HTHP) applications, including high operating temperature tolerance, no need for lubrication, low maintenance, and low cost. The coefficient of performance of ejector HTHPs could be improved by selecting binary fluids with unique thermodynamic properties. Although supersonic ejectors have been widely used in refrigeration systems, their application in spaces and water heating is limited. This study explores the theoretical potential of ejector HTHPs with a sink temperature of 100°C–130°C and a lift temperature of 10°C–30°C. Ejector HTHPs were evaluated with single-fluid ejectors (SFEs) and binary-fluid ejectors (BFEs). A comprehensive, geometry-free theoretical model of BFEs was developed to predict the theoretical maximum entrainment ratios. HFE7500 and R718 (water) were selected as working fluids for SFEs and BFEs. SFEs operating with R718 provided a higher coefficient of performance of ejector HTHPs than SFEs operating with HFE7500 and BFEs. This study preliminarily demonstrates the technical potential of ejector HTHP applications in recovering moderate-temperature heat sources

Wang, Pengtao↗

Linear analysis of the Kelvin–Helmholtz instability in relativistic magnetized symmetric flows

ABSTRACT We study the linear stability of a planar interface separating two fluids in relative motion, focusing on the symmetric configuration where the two fluids have the same properties (density, temperature, magnetic field strength, and direction). We consider the most general case with arbitrary sound speed cs, Alfvén speed vA, and magnetic field orientation. For the instability associated with the fast mode, we find that the lower bound of unstable shear velocities is set by the requirement that the projection of the velocity on to the fluid-frame wavevector is larger than the projection of the Alfvén speed on to the same direction, i.e. shear should overcome the effect of magnetic tension. In the frame where the two fluids move in opposite directions with equal speed v, the upper bound of unstable velocities corresponds to an effective relativistic Mach number $M_{\rm re}\equiv v/v_{\rm {f}\perp }\sqrt{(1-v_{\rm {f}\perp }^2)/(1-v^2)} \cos \theta =\sqrt{2}$, where $v_{\rm {f}\perp }=[v_{\rm {A}}^2+c_{\rm s}^2(1-v_{\rm {A}}^2)]^{1/2}$ is the fast speed assuming a magnetic field perpendicular to the wavevector (here, all velocities are in units of the speed of light), and θ is the laboratory-frame angle between the flow velocity and the wavevector projection on to the shear interface. Our results have implications for shear flows in the magnetospheres of neutron stars and black holes – both for single objects and for merging binaries – where the Alfvén speed may approach the speed of light.

Astronomy & Astrophysics↗

Collective excitations of a bound-in-the-continuum condensate

Spectra of low-lying elementary excitations are critical to characterize properties of bosonic quantum fluids. Usually these spectra are difficult to observe, due to low occupation of non-condensate states compared to the ground state. Recently, low-threshold Bose-Einstein condensation was realised in a symmetry-protected bound state in the continuum, at a saddle point, thanks to coupling of this electromagnetic resonance to semiconductor excitons. While it has opened the door to long-living polariton condensates, their intrinsic collective properties are still unexplored. Here we unveil the peculiar features of the Bogoliubov spectrum of excitations in this system. Thanks to the dark nature of the bound-in-the-continuum state, collective excitations lying directly above the condensate become observable in enhanced detail. We reveal interesting aspects, such as energy-flat parts of the dispersion characterized by two parallel stripes in photoluminescence pattern, pronounced linearization at non-zero momenta in one of the directions, and a strongly anisotropic velocity of sound.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Joint impedance and facies inversion of time-lapse seismic data for improving monitoring of CO 2 incidentally stored from CO 2 EOR

Time-lapse seismic monitoring is an effective and proven technology for mapping the distribution of CO 2 in a subsurface reservoir. When injected CO 2 displaces other reservoir fluids, porous-medium properties are changed and thus the seismic impedance changes, causing time-lapse seismic amplitude differences in the injection zones. The analysis and interpretation of images created from these amplitude differences can provide information about reservoir architecture and the CO 2 migration within the reservoir. Incorporating seismic inversion and rock physics into the interpretation of time-lapse seismic data can considerably improve the modeling and monitoring to detect and assess the location of CO 2 over time. The joint inversion method presented in this paper has an integral representation of the geology in the inversion algorithm using elastic facies, which provides information about the spatial distribution of the geologic heterogeneities controlling the movement of fluids in the reservoir. The method was successfully applied to time-lapse seismic data from a mature oil field undergoing CO 2 enhanced oil recovery. The estimated seismic acoustic impedances and facies reflect the characteristics of individual geologic facies and fluid conditions of the reservoir subject to CO 2 injection. Finally, the probabilities estimated by the joint impedance and facies inversion for the reservoir's litho-fluid facies can be used for forecasting CO 2 saturation and pressure changes within the target reservoir.

4D seismic↗

Finite Element Analysis (FEA) for Water-Foam Fracturing of Granite Rock

In addition to the foam data that were obtained from literature and that were collected from the current study, simulation data was also generated from finite element analysis (FEA) conducted in this study using COMSOL Multiphysics software. The FEA models were built to simulate the experiments conducted at Oak Ridge National Laboratory (ORNL) on cement and granite samples. In these FEA models, temperature was kept at ambient while the pressure profile resembled the loading conditions during the ORNL experiments, where pressure was either monotonically increased or applied cyclically. The cement material was used as a model material and was used to study Von Mises stress and tensile stress distribution for different bore hole length geometry using a parametric sweep with water as fracturing fluid using solid-fluid interaction module. For the granite material, FEA models were developed for stress analysis of cylindrical samples with water or foam fluids. The solid mechanics module in COMSOL was implemented to solve for Von Mises stress and tensile stress. The fluid-structure interaction module was implemented to solve for water-foam interaction on granite cylinder with addition of fluid-loading on structure, i.e., large deformation in solid mechanics with no impact on fluid deformation. Foam was considered as a pseudo single-phase compressible fluid for which material properties were calculated from water and gas (nitrogen) phases. The density of foam is calculated as a function of the densities of water and nitrogen, while viscosity is a function of temperature. Four types of FEA analyses were modelled: 1. Monotonic injection with water 2. Monotonic injection with foam 3. Cyclic injection with water 4. Cyclic injection with foam All the COMSOL files are converted to a zip file which is save in .mph.

15 GEOTHERMAL ENERGY↗

Online Alpha Monitoring of High Cs-137 Hanford and SRS High Level Waste with Tensioned Metastable Fluid Detectors

The Department of Energy’s Hanford and Savannah River Sites maintain millions of gallons of caustic supernate and salt high activity waste in their high-level waste (HLW) tank farm inventories. The Savannah River Site is currently treating this waste with a calixarene-based solvent extraction of Cs-137 to reduce these inventories. Hanford is employing an at-tank crystalline silico titanate (CST) based solid phase extraction methodology to reduce their liquid HLW inventories. Due to the high solubility of Cs-137 and the relative insolubility of the actinides in these caustic waste forms, the beta to alpha radioactivity ratio can often exceed six orders of magnitude in the feed solutions to these treatment processes. This unique characteristic leads to significant technical challenges in making rapid gross alpha measurements in the presence of the overwhelming beta, gamma, as well as dissolved sodium salt in these HLW matrices. Conventional radioanalytical techniques, such as liquid scintillation analysis or gas flow proportional counting require significant radiochemistry preparation prior to the radiometric measurements for gross alpha activity. These required pretreatments render these technologies untenable for rapid quantification of gross alpha activity that could be required to support on or at-line measurements ensuring a waste stream will meet regulatory requirements. The radiation measurement properties of Tensioned Metastable Fluid Detectors (TMFDs) have been studied by Purdue University’s Taleyarkhan research group for well over a decade. Fluids tensioned to the appropriate degree will rupture when struck by radiation, resulting in a measurable cavitation event. The negative pressure generating this tension can be adjusted by centrifugal rotation or by acoustic means in such a way that these cavitation events can be generated from alpha radiation but will not be generated by beta or gamma radiation. Purdue University and the Savannah River National Laboratory are currently collaborating to develop a gamma/beta blind, spectroscopic alpha measurement system based on the Tensioned Metastable Fluid Detector technology to provide a potential solution for performing rapid gross alpha measurements on these high gamma/beta sample matrices. Measurements using the Indirect Drive Acoustically Tensioned Metastable Fluid Detectors developed as part of this collaboration were performed with an alpha emitting radionuclide. Successful determination of gross alpha activity was observed, indicating a potential pathway for rapid gross alpha measurements in remote-handled shielded cells or in process situations requiring online alpha monitoring. Measurements using this system have been conducted on high beta activity solutions, demonstrating the beta blind capability of this system. Measurements are currently underway to test the system’s capability to measure gross alpha activity on Savannah River Site high level waste high Cs-137 samples that have been measured by the SRNL radiochemistry team. This work was supported by the DOE EM Technology Development program.

DiPrete, David [Savannah River National Laboratory↗

Temperature-Dependent Speciation of Ni(II) in Molten Chloride, Bromide, and Iodide Salts

Understanding the fundamental speciation and local structure of metal cations─such as nickel, which is present as a corrosion product─in molten salt media is essential for assessing their impact on the thermal and physical properties of salts used for advanced nuclear energy applications. In this study, we employed an integrated approach of combining X-ray absorption spectroscopy (XAS), UV–vis spectroscopy, and ab initio molecular dynamics (AIMD) to investigate the halide-dependent speciation of divalent nickel ions, Ni(II), in LiCl–KCl, LiBr–KBr, and LiI–KI eutectic salt mixtures from room temperature to 600 °C. Our findings show that both halide type and temperature significantly influence the local coordination of nickel in these molten salt systems. A notable aspect of Ni(II) speciation in LiI–KI salts is the formation of a polyiodide species. In conclusion, our findings provide critical insights into metal ion speciation in molten salt systems, enhancing our understanding of the factors that govern the physicochemical properties of these complex fluids.

36 MATERIALS SCIENCE↗

Rheological Properties of Enzymatically Hydrolyzed Corn Stover Pretreated via Deacetylation and Mechanical Refining

Lignocellulosic biomass is a feedstock for fuels and chemicals that does not compete with food resources and has less contaminants than refuse-derived biomass feedstocks. To convert lignocellulosics to biofuels or value-added products, multiple processing steps are typically necessary. One method of producing biofuels from lignocellulosic biomass utilizes a deacetylation and mechanical refining pretreatment and an enzymatic hydrolysis reaction to produce fermentable sugars from cellulose and hemicellulose. The rheological properties of biomass, such as yield stress and plastic viscosity, change during enzymatic hydrolysis and alter the energy requirements for pumping and mixing, an important consideration for the design of processing equipment. The dynamic changes in rheological properties that occur in a corn stover feedstock undergoing enzymatic hydrolysis are characterized in this work, and the influence on pressure losses in piping systems is estimated. Two rheometer geometries were fabricated with stereolithography 3D printing to reduce wall slip and sample ejection. The slurries have complex rheological behaviors that include shear-thinning behavior. Shear stress ramps were performed on samples at 20 and 50 degrees C using the custom geometries, and the Herschel-Bulkley model was fit to the data. The dynamic nature of the rheological properties is correlated with changes in the average fiber length at various extents of reaction, and the influence of solids concentration on the observed rheology and piping pressure losses is discussed.

09 BIOMASS FUELS↗

Perfectly Matched Layers and Characteristic Boundaries in Lattice Boltzmann: Accuracy vs Cost

Artificial boundary conditions (BCs) play a ubiquitous role in numerical simulations of transport phenomena in several diverse fields, such as fluid dynamics, electromagnetism, acoustics, geophysics, and many more. They are essential for accurately capturing the behavior of physical systems whenever the simulation domain is truncated for computational efficiency purposes. Ideally, an artificial BC would allow relevant information to enter or leave the computational domain without introducing artifacts or unphysical effects. Boundary conditions designed to control spurious wave reflections are referred to as nonreflective boundary conditions (NRBCs). Another approach is given by the perfectly matched layers (PMLs), in which the computational domain is extended with multiple dampening layers, where outgoing waves are absorbed exponentially in time. Here, in this work, the definition of PML is revised in the context of the lattice Boltzmann method. The impact of adopting different types of BCs at the edge of the dampening zone is evaluated and compared, in terms of both accuracy and computational costs. It is shown that for sufficiently large buffer zones, PMLs allow stable and accurate simulations even when using a simple zeroth-order extrapolation BC. Moreover, employing PMLs in combination with NRBCs potentially offers significant gains in accuracy at a modest computational overhead, provided the parameters of the BC are properly tuned to match the properties of the underlying fluid flow.

97 MATHEMATICS AND COMPUTING↗

Synthesis and anisotropic magnetism in quantum spin liquid candidates A YbSe 2 ( A = K and Rb)

The quantum spin liquid (QSL) state in rare-earth triangular lattices has attracted much attention recently due to its potential application in quantum computing and communication. Here, we report the single-crystal growth synthesis, crystal structure characterizations, and magnetic properties of AYbSe 2 (A = K and Rb) compounds. The x-ray diffraction analysis shows that AYbSe 2 (A = K and Rb) crystallizes in a trigonal space group, R-3m (No. 166) with Z = 3. AYbSe 2 possesses a two-dimensional (2D) Yb–Se–Yb layered structure formed by edged-shared YbSe 6 octahedra. The magnetic properties are highly anisotropic for both title compounds, and no long-range order is found down to 0.4 K, revealing the possible QSL ground state in these compounds. The isothermal magnetization exhibits a one-third magnetization plateau when the magnetic fields are applied in the ab-plane. Heat capacity is performed along both ab-plane and c axis and features the characteristic dome for triangular magnetic lattice compounds as a function of magnetic fields. Due to the change in the interlayer and intralayer distance of Yb 3+ , the dome shifts to low fields from KYbSe 2 to RbYbSe 2 . All these results indicate that the AYbSe 2 family presents unique frustrated magnetism close to the possible QSL and noncollinear spin states.

36 MATERIALS SCIENCE↗

Rheology characterization of ionic liquids under high pressure and high temperature

Ionic liquids (ILs) are liquid salts that exist at or below ambient temperatures and are composed of ion pairs. They offer promising alternatives to toxic, hazardous, highly flammable, and volatile solvents in various applications such as solution preparation, dispersion, gel formation, composites, and polymer melts. ILs possess unique and interesting characteristics, including excellent chemical and thermal stability and low vapor pressures. Understanding the rheological properties of ILs is essential to optimizing IL performance. This paper presents a comparative analysis of the rheological properties of two ionic liquids, Nsingle bondHexylpyridinium tetrafluoroborate (HPyBF4) and Nsingle bondHexylpyridinium bromide (HPyBr), under different shear rates, temperatures, and pressures. Rheological measurements were performed under varying controlled pressure and temperature conditions. The experimental investigation covered a pressure range of 689–12,411 kPa [100–1800 psi] and a temperature range from room temperature up to 522 kelvin (K) [480°F]. The primary objective is to explore and compare the flow behavior and viscoelastic characteristics of HPyBF4 and HPyBr under high-pressure and high-temperature conditions. The experimental data showed that HPyBF4 and HPyBr exhibited shear-thinning behavior, and pressure had an insignificant effect on rheology compared to the temperature effect. Under the same testing conditions, HPyBr showed higher shear stress and viscosity than HPyBF4. This research significantly contributes to the improved understanding of the rheological behavior of these specific ionic liquids and their suitability for diverse industrial and scientific applications, particularly in high-pressure and high-temperature environments.

15 GEOTHERMAL ENERGY↗