SEARCH · Search NASA

Results for “force control”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 469 records · Page 26

All sources NASA NTRS DOE OSTI

DOE OSTI2023

Voltage Support With PV Inverters in Low-Voltage Distribution Networks: An Overview

Large solar photovoltaic (PV) penetration using inverters in low voltage (LV) distribution networks may pose several challenges, such as reverse power flow and voltage rise situations. These challenges will eventually force grid operators to carry out grid reinforcement to ensure continued safe and reliable operations. However, smart inverters with reactive power control capability enable PV systems to support voltage quality in the distribution network better. Here, this paper gives an overview of the current state-of-the-art control strategies for handling voltage problems through PV inverters and other devices. In addition, the (control) technical issues of PV systems integrated into the LV distribution network are considered from a control point of view. By comparing the control issues of PV integration into the grid, the paper aims to help distribution system operators to expand the volume of PV generation in the distribution system in an efficient and safe manner. Additionally, it will help control engineers and researchers select proper control strategies for PV systems as well as other distributed renewable sources.

42 ENGINEERING↗

DOE OSTI2021

Improved Contacts for Tandem Cells with Enhanced Effciency Grown by D-HVPE

Dynamic hydride vapor phase epitaxy (D-HVPE) has provided a potential route towards lower cost material growth of III-V photovoltaic devices. A highly doped Ga0.5In0.5P emitter layer is grown in attempt to force passivation at the front of the cell to account for the absence of a suitable window layer in our reactor. Control over dopant diffusion in these materials is thus critical to achieving high-efficiency device performance. Here, we institute a two-step contact layer design that boosts tandem cell efficiencies by minimizing Se diffusion into the underlying emitter layer while still providing sufficiently low contact resistance. Short-circuit current in these devices improves to 8.9 mA/cm2 in an uncoated GaInP/GaAs cell, an improvement of 1.2 mA/cm2 over our previous best cell. Conversion efficiency with the new contact doping scheme is projected to be ~27% with the inclusion of a standard ARC layer.

27 ARPA - Advanced Research Projects Agency-Energy↗

DOE OSTI2022

7.2 kV Three-Port SiC Single-Stage Current-Source Solid-State Transformer With 90 kV Lightning Protection

This article proposes a multiport modular single-stage current-source solid-state transformer (SST) for applications like photovoltaic, energy storage integration, electric vehicle fast charging, data center, etc. The 7.2 kV 50 kVA current-source SST consists of five input-series output-parallel modules, each based on 3.3 kV SiC reverse-blocking MOSFET-plus-diode modules. The proposed SST has some unique features. First, compared to the voltage-source or matrix converter-based SSTs, the current-source SST has a unique advantage of single-stage isolated AC/DC or AC/AC conversion with an inductive DC link, but no medium-voltage (MV) AC experiments have been reported. This article for the first time demonstrates MV AC current-source SST up to 7.5 kV peak. Second, the multiport SST has a buffer port for active power decoupling (APD) or energy storage integration. The double-line-frequency power ripple from single-phase AC grid normally results in a large capacitor size in MV SSTs. The APD scheme is proposed in MV applications for the first time to enable a reduced DC link and the electrolytic capacitor-less SST with high reliability. Third, as a direct grid-connected converter without line-frequency transformer, insulation and protection are critical. A medium-frequency transformer design passes 55 kV basic-insulation level (BIL) and 60 kV high potential dielectrics withstand test with only 0.09% leakage inductance. Importantly, a lightning protection scheme is presented to protect the SST itself from 90 kV BIL impulse. Fourth, the proposed current-source SST topology is a modular soft-switching solid-state transformer (M-S4T) with full-range zero-voltage switching and controlled dv/dt for low electromagnetic interference. Furthermore, these concepts are verified in a three-port M-S4T prototype with forced oil cooling under single-module, stacked-module, steady-state, and dynamic operations.

14 SOLAR ENERGY↗

DOE OSTI2021

TEAMER: Biofouling Analysis for Wave Energy Piston Design - Load Cell Data

Biofouling and corrosion are a major concern for all ocean-deployed components, especially when mechanical motion is involved. Triton has developed the concept of a biofouling mitigation seal as part of the piston sealing assembly for the Triton Wave Energy Converter (TSI-WEC). This mitigation seal has the purpose of preventing the formation of a biofilm on the inside of the piston cylinder. It is hypothesized that the prevention of a biofilm will reduce the amount of macro-biofouling that can occur in the piston assembly. The mitigation seal can also reduce the wear on the main dynamic seal, helping to maintain smooth operation and water-tightness. The cylinder is made from a thermoset composite epoxy, which is resistant to corrosion. However, no studies have researched the material's performance with biofouling. Triton placed two prototype Power Take-Off (PTO) assemblies in a PNNL biofouling tank, one with a biofouling mitigation seal and one without, allowing for an evaluation of seal effectiveness at the prevention of biofouling. In actual WEC operation, wave action would react against the piston, which would drive the linear actuator and electric generator, providing electrical power. In the test setup, this was reversed; a linear actuator was powered to drive the piston in a consistent motion within the cylinder. There are two assemblies: one has a biofouling mitigation seal, the other (control) does not. The following data encompasses a 4 month test period, with load cells being used to monitor piston friction force. Results from this testing will be used to improve seal design and material selection, mitigating risk of premature failure during open water testing and evaluation. This project is part of the TEAMER RFTS 3 (request for technical support) program.

16 TIDAL AND WAVE POWER↗

DOE OSTI2022

kynema-fmb [SWR-23-07]

Kynema-FMB (FKA: Kynema) is an open-source performance portable flexible multibody (FMB) dynamics solver designed for time-domain simulations. While originally tailored for wind turbine structural dynamics, the formulation and implementation are those of a general flexible-multidbody dynamics solver that can readily be applied to a wide range of systems. Kynema was designed with a narrow focus, namely to provide a lightweight, fast, accurate FMD solver for coupling to computational-fluid-dynamics (CFD) codes, especially the CFD codes in the Kynema suite, for fluid-structure-interaction (FSI) simulations. Kynema-FMB is equipped to model systems that can be represented as a collection of beams and rigid bodies that are connected through constraints. Degrees of freedom are defined in the inertial/global frame of reference and include displacements and rotations (formally as rotation matrices, but stored as quaternions). The underlying formulation is built on a Lie-group time integrator designed for index-3 differential-algebraic equations, which is second-order accurate in time (Bruls et al., 2012). Beam models are based on geometrically exact beam theory and are discretized as high-order spectral finite elements similar to those in BeamDyn (Wang et al., 2017). The governing equations for a FMD system like a wind turbine constitute a highly nonlinear system of constrained partial-differential equations. Kynema-FMB uses analytical Jacobians in the nonlinear-system solves in each time step. Linear systems use sparse storage and several third-party sparse-linear-system solvers are enabled. Ill conditioning of linear systems is mitigated with preconditioning described in Bottasso et al, 2008. Kynema-FMB is integrated with a simple open-source controller (ROSCO). There is an application programming interface (API) for coupling to geometry-resolved CFD (like that in Sharma et al., 2023) and actuator-force CFD (like that in Kuhn et al., 2025). In the latter, for actuator-line models, Kynema-FMB includes an internal blade-element solver that depends on user-provided lookup tables for coefficients of lift and drag, i.e., aerodynamic polars. Kynema-FMB is written in C++ and leverages Kokkos and Kokkos-Kernels (KokkosEcosystem) as its performance portability layer enabling simulations on both CPU and GPU systems. The repository is equipped with extensive automated testing at the unit and regression/system levels. The following describes the high-level development objectives conceived for Kynema: *Kynema will follow modern software development best practices, including test-driven development (TDD), version control, hierarchical automated testing, and continuous integration (CI) for a robust development environment. *The core data structures are memory efficient and enable vectorization and parallelization at multiple levels. *Data structures are data-oriented to exploit methods for accelerated computing including high utilization of chip resources (e.g., single instruction multiple data (SIMD) instruction sets) and parallelization using GP-GPUs. *The computational algorithms incorporate robust open-source libraries for mathematical operations, resource allocation, and data management. *The API design considers multiple stakeholder needs and ensure integration with existing and future ecosystems for data science, machine learning, and AI. *Kynema-FMB is written in modern C++ and leverages Kokkos as its performance-portability library with inspiration from the kynema stack.

Sprague, MichaelA.↗

DOE OSTI2021

Optical Forces on an Oscillating Dipole Near VO2 Phase Transition

We investigate optical forces on oscillating dipoles close to a phase change vanadium dioxide (VO2) film, which exhibits a metal-insulator transition around 340 K and low thermal hysteresis. This configuration emulates the interaction between an illuminated nanosphere and an interface and we employ a classical description to capture its important aspects. We consider both electric and magnetic dipoles for two different configurations, namely with the dipole moments parallel and perpendicular to the VO2 film. By using Bruggeman theory to describe the effective optical response of the material, we show that the thermal hysteresis present in the VO2 transition clearly shows up in the behavior of optical forces. In the near-field regime, the force on both dipoles can change from attractive to repulsive just by heating (or cooling) the film for a selected frequency range. We also verified that the optical forces are comparable to the Casimir-Polder force in a similar system, revealing the possibility of modulating or even changing the sign of the resultant force on an illuminated nano-object due to the presence of a thermochromic material. We hope that this work contributes to set the grounds for alternative approaches to control light-matter interactions using phase-change materials.

36 MATERIALS SCIENCE↗

DOE OSTI2025

Steric Modulation of Protein‐Mediated Nanoparticle Assembly: Controlling Cluster Size, Polydispersity, and FRET Responses by Rebalancing Short‐ and Long‐Range Interactions

Understanding and manipulating protein-nanoparticle interactions is of broad interest to fields ranging from nanomedicine to the biological fabrication of functional hierarchical materials. This study investigates how steric forces introduced by a pegylated derivative of superfolder green fluorescent protein (sfGFP) that is monofunctional for silica binding modulate the delicate interplay of long-range (electrostatic and van der Waals) and short-range (protein-mediated) interactions in pH-responsive silica nanoparticle (SiNP) assembly by bifunctional silica-binding sfGFP. Increasing the length of the PEG segment and pre-incubating SiNPs with increasing concentrations of pegylated proteins enables precise control over cluster size within the 800–1450 nm range with a sixfold decrease in polydispersity index to a remarkable 0.1 endpoint. Weakening short-range attractive interactions via mutagenesis extends this control to clusters in the 50–250 nm range and reveals that the Förster resonance energy transfer (FRET) efficiency of clusters scales linearly with cluster diameter below 230 nm but increases only by 15% as clusters grow to 1450 nm. Furthermore, these findings enable the development of a system that provides an optical readout to dynamic changes in solution conditions enacted by a combination of pH adjustment and ion charge screening.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

DOE OSTI2025

Control and synchronization of rapid nanoscale DNA heat engine by local heating

To further activate devices based on DNA nanotechnology, we introduce an approach that notably enhances both the speed and force of DNA powered machines and artificial hinge machine. A microheater, with millisecond response, heats or recools DNA origami constructs, hybridizing or dehybridizing sticky ends. Because anything within 20 micrometers of the heater equilibrates to a programmed temperature change in milliseconds, sticky ends of a compound DNA origami machine can open and close synchronously and operate cooperatively, in phase, additively increasing the drive force compared to single pair of sticky ends DNA machine (the six-helix bundle DNA origami hinge machine). In our demonstrations, we fold and unfold two square origami with 10 pairs of complementary sticky ends to drive a bead on the end of a rod like origami to speeds exceeding 30 micrometers per second. Our device envisions the creation of complex, synchronized DNA machines.

Science & Technology - Other Topics↗

DOE OSTI2021

Texture Formation in Polycrystalline Thin Films of All‐Inorganic Lead Halide Perovskite

Abstract Controlling grain orientations within polycrystalline all‐inorganic halide perovskite solar cells can help increase conversion efficiencies toward their thermodynamic limits; however, the forces governing texture formation are ambiguous. Using synchrotron X‐ray diffraction, mesostructure formation within polycrystalline CsPbI 2.85 Br 0.15 powders as they cool from a high‐temperature cubic perovskite (α‐phase) is reported. Tetragonal distortions (β‐phase) trigger preferential crystallographic alignment within polycrystalline ensembles, a feature that is suggested here to be coordinated across multiple neighboring grains via interfacial forces that select for certain lattice distortions over others. External anisotropy is then imposed on polycrystalline thin films of orthorhombic (γ‐phase) CsPbI 3‐ x Br x perovskite via substrate clamping, revealing two fundamental uniaxial texture formations; i) I‐rich films possess orthorhombic‐like texture (<100> out‐of‐plane; <010> and <001> in‐plane), while ii) Br‐rich films form tetragonal‐like texture (<110> out‐of‐plane; <110> and <001> in‐plane). In contrast to relatively uninfluential factors like the choice of substrate, film thickness, and annealing temperature, Br incorporation modifies the γ‐CsPbI 3− x Br x crystal structure by reducing the orthorhombic lattice distortion (making it more tetragonal‐like) and governs the formation of the different, energetically favored textures within polycrystalline thin films.

Steele, Julian A.↗

DOE OSTI2024

Performance characterization and modeling of an oscillating surge wave energy converter

Abstract Testing wave energy converters in the ocean could be expensive and complex, which necessitates the use of numerical modeling. However, accurately modeling the response of wave energy converters with high-fidelity simulations can be computationally intensive in the design stage where different configurations must be considered. Reduced-order models based on simplified equations of motion can be very useful in the design, optimization, or control of wave energy converters. Given the complex dynamics of wave energy converters, accurate representation, and evaluation of relative contributions by different forces are required. This effort is concerned with a performance characterization of the hydrodynamic response of an oscillating surge wave energy converter that is based on a reduced-order model. A state-space model is used to represent the radiation damping term. Morison’s representation of unsteady forces is used to account for the nonlinear damping. Wave tank tests are performed to validate simulations. A free response simulation is used to determine the coefficients of the state-space model. Torque-forced simulations are used to identify the coefficients of the nonlinear damping term for different amplitudes and wave frequencies. The impact of varying these coefficients on the response is investigated. An assessment of the capability of the model in predicting the hydrodynamic response under irregular forcing is performed. The results show that the maximum error is 3% when compared with high-fidelity simulations. It is determined that the nonlinear damping is proportional to the torque amplitude and its effects are more pronounced as the amplitude of the flap oscillations increases.

16 TIDAL AND WAVE POWER↗

DOE OSTI2022

Dual-scale folding in cutting of commercially pure aluminum alloys

We examine a, hitherto, little-studied and curious machining chip morphology, with tell-tale signs of folding at two different length scales, that is common in cutting of certain ductile and highly strain-hardening metals like soft aluminum alloys, tantalum and niobium. This chip morphology does not appear in the usual catalogues of common chip types. The mechanics of formation of the “dual-scale folded chip” is studied in model material systems of commercially pure aluminum alloys (AA 1100 and AA 8040), that prominently exhibit this chip morphology. The flow, folding and associated plastic instabilities are investigated using micro/macro structure observations of the chip in a plane-strain cutting framework, with high-speed in situ imaging and image analysis of material flow; and force measurements. The smaller-scale folding is shown to develop in the primary deformation zone while the larger-scale folding occurs as the chip traverses the rake face of the tool. The resulting chip is composed of irregularly-spaced large folds, superimposed onto which are the quasi periodic small folds. The representative wavelengths of the two folds differ on average by an order of magnitude, 0.1 mm vs. 2 mm. The observations reveal a direct coupling between the material flow and chip morphology, and how specific attributes of the dual-scale folded chip arise from the flow mechanism. Plastic buckling is found to play a key role in the folding at both length scales. The small-scale folds are characteristic of a sinuous plastic flow mode, while the large-scale folding is characterized by buckling and stick-slip along the tool rake face, triggered by adhesive pinning of the chip to the tool. Important consequences of the dualscale folding are very large cutting forces, and force oscillations of large amplitude, despite the alloys being very soft, only ~ 25 HV. Here, the dual-scale folding is why many of these alloys are classified as “gummy” to machine. Since the dual-scale folded chip is associated with large cutting forces and poor surface quality, there is much to be gained by disrupting this flow type in practical machining applications. Methods for controlling the folding to improve machining performance with the gummy alloys are briefly discussed.

36 MATERIALS SCIENCE↗

DOE OSTI2022

Hybrid data-driven and model-informed online tool wear detection in milling machines

Precision machining tool wear is responsible for low product throughput and quality. Monitoring the tool wear online is vital to prevent degradation in machining quality. However, direct real-time tool wear measurement is not practical. This paper presents residual-based anomaly detection models, combining a hybrid model comprised of a physics-based model and a data-driven model (a decision tree or a neural network) to predict signals of interest (e.g., power or forces) under nominal conditions, followed by Page’s cumulative sum test for detecting tool wear on-line using the computer numerical control machine measurements. The most informative features are ranked using dynamic programming and its approximation variants from real-time measurements and machine settings, such as the width of cut, depth of cut, feed rate and spindle speed, that serve as inputs to the predictive models. The baseline nominal model is incrementally updated with experimental data via a gradient boosted adaptation model to generate the residuals that account for discrepancies between the actual machine data under normal conditions and the baseline nominal model predictions. The hybrid model is validated against 20 Mazak milling machine experimental tests and one Haas run-to-failure experiment. The proposed anomaly detector is applied to synthetic data from simulations of the physics-based model at different operating conditions, measurement noise levels, and tool wear levels, and the methods were able to achieve an overall 92% accuracy in data with 1% noise. The anomaly detection methods based on hybrid model reduced the false alarms of either the data-driven or physical-based models alone, and are found to be capable of good online detection of tool wear.

42 ENGINEERING↗

DOE OSTI2024

Aerosol-Cloud Interactions From Aviation Soot Emissions

Current models estimate global aviation contributes approximately 5% to the total anthropogenic climate forcing, with aerosol-cloud interactions having the greatest effect. However, radiative forcing estimates from aviation aerosol-cloud interactions remain undetermined. There is an expected significant increase in aircraft emissions with aviation demand expected to rise by over 4% per year. Soot may play an important role in the ice nucleation of aircraft-induced cirrus formation due to a high emission rate, but the ice nucleating properties are poorly constrained. Understanding the microphysical processes leading to atmospheric ice crystal formation is crucial for the reliable parameterization of aerosol-cloud interactions in climate models due to their impact on precipitation and cloud radiative properties. Ice nucleation of aircraft-emitted soot is potentially affected by particle morphology with condensation of supercooled water occurring in pores followed by ice nucleation. However, soot has heterogeneous properties and undergoes atmospheric aging and oxidation that could change surface properties and contribute to complex ice nucleation processes. Further, this review synthesizes current knowledge of ice nucleation catalyzed by aviation in the cirrus regime and its effects on global radiative forcing. Further research is required to determine the ice nucleation and microphysical processes of cirrus cloud formation from aviation emissions in both controlled laboratory and field investigations to inform models for more accurate climate predictions and to provide efficient mitigation strategies.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2020

Vaporizable endoskeletal droplets via tunable interfacial melting transitions

Liquid emulsion droplet evaporation is of importance for various sensing and imaging applications. The liquid-to-gas phase transformation is typically triggered thermally or acoustically by low–boiling point liquids, or by inclusion of solid structures that pin the vapor/liquid contact line to facilitate heterogeneous nucleation. However, these approaches lack precise tunability in vaporization behavior. Here, we describe a previously unused approach to control vaporization behavior through an endoskeleton that can melt and blend into the liquid core to either enhance or disrupt cohesive intermolecular forces. This effect is demonstrated using perfluoropentane (C 5 F 12 ) droplets encapsulating a fluorocarbon (FC) or hydrocarbon (HC) endoskeleton. FC skeletons inhibit vaporization, whereas HC skeletons trigger vaporization near the rotator melting transition. Our findings highlight the importance of skeletal interfacial mixing for initiating droplet vaporization. Tuning molecular interactions between the endoskeleton and droplet phase is generalizable for achieving emulsion or other secondary phase transitions, in emulsions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2022

Mitigation of Molten Salt Corrosion

The desire to increase the peak operational temperature of Concentrated Solar Power (CSP) plants to lower the levelized cost of CSP-derived electricity, has led to consideration of an earth-abundant, low-cost, thermally-stable, molten MgCl 2 -KCl-NaCl (40-40-20 mole%) salt for high-temperature thermal energy storage and heat transfer. The reliable, long-term containment of such a molten salt is challenging, however, due to the corrosion of cost-effective, high-temperature structural metal alloys in this salt (particularly if this salt is contaminated with oxygen-bearing species). This project has been focused on: i) development of equipment for evaluating corrosion kinetics in molten chlorides under well-controlled thermal, fluid dynamic, and atmospheric conditions at CSP-relevant temperatures, ii) use of such equipment to evaluate the influence of forced convection on corrosion kinetics (to determine the role of liquid phase diffusion on such corrosion), and iii) identifying a corrosion-mitigation strategy that can be industrially applied to achieve a projected annual corrosion ≤30 μm at CSP-relevant temperatures. The overall goals of this project have been: i) to understand the rate-limiting step(s) and associated kinetic laws governing corrosion in molten chlorides at 750°C, ii) to determine key operational parameters affecting the corrosion kinetics in molten salts, and iii) to demonstrate the tailoring of such parameters to appreciably reduce such corrosion.

14 SOLAR ENERGY↗

DOE OSTI2024

Optical and electronic functionality arising from controlled defect formation in nanoscale complex oxide lateral epitaxy

Epitaxial crystallization of complex oxides provides the means to create materials with precisely selected composition, strain, and orientation, thereby controlling their functionalities. Extending this control to nanoscale three-dimensional geometries can be accomplished via a three-dimensional analog of oxide solid-phase epitaxy, lateral epitaxial crystallization. The orientation of crystals within laterally crystallized SrTiO 3 systematically changes from the orientation of the SrTiO 3 substrate. This evolution occurs as a function of lateral crystallization distance, with a rate of approximately 50° μm -1 . The mechanism of the rotation is consistent with a steady-state stress of tens of megapascal over a 100–nanometer scale region near the moving amorphous/crystalline interface arising from the amorphous-crystalline density difference. Second harmonic generation and piezoelectric force microscopy reveal that the laterally crystallized SrTiO 3 is noncentrosymmetric and develops a switchable piezoelectric response at room temperature, illustrating the potential to use lateral crystallization to control the functionality of complex oxides.

36 MATERIALS SCIENCE↗

DOE OSTI2020

Business Case Analysis for Digital Safety-Related Instrumentation & Control System Modernizations

The commercial nuclear sector faces unprecedented financial challenges driven by low natural gas prices and subsidized renewables in a marketplace that does not reward carbon-free baseload capacity. These circumstances, along with increasingly antiquated labor-centric operating models and analog technology, have forced the premature closure of multiple nuclear facilities and placed a much larger population of nuclear power stations at risk. Nuclear plant economic survival in current and forecasted market conditions requires an efficient and technology-centric operating model that harvests the native efficiencies of advanced technology. This is analogous to transformations in nearly every other industry. In light of previous industry experience in modernizing safety Instrumentation and Controls (I&C) systems, nuclear utilities are reluctant to pursue these upgrades due to uncertainty in licensing and cost.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2021

Radial storage efficiency for CO 2 injection: Quantifying effectiveness of local flow control methods

Basin-scale geologic carbon sequestration will require hundreds of injection wells, each of which has costs related to property rights and regulatory requirements that correlate with the areal size of the carbon dioxide plume. These surface-footprint-related issues motivate maximizing storage efficiency radially away from each well. Radial storage efficiency (RSE) is defined here as the ratio of the volumetrically weighted carbon dioxide free-phase saturation within a given radius away from the injection well to the available pore space within that same radius. Maximizing RSE effectively minimizes the radial extent of the CO 2 plume. Optimizing RSE around individual wells requires local flow control injection strategies that can increase the uniformity of the filling of pore space around the well over the entire length of the perforated injection interval despite differences in local formation transmissivity. The goal of uniform filling of the storage reservoir starting from the well outward is to maximize carbon dioxide sweep and trapping locally outward from the well in all of the layers of the storage region before buoyancy forces predominate and drive carbon dioxide upward where it will spread laterally under lower-permeability layers. Example simulations of carbon dioxide injection into a layered storage system with and without local flow control are presented to show the advantage of uniformly filling all layers and how RSE can be used to quantify storage efficiency for the two different injection approaches. Published 2021. This article is a U.S. Government work and is in the public domain in the USA.

03 NATURAL GAS↗