Search NASA⌕ Search

SEARCH · Search NASA

Results for “dynamic systems analysis”

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.

878 records · Page 38
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Identification of potent inhibitors of JUN N-terminal kinases for treatment of endometriosis and associated pain

Endometriosis, defined as the ectopic growth of endometrial tissue outside of the uterine cavity, is an inflammatory and hormone-dependent disease that causes excruciating pelvic pain, infertility, and significantly decreases quality of life in affected patients. The JUN N-terminal kinases (JNKs) are a leading class of nonhormonal therapeutic targets that have been validated in preclinical models of endometriosis and in a Phase 1/2 clinical trial. Despite their therapeutic potential, JNK inhibitors with increased potency and specificity are needed to address the inflammatory pathology of endometriosis and to prevent disease progression. Leveraging a DNA-encoded chemical library collection of ~4 billion compounds, we identified lead inhibitor CDD-2428 and optimized derivatives, CDD-2728 and CDD-3013, with excellent binding affinity to JNK1-3 (K d = 0.12 to 3.7 nM), enhanced selectivity, metabolic stability, and cellular permeability. Crystallographic and biochemical studies confirmed that CDD-3013 exhibited superior kinase selectivity with improved efficacy compared to existing JNK inhibitors. In primary endometriosis cell models, CDD-2728 and CDD-3013 suppressed JNK-dependent inflammatory signaling, dampening pathways linked to pain, invasion, angiogenesis, and macrophage recruitment. In an endometriosis mouse model, both CDD-2728 and CDD-3013 reduced endometriotic lesion size, macrophage infiltration, and cellular proliferation, showing in vivo efficacy. When tested in a lipopolysaccharide-induced hyperalgesia model, CDD-2728 and CDD-3013 decreased markers of induced pain, as measured by changes in a dynamic weight bearing test and Grimace scores. These findings nominate CDD-2728 and CDD-3013 as potent, nonhormonal therapeutic candidates for endometriosis with broad anti-inflammatory and analgesic activity, addressing a critical unmet clinical need.

Madasu, Chandrashekhar [Department of Pathology an↗
DOE OSTI2025

Intermodal microwave-to-optical transduction using silicon-on-sapphire optomechanical ring resonator

Optomechanical and electro-optomechanical systems have emerged as one of the most promising approaches for quantum microwave-to-optical transduction to interconnect distributed quantum modalities for scaling the quantum systems. These systems use suspended structures to increase mode overlap and mitigate loss to achieve high efficiency. However, the suspended design’s poor heat dissipation under strong drive limits the ultimate efficiency. Here, we demonstrate an unsuspended optomechanical ring resonator (OMR) based on the silicon-on-sapphire (SOS) platform for microwave-to-optical frequency conversion. The OMR achieves a triply resonant optical-to-optical conversion with an enhanced coupling rateG b = 3.6 gigahertz per square-root milliwatt at a peak conversion efficiency of 1.2% with 3.6-milliwatt microwave drive power and a microwave-to-optical conversion efficiency of 1.5 × 10 −5 at 10-milliwatt optical drive power. Our results show that the unsuspended SOS platform, which mitigates the thermal effect and is compatible with superconducting qubits, is a promising platform for optomechanical circuitry and quantum transduction.

Science & Technology - Other Topics↗
DOE OSTI2025

Electricity Markets and Long-Duration Energy Storage: A Survey of Grid Services and Revenue Streams

Purpose of Review Long Duration Energy Storage (LDES) is increasingly viewed as a potential resource for providing grid services that enhance the stability and flexibility of electricity systems. While some LDES services are integrated into existing market frameworks, traditional mechanisms may not fully account for their operational characteristics, potentially leading to undervaluation. Within this context, this paper reviews the literature and industry practices to assess potential grid services for LDES, evaluates existing compensation mechanisms, and identifies challenges to full market integration. Recent Findings We first review existing literature and identify key grid services unique to LDES, including enhancing grid resilience during extreme weather events, enabling long-term energy shifting, and providing flexible and firm energy in systems with limited dispatchable resources. Here, we also review how LDES services are compensated in current market frameworks and the challenges associated with the full realization of LDES values. Additionally, we summarize market mechanisms for storage technologies across U.S. wholesale markets. We find that some markets are adjusting incentive structures, such as incorporating storage duration in capacity accreditation, to better align with system needs and LDES contributions to the grid. However, further refinements in capacity remuneration and dispatch timeframes may be needed for more effective realization of LDES value. Summary This review evaluates potential grid services for LDES, examines existing compensation mechanisms for LDES technologies, and identifies gaps between these mechanisms and LDES operational characteristics. The review concludes by outlining potential market enhancements for more effective LDES integration and articulating additional research needs to support its efficient participation in future power systems.

Flexible resources↗
DOE OSTI2026

Programmable Digital Devices used in Advanced Reactors

This paper introduces the concepts of common cause failure, diversity, and defense-in-depth used by the nuclear industry to analyze resilience in reactors. A survey of publicly traded and private companies building advanced reactors and their licensing status is presented. Safety and non-safety systems found in the NuScale Power design are summarized and the likely hardware and software categories used by those systems are enumerated. The importance of industry partners is highlighted. This paper also identifies an alternate path forward without industry partners to advance the knowledge needed to use artificial intelligence to analyze HBOMs and SBOMs to better understand reactor resiliency.

cybersecurity↗
DOE OSTI2026

Oversizing and Part-Load Problems

Oversizing, the common engineering practice of specifying devices with capacity exceeding the actual load requirement, is a widespread practice across virtually all building technologies end-use categories, including HVAC, electrical systems, lighting, appliances, and plug loads. This practice, driven by factors like design uncertainty, institutional pressures, and risk aversion, results in wasted capital investment, control difficulties, and excessive energy consumption due to inefficient part-load operation. Part-load operation, where devices run below maximum capacity, is the dominant operating mode in most energy systems and presents a complex design challenge. Solutions to match output to load fall into three broad categories: constraining the output, adjusting the device’s internal behavior, and linking output to energy storage or other waste-heat reuse applications. The energy implications of part-load are critical, as efficiency often drops sharply as load decreases across a wide diversity of devices. To quantify the extent of this problem, we derive a dimensionless Part-Load Metric (PLM) based on device efficiency and its frequency distribution of operating hours at various output levels. The PLM quantifies the deviation of a device's actual efficiency from its maximum design efficiency. This metric also serves as a measure of "capital inefficiency," enabling engineers to compare the impact of different part-load solutions and providing a unified framework for evaluating performance across various devices and systems.

Meier, Alan↗
DOE OSTI2026

Ecohydrology of urban environments

Urban watersheds include a gradient of highly managed to unmanaged ecosystems that vary in vegetation cover, built and natural drainage systems, and impervious cover. A distinct difference to undeveloped ecosystems is the dominance of material and energy import and export, the disruption of natural catenae, the covariation of soils, terrain and canopy conditions. In this chapter, we discuss the distinct characteristics of urban ecohydrology compared to non-urban systems. We also investigate emergent principles of unmanaged watersheds that underscore the promotion of ecohydrological resilience, and the need to abstract and adapt these principles to support the goal-oriented restoration of impaired urban ecosystems.

42 ENGINEERING↗
DOE OSTI2025

Utility Experience with Inverter Based Resource Impacts on Transmission Protection

High penetration of inverter-based resources (IBR) can adversely affect the transmission protection schemes in the area. With the proliferation of IBRs, utilities are finding out that conventional protection schemes are not adequately equipped to protect the electric systems because the existing system has been designed based on the fault current response of conventional rotating-machine-based generators. In several cases, the available literature does not provide any clear solution for the issues when the protection scheme does not operate properly near IBRs. This report identifies various protection challenges due to IBRs that industry is facing, from the utility perspective. It includes a broad review of all challenges that system protection has experienced so far with high penetration of IBRs. Key issues have been identified based on a survey of different utility protection engineers. Finally, fault event recordings are shown to demonstrate some actual cases of mis-operations of relays near IBRs.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2021

AI-Batt (Autonomous Identification of Battery Life Models) [SWR 21-36]

Autonomous Identification of Battery Life Models (AI-Batt) AI-Batt is a MATLAB code base for developing lifetime models for batteries from accelerated aging data. The code base provides many functions for processing, visualizing, and modeling battery aging data, making the data processing, exploration, and modeling workflow substantially faster. These tools are tailored for working with battery aging data sets, which usually consist of many separate time-series for each cell, with many test conditions and possible replicates at each condition, which makes it difficult to simply process or visualize the data set. Complex modeling tasks, such as cross-validation, sensitivity analysis, and uncertainty quantification have been implemented to enable thorough statistical investigation of model predictions. Additionally, several machine-learning algorithms are implemented to autonomously identify suitable models via symbolic regression. Data processing functions automatically cast data from the struct data type, which is commonly used to store experimental data, but is not an acceptable input for most algorithms, to the table data type, which can be easily used as input to any optimization algorithm. Also, the data can be separated into time-invariant and time-variant data tables, which is helpful for exploring the data set as well as developing separate models for time-variant and time-invariant aging mechanisms. For example, in aging tests with constant temperature, temperature is a time-invariant experimental condition. Visualization tools enable plotting of data, model fits, and model simulations possible with single-line function calls, empowering data exploration of complex data sets with both time-varying and time-invariant trends. Plots can be automatically generated for the whole data set, or separated by data group (groups of test replicates) or individual data series. Data points or data series can be automatically colored by the value of a variable with a variety of color maps, and model predictions can also be colored by the value of a fit statistic. Comparisons between data sets and the predictions/simulations of different models on the same data set can be easily plotted as well. Distributions of parameter values from bootstrap resampling can be plotted to visualize the reliability of parameter estimation, or determine any correlations between parameters. Modeling tools handle the complex task of creating and parsing symbolic equations for modeling battery lifetime. Equations are parsed to grab relevant data variables, parameter values, or specified sub-models for input into optimization, evaluation, or simulation functions. Models can be optimized locally (one set of parameters for each data series), bi-level (some parameters shared across the data set), or globally (single set of parameters for all data). Functions implementing symbolic regression algorithms help users to discover effective model equations, even in poorly sampled, high-dimensional data.

Smith, Kandler [National Renewable Energy Lab. (NR↗
DOE OSTI2024

Thermodynamically informed priors for uncertainty propagation in first-principles statistical mechanics

Here, this work demonstrates how first-principles statistical mechanics approaches within a Bayesian framework can quantify and propagate uncertainties to downstream thermodynamic calculations. To address the issue of Bayesian prior selection, knowledge of 0 K ground states in the material system of interest is incorporated into the prior. The effectiveness of this framework is shown by creating a phase diagram for the fcc zirconium nitride system, including confidence intervals on order-disorder transition temperatures.

Bayesian methods↗
DOE OSTI2026

DESI Strong Lens Foundry. III. Keck Spectroscopy for Strong Lenses Discovered Using Residual Neural Networks

We present spectroscopic data of strong lenses and their source galaxies using the Keck Near-Infrared Echellette Spectrometer (NIRES) and the Dark Energy Spectroscopic Instrument (DESI), providing redshifts necessary for nearly all strong-lensing applications with these systems, especially the extraction of physical parameters from lensing modeling. These strong lenses were found in the DESI Legacy Imaging Surveys using residual neural networks and followed up by our Hubble Space Telescope program, with all systems displaying unambiguous lensed arcs. With NIRES, we target eight lensed sources at redshifts difficult to measure in the optical range and determine the source redshifts for six, between z s = 1.675 and 3.332. DESI observed one of the remaining source redshifts, as well as an additional source redshift within the six systems. The two systems with nondetections by NIRES were observed for a considerably shorter 600 s at high airmass. Combining NIRES infrared spectroscopy with optical spectroscopy from our DESI Strong Lensing Secondary Target Program, these results provide the complete lens and source redshifts for six systems, a resource for refining automated strong lens searches in future deep- and wide-field imaging surveys and addressing a range of questions in astrophysics and cosmology.

Agarwal, Shrihan [University of Chicago, IL (Unite↗
DOE OSTI2026

Durability of ABC Team Wall Assemblies

This study evaluates the long-term hygrothermal durability of four advanced retrofit wall systems using a field test facility located in Hollywood, South Carolina, representative of a mixed-humid coastal climate (Zone 3). The research focuses on assessing the thermal and moisture performance of a baseline wall assembly retrofitted with (1) Tremco/Dryvit prefabricated panel systems (Revitalite and Fedderlite), (2) Reinforced Fiberglass Plastic (RFP) panels developed by Oak Ridge National Laboratory (ORNL), and (3) Vacuum Insulated Panels (VIP) integrated with EPS by Home Innovation Research Labs (HIRL) and ORNL. A baseline wall representing a typical uninsulated wood-frame construction was used for comparison.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2026

Fine-Grained Power and Energy Attribution on AMD GPU/APU-Based Exascale Nodes

Modern exascale GPU- and APU-based systems provide multiple power and energy sensors, but differences in scope, update rate, timing, and filtering complicate the attribution of short-lived accelerator activity. This paper presents a methodology to characterize and correct these effects on Cray EX systems with AMD Instinct MI250X GPUs (Frontier) and MI300A APUs (Portage). Using controlled square-wave workloads, we quantify update intervals, delay, aliasing, and variability across up to 512 GPUs and 480 APUs with on-chip (rocm-smi/amd-smi) and off-chip Cray Power Management sensors. We reconstruct power from cumulative energy counters to achieve faster response times, validate it against on-chip, off-chip, and node-level sensors, and integrate the resulting streams into a Score-P/PAPI-based tool for time-aligned, phase-level attribution. Applied to rocHPL, rocHPL-MxP, and HPG-MxP, the method separates energy savings due to reduced runtime from changes in power. Mixed precision reduces node energy on Frontier by 79% for rocHPL-MxP and 31% for HPG-MxP, with similar trends on Portage. These results provide portable guidance for sensor validation and power-aware optimization on current and future exascale systems.

Mcdaniel, Adam [ORNL] (ORCID:000000016926028X)↗
DOE OSTI2025

Erosion of high-Z refractory coatings for helicon plasma source window

Proto-MPEX (Materials Plasma Exposure eXperiment), a linear plasma device (LPD), using a high power radio frequency (RF) (⩾100 kW, 13.56 MHz) helicon plasma source, has suffered RF sheath-induced window erosion. The sputtered impurities from the window surface transport toward the downstream target affecting plasma-material interaction studies. The rectified sheath voltage formed was high enough to cause window erosion by light ions due to its low-Z components (e.g. Si 3 N 4 and AlN). Hence, we are proposing impurity mitigation strategies, which involve the application of a Faraday screen to lower the rectified sheath voltage and the application of high-Z refractory coatings on the existing window plasma-facing surface. In this work, we report the erosion of two different coatings, i.e., tantalum oxide (Ta 2 O 5 ) and hafnium oxide (HfO 2 ) on a silicon nitride (Si 3 N 4 ) window material under conditions of low-energy deuterium (D) ion impact, well below the Ta 2 O 5 sputtering energy threshold (i.e. 250 eV). The test samples were RF-biased and exposed to a high D-ion fluence (∼10 26 m −2 ) in Plasma Interaction with Surface Component Experimental Station (PISCES-A) LPD. The experimentally measured sputtering yields of the high-Z refractory coatings below the sputtering energy threshold of Ta 2 O 5 indicate an increased erosion due to the plasma impurities, especially due to oxygen. Improving the vacuum level could reduce the oxygen impurities and increase the lifespan of the coated helicon window for plasma operation. Post-surface analysis indicates no preferential erosion of oxygen or enrichment of the high-Z surface component. This is likely due to an effective energy transfer from the impurity ion for sputtering of the high-Z component in the coating. With the reduced rectified sheath voltage at the window surface and improved vacuum conditions, the proposed high-Z refractory coatings offer a promising solution for reducing window erosion in future MPEX plasma operations at ORNL.

RF bias↗
DOE OSTI2026

AT2025ulz and S250818k: Leveraging DESI Spectroscopy in the Hunt for a Kilonova Associated with a Subsolar-mass Gravitational-wave Candidate

On 2025 August 18, the LIGO–Virgo–KAGRA collaboration reported a subthreshold gravitational-wave candidate detection consistent with a subsolar-mass neutron star merger, denoted S250818k. An optical transient, AT2025ulz, was discovered within the localization region. AT2025ulz initially appeared to meet the expected behavior of kilonova emission, the telltale signature of a binary neutron star merger. The transient subsequently rebrightened after ∼5 days and was classified as a Type IIb supernova. In this work, we analyze the observations of its host galaxy obtained by the Dark Energy Spectroscopic Instrument (DESI). From the DESI spectrum, we obtain a secure redshift of z = 0.084840 ± 0.000006. If S250818k has an astrophysical origin, this places the transient within 2σ of the gravitational-wave distance and results in an integral overlap between the gravitational-wave alert and the transient location of log 10 $\mathcal{I} ≈ 3.9 - 4.2$. Our analysis of the host galaxy’s spectral energy distribution reveals a star-forming, dusty galaxy with stellar mass ∼10 10 M ⊙ , broadly consistent with the population of both short gamma-ray bursts and core-collapse supernova host galaxies. We also present our follow-up of DESI-selected candidate host galaxies using the Fraunhofer Telescope at the Wendelstein Observatory, and show the promise of DESI for associating or rejecting candidate electromagnetic counterparts to gravitational-wave alerts. These results emphasize the value of DESI’s extensive spectroscopic dataset in rapidly characterizing host galaxies, enabling spectroscopic host subtraction, and guiding targeted follow-up.

Hall, Xander J. [Carnegie Mellon University, Pitts↗
DOE OSTI2026

Design and Modeling of the Charge Readout of a SiMOS Quantum Dot with a Single Electron Transistor and CryoCMOS

Single electron spin qubits trapped in SiMOS quantum dots are a promising technology for scaling to thousand- or million-qubit systems due to their compatibility with mature CMOS manufacturing processes. A readout system that combines a single electron transistor with a custom cryogenic CMOS amplification and digitization chain offers key advantages by avoiding the use of bulky RF components or room-temperature interconnects. We present design techniques and simulation results for an optimized qubit-SET cryoCMOS interface, culminating in the design of the QNDR1 ASIC, the first cryogenic readout ASIC designed under the Quandarum project, which targets the development of a many-channel spin qubit based detector for use in high energy physics.

Quinn, Adam [Fermilab] (ORCID:0009000605056371)↗
DOE OSTI2026

Case Study of Integrating High-Temperature Heat Pump with LiBr-H2O Absorption Chiller for Data Center Liquid Cooling

Data centers (DCs) are physical infrastructures that support artificial intelligence workloads. The rapid growth of artificial intelligence is putting substantial pressure on the US power grid. Most of electricity consumed by IT equipment, accounting for 50%-60% of total DC power, ultimately becomes waste heat. This heat is dissipated by DC’s cooling facilities, accounting for an additional 30%-40% of total DC power. Recovering and repurposing this waste heat offers a significant opportunity to enhance energy efficiency and reduce operating costs of DCs. One potential pathway is converting heat to cold using thermal-driven absorption chillers, therefore, reducing the power consumption in DC cooling facilities. Existing studies mainly demonstrate the technical and economic feasibility of repurposing DC’s waste heat for cooling applications but provide limited technical details on how to integrate the thermal-driven absorption chillers with DC cooling systems. In addition, the low-grade waste heat available from DCs must be upgraded to higher temperatures suitable for absorption chillers. This paper presents a case study on integrating high-temperature heat pumps with a LiBr-H2O absorption chiller to use DC waste heat for cooling. A thermodynamic model of single-effect, LiBr-H2O absorption chiller and an empirical model of high-temperature heat pumps were built. The case study considers ASHRAE W17 liquid-cooled DC, with facility service water supplied at 17.0℃ and returned at 25.3℃. The thermal behaviors of absorption chiller components were predicted for the generation temperature ranging from 75.0℃ to 115.0℃. Based on the available waste heat in the integrated system, two waste heat recovery strategies were evaluated: a facility service water-based strategy and cooling water-based strategy. Results indicated that the cooling water-based strategy achieves higher Coefficient of Performance (COPs) than the facility service water-based strategy. The relatively low cooling COPs of single-effect LiBr-H2O absorption chillers could be offset by high heating COP of high temperature heat pumps. The maximum cooling COP of absorption chiller and the overall COP of integrated systems occur at lower generation temperatures, but these conditions also yield lower cooling capacities. In practice, system operation should balance the trade-off between the COP and cooling capacity

Wang, Pengtao [ORNL] (ORCID:0000000214713429)↗
DOE OSTI2025

Topological prethermal strong zero modes on superconducting processors

Abstract Symmetry-protected topological phases 1–4 cannot be described by any local order parameter and are beyond the conventional symmetry-breaking model 5 . They are characterized by topological boundary modes that remain stable under symmetry respecting perturbations 1–4,6–8 . In clean, gapped systems without disorder, the stability of these edge modes is restricted to the zero-temperature manifold; at finite temperatures, interactions with mobile thermal excitations lead to their decay 9–11 . Here we report the observation of a distinct type of topological edge mode 12–14 , which is protected by emergent symmetries and persists across the entire spectrum, in an array of 100 programmable superconducting qubits. Through digital quantum simulation of a one-dimensional disorder-free stabilizer Hamiltonian, we observe robust long-lived topological edge modes over up to 30 cycles for a wide range of initial states. We show that the interaction between these edge modes and bulk excitations can be suppressed by dimerizing the stabilizer strength, leading to an emergent U(1) × U(1) symmetry in the prethermal regime of the system. Furthermore, we exploit these topological edge modes as logical qubits and prepare a logical Bell state, which exhibits persistent coherence, despite the system being disorder-free and at finite temperature. Our results establish a viable digital simulation approach 15–18 to experimentally study topological matter at finite temperature and demonstrate a potential route to construct long-lived, robust boundary qubits in disorder-free systems.

Science & Technology - Other Topics↗
DOE OSTI2025

Development of Hydrogen Burner for FT4000® Aeroderivative Engine - Final Report

This report details an effort to develop a retrofittable fuel/air injector for the FT4000® aeroderivative gas turbine that enables use of hydrogen as a carbon-free fuel for efficient power generation. The FT4000 engine’s low-NOx combustor was developed by Pratt & Whitney and RTX Technology Research Center with core technology from the Pratt & Whitney PW4000 turbofan aircraft engine. The current FT4000 production engine operates on either natural gas or No. 2 fuel oil with water injection to achieve high thermal efficiency and low emissions. This engine is fielded by Mitsubishi Power Aero and delivers 70 MW of power with a simple-cycle efficiency of over 41% when operating with wet compression. The work reported here advances the technology readiness level of the FT4000 combustor for operation with hydrogen, starting with an experimental assessment of the current production hardware with increasing hydrogen content mixed with natural gas and ending with improved nozzle concepts for fully robust operation with 100% hydrogen. High-pressure single-sector combustor rig tests have been completed, demonstrating the ability for the dual fuel nozzle to operate an FT4000 combustor on 100% hydrogen with low nitrogen oxide (NOx) emissions. Metal temperature measurements and video images of the flame structure from zero to 100% hydrogen highlight opportunities to improve the fuel nozzle robustness for high hydrogen conditions. The design of new fuel/air mixer concepts to improve durability and operability with high hydrogen levels was also completed. A total of eleven new concepts were developed and analyzed, ranging from modifications to the bill-of-materials nozzle to fully clean sheet designs. The concepts were evaluated with non-reacting and reacting flow evaluations to assess the performance of the new hardware designs. Non-reacting tests included Phase Doppler Particle Analysis (PDPA) for droplet size and velocity, mechanical patternation for liquid water flux, and planar laser induced fluorescence (PLIF) with acetone-seeding for gaseous fuel/air mixing characterization. Five scaled candidate fuel nozzle designs, in addition to a scaled bill-of-materials nozzle, were then successfully evaluated in an atmospheric pressure burner rig. The nozzles were evaluated for performance with natural gas, hydrogen/natural gas blends, and pure hydrogen. For 100% hydrogen, the nozzles were evaluated with and without water injection. Optical and infrared imaging of the flame and fuel nozzle was captured. NOx emissions were sampled from fixed emissions probes. The results show measurable differences between the various designs, and the data was used to down-select the two most promising designs to advance to future full pressure rig testing. Results from this study have cleared the current production FT4000 engines with dual fuel nozzles to operate at baseload power on blends of hydrogen mixed with natural gas and water. Two promising nozzle designs have been developed to enable fully robust operation with up to 100% hydrogen. These nozzles require validation at full baseload operating conditions before they can be introduced to the field.

03 NATURAL GAS↗