Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Dynamic equivalent”

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 847 records · Page 47

Three-Dimensional, Transgenic Cell Models to Quantify Space Genotoxic Effects

The space environment contains radiation and chemical agents known to be mutagenic and carcinogenic to humans. Additionally, microgravity is a complicating factor that may modify or synergize induced genotoxic effects. Most in vitro models fail to use human cells (making risk extrapolation to humans more difficult), overlook the dynamic effect of tissue intercellular interactions on genotoxic damage, and lack the sensitivity required to measure low-dose effects. Currently a need exists for a model test system that simulates cellular interactions present in tissue, and can be used to quantify genotoxic damage induced by low levels of radiation and chemicals, and extrapolate assessed risk to humans. A state-of-the-art, three-dimensional, multicellular tissue equivalent cell culture model will be presented. It consists of mammalian cells genetically engineered to contain multiple copies of defined target genes for genotoxic assessment,. NASA-designed bioreactors were used to coculture mammalian cells into spheroids, The cells used were human mammary epithelial cells (H184135) and Stratagene's (Austin, Texas) Big Blue(TM) Rat 2 lambda fibroblasts. The fibroblasts were genetically engineered to contain -a high-density target gene for mutagenesis (60 copies of lacl/LacZ per cell). Tissue equivalent spheroids were routinely produced by inoculation of 2 to 7 X 10(exp 5) fibroblasts with Cytodex 3 beads (150 micrometers in diameter). at a 20:1 cell:bead ratio, into 50-ml HARV bioreactors (Synthecon, Inc.). Fibroblasts were cultured for 5 days, an equivalent number of epithelial cells added, and the fibroblast/epithelial cell coculture continued for 21 days. Three-dimensional spheroids with diameters ranging from 400 to 600 micrometers were obtained. Histological and immunohistochemical Characterization revealed i) both cell types present in the spheroids, with fibroblasts located primarily in the center, surrounded by epithelial cells; ii) synthesis of extracellular matrix; and iii,, mitotic cells located throughout the spheroids. Spheroidal integrity and cell viability were retained for the 30-day test period after removal of spheroids from the bioreactor. Potential utility of this three-dimensional, transgenic model for genotoxicity was initially assessed by exposure of spheroids to 0-2 Gy neon at dose rates of 0.3 to 1.5 Gy/min (National Institute of Radiological Sciences, Chiba, Japan). Quantification of mutation at the lacl gene revealed a linear dose response for mutation induction. Limited sequencing analysis of mutant clones revealed higher frequencies of deletions and multiple base sequence changes with increasing dose. These results suggest that our three-dimensional, transgenic model is applicable to a wide variety of studies involving the quantification, identification, and characterization of genotoxicity incurred in space and on Earth. This model uniquely allows investigation of the interaction of relevant factors, namely cell-to-cell interactions and the mechanistic interaction of microgravity with radiation insults and DNA repair. Using this three-dimensional model will allow us to obtain dual genotoxic information (i.e., mutation rate plus chromosome aberration data) from the same system so that one endpoint can be used to reference the other, thereby increasing the fidelity of the data set. Moreover, the tissue-equivalent nature of the three-dimensional model provides high confidence for relevance of risk assessment, i.e., the establishment of quality factors directly applicable to the microgravity environment.

Gonda, S. R.↗

Chemical Reactor Network Modeling of Ammonia Rich-Quench-Lean Combustion Using a Partially Stirred Reactor Approach

Ammonia is a promising alternative to hydrogen with high energy density and favorable storage and transport characteristics. However, low flammability and a propensity for high nitrogen oxide (NO x ) emissions make direct utilization challenging. Recently, two-stage rich-quench-lean (RQL) combustion strategies have shown promise in achieving low NO x emissions with ammonia. In this approach, the rich stage serves to oxidize a portion of the fuel while thermally decomposing as much of the remaining ammonia as possible, generating hydrogen. In the second (lean) stage, air is rapidly introduced, burning out the hydrogen and residual ammonia. Two-stage RQL combustion of ammonia has been investigated in the open literature both experimentally and numerically. In general, idealized chemical reactor network (CRN) models predict NO x concentrations below those of 2D/3D computational fluid dynamics models and experiments. The primary drivers of these discrepancies may be largely attributed to finite rate mixing nonadiabatic operation. The typical CRN model is comprised of a perfectly-stirred-reactor (PSR), followed by a plug-flow-reactor (PFR), meant to represent the flame, and postflame zones, respectively. In the two-stage RQL approach two PSR-PFR networks are arranged sequentially, corresponding to the rich and lean stages, with secondary air injection in between. In the authors' past work, this arrangement has demonstrated the significant sensitivity of exit NO x to the rich stage equivalence ratio, while the amount of secondary air injection was shown to be less critical. In this paper, the CRN model is extended to (1) include the impacts of heat loss and (2) utilize a partially-stirred-reactor (PaSR) approach to study the impacts of mixing on emissions performance. Varying amounts of heat loss are applied to the rich relaxation zone to understand emissions performance and changes to optimization of equivalence ratio and residence time. Premixed and nonpremixed configurations are considered in the rich stage PaSR, with varying degrees of mixing intensity to study the interaction between mixing, transport, and kinetic timescales. Critically, the impact of mixing between hot products and secondary air injection is studied to understand practical injector needs. Results show unburnt ammonia leaving the rich stage as a primary contributor to NO x emissions – driven both by increased heat loss and reduced mixing rates. Furthermore, heat losses have been shown to create conditions that are conducive to increased N 2 O formation in the lean stage. In conclusion, the results of this study will be considered in the context of developing optimized two-stage RQL combustors for ammonia.

Combustion↗

Altering the Supply of Shielding Gases to Fabricate Distinct Geometry in GMA Additive Manufacturing

Wire arc additive manufacturing (WAAM) is the process by which large, metallic structures are built, layer-by-layer, using a welding arc to melt wire feedstock. In this process, the proper selection of the shielding gas plays a vital role in the achievement of structurally acceptable part geometries and quality surface finishes. In this study, the authors used either a ternary mix (He, Ar and CO 2 ) or a binary mix (Ar and CO 2 ) of shielding gases to deposit wall geometries using an open loop-controlled WAAM system developed at Oak Ridge National Laboratory’s Manufacturing Demonstration Facility. The binary blend produced a wider and shorter geometry, while the ternary blend resulted in a narrower build that was more equivalent to the CAD geometry. The data indicated that the binary blend provided a higher oxygen concentration in the weld as compared to that of the ternary blend. The results imply that the arc characteristics and heat input had a significantly higher impact on the weld penetration than the surface tension effect of surface active elements. This was further verified by developing and applying a high-fidelity computational fluid dynamics (CFD) model of the thermophysical properties of gas mixtures. The results from the model showed that, while the influence of increased oxygen concentration on the surface tension for the binary blend led to a deeper penetration, the ternary blend gave rise to heat flux to the workpiece.

36 MATERIALS SCIENCE↗

OH PLIF Visualization of the UVa Supersonic Combustion Experiment: Configuration C

Non-intrusive hydroxyl radical (OH) planar laser-induced fluorescence (PLIF) measurements were obtained in configuration C of the University of Virginia supersonic combustion experiment. The combustion of hydrogen fuel injected through an unswept compression ramp into a supersonic cross-flow was imaged over a range of streamwise positions. Images were corrected for optical distortion, variations in the laser sheet profile, and different camera views. Results indicate an effect of fuel equivalence ratio on combustion zone shape and local turbulence length scale. The streamwise location of the reaction zone relative to the fuel injector was also found to be sensitive to the fuel equivalence ratio. The flow boundary conditions in the combustor section, which are sensitive to the fuel flow rate, are believed to have caused this effect. A combination of laser absorption and radiative trapping effects are proposed to have caused asymmetry observed in the images. The results complement previously published OH PLIF data obtained for configuration A along with other non-intrusive measurements to form a database for computational fluid dynamics (CFD) model validation.

McRae, Colin D.↗

Self-Aligned Focusing Schlieren and OH Planar Laser-Induced Fluorescence Flow Visualization in a Dual-Mode Scramjet

Ahigh-speed self-aligned focusing schlieren (SAFS) system was used to visualize density gradients in and around the cavity flameholder of the combustor section of the University of Virginia Supersonic Combustion Facility (UVASCF). Images with this system were acquired at a framing rate of 110 kHz with no fuel injection, with fuel injection but no flame, and for fuel injection with combustion corresponding to a global equivalence ratio of 𝜙 = 0.18. Images with an air throttle in operation to modify the shock train location with fuel injection and with flame were also acquired. Simultaneous OH planar laser-induced fluorescence (PLIF) images were also acquired at a framing rate of 20 Hz. Results obtained with both visualization techniques are compared to one another to highlight how SAFS can complement more advanced flow visualization techniques and resolve dynamic behavior that may not otherwise be captured. Both proper orthogonal decomposition (POD) and dynamic mode decomposition (DMD) analysis techniques are applied to the SAFS image sequences to identify coherent periodic structures for the runs with fuel injection and combustion.

Brett F Bathel↗

Effect of hollow fiber membrane packing on the performance of modules formed with fiber tows

Membrane processes are the preferred option for many gas separations. Hollow fiber membranes are commonly used in these processes. Commercial hollow fiber membrane modules are often fabricated from small groups of fibers (i.e., a tow) that are arranged into a larger cylindrical bundle. The existence of fiber tows can lead to packing non-uniformity. This nonuniformity can result in nonideal fluid distribution in the shell and be detrimental to module performance. In this work, full three-dimensional models of flow and mass transfer within fiber bundles formed from tows have been developed. These models account for the detailed fiber arrangement within modules and are solved to evaluate module separation performance in terms of stage cut and required membrane area as a function product gas purity for carbon capture applications. Results are presented illustrating the effects of varying fiber packing arrangements within a tow and varying tow packing arrangements. Further, the results suggest smaller tows are preferred over larger tows, at fixed overall fiber packing and lower intra- and inter-tow packings are preferred. To reduce the computational cost of simulating large fiber bundles, the use of an equivalent planar bundle is proposed (EPB). A simple procedure for determining the geometry of the EPB is described and validated through comparisons of calculated performance metrics with results obtained from full three-dimensional simulations.

42 ENGINEERING↗

Three Jahn-Teller States of Matter in Spin-Crossover System Mn(taa)

Three high-spin phases recently discovered in the spin-crossover system Mn(taa) are identified through analysis by a combination of first-principles calculations and Monte Carlo simulation as a low-temperature Jahn-Teller ordered (solid) phase, an intermediate-temperature dynamically correlated (liquid) phase, and an uncorrelated (gas) phase. In particular, the Jahn-Teller liquid phase arises from competition between mixing with low-spin impurities, which drive the disorder, and intermolecular strain interactions. Additionally, the latter are a key factor in both the spin-crossover phase transition and the magnetoelectric coupling. Jahn-Teller liquids may exist in other spin-crossover materials and materials that have multiple equivalent Jahn-Teller axes.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Testing of a one-bladed 30-meter-diameter rotor on the DOE/NASA Mod-O wind turbine

Tests were conducted on the DOE/NASA Mod-O 200-kW horizontal-axis wind turbine in a one-bladed rotor configuration. The objectives of the test were to evaluate the performance, loads, and dynamic characteristics of a one-bladed rotor, and then to compare these parameters with those of an aerodynamically similar two-bladed rotor configuration. Test operations showed that this intermediate-size (15.2-m radius) one-bladed rotor configuration can be operated successfully. Test results show that the one-bladed rotor had cyclic blade loads comparable to those of a two-bladed rotor. A moderate power penalty equivalent to a reduction in windspeed of 1 m/sec occurred with the one-bladed rotor when operated at a rotor speed 50 percent higher than that of the two-bladed rotor.

Ensworth, C. B. F., III↗

Flexible robot control: Modeling and experiments

Described here is a model and its use in experimental studies of flexible manipulators. The analytical model uses the equivalent of Rayleigh's method to approximate the displaced shape of a flexible link as the static elastic displacement which would occur under end rotations as applied at the joints. The generalized coordinates are thereby expressly compatible with joint motions and rotations in serial link manipulators, because the amplitude variables are simply the end rotations between the flexible link and the chord connecting the end points. The equations for the system dynamics are quite simple and can readily be formulated for the multi-link, three-dimensional case. When the flexible links possess mass and (polar moment of) inertia which are small compared to the concentrated mass and inertia at the joints, the analytical model is exact and displays the additional advantage of reduction in system dimension for the governing equations. Four series of pilot tests have been completed. Studies on a planar single-link system were conducted at Carnegie-Mellon University, and tests conducted at Toshiba Corporation on a planar two-link system were then incorporated into the study. A single link system under three-dimensional motion, displaying biaxial flexure, was then tested at Carnegie-Mellon.

Oppenheim, Irving J.↗

Testing for Random Limit Load Versus Static Limit Load

This document is an effort to report the basic test findings in an ongoing quest for understanding how random load factors should be applied to structural components in order to verify the strength of space flight hardware. A Spacelab experiment known as the Atmospheric Emission Photometric Imager (AEPI) was subjected to both an expected flight random environment and the associated Miles' equation equivalent static load. During each of these tests, the fiberglass pedestal was instrumented with 16 triaxial strain gauges around its base. Component strains and invariant stresses were compared. As seen previously in other hardware tests, the stress distribution from the random environment was an order of magnitude below the comparable static stresses. With a proposed data acquisition system, a strain database will be developed that will quantify an empirical relationship between dynamic and static limit stresses. This event will allow a more accurate estimate of launch environment effects on new technology structural components.

Lee, H. M.↗

Synthesizing Monitors for Safety Properties

The problem of testing a linear temporal logic (LTL) formula on a finite execution trace of events, generated by an executing program, occurs naturally in time analysis of software. An algorithm which takes a past time LTL formula and generates an efficient dynamic programming algorithm is presented. The generated algorithm tests whether the formula is satisfied by a finite trace of events given as input and runs in linear time, its constant depending on the size of the LTL formula. The memory needed is constant, also depending on the size of the formula. Further optimizations of the algorithm are suggested. Past time operators suitable for writing succinct specifications axe introduced and shown definitionally equivalent to the standard operators. This work is part of the PathExplorer project, the objective of which it is to construct a flexible framework for monitoring and analyzing program executions.

Havelund, Klaus↗

Thermal energy evolution and mechanical deformation of monocrystalline yttria-stabilized zirconia nanoparticles in aerosol deposition processes

Aerosol deposition (AD) is a coating process wherein aerosol particles are impacted on a target substrate. There are fundamental differences between the AD process (cold impact), where particle translational kinetic energy is high and thermal energy is low, and thermal spray deposition (thermal impact), where translational energy is lower but thermal energy is high. To better compare cold and thermal impact effects on particles, we carried out molecular dynamics simulations using yttria-stabilized zirconia (YSZ) nanoparticles on YSZ substrates as a model system. We performed cold impact simulations at 300 K with variable impact velocity in the 500 ms –1 –1500 ms –1 range to understand how increasing translational kinetic energy affects thermal energy and mechanical evolution. We then performed thermal impact simulations at variable temperature and impact velocity, but where the total kinetic energy of the nanoparticle was equivalent to that of a 300 K, 1000 impact. In cold impact, the temperature increases in YSZ nanoparticles at a rate of 10 13 – 10 14 K s –1 , and large temperature gradients result. Conversely, in thermal impact, nanoparticle temperatures remain uniform. Furthermore, the temperature gradients during cold impact coincide with plastic deformation in nanoparticles, while with larger thermal energies, plastic deformation is reduced.

42 ENGINEERING↗

Analytical study of a free-wing/free-trimmer concept

The free-wing/free-trimmer is a NASA-Conceived extension of the free-wing concept intended to permit the use of high-lift flaps. Wing pitching moments are balanced by a smaller, external surface attached by a boom or equivalent structure. The external trimmer is, itself, a miniature free wing, and pitch control of the wing-trimmer assembly is effected through a trailing-edge control tab on the trimmer surface. The longitudinal behavior of representative small free-wing/free-trimmer aircraft was analyzed. Aft-mounted trimmer surfaces are found to be superior to forward trimmers, although the permissible trimmer moment arm is limited, in both cases, by adverse dynamic effects. Aft-trimmer configurations provide excellent gust alleviation and meet fundamental stick-fixed stability criteria while exceeding the lift capabilities of pure free-wing configurations.

Porter, R. F.↗

Experimenting with the GMAO 4D Data Assimilation

The Global Modeling and Assimilation Office (GMAO) has been working to promote its prototype four-dimensional variational (4DVAR) system to a version that can be exercised at operationally desirable configurations. Beyond a general circulation model (GeM) and an analysis system, traditional 4DV AR requires availability of tangent linear (TL) and adjoint (AD) models of the corresponding GeM. The GMAO prototype 4DVAR uses the finite-volume-based GEOS GeM and the Grid-point Statistical Interpolation (GSI) system for the first two, and TL and AD models derived ITom an early version of the finite-volume hydrodynamics that is scientifically equivalent to the present GEOS nonlinear GeM but computationally rather outdated. Specifically, the TL and AD models hydrodynamics uses a simple (I-dimensional) latitudinal MPI domain decomposition, which has consequent low scalability and prevents the prototype 4DV AR ITom being used in realistic applications. In the near future, GMAO will be upgrading its operational GEOS GCM (and assimilation system) to use a cubed-sphere-based hydrodynamics. This versions of the dynamics scales to thousands of processes and has led to a decision to re-derive the TL and AD models for this more modern dynamics, thus taking advantage of a two-dimensional MPI decomposition and improved scalability properties. With the aid of the Transformation of Algorithms in FORTRAN (l'AF) automatic adjoint generation tool and some hand-coding, a version of the cubed-sphere-based TL and AD models, with a simplified vertical diffusion scheme, is now available, enabling multiple configurations of standard implementations of 4DV AR in GEOS. Concurrent to this development, collaboration with the National Centers for Environmental Prediction (NCEP) and the Earth System Research Laboratory (ESRL) has allowed GMAO to implement a hybrid-ensemble capability within the GEOS data assimilation system. Both 3Dand 4D-ensemble capabilities are presently available thus allowing GMAO to now evaluate the performance and benefit of various ensemble and variational assimilation strategies. This presentation will cover the most recent developments taking place at GMAO and show results from various comparisons from traditional techniques to more recent ensemble-based ones.

Todling, R.↗

Dynamic Mode Decomposition of Unsteady Pressure-Sensitive Paint Measurements for the NASA Unitary Plan Wind Tunnel Tests

This paper discusses the Dynamic Mode Decomposition (DMD) of the Unsteady Pressure-Sensitive Paint (uPSP) measurements, which were collected with four Phantom high-speed cameras at a constant sample frequency in the Ascent Transient Aerodynamics Test (ATAT) of the Space Launch System (SLS) Block 1 cargo vehicle with the Unitary Plan Wind Tunnel (UPWT) 11-by-11-foot Transonic Wind Tunnel in September 2019 at NASA Ames Research Center. The conventional DMD algorithm is based on the Singular Value Decomposition (SVD). For the data with zero mean, the DMD is equivalent to the Discrete Fourier Transform (DFT). Since the uPSP is mainly used to determine the unsteady property of the aerodynamic flow, the DMD of the uPSP measurements is implemented in two steps: (1) subtract the mean value from the uPSP measurement; (2) apply the Fast Fourier Transform (FFT) on the resulting data with zero mean. The DMD of the uPSP measurements with FFT has two advantages: (1) the FFT algorithm is well known for its computational efficiency, therefore, compared to the SVD-based DMD algorithm, the DMD with FFT reduces the computation time; (2) the DMD with FFT can be easily implemented in parallel processing. The DMD outputs were generated with the execution in parallel of a code in C, with libraries of FFTW for FFT and MPI/OpenMP for parallel processing, on the NASA Pleiades supercomputer. In this paper, the results of DMD of the uPSP measurements in the tests of Mach sweep runs of the SLS ATAT are presented, and the effectiveness of the DMD of the uPSP measurements in the diagnosis of the unsteady, aerodynamic phenomena is demonstrated. The work described in this paper is a part of NASA’s development of a new state-of-the-art uPSP capability in production wind tunnels. Funding for this research was provided by the NASA Aeroscience Evaluation and Test Capabilities Project.

Pressure-Sensitive Paint↗

Integrating data types to estimate spatial patterns of avian migration across the Western Hemisphere

For many avian species, spatial migration patterns remain largely undescribed, especially across hemispheric extents. Recent advancements in tracking technologies and high-resolution species distribution models (i.e., eBird Status and Trends products) provide new insights into migratory bird movements and offer a promising opportunity for integrating independent data sources to describe avian migration. Here, we present a three-stage modeling framework for estimating spatial patterns of avian migration. First, we integrate tracking and band re-encounter data to quantify migratory connectivity, defined as the relative proportions of individuals migrating between breeding and nonbreeding regions. Next, we use estimated connectivity proportions along with eBird occurrence probabilities to produce probabilistic least-cost path (LCP) indices. In a final step, we use generalized additive mixed models (GAMMs) both to evaluate the ability of LCP indices to accurately predict (i.e., as a covariate) observed locations derived from tracking and band re-encounter data sets versus pseudo-absence locations during migratory periods and to create a fully integrated (i.e., eBird occurrence, LCP, and tracking/band re-encounter data) spatial prediction index for mapping species-specific seasonal migrations. To illustrate this approach, we apply this framework to describe seasonal migrations of 12 bird species across the Western Hemisphere during pre- and postbreeding migratory periods (i.e., spring and fall, respectively). We found that including LCP indices with eBird occurrence in GAMMs generally improved the ability to accurately predict observed migratory locations compared to models with eBird occurrence alone. Using three performance metrics, the eBird + LCP model demonstrated equivalent or superior fit relative to the eBird-only model for 22 of 24 species–season GAMMs. In particular, the integrated index filled in spatial gaps for species with over-water movements and those that migrated over land where there were few eBird sightings and, thus, low predictive ability of eBird occurrence probabilities (e.g., Amazonian rainforest in South America). This methodology of combining individual-based seasonal movement data with temporally dynamic species distribution models provides a comprehensive approach to integrating multiple data types to describe broad-scale spatial patterns of animal movement. Further development and customization of this approach will continue to advance knowledge about the full annual cycle and conservation of migratory birds.

59 BASIC BIOLOGICAL SCIENCES↗

Gaussian integral method for void fraction

Here, a novel method, the Gaussian Integral Method (GIM), is presented for calculating void fractions in Computational Fluid Dynamics–Discrete Element Method (CFD-DEM) simulations. GIM is versatile and applicable to various grid types, including structured and unstructured polyhedral meshes, without requiring special boundary treatments. An optimization technique is introduced to make GIM independent of grid resolution and type. The method is validated against experimental data from a fluidized bed, demonstrating that GIM produces realistic simulations closely resembling experimental observations. Additionally, unstructured polyhedral grids using GIM outperform structured grids of equivalent resolution, yielding results more aligned with experimental data. The gradient of the void fraction is computed in the CFD solver and utilized in the DEM solver for precise estimation at particle locations. Overall, GIM provides an effective solution for void fraction calculations in particulate media simulations with complex geometries, enhancing the accuracy and applicability of CFD-DEM simulations for industrial processes.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

One Earth Energy Static and Dynamic Reservoir Modeling

This report presents the static and dynamic reservoir modeling conducted for the CarbonSAFE Phase III Illinois Storage Corridor project to assess the feasibility of commercial-scale CO 2 storage in the Mt. Simon Sandstone at the One Earth Energy (OEE) site in McLean County, Illinois. Three-dimensional geocellular models of the Mt. Simon storage complex were developed in Petrel ® by integrating petrophysical log data, core analyses, and seismic surveys from the OEE #1 stratigraphic test well and two nearby wells, with multiple model versions created as new data became available. Dynamic reservoir simulations, performed using Landmark's Nexus software, progressed through three phases (preliminary, sensitivity, and UIC Class VI permit studies) evaluating injection scenarios across varying rates, well orientations, permeability models, and multi-well configurations. Results demonstrate that commercial-scale storage is feasible: three injection wells spaced approximately one mile apart can store a total of 90 million tonnes of CO 2 over 20 years, producing a combined plume with an equivalent radius of 3.2 miles and a maximum pressure-front-defined Area of Review of 178 mi 2 at the end of injection that diminishes to 34 mi 2 after 50 years of post-injection monitoring. Sensitivity analyses indicate that a 20% change in porosity or permeability yields approximately a 7% change in AoR radius, and that perforating the high-permeability arkosic zone minimizes the pressure front compared to injection in the upper Mt. Simon Sandstone.

09 BIOMASS FUELS↗