Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “generator parameter estimation”

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 649 records · Page 36

Monte Carlo simulations of the diurnal variation in seismic detection rate of sporadic meteoroid impacts on the moon

The rate of detection of meteoroid impacts on the moon by the lunar seismic network shows a characteristic diurnal variation. Assuming that these meteoroids have a flux and a preimpact orbital distribution similar to that of fireballs observed by terrestrial camera networks, one can compute the expected diurnal variation for a given set of parameters that describe the seismic wave generation and transmission on the moon. An iterative process to match the theoretical variation with the observed one has led us to the following results: (1) the majority of the detected impact events occur within a closer range of the network than was believed earlier. This results in higher meteoroid flux estimates from lunar seismic data that agree with the terrestrially measured flux. (2) For meteoroid masses smaller than 1000 g, seismic amplitude is approximately proportional to the one-fifth power of the impact speed; for larger masses it is approximately proportional to the eight-fifth power, provided that the terrestrial meteor data used for analysis are not biased. (3) Seismic efficiency of meteoroids smaller than 1000 g is significantly less than that of large meteoroids. (4) Using orbits of fireballs that represent meteorites, we predict that the share of meteorites among the detected impacts is approximately 15 percent assuming that seismic efficiency of the high-density meteorites is the same as that of average meteoroids. A greatly increased seismic efficiency for these high-density objects is not likely.

Oberst, J.↗

A case study of the Thunderstorm Research International Project storm of July 11, 1978. I - Analysis of the data base

The data obtained at the time of the Thunderstorm Research International Project storm at the Kennedy Space Center on July 11, 1978 are analyzed in a model-independent manner. The data base included data from three Doppler radars, a lightning detection and ranging system and a network of 25 electric field mills, and rain gages. Electric field measurements were used to analyze the charge moments transferred by lightning flashes, and the data were fitted to Weibull distributions; these were used to estimate statistical parameters of the lightning for both intracloud and cloud-to-ground flashes and to estimate the fraction of the flashes which were below the observation threshold. The displacement and the conduction current densities were calculated from electric field measurements between flashes. These values were used to derive the magnitudes and the locations of dipole and monopole generators by least squares fitting the measured Maxwell current densities to the displacement-dominated equations.

Nisbet, John S.↗

HOD-dependent systematics for luminous red galaxies in the DESI 2024 BAO analysis

In this paper, we present the estimation of systematics related to the halo occupation distribution (HOD) modeling in the baryon acoustic oscillations (BAO) distance measurement of the Dark Energy Spectroscopic Instrument (DESI) 2024 analysis. This paper focuses on the study of HOD systematics for luminous red galaxies (LRG). We consider three different HOD models for LRGs, including the base 5-parameter vanilla model and two extensions to it, that we refer to as baseline and extended models. The baseline model is described by the 5 vanilla HOD parameters, an incompleteness factor and a velocity bias parameter, whereas the extended one also includes a galaxy assembly bias and a satellite profile parameter. We utilize the 25 dark matter simulations available in the AbacusSummit simulation suite at z=0.8 and generate mock catalogs for our different HOD models. To test the impact of the HOD modeling in the position of the BAO peak, we run BAO fits for all these sets of simulations and compare the best-fit BAO-scaling parameters α iso and α AP between every pair of HOD models. Furthermore, we do this for both Fourier and configuration spaces independently, using post-reconstruction measurements. We find a 3.3σ detection of HOD systematic for α AP in configuration space with an amplitude of 0.19%. For the other cases, we did not find a 3σ detection, and we decided to compute a conservative estimation of the systematic using the ensemble of shifts between all pairs of HOD models. By doing this, we quote a systematic with an amplitude of 0.07% in α iso for both Fourier and configuration spaces; and of 0.09% in α AP for Fourier space.

79 ASTRONOMY AND ASTROPHYSICS↗

Estimating Evaporative Fraction From Readily Obtainable Variables in Mangrove Forests of the Everglades, U.S.A.

A remote-sensing-based model to estimate evaporative fraction (EF) the ratio of latent heat (LE; energy equivalent of evapotranspiration -ET-) to total available energy from easily obtainable remotely-sensed and meteorological parameters is presented. This research specifically addresses the shortcomings of existing ET retrieval methods such as calibration requirements of extensive accurate in situ micro-meteorological and flux tower observations, or of a large set of coarse-resolution or model-derived input datasets. The trapezoid model is capable of generating spatially varying EF maps from standard products such as land surface temperature [T(sub s)] normalized difference vegetation index (NDVI)and daily maximum air temperature [T(sub a)]. The 2009 model results were validated at an eddy-covariance tower (Fluxnet ID: US-Skr) in the Everglades using T(sub s) and NDVI products from Landsat as well as the Moderate Resolution Imaging Spectroradiometer (MODIS) sensors. Results indicate that the model accuracy is within the range of instrument uncertainty, and is dependent on the spatial resolution and selection of end-members (i.e. wet/dry edge). The most accurate results were achieved with the T(sub s) from Landsat relative to the T(sub s) from the MODIS flown on the Terra and Aqua platforms due to the fine spatial resolution of Landsat (30 m). The bias, mean absolute percentage error and root mean square percentage error were as low as 2.9% (3.0%), 9.8% (13.3%), and 12.1% (16.1%) for Landsat-based (MODIS-based) EF estimates, respectively. Overall, this methodology shows promise for bridging the gap between temporally limited ET estimates at Landsat scales and more complex and difficult to constrain global ET remote-sensing models.

trapezoid↗

Predicting the intensity mapping signal for multi-J CO lines

We present a novel approach to estimating the intensity mapping signal of any CO rotational line emitted during the Epoch of Reionization (EoR). Our approach is based on large velocity gradient (LVG) modeling, a radiative transfer modeling technique that generates the full CO spectral line energy distribution (SLED) for a specified gas kinetic temperature, volume density, velocity gradient, molecular abundance, and column density. These parameters, which drive the physics of CO transitions and ultimately dictate the shape and amplitude of the CO SLED, can be linked to the global properties of the host galaxy, mainly the star formation rate (SFR) and the SFR surface density. By further employing an empirically derived SFR−M relation for high redshift galaxies, we can express the LVG parameters, and thus the specific intensity of any CO rotational transition, as functions of the host halo mass M and redshift z. Integrating over the range of halo masses expected to host CO-luminous galaxies, we predict a mean CO(1-0) brightness temperature ranging from ∼ 0.6 μK at z = 6 to ∼ 0.03 μK at z = 10 with brightness temperature fluctuations of Δ{sub CO}{sup 2} ∼ 0.1 and 0.005 μK respectively, at k = 0.1 Mpc{sup −1}. In this model, the CO emission signal remains strong for higher rotational levels at z = 6, with ( T{sub CO} ) ∼ 0.3 and 0.05 μK for the CO J = 6arrow5 and CO J = 10arrow9 transitions respectively. Including the effects of CO photodissociation in these molecular clouds, especially at low metallicities, results in the overall reduction in the amplitude of the CO signal, with the low- and high-J lines weakening by 2–20% and 10–45%, respectively, over the redshift range 4 < z < 10.

79 ASTRONOMY AND ASTROPHYSICS↗

NEON AOP Survey of Upper East River CO Watersheds: Waveform LiDAR Binary Data

The waveform Light Detection and Ranging (LiDAR) data in this package were generated through a National Ecological Observatory Network Airborne Observation Platform (NEON AOP) acquisition over watersheds of interest surrounding Crested Butte, Colorado. The remote sensing imagery acquired by the NEON AOP supports an interdisciplinary project on hydrology, biogeochemistry, and ecosystem functioning in a snow-dominated headwater environment. These waveform LiDAR data enable spatially continuous estimation of vegetation structure parameters to facilitate analyses of the major environmental drivers of structural and compositional variability. The package contains 97 compressed file archives in 7-zip (.7z) format, each corresponding to one acquisition flightpath. Within each .7z archive is a set of constituent files describing properties of the LiDAR waveforms, such as return intensity, geolocation, outgoing pulse and other behavior of the sensor and signals. Once downloaded, the files must first be unzipped using the widely distributed command-line software utility 7z, using the command '7z x \[filename\].7z \[target_directory\]'. All files within the .7z archives can be opened in IDL, MatLab, or the open-source R statistical computing environment. Further details about the data package are in the attached user guide (neon_aop_crbu_waveformlidar_userguide.pdf). This work was supported by the Watershed Function Science Focus Area at Lawrence Berkeley National Laboratory funded by the US Department of Energy, Office of Science, Biological and Environmental Research under Contract No. DE-AC02-05CH11231.

2018 NEON and 2025 CHESS Campaigns↗

Moisture Durability Assessment of Selected Well-insulated Wall Assemblies

This report presents the results from studying the hygrothermal performance of two well-insulated wall assemblies, both complying with and exceeding international building codes (IECC 2015 2014, IRC 2015). The hygrothermal performance of walls is affected by a large number of influential parameters (e.g., outdoor and indoor climates, workmanship, material properties). This study was based on a probabilistic risk assessment in which a number of these influential parameters were simulated with their natural variability. The purpose of this approach was to generate simulation results based on laboratory chamber measurements that represent a variety of performances and thus better mimic realistic conditions. In total, laboratory measurements and 6,000 simulations were completed for five different US climate zones. A mold growth indicator (MGI) was used to estimate the risk of mold which potentially can cause moisture durability problems in the selected wall assemblies. Analyzing the possible impact on the indoor climate due to mold was not part of this study. The following conclusions can be reached from analyzing the simulation results. In a hot-humid climate, a higher R-value increases the importance of the airtightness because interior wall materials are at lower temperatures. In a cold climate, indoor humidity levels increase with increased airtightness. Air leakage must be considered in a hygrothermal risk assessment, since air efficiently brings moisture into buildings from either the interior or exterior environment. The sensitivity analysis of this study identifies mitigation strategies. Again, it is important to remark that MGI is an indicator of mold, not an indicator of indoor air quality and that mold is the most conservative indicator for moisture durability issues.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

On three dimensional object recognition and pose-determination: An abstraction based approach

A method of computing reliable Gaussian and mean curvature sign-map descriptors from the polynomial approximation of surfaces was demonstrated. Such descriptors which are invariant under perspective variation are suitable for hypothesis generation. A means for determining the pose of constructed geometric forms whose algebraic surface descriptors are nonlinear in terms of their orienting parameters was developed. This was done by means of linear functions which are capable of approximating nonlinear forms and determining their parameters. It was shown that biquadratic surfaces are suitable companion linear forms for cylindrical approximation and parameter estimation. The estimates provided the initial parametric approximations necessary for a nonlinear regression stage to fine tune the estimates by fitting the actual nonlinear form to the data. A hypothesis-based split-merge algorithm for extraction and pose determination of cylinders and planes which merge smoothly into other surfaces was developed. It was shown that all split-merge algorithms are hypothesis-based. A finite-state algorithm for the extraction of the boundaries of run-length regions was developed. The computation takes advantage of the run list topology and boundary direction constraints implicit in the run-length encoding.

Quek, Kok How Francis↗

Sub-mm Scale Fiber Guided Deep/Vacuum Ultra-Violet Optical Source for Trapped Mercury Ion Clocks

We demonstrate the functionality of a mercury capillary lamp with a diameter in the sub-mm range and deep ultraviolet (DUV)/ vacuum ultraviolet (VUV) radiation delivery via an optical fiber integrated with the capillary. DUV spectrum control is observed by varying the fabrication parameters such as buffer gas type and pressure, capillary diameter, electrical resonator design, and temperature. We also show spectroscopic data of the 199Hg+ hyper-fine transition at 40.5GHz when applying the above fiber optical design. We present efforts toward micro-plasma generation in hollow-core photonic crystal fiber with related optical design and theoretical estimations. This new approach towards a more practical DUV optical interface could benefit trapped ion clock developments for future ultra-stable frequency reference and time-keeping applications.

lamps↗

Influence of Vegetation Canopy Architecture on Dynamics of Spatiotemporal Light Distribution and its Effect on Land-Biosphere-Atmosphere Fluxes of Energy, Carbon, and Water

The dynamic spatiotemporal distribution of light and micro climatic variables (wind speed, atmospheric CO 2 concentration, humidity, pressure, and air temperature) within the plant canopy as determined by the canopy architecture play a critical role in land-biosphere-atmosphere (LBA) photosynthetic carbon uptake, transpiratory water loss and energy partitioning. Vegetation models that include effects of sun and shade flecks, leaf temperature and energy balance partitioning between sensible and latent heat fluxes, and reduced turbulent mixing in dense vegetation canopies can accurately predict canopy performance, particularly its shaded leaves, which can constitute 80% of leaf area and contribute up to 50% of carbon update and 65% of water loss. However, current vegetation models with their bulk characterization of shaded leaves poorly represent the complex three dimensional (3D) canopy architecture are insufficient to capture the spatiotemporal dynamics of shaded leaves. Here, we propose to employ an x-ray computed tomography (XCT) imaging technique in combination with an automated, high throughput, image processing pipeline (3DArch, to be developed) to obtain an explicit 3D triangulated surface mesh of a plant canopy (Fig 1). 3DArch would provide a much needed technological jump in 3D LBA modeling where computer generated artificial 3D models suffer from lack of high resolution detailed canopy architecture parameters. The 3D canopy surface mesh output from 3DArch is coupled with an explicit 3D forward ray tracing model (fastTracer) to obtain a high resolution representation of the dynamic spatiotemporal light environment within the canopy over the course of a day (Fig 1). The results from the light attenuation model are then coupled with a process based multi-scale bio-physical canopy and leaf model of photosynthesis (PlantGro) to get fluxes of carbon, water and energy at the leaf and canopy scale, in addition to estimating whole plant carbon uptake, and water loss. Our proposed model framework would be the first of its kind to provide a breakthrough capability in the field of process-based plant modeling that bridges the key technological gap in shaded leaf modeling through explicit representation of 3D canopy architecture. Computer generated artificial 3D model canopy simulations on grasses show that the canopy light distribution does not follow an exponential behavior as assumed in most 1D vegetation models, but rather exhibits a heterogeneous behavior, with a small but significant portion of interior canopy leaves experiencing peak light levels due to sun flecks while the rest of the canopy receives substantially lower light levels. These results indicate a bias in traditional models that overestimates shaded leaf contribution towards whole canopy photosynthesis by up to 15% and underestimate transpiration by up to 12% with an associated increase in sensible heat fluxes since photosynthesis, transpiration and energy partitioning are tightly coupled processes. However, these test results need to be validated using the 3D architecture obtained from a real plant imaged in 3D using 3DArch, rather than an artificial model plant. Our proposed model framework would allow us to identify inefficiencies in canopy light distribution and prescribe optimal canopy architectures (leaf optical properties, leaf area, leaf angles, row spacing, planting density, etc.) to traditional and genetic engineering based crop breeders, for maximizing target functions such as light use efficiency, water use efficiency, drought tolerance, crop yield, crop resilience etc., with important applications in biofuel crops such as sorghum, and switchgrass. Our model can also be applied for designing optimal canopy architecture that increase surface albedo as a mitigation strategy for global warming and projections of crop yield for future food and water security under climate change.

59 BASIC BIOLOGICAL SCIENCES↗

Overcoming Anode Instability in Solid‐State Batteries through Control of the Lithium Metal Microstructure

Abstract Enabling the lithium metal anode (LMA) in solid‐state batteries (SSBs) is the key to developing high energy density battery technologies. However, maintaining a stable electrode–electrolyte interface presents a critical challenge to high cycling rate and prolonged cycle life. One such issue is the interfacial pore formation in LMA during stripping. To overcome this, either higher stack pressure or binary lithium alloy anodes are used. Herein, it is shown that fine‐grained ( d = 20 µm) polycrystalline LMA can avoid pore formation by exploiting the microstructural dependence of the creep rates. In a symmetric cell set‐up, i.e., LiǀLi 6.25 Al 0.25 La 3 Zr 2 O 12 (LLZO)ǀLi, fine‐grained LMA achieves > 11.0 mAh cm −2 compared to ≈ 3.6 mAh cm −2 for coarse‐grained LMA ( d = 295 µm) at 0.1 mA cm −2 and at moderate stress of 2.0 MPa. Smaller diffusion lengths (≈ 20 µm) and higher diffusivity pathway along dislocations ( D d ≈ 10 −7 cm 2 s −1 ), generated during cell fabrication, result in enhanced viscoplastic deformation in fine‐grained polycrystalline LMA. The electrochemical performances corroborate well with estimated creep rates. Thus, microstructural control of LMA can significantly reduce the required stack pressure during stripping. These results are particularly relevant for “anode‐free” SSBs wherein both the microstructure and the mechanical state of the lithium are critical parameters.

25 ENERGY STORAGE↗

Bayesian Parameter Estimation of the k-ω Shear Stress Transport Model for Accurate Simulations of Impinging-Jet Heat Transfer

Given the lack of fusion-relevant component test facilities, current estimates of the thermo-fluid performance of plasma-facing components are based for the most part on numerical simulations. A major source of uncertainty in these simulations is the semiempirical turbulence (closure) models for the Reynolds stresses appearing in the governing Reynolds-averaged Navier-Stokes equations, which involve a set of constants that depend upon the flow. The objective of this study is to evaluate Bayesian parameter estimation of turbulence closure constants in ANSYS Fluent to model heat transfer in impinging jets. The Bayesian statistical calibration produces a probability distribution for these constants from experimental data; the maximum a posteriori estimates are then taken to be the calibrated constants, or parameters. The turbulence model constants are calibrated using an experimental study of a submerged jet of air impinging on a flat heated surface at Reynolds numbers Re = O(10 4 ) and impingement distance in jet diameters H/d = 2. Numerical predictions using the calibrated model parameters are then compared with those generated using the default constants. Predictions obtained with model parameters calibrated on datasets of two different sizes are compared to evaluate the effect of the number of calibration samples. Lastly, the extrapolative ability of the calibrated model is examined by predictions at a Re beyond the calibration values.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Use of multiple-stage remote sensing techniques to develop forest stocking equations

Portions of a forest supporting density levels that affect natural resource yields must be quantified to estimate the operational feasibility of a proposed management system. Estimates of average parameters do not necessarily provide complete knowledge, particularly with frequently skewed forest parameters. Instead, another statistic, the portion of a forest supporting minimum density levels would be useful to set realistic limits to the implementation of management systems and to establish priorities for operational programs. Such a statistic can be obtained from solutions of forest stocking equations generated by applications of multiple-stage remote sensing techniques.

Bisson, H. R.↗

The area-time-integral technique to estimate convective rain volumes over areas applied to satellite data - A preliminary investigation

The use of the area-time-integral (ATI) technique, based only on satellite data, to estimate convective rain volume over a moving target is examined. The technique is based on the correlation between the radar echo area coverage integrated over the lifetime of the storm and the radar estimated rain volume. The processing of the GOES and radar data collected in 1981 is described. The radar and satellite parameters for six convective clusters from storm events occurring on June 12 and July 2, 1981 are analyzed and compared in terms of time steps and cluster lifetimes. Rain volume is calculated by first using the regression analysis to generate the regression equation used to obtain the ATI; the ATI versus rain volume relation is then employed to compute rain volume. The data reveal that the ATI technique using satellite data is applicable to the calculation of rain volume.

Doneaud, Andre A.↗

Dynamic, resilient virtual sensing system and shadow controller for cyber-attack neutralization

An industrial asset may have monitoring nodes (e.g., sensor or actuator nodes) that generate current monitoring node values. An abnormality detection and localization computer may receive the series of current monitoring node values and output an indication of at least one abnormal monitoring node that is currently being attacked or experiencing a fault. An actor-critic platform may tune a dynamic, resilient state estimator for a sensor node and output tuning parameters for a controller that improve operation of the industrial asset during the current attack or fault. The actor-critic platform may include, for example, a dynamic, resilient state estimator, an actor model, and a critic model. According to some embodiments, a value function of the critic model is updated for each action of the actor model and each action of the actor model is evaluated by the critic model to update a policy of the actor-critic platform.

Roychowdhury, Subhrajit↗

On second-harmonic generation in nonuniformly magnetized media

This paper presents a theoretical study of new effects of double-frequency optical signal generation in magnetic media with a nonuniform magnetization distribution. The study is based on the hydrodynamic approximation of the motion of conduction electrons in the field of an electromagnetic wave and equations for the mean electron spin. Within this approach, the mechanisms of the toroidal moment effect, which was experimentally discovered earlier in a system of magnetic particles with a vertical magnetization distribution as well as in a multilayered magnetic system (which is a collinear nonuniform magnetic system), are demonstrated. A new effect that appears due to the presence of equilibrium spin currents in a nonuniform magnetic system is studied in detail. This effect was predicted from the symmetry and can occur only in noncollinearly magnetized media. It is shown that this effect has a resonance nature with a resonance at a pump frequency equal to the plasma frequency of conduction electrons. Estimates of the susceptibility at the double frequency for the parameters of nickel and the typical scale of variation in the magnetization, which is equal to 10 nm, show that the susceptibility near the frequency of the plasma resonance at the chosen parameters is on the order of 10{sup –9} esu, which holds up a hope of the experimental detection of this effect.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

On hydrologic similarity: A dimensionless flood frequency model using a generalized geomorphologic unit hydrograph and partial area runoff generation

One of the shortcomings of the original theory of the geomorphologic unit hydrograph (GUH) is that it assumes that runoff is generated uniformly from the entire catchment area. It is now recognized that in many catchments much of the runoff during storm events is produced on partial areas which usually form on narrow bands along the stream network. A storm response model that includes runoff generation on partial areas by both Hortonian and Dunne mechanisms was recently developed by the authors. In this paper a methodology for integrating this partial area runoff generation model with the GUH-based runoff routing model is presented; this leads to a generalized GUH. The generalized GUH and the storm response model are then used to estimate physically based flood frequency distributions. In most previous work the initial moisture state of the catchment had been assumed to be constant for all the storms. In this paper we relax this assumption and allow the initial moisture conditions to vary between storms. The resulting flood frequency distributions are cast in a scaled dimensionless framework where issues such as catchment scale and similarity can be conveniently addressed. A number of experiments are performed to study the sensitivity of the flood frequency response to some of the 'similarity' parameters identified in this formulation. The results indicate that one of the most important components of the derived flood frequency model relates to the specification of processes within the runoff generation model; specifically the inclusion of both saturation excess and Horton infiltration excess runoff production mechanisms. The dominance of these mechanisms over different return periods of the flood frequency distribution can significantly affect the distributional shape and confidence limits about the distribution. Comparisons with observed flood distributions seem to indicate that such mixed runoff production mechanisms influence flood distribution shape. The sensitivity analysis also indicated that the incorporation of basin and rainfall storm scale also greatly influences the distributional shape of the flood frequency curve.

Sivapalan, Murugesu↗

Improving Distribution Resiliency with Microgrids and State and Parameter Estimation

Modern society relies on low-cost reliable electrical power, both to maintain industry, as well as provide basic social services to the populace. When major disturbances occur, such as Hurricane Katrina or Hurricane Sandy, the nation’s electrical infrastructure can experience significant outages. To help prevent the spread of these outages, as well as facilitating faster restoration after an outage, various aspects of improving the resiliency of the power system are needed. Two such approaches are breaking the system into smaller microgrid sections, and to have improved insight into the operations to detect failures or mis-operations before they become critical. Breaking the system into smaller sections of microgrid islands, power can be maintained in smaller areas where distribution generation and energy storage resources are still available, but bulk power generation is no longer connected. Additionally, microgrid systems can maintain service to local pockets of customers when there has been extensive damage to the local distribution system. However, microgrids are grid connected a majority of the time and implementing and operating a microgrid is much different than when islanded. This report discusses work conducted by the Pacific Northwest National Laboratory that developed improvements for simulation tools to capture the characteristics of microgrids and how they can be used to develop new operational strategies. These operational strategies reduce the cost of microgrid operation and increase the reliability and resilience of the nation’s electricity infrastructure. In addition to the ability to break the system into microgrids, improved observability into the state of the distribution grid can make the power system more resilient. State estimation on the transmission system already provides great insight into grid operations and detecting abnormal conditions by leveraging existing measurements. These transmission-level approaches are expanded to using advanced metering infrastructure and other distribution-level measurements to create a three-phase, unbalanced distribution state estimation approach. With distribution-level state estimation, the grid can be operated more efficiently, and outages or equipment failures can be caught faster, improving the overall resilience and reliability of the grid.

24 POWER TRANSMISSION AND DISTRIBUTION↗