Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “direct drive”

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 1,351 records · Page 75
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

On-the-fly response function generation method for composite coarse mesh

The hybrid stochastic deterministic transport code COMET, based on the incident response expansion theory, is used to model reactor cores with high fidelity and formidable computational speed. COMET models a reactor core using a library of incident flux response expansion coefficients that are pre computed for all the unique lattice cells (e.g., fuel assemblies, reflector blocks, etc.) in the core. In order to further improve its computational efficiency in pre-calculating the response library a new response function generation method is developed to compute the response functions for the composite coarse meshes made of a smaller set of unique lattices on the fly within the COMET's deterministic transport core sweep. The efficiency is achieved by eliminating a number of unique lattices that can be made up from the reduced set of unique meshes on the fly. The numerical process consists of the following steps. First, the boundary condition on composite coarse mesh boundaries is projected onto the expansion basis to compute the incident flux moments on external surfaces of all the basic (reduced set of unique) coarse meshes. Secondly, the deterministic sweeping solver in COMET is used to converge on the outgoing/incoming flux expansion moments crossing interfaces between the basic coarse meshes. Thirdly, the response functions for the composite coarse meshes are constructed as a superposition on the fly. The new response function generation method was tested on 88 composite coarse meshes consisting of CANDU fuel bundles and moderator blocks. It was found that response functions generated by the new method agree very well with those generated by direct Monte Carlo calculations. The average and maximum relative differences in the surface-to-surface response coefficients computed by the two methods are 0.10% and 0.20%, respectively. Similarly, the average and maximum relative differences in the response fission densities are 0.13% and 0.43%, respectively. These discrepancies are within one standard deviation of the stochastic uncertainties. The new method is five times faster than the original direct Monte Carlo method. The size of the response function library for the new method is five times smaller than that for the original method, leading to significantly less requirement for the computer hard drive space and memory. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2020

Empirical Validation of Multi-Zone HVAC System Model: Evaluation of Existing Infiltration Models used in Building Energy Simulation

Infiltration can have a significant impact on building loads. Studies have shown that infiltration can account for 15-40% of annual space conditioning needs in commercial buildings (Emmerich et al. 2019; Younes et al. 2012). The driving force of infiltration is the pressure difference across the building envelope caused by wind, the stack effect (known as buoyancy effect), and the operation of ventilation equipment. Wind pressure is governed by wind direction, speed, building shape, and other structures around the building. The stack effect is a function of the building height and air density differences of ambient air (Han, 2015). The effect of wind is dominant in low-rise residential buildings, and the stack effect is dominant in highrise buildings (ASHRAE 2017). In building energy simulation programs (e.g., EnergyPlus), various empirical infiltration models (e.g., the effective leakage area model, the flow coefficient model) are available to simulate infiltration rates. To help users in selecting a proper infiltration model for modeling of the two-story Flexible Research Platform (FRP), the team evaluates the existing infiltration models in EnergyPlus based on field measurements from the FRP. The blower door and tracer gas decay tests were performed in the FRP. The blower door test result was used to estimate input parameters required in the infiltration models. The actual infiltration rates were estimated with the tracer gas decay test results.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2025

Gain of Object/Slot Array On Infinite Cylinder Operating At Fundamental Cavity Modes

This report examines the coupling into a circumferential slot array on an infinite cylinder with a finite length interior cavity operating at the fundamental cavity modal frequencies. The directivity gain associated with the interior cavity field is examined, and in this frequency range (as well as lower frequencies), is found to be enhanced by the exterior short circuit current variation with angle driving the slots.

42 ENGINEERING↗
NASA NTRS1973

Analysis of the dynamics and frequency spectrum synthesis of an optical-mechanical scanning device

A two-coordinate optical-mechanical scanning device (OMSD), the operating unit of which is a scanning disk, with directional and focusing optics and a board, on which the data carrier is placed, is examined. The disk and board are kinematically connected by a transmission mechanism, consisting of a worm and complex gear drive and a tightening screw-nut with correcting device, and it is run by a synchronous type motor. The dynamic errors in the system depend, first, on irregularities in rotation of the disk, fluctuations in its axis and vibrations of the table in the plane parallel to the plane of the disk. The basic sources of the fluctuations referred to above are residual disbalance of the rotor and other rotating masses, the periodic component of the driving torque of the synchronous motor, variability in the resistance, kinematic errors in the drive and other things. The fluctuations can be transmitted to the operating units through the kinematic link as a flexural-torsional system, as well as through vibrations of the housing of the device.

Andryushkevichyus, A. I.↗
NASA NTRS1982

In-situ measurements of wave electric fields in the equatorial electrojet

Electric field wave measurements have been performed on two sounding rockets in the equatorial ionosphere. During a daytime flight from Chilca, Peru, intense electrostatic waves were detected on the upward-directed electron density gradient. During a nighttime flight from Kwajalein Atoll, similar waves were detected on a downward directed gradient. These results are in agreement with a gradient drift instability explanation of the generation of the waves. The wave amplitudes were as high as 5 mV/m, implying perturbation drifts comparable to the driving drift velocities. Power spectra from the turbulent region show a peak at long wavelengths, followed by a nearly flat spectral region before breaking into a power law form with negative index of 3.6-3.7 for wavelengths not greater than 30 m. Similarities between the spectra of the two flights suggest that the fundamental processes of the instabilities are the same in the day and nighttime conditions. The rocket data are consistent with radar results presented in a companion paper which show coherent, kilometer scale waves present in the electrojet.

Pfaff, R. F.↗
NASA NTRS1985

Space shuttle rudder/speedbrake actuation subsystem

The Rudder/Speedbrake (R/SB) Actuation Subsystem for use on the NASA Space Shuttle Orbiter is an electro-hydro-mechanical system which provides the control and positionary capability of the orbiter aero-dynamic primary flight control surface. The system is located in the vehicle's vertical stabilizer. The geared rotary actuators provide a power hinge feature of the split panel rudder. Actuation of both panels in the same direction provides conventional rudder control; actuating the panels differentially provides a speedbrake function intended to control both speed and pitch. The commands may be superimposed on one another. The system consists of one power drive unit which responds to quadredundant avionic signals to generate a rotary output, four geared rotary actuators, which develop rotary position and torque as outputs, and ten torque transmitting drive-shifts.

Naber, R. A.↗
NASA NTRS2012

High-Speed Transport of Fluid Drops and Solid Particles via Surface Acoustic Waves

A compact sampling tool mechanism that can operate at various temperatures, and transport and sieve particle sizes of powdered cuttings and soil grains with no moving parts, has been created using traveling surface acoustic waves (SAWs) that are emitted by an inter-digital transducer (IDT). The generated waves are driven at about 10 MHz, and it causes powder to move towards the IDT at high speed with different speeds for different sizes of particles, which enables these particles to be sieved. This design is based on the use of SAWs and their propelling effect on powder particles and fluids along the path of the waves. Generally, SAWs are elastic waves propagating in a shallow layer of about one wavelength beneath the surface of a solid substrate. To generate SAWs, a piezoelectric plate is used that is made of LiNbO3 crystal cut along the x-axis with rotation of 127.8 along the y-axis. On this plate are printed pairs of fingerlike electrodes in the form of a grating that are activated by subjecting the gap between the electrodes to electric field. This configuration of a surface wave transmitter is called IDT. The IDT that was used consists of 20 pairs of fingers with 0.4-mm spacing, a total length of 12.5 mm. The surface wave is produced by the nature of piezoelectric material to contract or expand when subjected to an electric field. Driving the IDT to generate wave at high amplitudes provides an actuation mechanism where the surface particles move elliptically, pulling powder particles on the surface toward the wavesource and pushing liquids in the opposite direction. This behavior allows the innovation to separate large particles and fluids that are mixed. Fluids are removed at speed (7.5 to 15 cm/s), enabling this innovation of acting as a bladeless wiper for raindrops. For the windshield design, the electrodes could be made transparent so that they do not disturb the driver or pilot. Multiple IDTs can be synchronized to transport water or powder over larger distances. To demonstrate the transporting action, a video camera was used to record the movement. The speed of particles was measured from the video images.

Bar-Cohen, Yoseph↗
NASA NTRS1997

Steady-State Pursuit Is Driven by Object Motion Rather Than the Vector Average of Local Motions

We have previously shown that humans can pursue the motion of objects whose trajectories can be recovered only by spatio-temporal integration of local motion signals. We now explore the integration rule used to derive the target-motion signal driving pursuit. We measured the pursuit response of 4 observers (2 naive) to the motion of a line-figure diamond viewed through two vertical bar apertures (0.2 cd/square m). The comers were always occluded so that only four line segments (93 cd/square m) were visible behind the occluding foreground (38 cd/square m). The diamond was flattened (40 & 140 degree vertex angles) such that vector averaging of the local normal motions and vertical integration (e.g. IOC) yield very I or different predictions, analogous to using a Type II plaid. The diamond moved along Lissajous-figure trajectories (Ax = Ay = 2 degrees; TFx = 0.8 Hz; TFy = 0.4 Hz). We presented only 1.25 cycles and used 6 different randomly interleaved initial relative phases to minimize the role of predictive strategies. Observers were instructed to track the diamond and reported that its motion was always coherent (unlike type II plaids). Saccade-free portions of the horizontal and vertical eye-position traces sampled at 240 Hz were fit by separate sinusoids. Pursuit gain with respect to the diamond averaged 0.7 across subjects and directions. The ratio of the mean vertical to horizontal amplitude of the pursuit response was 1.7 +/- 0.7 averaged across subjects (1SD). This is close to the prediction of 1.0 from vertical motion-integration rules, but far from 7.7 predicted by vector averaging and infinity predicted by segment- or terminator-tracking strategies. Because there is no retinal motion which directly corresponds to the diamond's motion, steady-state pursuit of our "virtual" diamond is not closed-loop in the traditional sense. Thus, accurate pursuit is unlikely to result simply from local retinal negative feedback. We conclude that the signal driving steady-state pursuit is not the vector average of local motion signals, but rather a more vertical estimate of object motion, derived in extrastriate cortical areas beyond V1, perhaps NIT or MST.

Stone, Leland S.↗
DOE OSTI2024

Fermi Surface Nesting Driving the RKKY Interaction in the Centrosymmetric Skyrmion Magnet Gd 2 PdSi 3

The magnetic skyrmions generated in a centrosymmetric crystal were recently first discovered in Gd 2 PdSi 3 . In light of this, we observe the electronic structure by angle-resolved photoemission spectroscopy (ARPES) and unveil its direct relationship with the magnetism in this compound. The Fermi surface and band dispersions are demonstrated to have a good agreement with the density functional theory (DFT) calculations carried out with careful consideration of the crystal superstructure. Most importantly, we find that the three-dimensional Fermi surface has extended nesting which matches well the q vector of the magnetic order detected by recent scattering measurements. Finally, the consistency we find among ARPES, DFT, and the scattering measurements suggests the Ruderman-Kittel-Kasuya-Yosida (RKKY) interaction involving itinerant electrons to be the formation mechanism of skyrmions in Gd 2 PdSi 3 .

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

A Zero-Emission Process for Direct Reduction of Iron by Hydrogen Plasma in a Rotary Kiln Reactor

This project’s goal was to demonstrate a hydrogen plasma (H-plasma)-rotary kiln process for reducing iron ore to iron as part of the steel manufacturing process. The H-plasma provides a greater thermodynamic driving force for reducing iron ores than thermal processes such as the DRI process, enabling lower reaction temperatures. We estimated that our process technology can reduce energy consumption by 45% compared to the blast furnace process and ~15% compared to the DRI process. Steel manufacturing produces about 1.8 tons of CO2/ton of steel with iron ore reduction accounting for about one-third of the CO2 produced in the overall manufacturing process. We estimated our process has the potential to reduce GHG emissions from ironmaking by 35% with today’s grid and by up to 88% with a future low-carbon grid while being cost competitive with the current blast furnace route. We demonstrated reduction of hematite and magnetite rich materials at temperatures from 600 to 800°C. We achieved 90-95% metallization on 100 gr samples in batch reduction experiments in the H-plasma rotary kiln furnace at 600-650°C. Attempts to perform tests in a continuous operation mode identified problems with the ore feed mechanism. We identified solutions but there was not time nor budget to correct these for this project

36 MATERIALS SCIENCE↗
DOE OSTI2025

A Zero-Emission Process for Direct Reduction of Iron by Hydrogen Plasma in a Rotary Kiln Reactor

This project’s goal was to demonstrate a hydrogen plasma (H-plasma)-rotary kiln process for reducing iron ore to iron as part of the steel manufacturing process. The H-plasma provides a greater thermodynamic driving force for reducing iron ores than thermal processes such as the DRI process, enabling lower reaction temperatures. We estimated that our process technology can reduce energy consumption by 45% compared to the blast furnace process and ~15% compared to the DRI process. Steel manufacturing produces about 1.8 tons of CO2/ton of steel with iron ore reduction accounting for about one-third of the CO2 produced in the overall manufacturing process. We estimated our process has the potential to reduce GHG emissions from ironmaking by 35% with today’s grid and by up to 88% with a future low-carbon grid while being cost competitive with the current blast furnace route. We demonstrated reduction of hematite and magnetite rich materials at temperatures from 600 to 800°C. We achieved 90-95% metallization on 100 gr samples in batch reduction experiments in the H-plasma rotary kiln furnace at 600-650°C. Attempts to perform tests in a continuous operation mode identified problems with the ore feed mechanism. We identified solutions but there was not time nor budget to correct these for this project

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2021

Crystalline materials for quantum computing: Semiconductor heterostructures and topological insulators exemplars

Abstract High-purity crystalline solid-state materials play an essential role in various technologies for quantum information processing, from qubits based on spins to topological states. New and improved crystalline materials emerge each year and continue to drive new results in experimental quantum science. This article summarizes the opportunities for a selected class of crystalline materials for qubit technologies based on spins and topological states and the challenges associated with their fabrication. We start by describing semiconductor heterostructures for spin qubits in gate-defined quantum dots and benchmark GaAs, Si, and Ge, the three platforms that demonstrated two-qubit logic. We then examine novel topologically nontrivial materials and structures that might be incorporated into superconducting devices to create topological qubits. We review topological insulator thin films and move onto topological crystalline materials, such as PbSnTe, and its integration with Josephson junctions. We discuss advances in novel and specialized fabrication and characterization techniques to enable these. We conclude by identifying the most promising directions where advances in these material systems will enable progress in qubit technology.

Scappucci, G. (ORCID:0000000325120079)↗
DOE OSTI2023

The BiDFET Device and Its Impact on Converters

The matrix converter topology for direct ac-to-ac conversion offers elimination of the bulky and unreliable D.C. link capacitors used in the popular Voltage-Source Inverter (VSI) with a front-end rectifier. The resulting more compact and higher efficiency implementation is a desirable solution for a wide variety of applications, such as photovoltaic energy generation, motor drives, and energy storage systems. The development of matrix converters has been hampered by the lack of availability of commercial bidirectional power switches with the ability to block high voltages in the first and third quadrants, carry on state current in both quadrants with low voltage drop, exhibit large forward-biased safe-operating-area (FBSOA), and low switching power losses. Consequently, many implementations have been tried in the past using discrete devices, as listed in Table I. Two implementations utilize commonly available asymmetric blocking IGBTs, while two cases make use of SiC power MOSFETs. They generally have a large parts count (4-6 individual switches) that occupies significant space in the converter where multiple bi-directional switches are required. They also have a high on-state voltage drop that degrades efficiency. As a result, one implementation utilizes symmetric blocking IGBT to achieve a low parts count (2), but its switching losses are unacceptably high.

42 ENGINEERING↗
NASA NTRS1999

Microgravity Transport Phenomena Experiment (MTPE) Overview

The Microgravity Transport Phenomena Experiment (MTPE) is a fluids experiment supported by the Fundamentals in Biotechnology program in association with the Human Exploration and Development of Space (BEDS) initiative. The MTP Experiment will investigate fluid transport phenomena both in ground based experiments and in the microgravity environment. Many fluid transport processes are affected by gravity. Osmotic flux kinetics in planar membrane systems have been shown to be influenced by gravimetric orientation, either through convective mixing caused by unstably stratified fluid layers, or through a stable fluid boundary layer structure that forms in association with the membrane. Coupled transport phenomena also show gravity related effects. Coefficients associated with coupled transport processes are defined in terms of a steady state condition. Buoyancy (gravity) driven convection interferes with the attainment of steady state, and the measurement of coupled processes. The MTP Experiment measures the kinetics of molecular migration that occurs in fluids, in response to the application of various driving potentials. Three separate driving potentials may be applied to the MTP Experiment fluids, either singly or in combination. The driving potentials include chemical potential, thermal potential, and electrical potential. Two separate fluid arrangements are used to study membrane mediated and bulk fluid transport phenomena. Transport processes of interest in membrane mediated systems include diffusion, osmosis, and streaming potential. Bulk fluid processes of interest include coupled phenomena such as the Soret Effect, Dufour Effect, Donnan Effect, and thermal diffusion potential. MTP Experiments are performed in the Microgravity Transport Apparatus (MTA), an instrument that has been developed specifically for precision measurement of transport processes. Experiment fluids are contained within the MTA fluid cells, designed to create a one dimensional flow geometry of constant cross sectional area, and to facilitate fluid filling and draining operations in microgravity. The fluid cells may be used singly for bulk solutions, or in a Stokes diaphragm configuration to investigate membrane mediated phenomena. Thermal and electrical driving potentials are applied to the experiment fluids through boundary plates located at the ends of the fluid cells. In the ground based instrument, two constant temperature baths circulate through reservoirs adjacent to the boundary plates, and establish the thermal environment within the fluid cells. The boundary plates also serve as electrodes for measurement and application of electrical potentials. The Fluid Manipulation System associated with the MTA is a computer controlled system that enables storage and transfer of experiment fluids during on orbit operations. The system is used to automatically initiate experiments and manipulate fluids by orchestrating pump and valve operations through scripted sequences. Unique technologies are incorporated in the MTA for measurement of fluid properties. Volumetric Flow Sensors have been developed for precision measurement of total fluid volume contained within the fluid cells over time. This data is most useful for measuring the kinetics of osmosis, where fluid is transported from one fluid cell to another through a semipermeable membrane. The MicroSensor Array has been designed to perform in situ measurement of several important fluid parameters, providing simultaneous measurement of solution composition at multiple locations within the experiment fluids. Micromachined sensors and interface electronics have been developed to measure temperature, electrical conductivity, pH, cation activity, and anion activity. The Profile Refractometer uses a laser optical system to directly image the fluid Index of Refraction profile that exists along the MTA fluid cell axis. A video system acquires images of the RI profile over time, and records the transport kinetics that occur upon application of chemical, thermal, or electrical driving potentials. Image processing algorithms have been developed to analyze the refractometer images on a pixel by pixel basis, calibrating and scaling the measured Index of Refraction profile to correlated solution properties of interest such as density, concentration, and temperature. Additional software has been developed to compile the processed images into a three dimensional matrix that contains fluid composition data as a function of experiment time and position in the fluid cell. These data are combined with data from the other sensor systems, and analyzed in the context of transport coefficients associated with the various transport phenomena. Analysis protocols have been developed to measure the transient kinetics, and steady state distribution of fluid components that occur in response to the applied driving potentials. The results are expressed in terms of effective transport coefficients. Experiments have been performed using a variety of solutes, and results generated are that are in agreement with published transport coefficient values.

Mason, Larry W.↗
DOE OSTI2021

Visually resolving the direct Z-scheme heterojunction in CdS@ZnIn 2 S 4 hollow cubes for photocatalytic evolution of H 2 and H 2 O 2 from pure water

We report the direct Z-scheme heterojunction has been recently emerging as an appealing architecture for photocatalysts design. Its efficiency depends on the interfacial and structural features of the photocatalysts. Herein, the two-dimensional ZnIn 2 S 4 nanosheets are grown on the surface of CdS hollow cubes to construct the CdS@ZnIn 2 S 4 hierarchical hollow photocatalysts with chemically bonded interface. The visualized measurements based on spatial-resolved surface photovoltage spectroscopy, combined with other spectroscopic and simulation investigations, clearly disclose that the CdS@ZnIn 2 S 4 hollow cubes constitute a highly efficient direct Z-scheme system. This accounts for the stoichiometric generation of H 2 and H 2 O 2 from pure water observed for the CdS@ZnIn 2 S 4 sulfide-only photocatalysts under visible light irradiation with an apparent quantum efficiency of 1.63 % at 400 nm. The present work demonstrates an effective protocol to achieve comprehensive insights into the charge transfer route at semiconductor heterojunction, and offers a viable way for constructing efficient sulfide-only photocatalysts for driving water splitting reaction.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2026

Chemomechanical origin of morphological disparity in lithium metal electrodeposition

The morphology of electrochemically deposited lithium (Li) critically governs the cycling stability and safety of Li metal batteries, yet the underlying controlling factors remain poorly understood. Even within the same coin cell, Li deposits can exhibit strikingly different morphologies, for example, sparse whiskers coexisting with particle-like deposits, indicting strong local variations in growth conditions. By combining cryogenic transmission electron microscopy with phase-field modeling, we here identify the root cause of this disparity. We reveal that the morphological divergence originates from the variations in structure and chemical composition of the initial SEI layer formed on the copper (Cu) substrate. These variations arise from substrate defects that induce localized electrochemical fields, driving selective dissociation of either solvent or salt, producing SEI layers with distinct chemomechanical properties. Solvent-derived organic SEI favors whisker growth, whereas salt-derived SEI promotes particle formation. This study establishes a direct mechanistic link between SEI chemomechanical properties and Li morphology, providing design principles for tailoring SEI layers to control Li deposition

energy storage↗
DOE OSTI2023

Direct Visualization of Charge Migration in Bilayer Tantalum Oxide Films by Multimodal Imaging

Inspired by biological neuromorphic computing, artificial neural networks based on crossbar arrays of bilayer tantalum oxide memristors have shown to be promising alternatives to conventional complementary metal-oxide-semiconductor (CMOS) architectures. In order to understand the driving mechanism in these oxide systems, tantalum oxide films are resistively switched by conductive atomic force microscopy (C-AFM), and subsequently imaged by kelvin probe force microscopy (KPFM) and spatially resolved time-of-flight secondary ion mass spectrometry (ToF-SIMS). These workflows enable induction and analysis of the resistive switching mechanism as well as control over the resistively switched region of the film. In this work it is shown that the resistive switching mechanism is driven by both current and electric field effects. Reversible oxygen motion is enabled by applying low (<1 V) electric fields, while high electric fields generate irreversible breakdown of the material (>1 V). Fully understanding oxygen motion and electrical effects in bilayer oxide memristor systems is a fundamental step toward the adoption of memristors as a neuromorphic computing technology.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Trends in best-in-class energy-efficient technologies for room air conditioners

Improving the efficiency of room air conditioners (RACs) could provide significant energy and associated emissions savings, particularly in emerging economies with hot climates where the cooling demand is expected to increase dramatically. To help accelerate efficiency improvements, this study identifies “best-in-class” high-efficiency RAC components and products. The findings show that manufacturers tend to minimize manufacturing costs by using RAC designs that are readily available or standardized to their production, and they share components across various models. High-efficiency RAC models use advanced compressor technologies optimized at a low frequency, large heat exchangers with thermodynamically effective materials and designs, highly efficient direct current fan motors, advanced metering devices, and smart sensors for temperature and humidity control. Recently RAC manufacturers have been improving seasonal efficiency – better reflecting part-load operation – especially via variable-speed (inverter) drives for compressor motors. Recent highest-efficiency RAC models use low global warming potential (GWP) refrigerants, having transitioned from conventional high-GWP refrigerants, in regions where RACs that use these low-GWP refrigerants are commercially available. Recently demonstrated innovative technologies show trends toward smart hybrid designs, evaporative cooling, and solid-state materials beyond the conventional vapor-compression technology. This information could help policymakers improve their RAC market-transformation programs to align with the most-efficient global technology.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗