Search NASA⌕ Search

SEARCH · Search NASA

Results for “transport limitations”

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 307 records · Page 17

Cycling phosphorus on the Archean Earth: Part I. Continental weathering and riverine transport of phosphorus

Phosphorus (P) is the key nutrient thought to limit primary productivity on geological timescales. Phosphate levels in Archean marine sediments are low, but quantification of the P cycle and how it changed through a billion years of recorded Archean history remain a challenge, hindering our understanding of the role played by P in biosphere/geosphere co-evolution on the early Earth. Here, we design kinetic and thermodynamic models to quantitatively assess one key component of the early P cycle – continental weathering – by considering the emergence and elevation of continents, as well as the evolution of climate, the atmosphere, and the absence of macroscopic vegetation during the Archean Eon. Our results suggest that the weathering rate of apatite, the major P-hosting mineral in the rocks, was at least five times higher in the Archean Eon than today, attributable to high levels of pCO2,g. Despite this, the weathering flux of P to the oceans was negligible in the early Archean Eon, increasing to a level comparable to or greater than the modern by the end of eon, a consequence of accelerating continental emergence. Furthermore, our thermodynamic calculations indicate high solubilities of primary and secondary P-hosting minerals in the acidic weathering fluids on land, linking to high Archean pCO2,g. Thus, weathering of P was both kinetically and thermodynamically favorable on the Archean Earth, and river water could transport high levels of dissolved P to the oceans, as also supported by the observed P-depletion in our new compilation of Archean paleosols. Lastly, we evaluated the relative rates of physical erosion and chemical weathering of silicates during the Archean Eon. The results suggest that continental weathering on the early and middle Archean Earth might have been transport-limited due to low erosion rates associated with limited subaerial emergence and low plateau elevations; by the late Archean, however, continental weathering would have transited to kinetically-limited state because of continental emergence, increased plateau elevation, and weakening weathering rates. Overall, our weathering calculations together with paleosol evidence indicate an increasing flux of bioavailable P to the oceans through time, associated with late Archean continental emergence, reaching levels comparable to or higher than modern values by the end of the eon. Increased P fluxes could have fueled increasing rates of primary production, including oxygenic photosynthesis, through time, facilitating the irreversible oxidation of the Earth’s atmosphere early in the Proterozoic Eon.

Phosphorus↗

Linear ion-scale microstability analysis of high and low-collisionality NSTX discharges and NSTX-U projections

Linear gyrokinetic simulations were conducted to investigate ion-gyroradius-scale micro-instability predictions for high-beta NSTX discharges and NSTX-U projections that span over an order of magnitude variation in collisionality. A complex mix of microtearing modes and hybrid trapped electron modes/kinetic ballooning modes (TEM/KBM) is predicted for all experimental or projected conditions. Ion temperature gradient (ITG) instabilities are typically stable in the NSTX discharges investigated, consistent with the observed neoclassical ion thermal transport. ITG thresholds inferred from the simulations are typically much higher than the experimental NSTX gradients, as well as the projected gradients in the NSTX-U scenario, which assumed ion temperatures limited by neoclassical transport only. The analysis suggests ITG instabilities are unlikely to contribute significant anomalous thermal losses in high-beta, lower collisionality NSTX-U scenarios. On the other hand, the NSTX experimental profiles and NSTX-U projections are predicted to be very close to the predicted onset of unstable KBM at most radii investigated. The proximity of the various discharges to the KBM instability threshold implies it may play an important role in setting profile shapes and limiting global energy confinement. Here, it remains to be understood and predicted how KBM contributes to multi-channel transport (thermal and particle transport, for both ions and electrons) in a way that is consistent with experimental inferences.

36 MATERIALS SCIENCE↗

DG-IMEX method for a two-moment model for radiation transport in the $\mathscr{O}$($v$/$c$) limit

Here, we consider neutral particle systems described by moments of a phase-space density and propose a realizability-preserving numerical method to evolve a spectral two-moment model for particles interacting with a background fluid moving with nonrelativistic velocities. The system of nonlinear moment equations, with special relativistic corrections to $\mathscr{O}$($v$/$c$), expresses a balance between phase-space advection and collisions and includes velocity-dependent terms that account for spatial advection, Doppler shift, and angular aberration. The model is conservative for the correct $\mathscr{O}$($v$/$c$) Eulerian-frame number density and is consistent, to $\mathscr{O}$($v$/$c$), with Eulerian-frame energy and momentum conservation. This model is closely related to the one promoted by Lowrie et al. and similar to models currently used to study transport phenomena in large-scale simulations of astrophysical environments. The proposed numerical method is designed to preserve moment realizability, which guarantees that the moments correspond to a nonnegative phase-space density. The realizability-preserving scheme consists of the following key components: (i) a strong stability-preserving implicit-explicit (IMEX) time-integration method; (ii) a discontinuous Galerkin (DG) phase-space discretization with carefully constructed numerical uxes; (iii) a realizability-preserving implicit collision update; and(iv) a realizability-enforcing limiter. In time integration, nonlinearity of the moment model necessitates solution of nonlinear equations, which we formulate as fixed-point problems and solve with tailored iterative solvers that preserve moment realizability with guaranteed global convergence. We also analyze the simultaneous Eulerian-frame number and energy conservation properties of the semi-discrete DG scheme and propose a "spectral redistribution" scheme that promotes Eulerian-frame energy conservation. Through numerical experiments, we demonstrate the accuracy and robustness of this DG-IMEX method and investigate its Eulerian-frame energy conservation properties.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

The Performance Constraints of Gain-Limited Propulsion Systems

Rapid transportation of human crews to destinations throughout the solar system will require propulsion systems having not only very high exhaust velocities (i.e., I(sub sp) greater or equal to 10(exp 4) to 10(exp 5) sec) but also extremely low mass-power ratios (i.e., alpha less than or equal to 10(exp -2) kg/kW). These criteria are difficult to meet with electric propulsion and other power-limited systems, but may be achievable with propulsion concepts that use onboard power to produce a net gain in energy via fusion or some other nuclear process. This paper compares the fundamental performance of these gain-limited systems with that of power-limited systems, and determines from a generic power balance the gains required for ambitious planetary missions ranging up to 100 AU. Results show that energy gain reduces the required effective mass-power ratio of the system, thus enabling shorter trip times than those of power-limited concepts.

Schmidt, G. R.↗

Predictive Contaminant Transport Simulation with P2R Model for the Composite Analysis Limited Source Sensitivity Case.

The Plateau to River (P2R) is a groundwater flow and contaminant fate and transport (F&T) simulation model used to support remedial activities conducted by the Central Plateau Cleanup Company at the Hanford Site in Washington State. Figure 1 illustrates the P2R Model extents, discretization, and boundary conditions. The P2R Model is utilized in the Composite Analysis (CA) for the Hanford Site as the computational engine for computing F&T predictions as described in CP-60406, Hanford Site Composite Analysis Technical Approach Description: Groundwater. The model simulates contaminants of concern within the saturated zone of the uppermost aquifer beneath the Central Plateau and downgradient to the Columbia River. CP-57037, Model Package Report for the Plateau to River Model Version 8.3 documents the current version of the P2R Model including a description of the conceptual site model, model development and calibration, and limitations to the model application. Simulations conducted to support the dose calculations required by the CA are documented in ECF-HANFORD-19-0119, Predictive Flow Simulation with the P2R Model for the Composite Analysis Base Case and ECF-HANFORD-19-0120, Contaminant Transport Simulation with the P2R Model for the Composite Analysis Base Case.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Role of poloidal E × B drift in divertor heat transport in DIII‐D

Abstract Simulations for DIII‐D high confinement mode plasmas with the multifluid code UEDGE show a strong role of poloidal E × B drifts on divertor heat transport, challenging the paradigm of conduction‐limited scrape‐off layer (SOL) transport. While simulations with reduced drift magnitude are well aligned with the assumption that electron heat conduction dominates the SOL heat transport, simulations with drifts predict that the poloidal convective E × B heat transport dominates over electron heat conduction in both attached and detached conditions. As poloidal E × B flow propagates across magnetic field lines, poloidal transport with shallow magnetic pitch angles can reach values that are of the same order as would be provided by sonic flows parallel to the field lines. These flows can lead to strong convection‐dominated divertor heat transport, increasing the poloidal volume of radiative power front, consistent with previous measurements at DIII‐D. Due to these convective flows, the Lengyel integral approach, assuming zero convective fraction, is expected to provide a pessimistic estimate for the radiative capability of impurities in the divertor. For the DIII‐D simulations shown here, the Lengyel integral approach underestimates the radiated power by a factor of 6, indicating that, for reliable DIII‐D divertor power exhaust predictions, full two‐dimensional (2D) calculations, including drifts, would be necessary.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Is It Worth It? - the Economics of Reusable Space Transportation

Over the past several decades billions of dollars have been invested by governments and private companies in the pursuit of lower cost access to space through earth-to-orbit (ETO) space transportation systems. Much of that investment has been focused on the development and operation of various forms of reusable transportation systems. From the Space Shuttle to current efforts by private commercial companies, the overarching belief of those making such investments has been that reusing system elements will be cheaper than utilizing expendable systems that involve throwing away costly engines, avionics, and other hardware with each flight. However, the view that reusable systems are ultimately a "better" approach to providing ETO transportation is not held universally by major stakeholders within the space transportation industry. While the technical feasibility of at least some degree of reusability has been demonstrated, there continues to be a sometimes lively debate over the merits and drawbacks of reusable versus expendable systems from an economic perspective. In summary, is it worth it? Based on our many years of direct involvement with the business aspects of several expendable and reusable transportation systems, it appears to us that much of the discussion surrounding reusability is hindered by a failure to clearly define and understand the financial and other metrics by which the financial "goodness" of a reusable or expandable approach is measured. As stakeholders, the different users and suppliers of space transportation have a varied set of criteria for determining the relative economic viability of alternative strategies, including reusability. Many different metrics have been used to measure the affordability of space transportation, such as dollars per payload pound (kilogram) to orbit, cost per flight, life cycle cost, net present value/internal rate of return, and many others. This paper will examine the key considerations that influence stakeholders as they make space transportation investment decisions, including primary metrics by which various stakeholders measure financial goodness and other factors that significantly shape decisions to invest in reusable or expendable systems. It must be noted at the outset that reusable systems take many forms and perform different transportation functions including, but not limited to, ETO payload delivery. The discussion in this paper is limited to the economics of ETO transportation systems.

Webb, Richard↗

Self-Sacrificial Template Synthesis of Fe-N-C Catalysts with Dense Active Sites Deposited on A Porous Carbon Network for High Performance in PEMFC

In this study, iron-nitrogen-carbon (Fe-N-C) single-atom catalysts are promising sustainable alternatives to the costly and scarce platinum (Pt) to catalyze the oxygen reduction reactions (ORR) at the cathode of proton exchange membrane fuel cells (PEMFCs). However, Fe-N-C cathodes for PEMFC are made thicker than Pt/C ones, in order to compensate for the lower intrinsic ORR activity and site density of Fe-N-C materials. The thick electrodes are bound with mass transport issues that limit their performance at high current densities, especially in H 2 /air PEMFCs. Practical Fe-N-C electrodes must combine high intrinsic ORR activity, high site density, and fast mass transport. Herein, it has achieved an improved combination of these properties with a Fe-N-C catalyst prepared via a two-step synthesis approach, constructing first a porous zinc-nitrogen-carbon (Zn-N-C) substrate, followed by transmetallating Zn by Fe via chemical vapor deposition. A cathode comprising this Fe-N-C catalyst has exhibited a maximum power density of 0.53 W cm -2 in H 2 /air PEMFC at 80 °C. The improved power density is associated with the hierarchical porosity of the Zn-N-C substrate of this work, which is achieved by epitaxial growth of ZIF-8 onto g-C 3 N 4 , leading to a micro-mesoporous substrate.

25 ENERGY STORAGE↗

Electrically Insulated Catalyst–Ionomer Anode Interfaces toward Durable Alkaline Membrane Electrolyzers

Anion-exchange-membrane water electrolysis (AEMWE) is an emerging technology for scalable hydrogen production. AEMWE has poor durability when operating without supporting electrolyte due to the oxidation of ionomers and membranes in contact with the anode oxygen evolution reaction (OER) catalyst. We report a new “passivated” anode architecture for AEMWE where the OER catalysts and ionomers are physically separated with a thin film amorphous oxide coating that is electrically insulating but conductive to hydroxide ions. We find that 2–3 nm of HfO x passivation layers show sufficient hydroxide ion transport to minimally limit the cell performance while suppressing ionomer degradation with both Ir (500 mA·cm –2 for 40 h) and CoO x (1.0 A·cm –2 for 100 h) model porous-transport-layer-supported catalysts in AEMWE. As a result, this interfacial engineering approach guides electrode design to improve the durability of AEMWE, particularly for systems operating with pure-water feed.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Tentative civil airworthiness flight criteria for powered-lift transports

Representatives of the U.S., British, French, and Canadian airworthiness authorities participated in a NASA/FAA program to formulate tentative civil airworthiness flight criteria for powered-lift transports. The ultimate limits of the flight envelope are defined by boundaries in the airspeed/path-angle plane. Angle of attack and airspeed margins applied to these ultimate limits provide protection against both atmospheric disturbances and disturbances resulting from pilot actions or system variability, but do not ensure maneuvering capability directly, as the 30% speed margin does for conventional transports. Separate criteria provide for direct demonstration of adequate capability for approach path control, flare and landing, and for go-around. Demonstration maneuvers are proposed, and appropriate abuses and failures are suggested. Taken together, these criteria should permit selection of appropriate operating points within the flight envelopes for the approach, landing, and go-around flight phases which are likely to be most critical for powered-lift aircraft.

Hynes, C. S.↗

OPTIMUS{sup TM} Packaging for Intermediate Level Wastes and Spent Fuel - 20129

NAC's OPTIMUS{sup TM} packaging is a family of Type B(U)F transportation packages designed for maximum flexibility and cost-efficiency to support shipment of a wide range of challenging wastes and materials. They are small, modular packaging designs meeting DOT weight limits for road transport. OPTIMUS{sup TM} contents include low level waste, mixed low-level waste, irradiated fuel waste and intermediate level waste. Contents can also include aerosol cans with compressed gas and liquefied gas propellants and standard DOT 3E lecture bottles for shipment to treatment facilities. OPTIMUS{sup TM} packaging provides a cost-effective solution for challenging and unique waste contents combined with design simplicity and operational flexibility meeting the requirements of both SSR-6 and 10 CFR 71 for Type B(U)F packaging. The modular design of the packaging components combined with the flexibility to add modified internals for enhanced packaging performance makes OPTIMUS{sup TM} the preferred solution for many radioactive material transportation needs. The OPTIMUS{sup TM} product line currently includes two packaging designs; OPTIMUS{sup TM}-L and OPTIMUS{sup TM}-H. OPTIMUS{sup TM}-L is a lightweight transportation packaging intended for low-activity wastes, whereas OPTIMUS{sup TM}-H is a larger, heavier package intended for high-activities wastes. Due to its small size and weight, six (6) to ten (10) OPTIMUS{sup TM}-L packages can be accommodated on a single legal-weight truck shipment, depending on the weight of the contents. For typical payloads, two (2) OPTIMUS{sup TM}-H packages can be accommodated on a single legal-weight truck shipment. The OPTIMUS{sup TM} packaging share the same cask containment vessel (CCV) in the OPTIMUS{sup TM}-H and OPTIMUS{sup TM}-L. The large cavity size of the CCV, which can accommodate a 416-liter (110-gallon) drum, combined with the small size, low weight and modularity of the OPTIMUS{sup TM} packaging provides unmatched flexibility. In this paper, NAC provides a technical overview of the OPTIMUS{sup TM} packaging and identifies the design features and technology advancements making the OPTIMUS{sup TM} a readily adaptable and flexible solution for packaging processing facilities, reactor and decommissioning wastes. Optimized shielding performance through the addition of internal shielding components and alternative shipping configurations are discussed. In addition, an overview of the first major deployment of the OPTIMUS{sup TM} packaging for the Whiteshell Laboratories Closure Project (WLCP) is discussed. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Radial–axial transport coordination enhances sugar translocation in the phloem vasculature of plants

Abstract Understanding mass transport of photosynthates in the phloem of plants is necessary for predicting plant carbon allocation, productivity, and responses to water and thermal stress. Several hypotheses about optimization of phloem structure and function and limitations of phloem transport under drought have been proposed and tested with models and anatomical data. However, the true impact of radial water exchange of phloem conduits with their surroundings on mass transport of photosynthates has not been addressed. Here, the physics of the Munch mechanism of sugar transport is re-evaluated to include local variations in viscosity resulting from the radial water exchange in two dimensions (axial and radial) using transient flow simulations. Model results show an increase in radial water exchange due to a decrease in sap viscosity leading to increased sugar front speed and axial mass transport across a wide range of phloem conduit lengths. This increase is around 40% for active loaders (e.g. crops) and around 20% for passive loaders (e.g. trees). Thus, sugar transport operates more efficiently than predicted by previous models that ignore these two effects. A faster front speed leads to higher phloem resiliency under drought because more sugar can be transported with a smaller pressure gradient.

58 GEOSCIENCES↗

Integration of the Back End of the Nuclear Fuel Cycle

Management of spent nuclear fuel and high-level radioactive waste consists of three main phases – storage, transportation, and disposal – commonly referred to as the back end of the nuclear fuel cycle. Current practice for commercial spent nuclear fuel management in the United States (US) includes temporary storage of spent fuel in both pools and dry storage systems at operating or shutdown nuclear power plants. Storage pools are filling to their operational capacity, and management of the approximately 2,200 metric tons of spent fuel newly discharged each year requires transferring older and cooler spent fuel from pools into dry storage. Unless a repository becomes available that can accept spent fuel for permanent disposal, projections indicate that the US will have approximately 136,000 metric tons of spent fuel in dry storage systems by mid-century, when the last plants in the current reactor fleet are decommissioned. Current designs for dry storage systems rely on large multi-assembly canisters, the most common of which are so-called “dual-purpose canisters” (DPCs). DPCs are certified for both storage and transportation, but are not designed or licensed for permanent disposal. The large capacity (greater number of spent fuel assemblies) of these canisters can lead to higher canister temperatures, which can delay transportation and/or complicate disposal. This current management practice, in which the utilities continue loading an ever-increasing inventory of larger DPCs, does not emphasize integration among storage, transportation, and disposal. This lack of integration does not cause safety issues, but it does lead to a suboptimal system that increases costs, complicates storage and transportation operations, and limits options for permanent disposal. This paper describes strategies for improving integration of management practices in the US across the entire back end of the nuclear fuel cycle. The complex interactions between storage, transportation, and disposal make a single optimal solution unlikely. However, efforts to integrate various phases of nuclear waste management can have the greatest impact if they begin promptly and continue to evolve throughout the remaining life of the current fuel cycle. A key factor that influences the path forward for integration of nuclear waste management practices is the identification of the timing and location for a repository. The most cost-effective path forward would be to open a repository by mid-century with the capability to directly dispose of DPCs without repackaging the spent fuel into disposalready canisters. Options that involve repackaging of spent fuel from DPCs into disposalready canisters or that delay the repository opening significantly beyond mid-century could add 10s of billions of dollars to the total system life cycle cost.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

The systems improved numerical fluids analysis code (SINFAC)

A computer program has been developed to model the steady-state performance of single and two-phase heat transport systems. The development effort had the intent of making the program general, relatively easy to implement and easy to integrate into currently available software. The primary use of the program is for, but not limited to, heat transport system for spacecraft. In order to show the program has achieved its intended goals, the program was first applied to a small loop and a quick parametric study done. Next a verification study was done in which the temperature data from a thermal loop flown on the space shuttle was compared to modeler predictions.

Schweickart, Russell B.↗

All-Inorganic Halide Perovskites as Potential Thermoelectric Materials: Dynamic Cation off-Centering Induces Ultralow Thermal Conductivity

Halide perovskites are anticipated to impact next generation high performance solar cells because of their extraordinary charge transport and optoelectronic properties. However, their thermal transport behavior has received limited attention. In this work, we studied the thermal transport and thermoelectric properties of the CsSnBr 3–x I x perovskites. We find a strong correlation between lattice dynamics and an ultralow thermal conductivity for series CsSnBr 3–x I x reaching 0.32 Wm –1 K –1 at 550 K. The CsSnBr 3–x I x also possess a decent Seebeck coefficient and controllable electrical transport properties. The crystallography data and theoretical calculations suggest the Cs atom deviates from its ideal cuboctahedral geometry imposed by the perovskite cage and behaves as a heavy atom rattling oscillator. This off-center tendency of Cs, together with the distortion of SnX 6 (X = Br or I) octahedra, produces a highly dynamic and disordered structure in CsSnBr 3–x I x , which gives rise to a very low Debye temperature and phonon velocity. Moreover, the low temperature heat capacity data suggests strong coupling between the low frequency optical phonons and heat carrying acoustical phonons. In conclusion, this induces strong phonon resonance scattering that induces the ultralow lattice thermal conductivity of CsSnBr 3–x I x .

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Study of the oxygen transport through Ag (110), Ag (poly), and Ag 2.0 Zr

The transport of oxygen through high-purity membranes of Ag (110), Ag (poly), Ag (nano), and Ag 2.0 Zr has been studied by an ultrahigh vacuum permeation method over the temperature range of 400-800 C. The data show that there are substantial deviations from ordinary diffusion-controlled transport. A surface limitation has been confirmed by glow-discharge studies where the upstream O2 supply has been partially converted to atoms, which, for the same temperature and pressure, gave rise to over an order of magnitude increase in transport flux. Further, the addition of 2.0 wt percent Zr to the Ag has provided increased dissociative adsorption rates, which, in turn, increased the transport flux by a factor of 2. It was also observed that below a temperature of 630 C, the diffusivity exhibits an increase in activation energy of over 4 kcal/mol, which has been attributed to trapping of the atomic oxygen and/or kinetic barriers at the surface and subsurface of the vacuum interface. Above 630 C, the activation barrier decreases to the accepted value of about 11 kcal/mol for Ag (poly), consistent with zero concentration at the vacuum interface.

Outlaw, R. A.↗

Integrated Multiscale Modeling for Design of Robust 3D Solid-State Lithium Batteries - FY21 Annual Report

This project is developing a multiscale, multi-physics modeling framework for probing the effects of materials microstructure and device architecture on ion transport within 3D ceramic solid-state battery materials, with the goal of enhancing performance and reliability. The project has three primary objectives: (1) integrate multi physics and multiscale model components; (2) understand interface- and microstructure-derived limitations on ion transport; and (3) derive key structure-performance relations for enabling future optimization.

25 ENERGY STORAGE↗

Direct Observations of Solute Dispersion in Rocks With Distinct Degree of Sub‐Micron Porosity

Abstract The transport of chemical species in rocks is affected by their structural heterogeneity to yield a wide spectrum of local solute concentrations. To quantify such imperfect mixing, advanced methodologies are needed that augment the traditional breakthrough curve analysis by probing solute concentration within the fluids locally. Here, we demonstrate the application of asynchronous, multimodality imaging by X‐ray computed tomography (XCT) and positron emission tomography (PET) to the study of passive tracer experiments in laboratory rock cores. The four‐dimensional concentration maps measured by PET reveal specific signatures of the transport process, which we have quantified using fundamental measures of mixing and spreading. We observe that the extent of solute spreading correlate strongly with the strength of subcore‐scale porosity heterogeneity measured by XCT, while dilution is enhanced in rocks containing substantial sub‐micron porosity. We observe that the analysis of different metrics is necessary, as they can differ in their sensitivity to the strength and forms of heterogeneity. The multimodality imaging approach is uniquely suited to probe the fundamental difference between spreading and mixing in heterogeneous media. We propose that when multi‐dimensional data is available, mixing and spreading can be independently quantified using the same metric. We also demonstrate that one‐dimensional transport models have limited predictive ability toward the internal evolution of the solute concentration, when the model is solely calibrated against the effluent breakthrough curves. The data set generated in this study can be used to build realistic digital rock models and to benchmark transport simulations that account deterministically for rock property heterogeneity.

Kurotori, Takeshi [Department of Chemical Engineer↗