Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Normal Distribution”

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,207 records · Page 67

Modeling The Effects of Shock Pressure and Pore Morphology on Hot Spot Mechanisms in HMX

Here, we investigate the effects of shock pressure and pore morphology on the formation and growth of hot spots in HMX (octahydro-1,3,5,7-tetranitro-1,3,5,7-tetrazocine). Both non-reactive and reactive ALE3D simulations are used in these studies. Our non-reactive simulations show a viscous-dominated pore collapse mode at lower shock pressures (2–10 GPa) with shear band formation and a hydrodynamic-dominated mode at higher shock pressures (20-40 GPa) due to bulk melting. When normalized by bulk shock heating, viscous-dominated pore collapse modes are more efficient at generating hot spots. Pore morphology influences the post-collapse temperature distributions and reaction rate for a fixed pore area and shock pressure. We find that multiple surface pores at the binder-grain interface tend to react the fastest. Due to their upstream location in the HMX grain, the surface pores collapse sooner than interior pores; thus, the extent of reaction will generally favor these morphologies because they have more time to grow. In general, multiple smaller hot spots tend to react faster than a single larger hot spot because they accelerate one another's burning. The rank order of morphology effects, however, is not the same for non-reactive and reactive simulations. For example, while multiple surface pores produce the highest reaction rates they do not produce the highest (non-reactive) hot spot temperatures. In conclusion, our numerical studies provide insights on hot spot mechanisms in lieu of direct measurements and can be used to develop advanced shock initiation models.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Electrospray Generated from the Tip-Sealed Fine Glass Capillary Inserted with an Acupuncture Needle Electrode

In electrospray, excess charges are supplied to a sample solution by the occurrence of electrochemical reactions. Recently, different versions of electrospray, e.g., dielectric barrier electrospray ionization, inductive desorption electrospray ionization, and electrostatic-ionization driven by dielectric polarization, have been reported in which the sample solution was not in direct contact with the metal electrode but separated by dielectric materials. The objective of the current work is to elucidate the mechanism of dielectric barrier electrospray. A sealed borosilicate glass capillary inserted with a fine acupuncture needle was used as a probe. A sample solution (~ 400 nL) was captured on the glass capillary tip and a positive high voltage (HV) pulse (+ 4.5 kV) was applied to the internal metal electrode. Mass spectra were measured as a function of the HV pulse width from μs to 10 s. Ions started to be detected with the pulse width of ~ 5 ms. The ion intensities increased slowly with time and reached a plateau in a few seconds. The charge distribution of cytochrome c [M + nH]{sup n+} shifted to higher n values from a few ms to seconds. In addition to cone-jet mode normal electrospray that lasted until all the liquid sample was depleted from the glass tip, the polarization-induced electrospray ionization was observed at the early stage of the HV application. .

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Mitochondrial nucleoid clusters protect newly synthesized mtDNA during Doxorubicin- and Ethidium Bromide-induced mitochondrial stress

Mitochondrial DNA (mtDNA) is compacted in ribonucleoprotein complexes called nucleoids, which can divide or move within the mitochondrial network. Mitochondrial nucleoids are able to aggregate into clusters upon reaction with intercalators such as the mtDNA depletion agent Ethidium Bromide (EB) or anticancer drug Doxorobicin (DXR). However, the exact mechanism of nucleoid clusters formation remains unknown. Resolving these processes may help to elucidate the mechanisms of DXR-induced cardiotoxicity. Therefore, we addressed the role of two key nucleoid proteins; mitochondrial transcription factor A (TFAM) and mitochondrial single-stranded binding protein (mtSSB); in the formation of mitochondrial nucleoid clusters during the action of intercalators. We found that both intercalators cause numerous aberrations due to perturbing their native status. By blocking mtDNA replication, both agents also prevented mtDNA association with TFAM, consequently causing nucleoid aggregation into large nucleoid clusters enriched with TFAM, co-existing with the normal nucleoid population. In the later stages of intercalation (> 48 h), TFAM levels were reduced to 25%. In contrast, mtSSB was released from mtDNA and freely distributed within the mitochondrial network. Nucleoid clusters mostly contained nucleoids with newly replicated mtDNA, however the nucleoid population which was not in replication mode remained outside the clusters. Moreover, the nucleoid clusters were enriched with p53, an anti-oncogenic gatekeeper. We suggest that mitochondrial nucleoid clustering is a mechanism for protecting nucleoids with newly replicated DNA against intercalators mediating genotoxic stress. These results provide new insight into the common mitochondrial response to mtDNA stress and can be implied also on DXR-induced mitochondrial cytotoxicity. - Highlights: • The mechanism for mitochondrial nucleoid clustering is proposed. • DNA intercalators (Doxorubicin or Ethidium Bromide) prevent TFAM binding to mtDNA. • Replicating nucleoids are less prone to DNA intercalator and preserve more TFAM. • Nucleoid clusters mostly contain nucleoids with newly replicated mtDNA. • Recently replicated nucleoids are protected in clusters by increased TFAM and p53.

60 APPLIED LIFE SCIENCES↗

Development of low angle grain boundaries in lightly deformed superconducting niobium and their influence on hydride distribution and flux perturbation

This study shows that low angle grain boundaries (LAGBs) can be created by small 5% strains in high purity (RRR ≥ 200) SRF-grade single crystalline niobium (Nb) and that these boundaries act as hydrogen traps as indicated by the distribution of niobium hydrides (Nb 1-x H x ). Nb 1-x H x is detrimental to superconducting radio frequency (SRF) Nb cavities due to its normal conducting properties at cavity operating temperatures. By designing a single crystal tensile sample extracted from a large grain (>5 cm) Nb ingot slice for preferred slip on one slip plane, LAGBs and dense dislocation boundaries developed. With chemical surface treatments following standard SRF cavity fabrication practice, Nb1-xHx phases were densely precipitated at the LAGBs upon cryogenic cooling (8-10 K/min). Micro-crystallographic analysis confirmed heterogeneous hydride precipitation, which included significant hydrogen atom accumulation in LAGBs. Magneto-optical imaging (MOI) analysis showed that these sites can then act as sites for both premature flux penetration and eventually flux trapping. However, this hydrogen related degradation at LAGBs did not completely disappear even after a 800 °C/2hrs anneal typically used for hydrogen removal in SRF Nb cavities. These findings suggest that hydride precipitation at a LAGB is facilitated by a non-equilibrium concentration of vacancy-hydrogen (H) complexes aided by mechanical deformation and the hydride phase interferes with the recovery process under 800°C annealing.

36 MATERIALS SCIENCE↗

Residual strain orientation in rolled titanium determined with synchrotron X-ray Laue microdiffraction

Previously, synchrotron X-ray Laue microdiffraction has been used to measure the magnitudes of residual strain in materials. Recently the method was advanced to determine the orientation of the strain ellipsoid and applied to naturally deformed quartzites; however, the deformation history of these quartzites is ambiguous due to their natural origin. Here, in this study, synchrotron X-ray Laue microdiffraction (µXRD) is used to measure the residual strain for the first time in a sample with known stress history, rolled titanium. A deviatoric strain tensor is calculated from each Laue diffraction image collected with two µXRD scans of a rolled titanium sheet in different sample orientations. The principal strain axes are calculated using an eigen decomposition of the deviatoric strain tensors. The results show that the principal axis of compression is aligned with the normal direction of the titanium sheet, and the principal axis of extension is aligned with the rolling direction. Pole figures are used to represent the 3D distribution of residual strain axes.

36 MATERIALS SCIENCE↗

Technical Note: Study of the electron transport parameters used in PENELOPE for the Monte Carlo simulation of Linac targets

Purpose: The Monte Carlo simulation of electron transport in Linac targets using the condensed history technique is known to be problematic owing to a potential dependence of absorbed dose distributions on the electron step length. In the PENELOPE code, the step length is partially determined by the transport parameters C1 and C2. The authors have investigated the effect on the absorbed dose distribution of the values given to these parameters in the target. Methods: A monoenergetic 6.26 MeV electron pencil beam from a point source was simulated impinging normally on a cylindrical tungsten target. Electrons leaving the tungsten were discarded. Radial absorbed dose profiles were obtained at 1.5 cm of depth in a water phantom located at 100 cm for values of C1 and C2 in the target both equal to 0.1, 0.01, or 0.001. A detailed simulation case was also considered and taken as the reference. Additionally, lateral dose profiles were estimated and compared with experimental measurements for a 6 MV photon beam of a Varian Clinac 2100 for the cases of C1 and C2 both set to 0.1 or 0.001 in the target. Results: On the central axis, the dose obtained for the case C1 = C2 = 0.1 shows a deviation of (17.2% ± 1.2%) with respect to the detailed simulation. This difference decreases to (3.7% ± 1.2%) for the case C1 = C2 = 0.01. The case C1 = C2 = 0.001 produces a radial dose profile that is equivalent to that of the detailed simulation within the reached statistical uncertainty of 1%. The effect is also appreciable in the crossline dose profiles estimated for the realistic geometry of the Linac. In another simulation, it was shown that the error made by choosing inappropriate transport parameters can be masked by tuning the energy and focal spot size of the initial beam. Conclusions: The use of large path lengths for the condensed simulation of electrons in a Linac target with PENELOPE conducts to deviations of the dose in the patient or phantom. Based on the results obtained in this work, values of C1 and C2 larger than 0.001 should not be used in Linac targets without further investigation.

07 ISOTOPE AND RADIATION SOURCES↗

High-Temperature Gas-Cooled Pebble-Bed Reactors Running In And Transient Modeling Capabilities Demonstration

This study presents a comprehensive benchmarking and verification effort of several thermal-hydraulic and multiphysics capabilities for high-temperature gas-cooled reactor (HTGR) applications. The first part of this effort focuses on the running-in verification of Griffin's multiphysics capabilities, specifically for simulating the evolution of Pebble Bed reactor cores from startup to equilibrium. In the absence of validation data, code-to-code comparisons are conducted with Kugelpy, showing good agreement for key quantities like maximum power density and fresh core k-eigenvalue predictions. However, discrepancies in equilibrium core predictions suggest potential issues with cross sections, underscoring the need for further refinement and evaluation. The HTTF system analysis code benchmark involves RELAP5-3D, SAM, and GAMMA+ to assess their predictive capabilities for HTTF behavior under both normal operation and pressurized conduction cooldown (PCC) transient conditions. While there is good agreement in predicting major parameters such as coolant temperature, solid temperature, and flow distribution, discrepancies in transient behavior highlight differences in modeling approaches, nodalizations, and heat transfer models. The HTTF lower plenum CFD benchmark employs nekRS to simulate flow mixing phenomena, successfully capturing relevant flow physics and demonstrating mesh independence in complex geometries. Preliminary results suggest a relatively uniform temperature field but significant unsteadiness in the flow, requiring time-averaging analyses. The GPBR200 system analysis code benchmark uses SAM's core channel and porous media models, incorporating an RCCS loop for decay heat removal. During steady-state and transient conditions, including protected de-pressurized and pressurized loss of forced cooling (DLOFC and PLOFC), both models show good agreement in predicting temperature profiles and key parameters. Notably, while the core channel model underpredicts convective heat transfer effects, both models maintain temperatures well below the TRISO fuel safety limit. These benchmarking efforts collectively enhance the predictive capabilities of the tools used in HTGR design and safety analysis, guiding developments to improve their accuracy and applicability.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

A numerical simulation of CCN impacts on weather modification efficiency

Aerosols affect development of clouds and precipitation by serving as cloud condensation nuclei (CCN) and ice nuclei (IN). Considering the dramatically changing ambient aerosol concentration, it is important to examine the potential “side effect” of aerosol pollution on precipitation enhancement by weather modification. In this study, the cloud seeding was performed on a precipitation event in Beijing in the summer of 2008, which is simulated by the NSSL two-moment cloud scheme of the Weather Research and Forecasting (WRF) model. Sensitivity tests were conducted by modifying the ambient aerosol concentration and the ice crystal seeding amount to investigate the cloud seeding efficacy in different CCN concentration scenarios. There was a slight difference in the precipitation distribution between the simulations with two ambient CCN concentrations: the northern precipitation center in polluted scenario was weaker and the southern center was stronger. Compared with normal CCN scenario, the cloud liquid water mass and ice crystal mass in the severe pollution scenario is larger, and the total contents of snow and graupel were not sensitive to the CCN concentration. With the same amount of man-made ice crystals seeding, the precipitation enhancement was quite different under different CCN conditions. The higher the CCN concentration usually leads to stronger precipitation suppression. As CCN concentration increase, the deposition growth of snow, auto-conversion and accretion of ice crystals to snow were weakened, as well as the conversion of melting snow and graupel into rainwater.

54 ENVIRONMENTAL SCIENCES↗

A Study of the Properties and Dynamics of the Disk of Satellites in a Milky-Way-like Galaxy System

The dynamics of the satellite systems of Milky-Way-like galaxies offer a useful means by which to study the galaxy formation process in the cosmological context. It has been suggested that the currently observed anisotropic distribution of the satellites in such galaxy systems is inconsistent with the concordance ΛCDM cosmology model on the galactic scale if the observed satellites are random samples of the dark matter (DM) sub-halos that are nearly isotropically distributed around the central galaxy. In this study, we present original high-resolution zoom-in studies of central galaxies and satellite systems based upon initial conditions for the DM distribution from the Aquarius simulations but with substantial high-resolution baryon physics added. We find that the galaxy most like the Milky Way in this study does indeed contain a disk of satellites (DOS). Although one galaxy DOS system does not answer the question of how common such disks are, it does allow the opportunity to explore the properties and dynamics of the DOS system. Our investigation centers on the spatial arrangement (distances, angles, etc.) of satellites in this Milky-Way-like galaxy system with a specific emphasis on identifying and analyzing the disk-like structure along with its dynamical and morphological properties. Among the conclusions from this study, we find that the satellites and DM sub-halos in the galaxy simulations are anisotropically distributed. The dynamical properties of the satellites, however, indicate that the direction of the angular momentum vector of the whole satellite system is different from the normal direction of the fitted DOS and from the normal direction of the velocity dispersion of the system. Hence, the fitted DOS appears to be comprised of infalling sub-halos and is not a rotationally supported system.

Zhao, Xinghai↗

Design and Commissioning of a Magnetic Field Scanning System for SRF Cavities

Trapped magnetic vortices are one of the leading sources of residual losses in SRF cavities. Mechanisms of flux pinning depend on the materials treatment and cool-down conditions. A magnetic field scanning system using flux-gate magnetometers and Hall probes has been designed and built to allow measuring the local magnetic field of trapped vortices normal to the outer surface of 1.3 GHz single-cell SRF cavities at cryogenic temperatures. Such system will allow inferring the key information about the distribution and magnitude of trapped flux in the SRF cavities for different material, surface preparations and cool-down conditions.

Parajuli, I. P.↗

Profiling of adrenocorticotropic hormone and arginine vasopressin in human pituitary gland and tumor thin tissue sections using droplet-based liquid-microjunction surface-sampling-HPLC–ESI-MS–MS

Described here are the results from the profiling of the proteins arginine vasopressin (AVP) and adrenocorticotropic hormone (ACTH) from normal human pituitary gland and pituitary adenoma tissue sections using a fully automated droplet-based liquid microjunction surface sampling-HPLC-ESI-MS/MS system for spatially resolved sampling, HPLC separation, and mass spectral detection. Excellent correlation was found between the protein distribution data obtained with this droplet-based liquid microjunction surface sampling-HPLC-ESI-MS/MS system and those data obtained with matrix assisted laser desorption ionization (MALDI) chemical imaging analyses of serial sections of the same tissue. The protein distributions correlated with the visible anatomic pattern of the pituitary gland. AVP was most abundant in the posterior pituitary gland region (neurohypophysis) and ATCH was dominant in the anterior pituitary gland region (adenohypophysis). The relative amounts of AVP and ACTH sampled from a series of ACTH secreting and non-secreting pituitary adenomas correlated with histopathological evaluation. ACTH was readily detected at significantly higher levels in regions of ACTH secreting adenomas and in normal anterior adenohypophysis compared to non-secreting adenoma and neurohypophysis. AVP was mostly detected in normal neurohypophysis as anticipated. This work demonstrates that a fully automated droplet-based liquid microjunction surface sampling system coupled to HPLC-ESI-MS/MS can be readily used for spatially resolved sampling, separation, detection, and semi-quantitation of physiologically-relevant peptide and protein hormones, such as AVP and ACTH, directly from human tissue. In addition, the relative simplicity, rapidity and specificity of the current methodology support the potential of this basic technology with further advancement for assisting surgical decision-making.

59 BASIC BIOLOGICAL SCIENCES↗

Explanation and Application of the Evolving Contact Traction Fields in Shallow Foundation Systems

The present paper provides a qualitative discussion of the evolution of contact traction fields beneath rigid shallow foundations resting on granular materials. A phenomenological similarity is recognized in the measured contact traction fields of rigid footings and at the bases of sandpiles. This observation leads to the hypothesis that the stress distributions are brought about by the same physical phenomena, namely the development of arching effects through force chains and mobilized intergranular friction. A set of semi-empirical equations are suggested for the normal and tangential components of this contact traction based on past experimental measurements and phenomenological assumptions of frictional behaviors at the foundation system scale. These equations are then applied to the prescribed boundary conditions for the analysis of the settlement, resistance, and stress fields in supporting granular materials beneath the footing. A parametric sensitivity study is performed on the proposed modelling method, highlighting solutions to the boundary-value problems in an isotropic, homogeneous elastic half-space.

42 ENGINEERING↗

Progress on the National Solar Radiation Data Base (NSRDB): A New DNI Computation

This study introduces a new technique to compute direct normal irradiance (DNI) for improving the National Solar Radiation Data Base (NSRDB). A finite-surface integration algorithm is developed to compute solar radiation in differential solid angles and efficiently infer its contribution to a surface perpendicular to the solar direction. A lookup table of cloud bi-directional transmittance distribution function (BTDF) is developed by use of the discrete ordinates radiative transfer (DISORT) model for possible solar and observing directions and various cloud optical and microphysical properties. In each solar incident direction, DNI is given by the cloud BTDFs from approximately 200 differential solid angles. The simulated DNI is calibrated and evaluated using surface observations at the National Renewable Energy Laboratory's (NREL's) Solar Energy Research Laboratory (SRRL) and the Atmospheric Radiation Measurement (ARM) Southern Great Plains (SGP) facility.

41 EE - Solar Energy Technologies Office (EE-4S)↗

Triggering the Formation of Direct Collapse Black Holes by Their Congeners

Direct collapse black holes (DCBHs) are excellent candidates for seeds of supermassive black holes observed at z ≳ 6. The formation of a DCBH requires a strong external radiation field to suppress H{sub 2} formation and cooling in a collapsing gas cloud. Such a strong field is not easily achieved by first stars or normal star-forming galaxies. Here we investigate a scenario in which a previously formed DCBH can provide the necessary radiation field for the formation of additional ones. Using a one-zone model and simulated DCBH Spectral Energy Distributions (SEDs) filtered through absorbing gas initially having column density N {sub H}, we derive the critical field intensity, J{sub LW}{sup crit}, to suppress H{sub 2} formation and cooling. For the SED model with N{sub H}=1.3×10{sup 25} cm{sup −2}, 8.0 × 10{sup 24} cm{sup −2}, and 5.0 × 10{sup 24} cm{sup −2}, we obtain J{sub LW}{sup crit}≈22, 35, and 54, all much smaller than the critical field intensity for normal star-forming galaxies J{sub LW}{sup crit}≳1000 X-ray photons from previously formed DCBHs build up a high-z X-ray background (XRB) that may boost the J{sub LW}{sup crit}. However, we find that in the three SED models, J{sub LW}{sup crit} only increases to ≈80, 170, and 390, even when ρ-dot {sub ∙} reaches the maximum value allowed by the present-day XRB level (0.22, 0.034, 0.006 M {sub ⊙} yr{sup −1} Mpc{sup −3}), which is still much smaller than the galactic value. Although considering the XRB from first galaxies may further increase J{sub LW}{sup crit}, we conclude that our investigation supports a scenario in which DCBHs may be more abundant than predicted by models only including galaxies as external radiation sources.

79 ASTRONOMY AND ASTROPHYSICS↗

Biparametric analyses of charge trapping in Cd{sub 0.9}Zn{sub 0.1}Te based virtual Frisch grid detectors

Detector grade Cd{sub 0.9}Zn{sub 0.1}Te (CZT) single crystals were grown from zone refined Cd, Zn, and Te ({approx}7 N) precursor materials, using a tellurium solvent method. Detectors with virtual Frisch grid configuration were fabricated using these crystals. I-V measurements revealed low leakage currents at room-temperature, {approx}11 nA for one such detector D1 and {approx}8 nA for another detector D2 at 1100 V. The spectroscopic performances of the two CZT virtual Frisch grid detectors have been evaluated and compared for high energy gamma ray detection. Detector D1 showed a well-resolved pulse-height spectrum with an energy resolution of {approx}1.6% for the 662 keV gamma rays. Detector D2 also showed a distinct 662 keV peak but with a broader pulse-height distribution. A digital biparametric correlation study of the depth of interaction and energy deposited by the 662 keV gamma rays was carried out. A different kind of correlation pattern from that observed normally for hole trapping was noticed in the case of detector D2. Correlation of results from thermally stimulated current measurement studies suggested that the anomalous biparametric correlation pattern was due to the trapping of holes but modified by the virtual Frisch grid effect. The results also suggested that the effect of electron trapping could not be ruled out either. Finally, a digital correction scheme was applied to recover the {sup 137}Cs spectrum from the effect of charge loss.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Spatial fractionation of the dose using neon and heavier ions: A Monte Carlo study

Purpose: This work explores a new radiation therapy approach which might trigger a renewed use of neon and heavier ions to treat cancers. These ions were shown to be extremely efficient in radioresistant tumor killing. Unfortunately, the efficient region also extends into the normal tissue in front of the tumor. The strategy the authors propose is to profit from the well-established sparing effect of thin spatially fractionated beams, so that the impact on normal tissues might be minimized while a high tumor control is achieved. The main goal of this work is to provide a proof of concept of this new approach. With that aim, a dosimetric study was carried out as a first step to evaluate the interest of further explorations of this avenue. Methods: The GATE/GEANT4 v.6.1 Monte Carlo simulation platform was employed to simulate arrays of rectangular minibeams (700 μm × 2 cm) of four ions (Ne, Si, Ar, and Fe). The irradiations were performed with a 2 cm-long spread-out Bragg peak centered at 7 cm-depth. Dose distributions in a water phantom were scored considering two minibeams center-to-center distances: 1400 and 3500 μm. Peak and valley doses, peak-to-valley dose ratios (PVDRs), beam penumbras, and relative contribution of nuclear fragments and electromagnetic processes were assessed as figures of merit. In addition, the type and proportion of the secondary nuclear fragments were evaluated in both peak and valley regions. Results: Extremely high PVDR values (>100) and low valley doses were obtained. The higher the atomic number (Z) of the primary ion is, the lower the valleys and the narrower the penumbras. Although the yield of secondary nuclear products increases with Z, the actual dose being deposited by the secondary nuclear fragments in the valleys starts to be the dominant contribution at deeper points, helping in the sparing of proximal normal tissues. Additionally, a wider center-to-center distance leads to a minimized contribution of heavier secondary fragments in valleys. Conclusions: The computed dose distributions suggest that a spatial fractionation of the dose combined to the use of submillimetric field sizes might allow profiting from the high efficiency of neon and heavier ions for the treatment of radioresistant tumors, while preserving normal tissues. The authors’ results support the further exploration of this avenue. Next steps include the realization of biological experiment to confirm the shifting of normal tissue complication probability curves.

60 APPLIED LIFE SCIENCES↗

Virginia Solar Pathways Project: Economic Study of Utility-Administered Solar Programs: Soft Costs, Community Solar, and Tax Normalization Considerations

This report presents economic considerations for solar development in support of the Virginia Solar Pathways Project (VSPP), an effort funded by the U.S. Department of Energy (DOE) SunShot Initiative that seeks to develop a collaborative utility-administered solar strategy for the Commonwealth of Virginia. The results presented are intended to be considered alongside the results of other studies conducted under the VSPP that evaluate the impacts of solar energy on the electric distribution, transmission, and generation systems in Virginia.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗