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,441 records · Page 80

A Novel Decentralized PWM Interleaving Technique for Ripple Minimization in Series-Stacked DC-DC Converters

Cascaded dc-dc converters are commonly used in applications where distributed energy sources or loads are connected to elevated voltage levels for power transfer. In such systems, it is advantageous to minimize the ripple on the bus current and voltage by proper phase shifting of the pulse-width modulation (PWM) pulses among the converters via a method known as interleaving. Existing approaches use either a centralized controller or separate communication lines among the stacked converters to control their relative PWM switch transitions. The key drawbacks are that these methods entail significant wiring, the central controller acts as a single point of failure, and implementation on very large numbers of units is impractical. In this paper, we introduce a decentralized interleaving control (DIC) strategy that acts on local current measurements at every converter and achieves communication-free PWM interleaving among the series-stacked converters. The proposed controller is simple in structure and is shown to converge asymptotically to the interleaved state irrespective of clock drifts among the digital signal processors. Experimental results are provided for a system of five series-connected converters showing a 10× reduction in the current ripple compared to normal operation.

cascaded converter↗

Natural Gas Community of the Future (Final Report)

Nicor Gas, a subsidiary of Southern Company Gas and its parent company Southern Company, is developing the concept of the Carbon-neutral, Affordable, Resilient, and Equitable Community ("CARE Community") to address the challenges of accessibility and affordability of grid-interactive efficient buildings in historically disadvantaged communities. The CARE Community is a new construction net-zero energy and carbon neutral residential neighborhood consisting of 50 homes with electricity and natural gas services in South Suburban Chicago, built as low-income affordable housing. The National Renewable Energy Laboratory (NREL) is assisting Southern Company to explore how energy efficiency, renewable technologies, and the reliability and resilience of natural gas can be combined and optimized in an islandable energy system as part of a broader effort to understand the role natural gas plays in future communities integrated with renewable energy and improved energy efficiency. The CARE Community will demonstrate how energy efficiency, distributed energy resources (DERs), and advanced controls can be combined with existing natural gas infrastructure to serve historically disadvantaged communities in a cold climate with low energy cost for affordable housing programs and enhanced resilience to withstand extreme weather. Homes will be a part of an islandable energy system, integrating multiple DER technologies such as solar panels, natural gas fuel cells, reciprocating natural gas engines, and battery storage. The CARE Community will provide a replicable template that affordable housing developers and energy companies can adopt throughout the United States. In this report, we present the design and modeling of the community and potential approaches for integrating DERs. We analyze and discuss annual simulations under normal operating conditions and week-long simulations under electric power grid outage conditions.

03 NATURAL GAS↗

Modeling the Effects of Artificial Drainage on Agriculture-dominated Watersheds using a Fully Distributed Integrated Hydrology Model: Datasets, scripts, model files

This model-data archive supports the research paper that demonstrates the integration of agricultural drainage features—specifically, narrow engineered ditches and tile drains—into a fully distributed, basin-scale integrated surface-subsurface hydrology model (ISSHM), Amanzi-ATS. The model employs innovative computational meshes aligned with agricultural ditches and incorporates the physically based Hooghoudt's drainage equation to simulate tile drainage, offering a novel strategy that enhances the accuracy of hydrological simulations.The archived dataset includes input parameters, model configurations, and select simulation outputs for the Amanzi-ATS model that successfully captured the streamflow patterns in the Portage River Watershed as validated by USGS gauge readings. Jupyter notebook for the preparation of model inputs and post-processing of outputs are also included. The model's predictive performance achieved a normalized Kling-Gupta Efficiency (KGE) of 0.81, surpassing SWAT without the necessity for site-specific calibration.The Amanzi-ATS model presented in this modeL-data archive allows for numerical experiments to explore the shifts in the flow structure under different drainage scenarios. As a tool for advancing the understanding of distributed hydrological responses and nutrient cycling, this archived model provides valuable insights for researchers, modelers, and decision-makers involved in watershed management and environmental modeling.The Watershed Workflow package is implemented in Python3. The Jupyter notebooks can be executed through multiple open-source tools, for example, Anaconda Jupyter Lab, VS Studio Code, etc. Other data files include CSV and HDF5 files, which can be read through Python scripts. The input files for the ATS model, open-source integrated hydrology, and transport model, are in XML format and can be edited in any commonly used text editors.

54 ENVIRONMENTAL SCIENCES↗

An innovative and alternative approach toward gear fabrication

In this paper, we present the first demonstration of single-step gear manufacturing via an innovative in-situ friction stir forging approach. This novel process is a natural extension to the friction stir processing technique and relies on shear and normal stresses to form the complex shapes. Friction stir welding pin tool with shoulder is plunged at high revolutions per minute (rpm) into the sample to heat it to the desired temperature and enable plastic flow. Then, the tool is forged at high plunge rate to form the desired structurally-sound and defect free complex shapes. To better visualize the process the experiments are complemented by a 3D thermomechanically coupled, smoothed particle hydrodynamics (SPH) model to understand the material flow pattern both radially and through thickness in the gear teeth, plastic strain distribution, and temperature profile. We have demonstrated application of I-FSF to inherently poor formable AZ31 Mg, AA 7075-T6, and AA5083–10 vol% TiB 2 composite and fabricated 14-tooth spur gears from these materials with minimal farther machining. In conclusion, microstructure of the forged gears exhibited a directional material flow pattern due to the normal forging action and refinement of grains to 6–13 μm through dynamically recrystallization across the radial, through-thickness direction, which is in good agreement with the corresponding SPH simulation results.

42 ENGINEERING↗

ARE SOME MILKY WAY GLOBULAR CLUSTERS HOSTED BY UNDISCOVERED GALAXIES?

The confirmation of a globular cluster (GC) in the recently discovered ultrafaint galaxy Eridanus II (Eri II) motivated us to examine the question posed in the title. After estimating the halo mass of Eri II using a published stellar mass—halo mass relation, the one GC in this galaxy supports extending the relationship between the number of GCs hosted by a galaxy and the galaxy’s total mass about two orders of magnitude in stellar mass below the previous limit. For this empirically determined specific frequency of between 0.06 and 0.39 GCs per 10{sup 9} M {sub ⊙} of total mass, the surviving Milky Way (MW) subhalos with masses smaller than 10{sup 10} M {sub ⊙} could host as many as 5–31 GCs, broadly consistent with the actual population of outer halo MW GCs, although matching the radial distribution in detail remains a challenge. Using a subhalo mass function from published high-resolution numerical simulations and a Poissonian model for populating those halos with the aforementioned empirically constrained frequency, we find that about 90% of these GCs lie in lower-mass subhalos than that of Eri II. From what we know about the stellar mass–halo mass function, the subhalo mass function, and the mass-normalized GC specific frequency, we conclude that some of the MW’s outer halo GCs are likely to be hosted by undetected subhalos with extremely modest stellar populations.

79 ASTRONOMY AND ASTROPHYSICS↗

Modeling the Effects of Artificial Drainage on Agriculture-Dominated Watersheds Using a Fully Distributed Integrated Hydrology Model

In agriculture-dominated watersheds where natural drainage is poor, agricultural ditches (narrow engineered channels) and tile drains (perforated pipes) are widely employed to enhance surface and subsurface drainage, respectively. Despite their relatively small scale, these features exert substantial control over the hydro-biogeochemical function of watersheds and their effects need to be represented in the models. We introduce a novel strategy to incorporate the effects of artificial agricultural drainage into a fully distributed basin-scale integrated surface-subsurface hydrology models. In our approach, narrow agriculture ditches for surface drainage are resolved efficiently using ditch-aligned computational meshes that are hydrologically conditioned to ensure connectivity in the stream/ditch network. For tile drainage in the subsurface, we use the physically based Hooghoudt's drainage equation as a subgrid model and route the water drained through tiles to the nearest ditch. Without site-specific calibration, this model reproduced observed streamflow in the Portage River Watershed (>1,000 km 2 ) as recorded by a USGS gauge with good accuracy (normalized KGE = 0.81) and outperformed a calibrated SWAT model (normalized KGE = 0.68). Numerical experiments confirm that artificial drainage reduces surface inundations and effectively controls the water table. At the watershed scale, artificial drainage increases baseflow but has little effect on watershed discharges above the 90th percentile. The strong physical underpinnings and reduced need for calibration allow us to study the impacts of artificial drainage on distributed hydrological response in terms of fluxes and states and provide a platform for investigating watershed-scale nutrient transport.

54 ENVIRONMENTAL SCIENCES↗

TH-A-9A-03: Dosimetric Effect of Rotational Errors for Lung Stereotactic Body Radiotherapy

Purpose: To evaluate the dosimetric effects on target volume and organs at risk (OARs) due to roll rotational errors in treatment setup of stereotactic body radiation therapy (SBRT) for lung cancer. Methods: There were a total of 23 volumetric modulated arc therapy (VMAT) plans for lung SBRT examined in this retrospective study. Each CT image of VMAT plans was intentionally rotated by ±1°, ±2°, and ±3° to simulate roll rotational setup errors. The axis of rotation was set at the center of T-spine. The target volume and OARs in the rotated CT images were re-defined by deformable registration of original contours. The dose distributions on each set of rotated images were re-calculated to cover the planning target volume (PTV) with the prescription dose before and after the couch translational correction. The dose-volumetric changes of PTVs and spinal cords were analyzed. Results: The differences in D95% of PTVs by −3°, −2°, −1°, 1°, 2°, and 3° roll rotations before the couch translational correction were on average −11.3±11.4%, −5.46±7.24%, −1.11±1.38% −3.34±3.97%, −9.64±10.3%, and −16.3±14.7%, respectively. After the couch translational correction, those values were −0.195±0.544%, −0.159±0.391%, −0.188±0.262%, −0.310±0.270%, −0.407±0.331%, and −0.433±0.401%, respectively. The maximum dose difference of spinal cord among the 23 plans even after the couch translational correction was 25.9% at −3° rotation. Conclusions: Roll rotational setup errors in lung SBRT significantly influenced the coverage of target volume using VMAT technique. This could be in part compensated by the translational couch correction. However, in spite of the translational correction, the delivered doses to the spinal cord could be more than the calculated doses. Therefore if rotational setup errors exist during lung SBRT using VMAT technique, the rotational correction would rather be considered to prevent over-irradiation of normal tissues than the translational correction.

60 APPLIED LIFE SCIENCES↗

SU-E-T-73: Commissioning of a Treatment Planning System for Proton Spot Scanning

Purpose: A treatment planning system (TPS) was commissioned for clinical use with a fixed beam line proton delivery system. An outline of the data collection, modeling, and verification is provided. Methods: Beam data modeling for proton spot scanning in CMS Xio TPS requires the following measurements: (i) integral depth dose curves (IDDCs); (ii) absolute dose calibration; and (iii) beam spot characteristics. The IDDCs for 18 proton energies were measured using an integrating detector in a single spot field in a water phantom. Absolute scaling of the IDDCs were performed based on ion chamber measurements in mono-energetic 10×10 cm{sup 2} fields in water. Beam spot shapes were measured in air using a flat panel scintillator detector at multiple planes. For beam model verification, more than 45 uniform dose phantom and patient plans were generated. These plans were used to measure range, point dose, and longitudinal and lateral profiles. Tolerances employed for verification are: point dose and longitudinal profiles, ±2%; range, ±1 mm; FWHM for lateral profiles, ±2 mm; and patient plan dose distribution, gamma index of >90% at 3%/3 mm criteria. Results: More than 97% of the point dose measurements out of 115 were within +/-2% with maximum deviation of 3%. 98% of the ranges measured were within 1 mm with maximum deviation of 1.4mm. The normalized depth doses were within 2% at all depths. The maximum error in FWHM of lateral profiles was found to be less than 2mm. For 5 patient plans representing different anatomic sites, a total of 38 planes for 12 beams were analyzed for gamma index with average value of 99% and minimum of 94%. Conclusions: The planning system is successfully commissioned and can be safely deployed for clinical use. Measurements of IDDCs on user beam are highly recommended instead of using standard beam IDDCs.

60 APPLIED LIFE SCIENCES↗

Time domain-induced polarization geophysical data collected in the Rifle Floodplain, Colorado, USA

The IP data set was collected using a time domain-induced polarization (TDIP) method along 65 profiles of various lengths over the floodplain. Surface time domain?induced polarization (TDIP) data was collected to create the 3?D images of the complex electrical resistivity, in terms of magnitude and phase, which are associated with mineral precipitation and other lithological properties. The TDIP data was then used to estimate the spatial distribution of naturally reduced zones, which are known to be biogeochemical hotspots in floodplains having increased uranium and other redox-sensitive metal concentrations as well as high carbon contents. In the TDIP method, the transient decay of voltage is measured after current shut-off, typically in the form of an integral of decay curves over a predefined time window (so-called integral chargeability). TDIP measurements at the site were collected using the Syscal Iris Pro Switch equipment with a square-wave current injection, 50% duty cycle, and a pulse length of 2 s. The integral chargeability measurements were carried out using 20 windows during voltage decay between 240 and 1840 ms after current shut-off. Tomographic measurements were collected by deploying stainless steel electrodes with an electrode separation of 1.8 m and using a dipole-dipole ‘‘skip-2’’ and ‘‘skip-3’’ measuring protocol (i.e., for a dipoles length of 5.4 and 7.2 m, respectively). The sequence of dipole-dipole measurements was carefully arranged to (1) minimize unwanted electromagnetic coupling effects in the data, avoiding potential measurements with electrodes located inside the current dipole, (2) prevent voltage measurements using electrodes, which might be polarized due to previous current injection, and (3) increase the signal-to-noise ratio for an intended exploration depth of 8 m, i.e., the bottom of the aquifer. All measurements were collected as normal and reciprocal pairs for estimation of the data error. The IP measurements were collected with symmetric arrays (i.e., the measuring equipment placed at the center of the electrode array) with a maximum of 36 electrodes, considering that longer profiles revealed a significant increase in the normal-reciprocal misfit for the measurements of the decay curve, probably due to greater impact of electromagnetic coupling on the data. The TDIP data sets were inverted in a twodimensional domain along each transect using CRTomo, which is a smoothness-constraint inversion code based on a finite element algorithm. The resistivity and phase shift values at each pixel were then assigned at the corresponding point within the 3-D floodplain domain (the black rectangle Figure 1b) and used in the 3-D estimation. The TDIP inversion results provided the distribution of the complex resistivity, expressed in terms of its magnitude and phase-shift. The analysis of the normal-reciprocal misfit was used to estimate the relative error. We removed the measurement with smallest voltage difference (<2 mV), representing about 2% variations in the data. For the inversion of TDIP measurements, chargeability values were linearly converted to frequency domain phase values (at the fundamental frequency of 0.125 Hz), by assuming a constant-phase response. This approach has been demonstrated to provide consistent results in previous studies

54 ENVIRONMENTAL SCIENCES↗

Determination of the top-quark pole mass using tt¯ + 1-jet events collected with the ATLAS experiment in 7TeV pp collisions

In this study, the normalized differential cross section for top-quark pair production in association with at least one jet is studied as a function of the inverse of the invariant mass of the tt¯ + 1-jet system. This distribution can be used for a precise determination of the top-quark mass since gluon radiation depends on the mass of the quarks. The experimental analysis is based on proton-proton collision data collected by the ATLAS detector at the LHC with a centre-of-mass energy of 7 TeV corresponding to an integrated luminosity of 4.6 fb –1 . The selected events were identified using the lepton+jets top-quark-pair decay channel, where lepton refers to either an electron or a muon. The observed distribution is compared to a theoretical prediction at next-to-leading-order accuracy in quantum chromodynamics using the pole-mass scheme. With this method, the measured value of the top-quark pole mass, m pole t , is: m pole t = 173.7 ± 1.5(stat.) ± 1.4(syst.) +1.0 –0.5 (theory) GeV.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Very low 1/f barrier noise in sputtered MgO magnetic tunnel junctions with high tunneling magnetoresistance

Low frequency 1/f barrier noise has been investigated in sputtered MgO magnetic tunnel junctions (MTJs) with a tunneling magnetoresistance ratio of up to 330% at room temperature. The lowest normalized noise parameter {alpha} of the tunnel barrier reaches 2.5 Multiplication-Sign 10{sup -12}-2.1 Multiplication-Sign 10{sup -11} {mu}m{sup 2}, which is comparable to that found in MTJs with the MgO barrier grown by MBE or electron-beam evaporation. This normalized barrier noise is almost bias independent in the voltage range of up to {+-}1.2 V. The low noise level and high voltage stability may reflect the high quality of the sputtered MgO with a uniform distribution of defects in the MgO layer.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Short-lived species detection of nitrous acid by external-cavity quantum cascade laser based quartz-enhanced photoacoustic absorption spectroscopy

Spectroscopic detection of short-lived gaseous nitrous acid (HONO) at 1254.85 cm{sup −1} was realized by off-beam coupled quartz-enhanced photoacoustic spectroscopy (QEPAS) in conjunction with an external cavity quantum cascade lasers (EC-QCL). High sensitivity monitoring of HONO was performed within a very small gas-sample volume (of ∼40 mm{sup 3}) allowing a significant reduction (of about 4 orders of magnitude) of air sampling residence time which is highly desired for accurate quantification of chemically reactive short-lived species. Calibration of the developed QEPAS-based HONO sensor was carried out by means of lab-generated HONO samples whose concentrations were determined by direct absorption spectroscopy involving a ∼109.5 m multipass cell and a distributed feedback QCL. A minimum detection limit (MDL) of 66 ppbv (1 σ) HONO was achieved at 70 mbar using a laser output power of 50 mW and 1 s integration time, which corresponded to a normalized noise equivalent absorption coefficient of 3.6 × 10{sup −8 }cm{sup −1} W/Hz{sup 1/2}. This MDL was down to 7 ppbv at the optimal integration time of 150 s. The corresponding 1σ minimum detected absorption coefficient is ∼1.1 × 10{sup −7 }cm{sup −1} (MDL ∼ 3 ppbv) in 1 s and ∼1.1 × 10{sup −8 }cm{sup −1} (MDL ∼ 330 pptv) in 150 s, respectively, with 1 W laser power.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

THE TWO REGIMES OF PHOTOSPHERIC MOTIONS IN {alpha} HYDRA

High-resolution spectroscopic observations of {alpha} Hya were acquired between 2003 and 2010. Analysis of line shifts, differential shifts, line widths, and line bisectors points to two regimes of velocity fields in the photosphere of {alpha} Hya: (1) normal granulation embedded in (2) large convection cells. Variations occur on a wide range of timescales, from several years on down. Radial velocity variations, which are irregular and span 786 m s{sup -1}, have a distribution consistent with a true mean rise velocity of the large cells of {approx}725 m s{sup -1} and a dispersion of {approx}220 m s{sup -1}. The distribution of granulation velocities, as measured from the widths of spectral lines, shows only small variations, consistent with the two regime concepts. On the multi-year timescale, radial velocity changes, small temperature variations ({approx}10 K), and small line-width variations ({approx}<0.8%) track each other, possibly with phase shifts. The granulation velocity gradient for {alpha} Hya is about half as large as the Sun's and no variation with time was seen, implying that any variation in velocity gradient from one large cell to the next must be less than a few percent. The asymmetry in the granulation velocity distribution, as specified in the flux deficit, is smaller than expected for {alpha} Hya's position in the HR diagram and appears to be variable.

79 ASTRONOMY AND ASTROPHYSICS↗

Simultaneous inference of plate boundary stresses and mantle rheology using adjoints: large-scale 2-D models

SUMMARY Plate motions are a primary surface constraint on plate and mantle dynamics and rheology, plate boundary stresses and the occurrence of great earthquakes. Within an optimization method, we use plate motion data to better constrain uncertain mantle parameters. For the optimization problem characterizing the maximum a posteriori rheological parameters we derive gradients using adjoints and expressions to approximate the posterior distributions for stresses within plate boundaries. We apply these methods to a 2-D cross section from the western to eastern Pacific, with temperature distributions and fault zone geometries developed primarily from seismic and plate motion data. We find that the best-fitting stress exponent, n, is about 2.8 and the yield stress about 100 MPa or less. The normal stress on the interplate fault zones is about 100 MPa and the shear stresses about 10 MPa or less.

Rudi, Johann (ORCID:0000000265639265)↗

Gluon pseudo-distributions at short distances

We present the results that are necessary in the ongoing lattice calculations of the gluon parton distribution functions (PDFs) within the pseudo-PDF approach. We identify the two-gluon correlator functions that contain the invariant amplitude determining the gluon PDF in the light-cone z^2 \to 0 z 2 → 0 limit, and perform one-loop calculations in the coordinate representation in an explicitly gauge-invariant form. Ultraviolet (UV) terms, which contain \ln (-z^2) ln ( − z 2 ) -dependence cancel in the reduced Ioffe-time distribution (ITD), and we obtain the matching relation between the reduced ITD and the light-cone ITD. Using a kernel form, we get a direct connection between lattice data for the reduced ITD and the normalized gluon PDF.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Sensitivity and Uncertainty of the IFR-1 BISON Benchmark

The fuel performance code BISON is being used to evaluate metallic fuel for a new fast-spectrum test reactor called the Versatile Test Reactor, which is being considered by the US Department of Energy. To quantify the accuracy of BISON predictions, researchers at Oak Ridge National Laboratory have been developing a series of benchmarks based on legacy metallic fuel experiments. As part of this effort, the sensitivity of BISON predictions to variations in model inputs and the uncertainties associated with BISON predictions must be established. This report summarizes efforts to perform a comprehensive sensitivity analysis (SA) and uncertainty quantification (UQ) on a benchmark based on the IFR-1 experiment. For the SA, at least one input was chosen from every BISON model and physics module used in the benchmark. The inputs were varied individually in a series of BISON simulations. The resulting variations in benchmark predictions were normalized to calculate sensitivities. The strongest sensitivities were identified and used to inform input selections for the UQ. The UQ was performed using the Monte Carlo UQ method. A literature review was conducted to estimate uncertainty distributions for the selected inputs, and values were sampled randomly from each distribution in a series of BISON simulations. Variations in the benchmark predictions were used to estimate uncertainty distributions and confidence intervals. It was found that nearly 100% of benchmark predictions matched the corresponding legacy values within the confidence intervals. However, this is at least partially because of the wide confidence intervals associated with the benchmark predictions. The uncertainty contributions of assumptions in the benchmark, experimental uncertainties, and BISON models were quantified. Some analysis was performed to identify inputs that contributed to the uncertainties. Finally, recommendations are made for future benchmark development and future BISON development.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Estimating Total Methane Emissions from the Denver-Julesburg Basin Using Bottom-Up Approaches

Methane is a powerful greenhouse gas with a 25 times higher 100-year warming potential than carbon dioxide and is a target for mitigation to achieve climate goals. To control and curb methane emissions, estimates are required from the sources and sectors which are typically generated using bottom-up methods. However, recent studies have shown that national and international bottom-up approaches can significantly underestimate emissions. In this study, we present three bottom-up approaches used to estimate methane emissions from all emission sectors in the Denver-Julesburg basin, CO, USA. Our data show emissions generated from all three methods are lower than historic measurements. A Tier 1/2 approach using IPCC emission factors estimated 2022 methane emissions of 358 Gg (0.8% of produced methane lost by the energy sector), while a Tier 3 EPA-based approach estimated emissions of 269 Gg (0.2%). Using emission factors informed by contemporary and region-specific measurement studies, emissions of 212 Gg (0.2%) were calculated. The largest difference in emissions estimates were a result of using the Mechanistic Air Emissions Simulator (MAES) for the production and transport of oil and gas in the DJ basin. The MAES accounts for changes to regulatory practice in the DJ basin, which include comprehensive requirements for compressors, pneumatics, equipment leaks, and fugitive emissions, which were implemented to reduce emissions starting in 2014. The measurement revealed that normalized gas loss is predicted to have been reduced by a factor of 20 when compared to 10-year-old normalization loss measurements and a factor of 10 less than a nearby oil and production area (Delaware basin, TX); however, we suggest that more measurements should be made to ensure that the long-tail emission distribution has been captured by the modeling. This study suggests that regulations implemented by the Colorado Department of Public Health and Environment could have reduced emissions by a factor of 20, but contemporary regional measurements should be made to ensure these bottom-up calculations are realistic.

03 NATURAL GAS↗