Search NASA⌕ Search

SEARCH · Search NASA

Results for “control volume methods”

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 379 records · Page 21
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

System and method for in situ volumetric sensing of 3D cure state of resin being used in an additive manufacturing system

An in situ cure monitoring control system for use with a volumetric three dimensional (3D) printing system, wherein the volumetric 3D system includes a container defining a build volume, the container includes a photocurable resist used for making a 3D part, and the resist is responsive to an optical curing beam from a light source, which is controlled by a controller, and which is passed through the resist. The cure monitoring control system includes an optical signal source which generates optical signals having a wavelength selected in relation to a characteristic of the resist, and directed to pass through the build volume. A detector detects the optical signals and generates output signals in accordance therewith. Software monitors the output signals and uses the output signals to modify the curing beam to help optimize curing of the resist.

Panas, Robert Matthew↗
DOE OSTI2018

Controlling anisotropy in stereolithographically printed polymers

The ability to print three-dimensional objects was first developed in the 1980s and was originally strictly limited to polymeric materials. Through most of the intervening years, the approach has been thought of as a rapid prototyping method. This allowed low volume, high fidelity structures to be quickly fabricated to test things like fit and finish. However, more recently the term additive manufacturing has entered usage to represent the same methods and implies that the field is transitioning to creating finished parts. When we begin to think about these techniques as a real manufacturing approach, then the materials properties achieved during manufacture become much more important. This is true for all materials, as additive manufacturing has grown to encompass structures formed from metals and ceramics, however the bulk of work remains centered on polymers.

36 MATERIALS SCIENCE↗
NASA NTRS2003

Baroreflex-mediated heart rate and vascular resistance responses 24 h after maximal exercise

INTRODUCTION: Plasma volume, heart rate (HR) variability, and stimulus-response relationships for baroreflex control of forearm vascular resistance (FVR) and HR were studied in eight healthy men after and without performing a bout of maximal exercise to test the hypotheses that acute expansion of plasma volume is associated with 1) reduction in baroreflex-mediated HR response, and 2) altered operational range for central venous pressure (CVP). METHODS: The relationship between stimulus (DeltaCVP) and vasoconstrictive reflex response (DeltaFVR) during unloading of cardiopulmonary baroreceptors was assessed with lower-body negative pressure (LBNP, 0, -5, -10, -15, -20 mm Hg). The relationship between stimulus (Deltamean arterial pressure (MAP)) and cardiac reflex response (DeltaHR) during loading of arterial baroreceptors was assessed with steady-state infusion of phenylephrine (PE) designed to increase MAP by 15 mm Hg alone and during application of LBNP (PE+LBNP) and neck pressure (PE+LBNP+NP). Measurements of vascular volume and autonomic baroreflex responses were conducted on two different test days, each separated by at least 1 wk. On one day, baroreflex response was tested 24 h after graded cycle exercise to volitional exhaustion. On another day, measurement of baroreflex response was repeated with no exercise (control). The order of exercise and control treatments was counterbalanced. RESULTS: Baseline CVP was elevated (P = 0.04) from a control value of 10.5 +/- 0.4 to 12.3 +/- 0.4 mm Hg 24 h after exercise. Average DeltaFVR/DeltaCVP during LBNP was not different (P = 0.942) between the exercise (-1.35 +/- 0.32 pru x mm Hg-1) and control (-1.32 +/- 0.36 pru x mm Hg-1) conditions. However, maximal exercise caused a shift along the reflex response relationship to a higher CVP and lower FVR. HR baroreflex response (DeltaHR/DeltaMAP) to PE+LBNP+NP was lower (P = 0.015) after maximal exercise (-0.43 +/- 0.15 beats x min-1 x mm Hg-1) compared with the control condition (-0.83 +/- 0.14 beats x min-1 x mm Hg-1). CONCLUSION: Expansion of vascular volume after acute exercise is associated with altered operational range for CVP and reduced HR response to arterial baroreceptor stimulation.

Non-NASA Center↗
DOE OSTI2022

Measuring representative volume elements from high-resolution grain-scale strain fields

Most crystalline materials present a highly heterogeneous response at the microscale, which can be affected by both internal factors (such as microstructural parameters) and external factors (such as loading). Relating microscale inhomogeneities to the macroscale response of a material requires the use of homogenisation techniques, usually based on the concept of a representative volume element (RVE)—the smallest volume of material that represents the global average response. Here, in this work, we present a new and robust experimental method of measuring the size of a strain-based RVE from high-resolution grain-scale strain fields obtained using digital image correlation (DIC). The proposed method is based on the statistical (stereological) nature of the RVE, which has been widely adopted in numerical studies, and involves dividing a strain field into randomly selected regions of varying sizes and statistically analysing the distributions of average strains within them. To validate the new method, we generate a large number of synthetic strain fields from a fractional Gaussian noise algorithm. The proposed stereological method is shown to be capable of producing reliable RVE measurements from a very large range of possible microscale strain fields while at the same time being robust in that it can produce RVE measurement results even in cases where other existing methods may be unable to do so. The proposed method has a low field-of-view requirement, only needing a field-of-view about 1.2 times as large as the RVE to produce reliable measurements. In addition, the stereological method offers significant flexibility since its statistical nature allows for control over how strict the RVE measurement should be in each case.

36 MATERIALS SCIENCE↗
NASA NTRS1997

System identification of closed-loop cardiovascular control: effects of posture and autonomic blockade

We applied system identification to the analysis of fluctuations in heart rate (HR), arterial blood pressure (ABP), and instantaneous lung volume (ILV) to characterize quantitatively the physiological mechanisms responsible for the couplings between these variables. We characterized two autonomically mediated coupling mechanisms [the heart rate baroreflex (HR baroreflex) and respiratory sinus arrhythmia (ILV-HR)] and two mechanically mediated coupling mechanisms [the blood pressure wavelet generated with each cardiac contraction (circulatory mechanics) and the direct mechanical effects of respiration on blood pressure (ILV-->ABP)]. We evaluated the method in humans studied in the supine and standing postures under control conditions and under conditions of beta-sympathetic and parasympathetic pharmacological blockades. Combined beta-sympathetic and parasympathetic blockade abolished the autonomically mediated couplings while preserving the mechanically mediated coupling. Selective autonomic blockade and postural changes also altered the couplings in a manner consistent with known physiological mechanisms. System identification is an "inverse-modeling" technique that provides a means for creating a closed-loop model of cardiovascular regulation for an individual subject without altering the underlying physiological control mechanisms.

NASA Discipline Number 18-10↗
DOE OSTI2022

Comprehensive single-shot biophysical cytometry using simultaneous quantitative phase imaging and Brillouin spectroscopy

Abstract Single-cell analysis, or cytometry, is a ubiquitous tool in the biomedical sciences. Whereas most cytometers use fluorescent probes to ascertain the presence or absence of targeted molecules, biophysical parameters such as the cell density, refractive index, and viscosity are difficult to obtain. In this work, we combine two complementary techniques—quantitative phase imaging and Brillouin spectroscopy—into a label-free image cytometry platform capable of measuring more than a dozen biophysical properties of individual cells simultaneously. Using a geometric simplification linked to freshly plated cells, we can acquire the cellular diameter, volume, refractive index, mass density, non-aqueous mass, fluid volume, dry volume, the fractional water content of cells, both by mass and by volume, the Brillouin shift, Brillouin linewidth, longitudinal modulus, longitudinal viscosity, the loss modulus, and the loss tangent, all from a single acquisition, and with no assumptions of underlying parameters. Our methods are validated across three cell populations, including a control population of CHO-K1 cells, cells exposed to tubulin-disrupting nocodazole, and cells under hypoosmotic shock. Our system will unlock new avenues of research in biophysics, cell biology, and medicine.

Steelman, Zachary A.↗
DOE OSTI2016

High density plasmas and new diagnostics: An overview (invited)

One of the limiting factors for the full understanding of Electron Cyclotron Resonance Ion Sources (ECRISs) fundamental mechanisms consists of few types of diagnostic tools so far available for such compact machines. Microwave-to-plasma coupling optimisation, new methods of density overboost provided by plasma wave generation, and magnetostatic field tailoring for generating a proper electron energy distribution function, suitable for optimal ion beams formation, require diagnostic tools spanning across the entire electromagnetic spectrum from microwave interferometry to X-ray spectroscopy; these methods are going to be implemented including high resolution and spatially resolved X-ray spectroscopy made by quasi-optical methods (pin-hole cameras). The ion confinement optimisation also requires a complete control of cold electrons displacement, which can be performed by optical emission spectroscopy. Several diagnostic tools have been recently developed at INFN-LNS, including “volume-integrated” X-ray spectroscopy in low energy domain (2-30 keV, by using silicon drift detectors) or high energy regime (>30 keV, by using high purity germanium detectors). For the direct detection of the spatially resolved spectral distribution of X-rays produced by the electronic motion, a “pin-hole camera” has been developed also taking profit from previous experiences in the ECRIS field. The paper will give an overview of INFN-LNS strategy in terms of new microwave-to-plasma coupling schemes and advanced diagnostics supporting the design of new ion sources and for optimizing the performances of the existing ones, with the goal of a microwave-absorption oriented design of future machines.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2024

A modern concept of Lagrangian hydrodynamics

Here, we offer a modern interpretation of Lagrangian hydrodynamics as employed in Lagrangian simulations of compressible fluid flow. Our main result is to show that artificial viscosity, traditionally viewed as a numerical artifice to control unphysical oscillations in flows with shocks, actually represents a physical process and is necessary to derive accurate simulations in any compressible flow. We begin by reviewing the origins of two numerical devices, artificial viscosity and finite-volume methods. We proceed to construct a mathematical (PDE) model that incorporates those numerics and in which a new length scale, the observer, arises representing the discretization. Associated with that length scale, there are new inviscid fluxes that are the artificial viscosity as first formulated by Richtmyer and an artificial heat flux postulated by Noh but typically not included in Lagrangian codes. We discuss the connection of our results to bivelocity hydrodynamics. We conclude with some speculation as to the direction of future developments in multidimensional Lagrangian codes as computers get faster and have larger memories.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1992

Active control of asymmetric vortical flows around cones using injection and heating

The effectiveness of certain active-control methods for asymmetric flows around circular cones is investigated by using computational solution of the unsteady, compressible full Navier-Stokes equations. Two main methods of active control which include flow injection and surface heating are used. For the flow-injection-control method, flow injection is used either in the normal direction to the surface or in the tangential direction to the surface. For the surface-heating-control method, the temperature of the cone surface is increased. The effectiveness of a hybrid method of flow control which combines normal injection with surface heating has also been studied. The Navier-Stokes equations, subjected to various surface boundary conditions, are solved by using an implicit, upwind, flux-difference splitting, finite-volume scheme for locally-conical flow solutions.

Kandil, Osama A.↗
NASA NTRS1992

Design of an essentially non-oscillatory reconstruction procedure in finite-element type meshes

An essentially non oscillatory reconstruction for functions defined on finite element type meshes is designed. Two related problems are studied: the interpolation of possibly unsmooth multivariate functions on arbitary meshes and the reconstruction of a function from its averages in the control volumes surrounding the nodes of the mesh. Concerning the first problem, the behavior of the highest coefficients of two polynomial interpolations of a function that may admit discontinuities of locally regular curves is studied: the Lagrange interpolation and an approximation such that the mean of the polynomial on any control volume is equal to that of the function to be approximated. This enables the best stencil for the approximation to be chosen. The choice of the smallest possible number of stencils is addressed. Concerning the reconstruction problem, two methods were studied: one based on an adaptation of the so called reconstruction via deconvolution method to irregular meshes and one that lies on the approximation on the mean as defined above. The first method is conservative up to a quadrature formula and the second one is exactly conservative. The two methods have the expected order of accuracy, but the second one is much less expensive than the first one. Some numerical examples are given which demonstrate the efficiency of the reconstruction.

Abgrall, Remi↗
NASA NTRS2006

Utilization of Porous Media for Condensing Heat Exchangers

The use of porous media as a mean of separating liquid condensate from the air stream in condensing heat exchangers has been explored in the past inside small plant growth chambers and in the Apollo Command Module. Both applications used a cooled porous media made of sintered stainless steel to cool and separate condensation from the air stream. However, the main issues with the utilization of porous media in the past have been the deterioration of the porous media over long duration, such as clogging and changes in surface wetting characteristics. In addition, for long duration usage, biofilm growth from microorganisms on the porous medial would also be an issue. In developing Porous Media Condensing Heat Exchangers (PMCHX) for future space applications, different porous materials and microbial growth control methods will need to be explored. This paper explores the work performed at JSC and GRC to evaluate different porous materials and microbial control methods to support the development of a Porous Media Condensing Heat Exchanger. It outlines the basic principles for designing a PMCHX and issues that were encountered and ways to resolve those issues. The PMCHX has potential of mass, volume, and power savings over current CHX and water separator technology and would be beneficial for long duration space missions.

Tuan, George C.↗
NASA NTRS2010

NASA Tech Briefs, November 2010

Topics covered include: Portable Handheld Optical Window Inspection Device; Salience Assignment for Multiple-Instance Data and Its Application to Crop Yield Prediction; Speech Acquisition and Automatic Speech Recognition for Integrated Spacesuit Audio Systems ; Predicting Long-Range Traversability from Short-Range Stereo-Derived Geometry; Browser-Based Application for Telemetry Monitoring of Robotic Assets; Miniature Low-Noise G-Band I-Q Receiver; Methods of Using a Magnetic Field Response Sensor Within Closed, Electrically Conductive Containers; Differential Resonant Ring YIG Tuned Oscillator; Microfabricated Segmented-Involute-Foil Regenerator for Stirling Engines; Reducing Seal Adhesion in Low Impact Docking Systems; Optimal Flow Control Design; Corrosion-Resistant Container for Molten-Material Processing; Reusable Hot-Wire Cable Cutter; Deployment of a Curved Truss; High-Volume Airborne Fluids Handling Technologies to Fight Wildfires; Modeling of Alkane Oxidation Using Constituents and Species; Fabrication of Lanthanum Telluride 14-1-11 Zintl High-Temperature Thermoelectric Couple; A Computer Model for Analyzing Volatile Removal Assembly; Analysis of Nozzle Jet Plume Effects on Sonic Boom Signature; Optical Sidebands Multiplier; Single Spatial-Mode Room-Temperature-Operated 3.0 to 3.4 micrometer Diode Lasers; Self-Nulling Beam Combiner Using No External Phase Inverter; Portable Dew Point Mass Spectrometry System for Real-Time Gas and Moisture Analysis; Maximum Likelihood Time-of-Arrival Estimation of Optical Pulses via Photon-Counting Photodetectors; Handheld White Light Interferometer for Measuring Defect Depth in Windows; Decomposition Algorithm for Global Reachability on a Time-Varying Graph; Autonomous GN and C for Spacecraft Exploration of Comets and Asteroids; Efficient Web Services Policy Combination; Using CTX Image Features to Predict HiRISE-Equivalent Rock Density; Isolation of the Paenibacillus phoenicis, a Spore-Forming Bacterium; Monolithically Integrated, Mechanically Resilient Carbon-Based Probes for Scanning Probe Microscopy; Cell Radiation Experiment System; Process to Produce Iron Nanoparticle Lunar Dust Simulant Composite; Inversion Method for Early Detection of ARES-1 Case Breach Failure; Use of ILTV Control Laws for LaNCETS Flight Research;and Evaluating Descent and Ascent Trajectories Near Non-Spherical Bodies.

Source record↗
NASA NTRS1979

Dynamic loads analysis system (DYLOFLEX) summary. Volume 1: Engineering formulation

The DYLOFLEX computer program system expands the aeroelastic cycle from that in the FLEXSTAB computer program system to include dynamic loads analyses involving active controls. Two aerodynamic options exist within DYLOFLEX. The analyst can formulate the problem with unsteady aerodynamics calculated using the doublet lattice method or with quasi-steady aerodynamics formulated from either FLEXSTAB or doublet lattice steady state aerodynamics with unsteady effects approximated by indicial lift growth functions. The equations of motion are formulated assuming straight and level flight and small motions. Loads are calculated using the force summation technique. DYLOFLEX consists of nine standalone programs which can be linked with each other by magnetic files used to transmit the required data between programs.

Miller, R. D.↗
DOE OSTI2015

High performance anode for advanced Li batteries

The overall objective of this Phase I SBIR effort was to advance the manufacturing technology for ASI’s Si-CNF high-performance anode by creating a framework for large volume production and utilization of low-cost Si-coated carbon nanofibers (Si-CNF) for the battery industry. This project explores the use of nano-structured silicon which is deposited on a nano-scale carbon filament to achieve the benefits of high cycle life and high charge capacity without the consequent fading of, or failure in the capacity resulting from stress-induced fracturing of the Si particles and de-coupling from the electrode. ASI’s patented coating process distinguishes itself from others, in that it is highly reproducible, readily scalable and results in a Si-CNF composite structure containing 25-30% silicon, with a compositionally graded interface at the Si-CNF interface that significantly improve cycling stability and enhances adhesion of silicon to the carbon fiber support. In Phase I, the team demonstrated the production of the Si-CNF anode material can successfully be transitioned from a static bench-scale reactor into a fluidized bed reactor. In addition, ASI made significant progress in the development of low cost, quick testing methods which can be performed on silicon coated CNFs as a means of quality control. To date, weight change, density, and cycling performance were the key metrics used to validate the high performance anode material. Under this effort, ASI made strides to establish a quality control protocol for the large volume production of Si-CNFs and has identified several key technical thrusts for future work. Using the results of this Phase I effort as a foundation, ASI has defined a path forward to commercialize and deliver high volume and low-cost production of SI-CNF material for anodes in Li-ion batteries.

25 ENERGY STORAGE↗
NASA NTRSDate not supplied

Development and Validation of a Slosh-Based Ullage Collapse Model

Ullage collapse modeling is an essential tool for the success of spacecraft powered by cryogenic propellants. The degree of ullage collapse can have a major impact on the design of cryogenic propellant tanks as well as the mission’s concept of operations. Additional commodities (e.g., helium gas for pressure control) and their quantities are dependent on the amount of heat transfer and phase change of propellants expected during flight. Currently, high-fidelity computational fluid dynamics (CFD) with heat transfer and phase change is the most accurate method for predicting ullage collapse. CFD simulations of this kind, however, can become very costly for problems requiring large domains, long simulation times, or significant liquid motion. Consequently, there is a need for a reduced order model to provide fast responses which can be used in engineering design decisions. The proposed numerical model is a nodal approach, augmented by outputs from CFD analyses that are generated rapidly and reliably. These outputs include transient values for liquid-vapor interface area, interface velocity, and liquid position using a volume of fluid (VoF) method. There are several benefits of combining a nodal approach with slosh predictions from CFD simulations: leveraging of the Loci/STREAM-VoF solver that has been validated in high-g and low-g environments, rapid modeling of VoF with isothermal propellant simplifications and relatively coarse meshes, and improved model predictions from the nodal tool. The nodal approach, informed by the slosh predictions from CFD, assumes uniform conditions in the liquid and gas domains. Modifications to account for temperature stratification in the liquid and gas may be implemented as a future improvement. Thermodynamic changes within the tank are solved at user-defined timesteps. First, impacts due to mass entering or leaving the control volume are captured. Second, heat transfer impacts are accounted for by tracking the volume of the gas which has been contacted by liquid in that timestep and rebalancing the internal energy. Finally, phase change is approximated using the area and velocity of the interface. The tank conditions are re-computed by conserving mass and energy. Calibration of the model to several test flights and experiments are required to baseline the predictive power and uncertainty. The final product will enable designers and engineers to create quick-turnaround predictions for ullage collapse and inform design feasibility.

Computational Fluid Dynamics↗
DOE OSTI2021

Modelling the complex evaporated gas flow and its impact on particle spattering during laser powder bed fusion

The additive manufacturing (AM) of metals is becoming an increasingly important production process with the potential to replace traditional techniques such as casting. Laser Powder Bed Fusion (LPBF) is used in many applications to print metal parts from powder. The metal powder is heated locally with sufficient laser radiation that the liquid melt easily reaches its boiling temperature, which leads to a metallic vapour jet that can entrain both powder bed particles and molten droplets. The small size of laser-matter interaction site makes a detailed experimental analysis of the process challenging. Synchrotron X-ray imaging experiments are one of the few methods which can capture the dynamic melting and solidification processes. Comparing such experiments with computer simulations of the process is an important approach in order to better understand the manufacturing process and to analyse the influence of process parameters on the evaporated gas jet and the subsequent impact on particle ejection, leading to potentially reduced AM component quality. The melting and solidification of the metal powder is simulated using an Eulerian multiphase approach based on a control volume discretization of powder bed and substrate and a volume of liquid separation from melt and gas phase. The gas phase modelled as an ideal gas reaches velocities up to 100 m/s. Lagrangian particle tracking in the simulation demonstrates that the velocity fields calculated by the Eulerian multi-phase approach in combination with a standard drag-force model lead to particle accelerations in good agreement with those measured experimentally. In order to avoid numerical laborious Lagrangian calculations, a direct method to compare an Eulerian multiphase simulation with synchrotron X-ray experiments was introduced and validated. This approach is used to analyse the influence of process parameters including laser power and laser speed on the maximal acceleration of particles from the melt pool area. While the particle acceleration increases linearly with line energy in the conduction mode, a linear decrease of the acceleration with increasing line energy can be found in the transition mode before the accel-eration increases again with line energy in the keyhole mode.

36 MATERIALS SCIENCE↗
DOE OSTI2018

Improved electrochemical cycling stability of intercalation battery electrodes via control of material morphology

Using a model tunnel manganese oxide with the todorokite crystal structure (T-MnO 2 ), we reveal that controlling the morphology of the active material can improve the cycling stability of intercalation battery electrodes. The T-MnO 2 structure is built from tunnels that provide spacious 1D diffusion channels for charge-carrying ions. Taking advantage of the unique ability to synthesize T-MnO 2 in the form of both highly crystalline two-dimensional (2D) nanoplatelets and one-dimensional (1D) nanowires through a facile hydrothermal growth method, we investigated the effect of nanoscale particle dimensions on reversible battery cycling. Insertion of ions into the tunnels results in anisotropic expansion of the structure, making T-MnO 2 with different morphologies an excellent model platform to understand how intercalation-induced volume change, typically leading to the deterioration of the electrode performance over extended cycling, can be controlled through synthesis of targeted morphologies. T-MnO 2 nanowires showed not only significantly improved capacity retention but also substantially higher specific capacity than the T-MnO 2 nanoplatelets. The enhanced electrochemical properties of the nanowire electrodes could be attributed to the larger surface-to-volume ratio than that of nanoplatelets, resulting in higher contact area with electrolyte for the nanowires. Moreover, due to the smaller cross-sectional area of the nanowires, volume expansion and contraction perpendicular to the structural tunnels induced by reversible ion intercalation occurs in a more facile fashion. In this research we demonstrate that chemically controlling morphology and producing particles with nanostructure dimensionality replicating that of atomic structure (i.e., 1D morphology and 1D structure) makes it possible to enhance material performance.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2020

An experimental system and procedure of unsteady-state relative permeability test for gas hydrate-bearing sediments

Reliable estimations of the relative permeability of gas and water in hydrate-bearing sediments (HBS) and the dependency of the relative permeability on hydrate saturation are critical to predict the productivity of a hydrate reservoir. Yet, this remains poorly estimated owing to lack of experimental data associated with difficulties in conducting multiphase flow experiments in HBS. Recognizing the experimental challenges, this study intends to develop and validate a new experimental system and procedure of unsteady-state relative permeability test that can generate reliable and reproducible flow measurements in HBS. Furthermore, gas hydrate is considered as a part of solid matrix in the sediment, so one of the challenges is to maintain a constant hydrate saturation, which is achieved in this experimental study using tight pressure-temperature (P-T) control near the hydrate stability boundary. The measured differential pressure across the specimen, methane injection flow rate, and volume of displaced brine are used to calculate the relative permeability by adopting a conventional Buckley-Leverett theory-based interpretation method. Residual brine saturation calculated for the hydrate-bearing specimen is higher than that of hydrate-free specimen, presumably due to decrease in pore size, increase in heterogeneity of solid matrix, and increase in size distribution of solid matrix and pore in the presence of hydrates. Further studies are necessary to represent the results of the unsteady-state flow experiment in HBS with a gas hydrate-dependent relative permeability model.

03 NATURAL GAS↗