Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “System Testing”

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 253 records · Page 14

Photovoltaic Analysis and Response Support (PARS) Platform for Solar Situational Awareness and Resiliency Services

The project's primary objective is to develop a digital-twin based Photovoltaic (PV) Analysis and Response Support (PARS) platform, which aims to provide real-time situational awareness and optimal response plans. This platform is designed to enhance the performance of hybrid PV systems, making them competitive with or even superior to conventional generation resources. The PARS platform enabled the project team to develop and evaluate an extensive suite of grid support functionalities for the hybrid PV systems to enhance grid performance, across key areas including visibility, dispatchability, security, resilience, and reliability. Given the global push toward achieving 100% clean energy by 2035, there is a significant increase in the integration of inverter-based resources (IBRs) throughout the energy grid. Effectively managing the inherent variability and uncertainty associated with IBRs is crucial for ensuring cost-effectiveness, reliability, and security in both the main grid and islanded microgrids. Constrained to a limited array of IEEE test systems or standard feeder models, traditional IBR modeling struggles to assimilate new field data, accurately reflect system dynamics, and adapt to the evolving energy landscape. In our project, we embraced a Digital Twin (DT) strategy for crafting the PARS platform. A digital twin acts as a precise virtual counterpart of a physical system, built on historical data and continuously honed with real-time insights. This enables the high-fidelity DT to accurately mirror current system operations and forecast future scenarios. Consequently, the PARS platform becomes an ideal environment for testing and refining monitoring, control, power, and energy management algorithms designed to boost hybrid PV system performance. The defining feature of the PARS platform, distinguishing it from other advanced simulation tools, is its exceptional adaptability. This is achieved by employing actual network topologies and utilizing real-time field data for fine-tuning and calibration, ensuring a close emulation of real-world conditions. The project deliverables include: 1) High-fidelity IBR models and tools for real-time parameterization, utilizing real-time field measurements to refine IBR models for enhanced accuracy and performance; 2) Grid-forming and Grid-following capabilities to deliver resilience services, including blackstart, voltage and frequency support, cold-load pick-up, power reserves, and three-phase load balancing across grid-connected and microgrid settings; 3) Machine learning-based forecasting tools and methods for generating synthetic data and topologies, creating diverse and realistic simulation environments for evaluating varied operational scenarios; 4) Advanced microgrid power and energy management algorithms for optimizing the integration and operation of PV, storage, and demand response resources within both feeder and community scales. The power grid data sets are provided by four utility companies in North Carolina and the New York Power Administration. Acting as industry advisors, our industry partners communicated stakeholder needs and regulatory standards to the research teams, aiding technology transfer by incorporating the developed methodologies into their daily operations. This collaboration ensures that the PARS platform, functioning as a power system digital twin, enhances our understanding of IBR dynamic behaviors and enables the development and evaluation of IBR control functions that match or exceed the capabilities of conventional synchronous generators.

14 SOLAR ENERGY↗

A Transactive Approach for Service Restoration Utilizing Customer Load Flexibility and Grid-Edge Resources: Preprint

This paper develops a transactive energy system model to restore electricity to customers in an isolated distribution system after an outage. The model engages a variety of customer types -- prosumers, flexible loads, critical/noncritical customers, and distributed generators -- as active participants in the restoration process. Unlike many existing transactive approaches, the proposed model is developed for service restoration and accounts for various customer types and their \textit{autonomy} and \textit{privacy} through an iterative approach to determine the optimal market price, while maintaining system-level power flow and voltage constraints. The advantages of the proposed approach are numerically validated on a modified IEEE 123-bus test system.

distributed energy resources↗

A Transactive Approach for Service Restoration Utilizing Customer Load Flexibility and Grid-Edge Resources

This paper develops a transactive energy system model to restore electricity to customers in an isolated distribution system after an outage. The model engages a variety of customer types -- prosumers, flexible loads, critical/noncritical customers, and distributed generators -- as active participants in the restoration process. Unlike many existing transactive approaches, the proposed model is developed for service restoration and accounts for various customer types and their autonomy and privacy through an iterative approach to determine the optimal market price, while maintaining system-level power flow and voltage constraints. The advantages of the proposed approach are numerically validated on a modified IEEE 123-bus test system.

distributed energy resources↗

High-Frequency Signature-Based Fault Detection for Future MV Distribution Grids: Preprint

Increased deployment of Distributed Energy Resources (DERs) and microgrids are significantly changing the operation of distribution grid. Increasing penetration of power electronics based DERs also raise several operational problems such as protection coordination issues, voltage raise, recloser synchronization, stability, and frequency control. Legacy systems are designed for unidirectional power flow. Bidirectional power flows, and varying short circuit levels are some of the major challenges in implementing these legacy protection schemes for future power grid with high penetration of power electronics based DERs. Inverter based DER’s can contribute to a maximum of 2 p.u. fault current. Interconnection standards might require even lower fault current contribution. If these protection related challenges are not addressed, it will result in limiting the penetration of DERs in distribution feeders. In this paper, the feasibility of high frequency signature generated by transients is evaluated to identify and locate faults in a distribution system. The paper will discuss the need for an accurate line modeling approach to accurately capture the high frequency signatures. These high frequency signatures are naturally occurring electromagnetic transients for events such as switching events in the network, faults events, etc. This is different than that implemented in transmission systems as the proposed approach utilizes advanced signal processing techniques to retrieve valuable information from the high frequency signatures generated during faults and does not just rely on time delay measurements between reflections of the high frequency signals. In the full paper, information of a modified IEEE distribution feeder test system to accurately model the high frequency signatures will be provided. The distribution lines will be modeled using universal line modeling (ULM) technique in an electromagnetic transients program (EMTP) modeling platform. Different test cases will be simulated and these test cases will include different fault types, multiple fault locations, capacitor bank switching, etc. The feeder system modeled for these simulations will also include power electronics based DERs. Additionally, experimental results demonstrating the high frequency signatures generated due to a transient event in a medium voltage (13.8 kV) distribution line will also be presented in the full paper. Furthermore, the analysis combined with simulation and experimental results in the full paper will help demonstrate the feasibility and efficacy of the proposed protection scheme. The developed protection scheme will not just help solve problems generated by bidirectional current flow, varying short-circuit current capacity, but also help in proliferation of DERs in distribution systems.

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

Tank Side Cesium Removal for Hanford Tank Waste - Actual Waste Demonstrations - 20467

The US DOE and its contractors are evaluating several options for treating Hanford tank supernate streams to provide feed for the Low-Activity Waste Facility at the Hanford Waste Treatment and Immobilization Plant. At this time, the intent is to deploy a Tank-Side Cesium Removal system. The key unit operations for this process are filtration and cesium removal by ion exchange. To assess the anticipated performance of the system, tests with actual waste were performed with samples from Hanford tanks AW-102 and AP-107 to provide basis data regarding the performance of the system in treating these materials. The decanted supernate waste samples from tanks AW-102 and AP-107 were received directly by the Radiochemical Processing Laboratory at Pacific Northwest National Laboratory. These samples were then filtered to provide data on dead-end filtration with a Mott{sup TM} filter media grade 5. Subsequently, the samples were run through 10-milliliter ion exchange columns containing crystalline silicotitanate to assess the dynamic ion exchange performance of these feeds. Test results demonstrated that the planned TSCR operations should meet process performance expectations. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Caution on Using Tetrahydrofuran for Processing Crystalline Silica Samples From Engineered Stone for XRD Analysis

Abstract We conducted laboratory experiments to investigate a suspected effect of tetrahydrofuran (THF) on quantifying crystalline silica in samples collected from working with engineered stone when THF is used to process samples prior to the X-ray diffraction (XRD) analysis. Two groups of samples from grinding either engineered stone or granite were simultaneously taken from a laboratory testing system, with one group of samples using THF for processing and another group using muffle furnace for ashing. For each stone type, we also tested four levels of respirable dust loading on the samples by varying the grinding time from 1 to 8 min. Statistical analysis of the experimental results on crystalline silica contents of the two groups of samples showed that the difference between the two methods was not significant (P ≥ 0.05) for the granite at all four levels of respirable dust loading and for the engineered stone at the two levels of respirable dust loading greater than 0.5 mg. However, the crystalline silica content from using THF processing was significantly lower (P = 0.001) than that from using muffle furnace ashing for engineered stone when the respirable dust loading levels were less than 0.5 mg. For the engineered stone dust samples with grinding times of 1 and 2 min, the average respirable dust loading was about 0.19 and 0.34 mg, respectively; while the crystalline silica content from using THF processing was 30.9 and 21.5% lower than that from using muffle furnace ashing, respectively. Since most full-shift samples from field assessments in this industry are expected to have respirable dust loading less than 0.5 mg, muffle furnace or radio frequency plasma ashing should be specified as the preferred sample processing method instead of the THF processing method for quantification of crystalline silica when engineered stone is expected to present to avoid artificially reduced silica content values, which are likely caused by the reactions between THF and the resins in engineered stone.

Public, Environmental & Occupational Health↗

Electro-Thermo-Mechanical (ETM) Study on Submarine Dynamic Power Cables (SDPC) for Offshore Wind Transmission

Understanding the failure mechanisms of submarine dynamic power cables (SDPC) is critical for innovative design to meet the 2035 cost reduction target of U.S. DOE Floating Offshore Wind Shot. This is important because the current design suffers a significant failure rate in the field. This project conducted a systematic electro-thermo-mechanical (ETM) experimental study on the power cores extracted from a 15kV power cable with three cores of copper conductor, ethylene propylene rubber (EPR) insulation, and continuously corrugated welded aluminum armor (CCWA). An ETM testing system was developed by integrating a high voltage (HV) amplifier, two ceramic heaters, and a rod-plate transverse compression setup into a mechanical testing machine. Increasing temperature from room temperature (RT, 22 degree C) to 90 degree C resulted in the 67% decrease in the failure mechanical load as defined by the dielectric breakdown. Under creep mechanical loading, the dielectric breakdown time was decreased by 70% for a given mechanical load when the specimen temperature increased from RT to 90oC. The failure strain in both monotonic and creep mechanical loading modes was related to the maximum mechanical load applied. Although the dielectric breakdown occurred in the ETM test, the post-test measurement revealed an impressive recovery of electric resistance. The failure analysis on cross section of tested specimens indicated a sizable gap across insulation layer near the area between copper conductor and loading rod, which apparently resulted from online dielectric failure and offline elastic recovery of components including conductor and insulation layer.

17 WIND ENERGY↗

DistOPF: Advanced Solutions for Distribution Optimal Power Flow Analysis - DistOPF v0.2 Documentation

To achieve an affordable and reliable energy system, research on power distribution system is often focused on integration of distributed generators, energy storage solution, EV charging, smart meters, and other advanced assets that may benefit from or require more advanced control and optimization techniques. Despite this focus on advanced distribution system topics, early researchers and grid scientists often start from scratch when developing optimization programs for power distribution systems. This report introduces DistOPF, a Python package that consolidates years of research into a versatile and modular tool. DistOPF provides researchers with essential capabilities to solve distribution system Optimal Power Flow (OPF) problems using standard network models. Additionally, it offers a platform to benchmark both new and existing algorithms against established test systems.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Optimal Sizing of Movable Energy Resources for Enhanced Resilience in Distribution Systems: A Techno-Economic Analysis

This article introduces a techno-economic analysis aimed at identifying the optimal total size of movable energy resources (MERs) to enhance the resilience of electric power supply. The core focus of this approach is to determine the total size of MERs required within the distribution network to expedite restoration after extreme events. Leveraging distribution line fragility curves, the proposed methodology generates numerous line outage scenarios, with scenario reduction techniques employed to minimize computational burden. For each reduced multiple line outage scenario, a systematic reconfiguration of the distribution network, represented as a graph, is executed using tie-switches within the system. To evaluate each locational combination of MERs for a specific number of these resources, the expected load curtailment (ELC) is calculated by summing the load curtailment within microgrids formed due to multiple line outages. This process is repeated for all possible locational combinations of MERs to determine minimal ELC for each MER total size. For every MER total size, the minimal ELCs are determined. Finally, a techno-economic analysis is performed using power outage cost and investment cost of MERs to pinpoint an optimal total size of MERs for the distribution system. To demonstrate the effectiveness of the proposed approach, case studies are conducted on the 33-node and the modified IEEE 123-node distribution test systems.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Study on Foam Fracturing of Granite for the Development of Enhanced Geothermal Systems

Foam fracturing has been considered as a potential technical approach to address water concerns with hydraulic stimulation in the development of enhanced geothermal systems (EGS). This paper reports the recent progress of a project sponsored by the U.S. DOE GTO Waterless Stimulation Initiative. The project investigates the foam fracturing as an alternative approach to the hydraulic fracturing. The paper describes the development of foam fracturing testing system at Oak Ridge National Laboratory (ORNL). The system can be used to perform the hydraulic fracturing of rock-like materials with both water and foamed fluids under pressure to 6,000 psi (41.4 MPa). The system features the sections of foam generation and foam injection and has the capability of testing materials in both monotonic and cyclic injection modes.Experimental results of foam fracturing on Charcoal Black granite will be reported. Cylindrical granite specimens with a blind hole were studied along with water and nitrogen-gas-in-water foam as the fracturing fluids. The effects of fracturing fluids and injection modes on the breakdown pressure and failure responses of the material were examined. The observations from the experimental work and the implications to the EGS application will be presented and discussed.

Wang, Hong↗

GRIDAPPSD/distopf (33583-E)

DistOPF is an open-source Python package providing a three-phase, asymmetric optimal power flow (OPF) tool specifically designed for distribution systems. The key inventive features include: - Asymmetrical 3-phase OPF modeling for distribution systems with unbalanced phases - Comprehensive control optimization supporting both active (P) and reactive (Q) power control variables - Built-in visualization and validation tools - Standard test system benchmarking platform for algorithm development and comparison - Modular CSV-based input system using Pandas DataFrames for flexible model specification - Standard power distribution model importer enabling direct conversion from CIM and OpenDSS format to optimization-ready models - Multiple solve interface compatibility (PYOMO, CVXPY, SciPy) with automatic solver selection based on problem type

Gray, Nathan [Pacific Northwest National Laborator↗

A cell-based Papain-like Protease (PLpro) activity assay for rapid detection of active SARS-CoV-2 infections and antivirals

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) and its variants are a continuous threat to human life. An urgent need remains for simple and fast tests that reliably detect active infections with SARS-CoV-2 and its variants in the early stage of infection. Here we introduce a simple and rapid activity-based diagnostic (ABDx) test that identifies SARS-CoV-2 infections by measuring the activity of a viral enzyme, Papain-Like protease (PLpro). The test system consists of a peptide that fluoresces when cleaved by SARS PLpro that is active in crude, unprocessed lysates from human tongue scrapes and saliva. Test results are obtained in 30 minutes or less using widely available fluorescence plate readers, or a battery-operated portable instrument for on-site testing. Proof-of-concept was obtained in a study on clinical specimens collected from patients with COVID-19 like symptoms who tested positive (n = 10) or negative (n = 10) with LIAT RT-PCR using nasal mid turbinate swabs. When saliva from these patients was tested with in-house endpoint RT-PCR, 17 were positive and only 5 specimens were negative, of which 2 became positive when tested 5 days later. PLpro activity correlated in 17 of these cases (3 out of 3 negatives and 14 out of 16 positives, with one invalid specimen). Despite the small number of samples, the agreement was significant (p value = 0.01). Two false negatives were detected, one from a sample with a late Ct value of 35 in diagnostic RT-PCR, indicating that an active infection was no longer present. The PLpro assay is easily scalable and expected to detect all viable SARS-CoV-2 variants, making it attractive as a screening and surveillance tool. Additionally, we show feasibility of the platform as a new homogeneous phenotypic assay for rapid screening of SARS-CoV-2 antiviral drugs and neutralizing antibodies.

60 APPLIED LIFE SCIENCES↗

A Novel Catalytic Membrane Reactor for DME Synthesis from Renewable Resources

Production of liquid fuels or chemicals from CO 2 (captured from the air or flue gases) and renewable hydrogen presents a new approach to producing clean fuels domestically. While significant progress has been made in the area of renewable electricity generation from solar and wind, a large gap remains with respect to the production of renewable liquid fuels/chemicals. Other processes for producing liquid fuels/chemicals from renewable electricity are constrained by thermodynamic limitations, making them prohibitively expensive and impractical. The team is overcoming these limitations and developing catalytic membrane reactor processes with high yields and low energy penalties. Supported by the Advanced Research Projects Agency-Energy (ARPA-E) of the US Department of Energy (DOE), GTI Energy and partners have been developing a technology for the production of renewable dimethyl ether (DME) from carbon dioxide (CO 2 ) and renewable hydrogen (H 2 ) using a novel catalytic membrane reactor and demonstration of this system at a scale of 1 kg/day. DME is a clean-burning, non-toxic fuel with a high cetane value (55-60), making it an excellent diesel alternative. DME can be stored as a liquid under moderate pressure, eliminating the need for the high-pressure containers used for CNG or cryogenics, as in the case of LNG. DME is also approved as a renewable fuel under the U.S. Environmental Protection Agency’s Renewable Fuels Standard (RFS), making it eligible for Renewable Identification Numbers (RINs) credits. By producing DME through the catalytic conversion of captured CO 2 and renewable H 2 , this process will produce renewable liquid transportation fuel and a means of large-scale utilization of captured CO 2 . In the DME synthesis process, CO 2 and H 2 are fed to a hollow fiber catalytic membrane reactor at 300-600 psig that contains a bi-functional catalyst that combines two reactions, methanol synthesis (CO 2 + 3H 2 → CH 3 OH + H 2 O) and methanol dehydration (2CH 3 OH → CH 3 OCH 3 + H 2 O), into a one-step process to produce DME. The bifunctional catalyst converts methanol to DME, enabling higher overall CO 2 conversion. A Cu/ZnO/ZrO 2 /Al 2 O 3 (CZZA) catalyst is used for methanol synthesis and is coupled with a zeolite catalyst H-ZSM-5 for dehydration. This one-step process intensifies a process that would otherwise require multiple reaction steps. However, combining these two reactions results in increased water production which inhibits catalytic activity. Here, the Na + -gated, water-transport membrane (Science, vol. 367, pp. 667, 2020), removes water in situ, shifting the thermodynamic equilibrium towards product formation while decreasing kinetic inhibition from water adsorption onto the catalyst surface. The Na + gated, water-transport nanochannel membrane showed H 2 O/CO 2 selectivity of 560 at 250 °C and 300 psig for H 2 O/CO 2 /CO/H 2 /MeOH gas mixtures. The selectivities of H 2 O/H 2 , H 2 O/CO, and H 2 O/MeOH were 190, 170, and 80, respectively. In a laboratory-scale membrane reactor, DME synthesis testing using this membrane, a DME production rate of 440 g DME /kg cat /h was achieved at 260 °C and 550 psig. Compared to the packed bed reactor, the CO 2 conversion and DME production rate in the membrane reactor were 80% and three times higher, respectively. A prototype test system (1 kg/day) was designed, constructed, and tested. A DME production rate of 1.31 kg/day and a DME productivity of 360 g/h/kg were achieved in the prototype membrane reactor. Good stability was demonstrated during 150-h continuous operation and multiple startups/shutdowns tests.

10 SYNTHETIC FUELS↗

Optimal Gas-Electric Energy System Decarbonization Planning

As energy utilities implement climate change mitigation policies, system planners require strategies for achieving affordable emissions reductions. Coordinated planning of electric power and natural gas systems will allow synergistic investments to address cross-sector operational constraints, competing uses for net-zero emissions fuels, and shifts in energy demands across energy carriers. In this study, we develop a novel optimization program that finds the cost-minimizing mix of infrastructure expansion or reduction across gas and electric systems to satisfy sector-specific emissions constraints. Alongside energy supply resources, our framework allows for central-planning of end-use equipment stocks to allow switching between gas and electric appliances upon failure or premature replacement. The proposed model is used to simulate case study scenarios for a benchmark 24-pipe gas network coupled to a 24-node power system test network. We find that electrification of greater than 80% of core gas demands is a component of the least-cost solution for modeled energy systems. Despite this substitution, the gas system is maintained to service difficult-to-electrify customers and to deliver net-zero emissions gas to electricity generators in times of peak electricity demand. Restricting electrification of gas appliances increases reliance on power- to-gas technologies and increases annual costs by 15% in 2040. Neglecting constraints on pipeline blending of hydrogen can produce a misleading result that relies on hydrogen blend fractions of greater than 50%. In all cases, we find the average costs of delivered gas increase nearly 5-fold across the decarbonization transition, highlighting the importance of future work investigating cost-allocation strategies for ensuring an equitable energy transition.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

The LSBmax algorithm for boosting resilience of electric grids post (N‐2) contingencies

Abstract A computationally improved algorithm is presented to find the best transmission switching (TS) candidate for boosting resilience of electricity grids subject to ( N ‐2) contingencies. Here, resilience is computed as the reduction in load shed after the above‐mentioned ( N‐ ) contingencies. TS is a planned line outage, and past research shows that changing the transmission system's topology changes the power flow and removes post contingency violations. Finding the best TS candidate in a computationally suitable time for effectively boosting resilience is a challenge. The best TS candidate is found using a novel heuristic method by decreasing the search space based on proximity to the bus with the maximum load shedding (LSB). The LSB algorithm is faster than existing algorithms in the literature; and, it is compatible with both the AC and DC optimal power flow formulations. To validate the authors' claims of speedup and accuracy, two metrics are used to analyze the results from the IEEE 39‐bus and 118‐bus systems. Finally, the inherent parallelism of the LSB algorithm is leveraged on a high‐performance computing platform and applied to the large‐scale Polish 2383‐bus test system to validate scalability in both size and speedup in computation time.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Axisymmetric hydrodynamics in numerical relativity using a multipatch method

In this paper, we describe a method of implementing the axisymmetric evolution of general-relativistic hydrodynamics and magnetohydrodynamics through modification of a multipatch grid scheme. In order to ease the computational requirements required to evolve the post-merger phase of systems involving binary compact massive objects in numerical relativity, it is often beneficial to take advantage of these system's tendency to rapidly settle into states that are nearly axisymmetric, allowing for 2D evolution of secular timescales. We implement this scheme in the spectral Einstein code and show the results of application of this method to four test systems including viscosity, magnetic fields, and neutrino radiation transport. Our results show that this method can be used to quickly allow already existing 3D infrastructure that makes use of local coordinate system transformations to be made to run in axisymmetric 2D with the flexible grid creation capabilities of multipatch methods. Our code tests include a simple model of a binary neutron star postmerger remnant, for which we confirm the formation of a massive torus which is a promising source of post-merger ejecta.

79 ASTRONOMY AND ASTROPHYSICS↗

Establishing the accuracy of position–specific carbon isotope analysis of propane by GC–pyrolysis–GC–IRMS

Rationale: Position-specific (PS) δ 13 C values of propane have proven their ability to provide valuable information on the evolution history of natural gases. Two major approaches to measure PS δ 13 C values of propane are isotopic 13 C nuclear magnetic resonance (NMR) and gas chromatography-pyrolysis-gas chromatography-isotope ratio mass spectrometry (GC-Py-GC-IRMS). Measurement accuracy of the isotopic 13 C NMR has been verified, but the requirements of large sample size and long experimental time limit its applications. GC-Py-GC-IRMS is a more versatile method with a small sample size, but its accuracy has not been demonstrated. Methods: We measured the PS δ 13 C values of propane from nine natural gases using both 13 C NMR and GC-Py-GC-IRMS, then evaluated the accuracy of the GC-Py-GC-IRMS method. Results: The results show that large carbon isotope fractionations occurred for both terminal and central carbons within propane during pyrolysis. The isotope fractionations during the pyrolysis are reproducible at optimum conditions, but vary between the two GC-Py-GC-IRMS systems tested, affected by experimental conditions (e.g., pyrolysis temperature, flow rate, and reactor conditions). Conclusions: It is necessary to evaluate and calibrate each GC-Py-GC-IRMS system using propane gases with accurately determined PS δ 13 C values. Finally, this study also highlights a need for PS isotope standards for propane and other molecules (e.g., butane and acetic acid).

03 NATURAL GAS↗