Search NASA⌕ Search

SEARCH · Search NASA

Results for “FLUID MECHANICS”

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 919 records · Page 51
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Time reversal symmetry breaking and odd viscosity in active fluids: Green–Kubo and NEMD results

Active fluids, which are driven at the microscale by non-conservative forces, are known to exhibit novel transport phenomena due to the breaking of time reversal symmetry. Recently, Epstein and Mandadapu [arXiv:1907.10041 (2019)] obtained Green-Kubo relations for the full set of viscous coefficients governing isotropic chiral active fluids, including the so-called odd viscosity, invoking Onsager's regression hypothesis for the decay of fluctuations in active non-equilibrium steady states. In this Communication, we test these Green-Kubo relations using molecular dynamics simulations of a canonical model system consisting of actively torqued dumbbells. We find the resulting odd and shear viscosity values from the Green-Kubo relations to be in good agreement with values measured independently through non-equilibrium molecular dynamics flow simulations. This provides a test of the Green-Kubo relations and lends support to the application of the Onsager regression hypothesis in relation to viscous behaviors of active matter systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS2023

Fitting Leak Test Report: Ground-Based Cryogenic Leak Test of Fittings for Cryogenic Fluid Management

EXECUTIVE SUMMARY Mechanically connected joints used in cryogenic fluid lines as part of space flight elements need to survive launch vibrations and remain leak-free to minimize the loss of on-board commodity and hazardous gas accumulation. In 2020, a cryogenic test apparatus was developed which can evaluate the leak performance of pressurized threaded fluid fittings. The fittings were mounted in the TVAC and cooled to cryogenic test temperature and pressurized with helium while the leak rate was measured using a calibrated GHe leak detector. The test articles for the initial proof of concept testing were ¼ and 1 inch Swagelok VCR fittings with three different types of seal rings copper, nickel, and Ni. Each fitting configuration (size/seal material) was subjected to two consecutive cryogenic thermal cycles, followed by exposure to a launch vibration profile at ambient temperature, after which two additional TVAC cycle tests were performed. The testing reported here is a continuation of the 2020 tests with a statistically significant large number of samples and test runs. Three Swagelok VCR fitting sizes were tested ¼, ½ and 1 inch, and five (5) samples of each fitting size, each sample was tested with SST and Ni seal rings (Total of 30 unique test articles). Each test article was subjected to four (4) thermal cycles. Half of these cycles were performed before vibration testing and half were performed after vibration testing. The vibration testing was performed to evaluate the ability of the fittings to survive launch-type vibration profiles and remain leak-free. Leak checking of each fitting was completed at temperatures between 20K – 30K. The test procedure in Section 8.0 was designed to facilitate a qualification test program by allowing a higher test throughput rate coupled with repeatable test profiles. Results were very positive and show that out of the 30 samples they all passed with leak rates a factor of 2-3 lower than the established 10-6sccs GHe leak threshold. The result showed the Ni seals had lower leak rate, but the SST was more rugged. There were two deviations where damage to the Ni seal ring during assembly resulted in a leaky fitting, this is discussed in Section 10.7 Test Deviations. These fittings show great promise for space flight use and further testing is recommended to fully qualify the fittings per the ASTM F1387-19 and/or other relevant NASA specifications. The test equipment hardware and software capability developed for this testing is generic and not restricted to VCR fittings. It can be employed to evaluate/qualify the leak performance of other types of fittings and a wide range of other cryogenic fluid components such as valves, gages, connectors, etc.

Cryogenic↗
DOE OSTI2023

Understanding Subsurface Fracture Evolution Dynamics Using Time‐Lapse Full Waveform Inversion of Continuous Active‐Source Seismic Monitoring Data

Abstract Predicting the behavior, geometry, and flow properties of subsurface fractures remains a challenging problem. Seismic models that can characterize fractures usually suffer from low spatiotemporal resolution. Here, we develop a correlative double‐difference time‐lapse full waveform inversion of continuous active source seismic monitoring data for determining high‐spatiotemporal‐resolution time‐lapse V p models of in‐situ fracture evolution at a shallow contamination site in Wyoming, USA. Assisted by rock physics modeling, we find that (a) rapidly increasing pore pressure initializes and grows the fracture, increasing the porosity slightly (from ∼13.7% to ∼14.6%) in the tight clay formation, thus decreasing V p (∼50 m/s); (b) the fluid injection continues decreasing V p , likely through the introduction of gas bubbles in the injectate; and (c) final V p reductions reach over ∼150 m/s due to a posited ∼4.5% gas saturation. Our results demonstrate that high‐resolution V p changes are indicative of mechanical and fluid changes within the fracture zone during hydrofracturing.

58 GEOSCIENCES↗
NASA NTRS1983

Drinking and water balance during exercise and heat acclimation

The interactions between fluid intake and balance, and plasma ion, osmotic, and endocrine responses during dehydration produced by exercise in cool and warm environments during acclimation are explored. Two groups of five male subjects performed 8 days of ergometer exercise in hot and thermoneutral conditions, respectively. The exercise trials lasted 2 hr each. Monitoring was carried out on the PV, osmotic, sodium, and endocrine concentrations, voluntary fluid intake, fluid balances, and fluid deficits. A negative correlation was observed between the plasma sodium and osmolality during acclimation. The presence of hypervolemia during acclimation is suggested as a cause of drinking, while the vasopressin concentration was not found to be a significant factor stimulating drinking. Finally, the predominant mechanism in fluid intake during exercise and heat exposure is concluded to be the renin-angiotensin II system in the presence of reductions in total body water and extracellular plasma volumes.

Greenleaf, J. E.↗
NASA NTRS2013

Atomic Force Microscope Mediated Chromatography

The atomic force microscope (AFM) is used to inject a sample, provide shear-driven liquid flow over a functionalized substrate, and detect separated components. This is demonstrated using lipophilic dyes and normal phase chromatography. A significant reduction in both size and separation time scales is achieved with a 25-micron-length column scale, and one-second separation times. The approach has general applications to trace chemical and microfluidic analysis. The AFM is now a common tool for ultra-microscopy and nanotechnology. It has also been demonstrated to provide a number of microfluidic functions necessary for miniaturized chromatography. These include injection of sub-femtoliter samples, fluidic switching, and sheardriven pumping. The AFM probe tip can be used to selectively remove surface layers for subsequent microchemical analysis using infrared and tip-enhanced Raman spectroscopy. With its ability to image individual atoms, the AFM is a remarkably sensitive detector that can be used to detect separated components. These diverse functional components of microfluidic manipulation have been combined in this work to demonstrate AFM mediated chromatography. AFM mediated chromatography uses channel-less, shear-driven pumping. This is demonstrated with a thin, aluminum oxide substrate and a non-polar solvent system to separate a mixture of lipophilic dyes. In conventional chromatographic terms, this is analogous to thin-layer chromatography using normal phase alumina substrate with sheardriven pumping provided by the AFM tip-cantilever mechanism. The AFM detection of separated components is accomplished by exploiting the variation in the localized friction of the separated components. The AFM tip-cantilever provides the mechanism for producing shear-induced flows and rapid pumping. Shear-driven chromatography (SDC) is a relatively new concept that overcomes the speed and miniaturization limitations of conventional liquid chromatography. SDC is based on a sliding plate system, consisting of two flat surfaces, one of which has a recessed channel. A fluid flow is produced by axially sliding one plate past another, where the fluid has mechanical shear forces imposed at each point along the channel length. The shear-induced flow rates are very reproducible, and do not have pressure or voltage gradient limitations. SDC opens up a new range of enhanced separation kinetics by permitting the sample confinement with submicron dimensions. Small, highly confined liquid is advantageous for chromatographic separation because the separation rate is known to scale according to the square of the confined sample diameter. In addition, because shear-driven flows are not limited by fluid velocity, shear-driven liquid chromatography may provide up to 100,000 plate efficiency.

Anderson, Mark S.↗
NASA NTRS2009

Simulating Bone Loss in Microgravity Using Mathematical Formulations of Bone Remodeling

Most mathematical models of bone remodeling are used to simulate a specific bone disease, by disrupting the steady state or balance in the normal remodeling process, and to simulate a therapeutic strategy. In this work, the ability of a mathematical model of bone remodeling to simulate bone loss as a function of time under the conditions of microgravity is investigated. The model is formed by combining a previously developed set of biochemical, cellular dynamics, and mechanical stimulus equations in the literature with two newly proposed equations; one governing the rate of change of the area of cortical bone tissue in a cross section of a cylindrical section of bone and one governing the rate of change of calcium in the bone fluid. The mechanical stimulus comes from a simple model of stress due to a compressive force on a cylindrical section of bone which can be reduced to zero to mimic the effects of skeletal unloading in microgravity. The complete set of equations formed is a system of first order ordinary differential equations. The results of selected simulations are displayed and discussed. Limitations and deficiencies of the model are also discussed as well as suggestions for further research.

Pennline, James A.↗
DOE OSTI2020

Review of Tools for Modeling Core Radial Expansion in Liquid Metal-Cooled Fast Reactors

Core radial expansion in liquid-metal cooled fast reactor systems is a well-known phenomenon that produces strong reactivity feedback effects. It is crucial from a safety standpoint to design a reactor with the proper core restraint system to guide expansion of the fuel assemblies into a formation that produces negative reactivity feedback in accident conditions. However, the modeling of core radial expansion and its associated reactivity feedback is extremely challenging. Liquid metal-cooled fast reactor designers (both industry and commercial) have pointed to improved modeling of core radial expansion as a key modeling challenge for these types of reactors. The Department of Energy Nuclear Energy Advanced Modeling and Simulation (DOE-NEAMS) program has an important role to play in helping to develop a robust, predictive tool for this phenomenon. A literature review has been performed to identify current simulation capabilities applicable to the prediction of core radial expansion in liquid metal-cooled fast reactors and the resultant reactivity feedback effects. An overview as well as a detailed description of the multi-physics phenomena underlying core radial expansion is given in the first two chapters. The objective of this report is to identify existing software that has been or could be applied to the core radial expansion problem, and then to identify possible paths to develop a robust, coupled Multiphysics tool for advanced analysis, in which all three physics (radiation transport, structural mechanics, thermal fluids) are taken into account considering feedback effects from each physics to the other. Several computer codes and workflows have been developed in the United States and abroad to perform individual aspects of radial expansion calculations or to perform loosely coupled analysis. No code system currently exists which tightly and robustly couples the neutronics, thermal hydraulics, and thermal mechanical physics with enough detail to fully resolve the complex core radial expansion reactivity feedback effects. The future availability of such a code system is of vital importance to fully understanding the reactivity feedback effects that occur due to radial expansion, and consequently to optimizing the design of the core restraint system. A potential code development path forward using MOOSE-based tools is suggested in order to leverage the natural tight coupling and robustness of MOOSE-based applications. Additionally, potential gaps for modeling and common assumptions are listed which require further investigation to ensure accuracy for prediction of deformations.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2020

Temporal dynamics of the alignment of the turbulent stress and strain rate

We report that the flux of energy between scales in turbulence can be cast as the result of the interaction between a turbulent stress and a rate of strain. Efficient energy transfer requires that these two tensors be oriented properly relative to each other, but previous work has shown that the instantaneous alignment between them is poor. Here, we consider the temporal dynamics of this alignment in a direct numerical simulation of isotropic turbulence. We show that the orientation of the stress lags behind that of the strain rate, both at a single location and along trajectories. However, the timescale of the reorientation of the stress in these cases does not follow the expected dynamical scaling with length scale. To capture the proper dynamical scaling, we reformulate the energy flux between scales using the right Cauchy-Green strain tensor and the second Piola-Kirchhoff stress tensor. Our results highlight the key role played by the deformation of fluid elements in the physics of the energy cascade, and suggest that their irreversible deformation is a Lagrangian manifestation of the cascade's broken time-reversal symmetry.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Investigating Porous Media for Relief Printing Using Micro-Architected Materials

Advances in printed electronics are predicated on the integration of sophisticated printing technologies with functional materials. Although scalable manufacturing methods, such as letterpress and flexographic printing, have significant history in graphic arts printing, functional applications require sophisticated control and understanding of nanoscale transfer of fluid inks. In this paper, a versatile platform is introduced to study and engineer printing forms, exploiting a microscale additive manufacturing process to design micro-architected materials with controllable porosity and deformation. Building on this technology, controlled ink transfer for submicron functional films is demonstrated. The design freedom and high-resolution 3D control afforded by this method provide a rich framework for studying mechanics of fluid transfer for advanced manufacturing processes.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Implosion dynamics of triple-nozzle gas-puff z pinches on COBRA

Experiments on the 1-MA, 220-ns COBRA generator at Cornell University were conducted to provide detailed measurements of structured cylindrical gas-puff z pinches. In the experiments, a 7 cm diameter triple-nozzle gas valve assembly with concentric outer and inner annular nozzles and a central gas jet initialize the z-pinch load with various working gases, radial density profiles, and externally applied axial magnetic fields. Planar laser-induced fluorescence provides a measure of the initial neutral gas density of the load, while three-frame laser shearing interferometry and multi-frame extreme ultraviolet (XUV) cameras reveal the formation and propagation of a magneto-Rayleigh–Taylor (MRT) unstable shock layer. Implosion trajectories are compared to simple, experimentally informed models and found to be in good agreement. Differences in the structure of the accelerating plasma sheath and evolution of the MRT instability are observed for different gas species and axial magnetic field strengths, correlating with differences in pinch uniformity and x-ray emission. Here, the average instability growth is compared to linear MRT theory predictions using the instantaneous acceleration of the best-fit implosion models and characteristic instability wavelength, with the effective Atwood number and seed perturbation size as fit parameters. For high density argon center jets, ionization prior to the arrival of the imploding plasma sheath suggests a heating mechanism consistent with photoionization by XUV self-emission.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

A Thermo‐Flow‐Mechanics‐Fracture Model Coupling a Phase‐Field Interface Approach and Thermo‐Fluid‐Structure Interaction

This work proposes a novel approach for coupling non-isothermal fluid dynamics with fracture mechanics to capture thermal effects within fluid-filled fractures accurately. This method addresses critical aspects of calculating fracture width in enhanced geothermal systems, where the temperature effects of fractures are crucial. The proposed algorithm features an iterative coupling between an interface-capturing phase-field fracture method and interface-tracking thermo-fluid-structure interaction using arbitrary Lagrangian–Eulerian coordinates. We use a phase-field approach to represent fractures and reconstruct the geometry to frame a thermo-fluid-structure interaction problem, resulting in pressure and temperature fields that drive fracture propagation. We developed a novel phase-field interface model accounting for thermal effects, enabling the coupling of quantities specific to the fluid-filled fracture with the phase-field model through the interface between the fracture and the intact solid domain. We provide several numerical examples to demonstrate the capabilities of the proposed algorithm. In particular, we analyze mesh convergence of our phase-field interface model, investigate the effects of temperature on crack width and volume in a static regime, and highlight the method's potential for modeling slowly propagating fractures.

fracture↗
DOE OSTI2021

Modelling a compressible packed bed flow-through washing and deacetylation reactor for corn stover pretreatment

The deacetylation process is a key step in the lignocellulosic biomass pretreatment process that currently is done in batch reactors. A flow-through diffuser type reactor is a promising approach to make the process more efficient and economically favorable. A washing step before deacetylation can further reduce chemical usage and increase the process efficiency. However, the highly compacted biomass particles can cause significant pressure-drop for fluid flow through the bed, which in the worst case can cause self-enhancing pressure-drop leading to reactor blocking. Detailed understanding of the fundamentals of flow through a compressible packed bed is necessary to scale up the washing and deacetylation process. This work closely examines the classical Kozeny Carman equation and proposes modifications to accurately predict the pressure-drop at different compaction states. Here, a compressible packed bed model is further developed based on the interaction between solid mechanics and fluid dynamics. Experimental measurements of water flow through a packed bed of corn stover show good agreement with the developed model. Corn stover washing process scale up simulations show that smaller reactor diameter, higher friction lining material, lower process throughput, and higher process temperature is favorable in stabilizing the bed from extensive compression when the bed height is fixed.

09 BIOMASS FUELS↗
DOE OSTI2021

Modeling the Nonlinear Rheology of Polymer Additive Manufacturing

This report summarizes molecular and continuum simulation studies focused on developing physics - based predictive models for the evolution of polymer molecular order during the nonlinear processing flows of additive manufacturing. Our molecular simulations of polymer elongation flows identified novel mechanisms of fluid dissipation for various polymer architectures that might be harnessed to enhance material processability. In order to predict the complex thermal and flow history of polymer realistic additive manufacturing processes, we have developed and deployed a high - performance mesh - free hydrodynamics module in Sandia's LAMMPS software. This module called RHEO – short for Reproducing Hydrodynamics and Elastic Objects – hybridizes an updated - Lagrange reproducing - kernel method for complex fluids with a bonded particle method (BPM) to capture solidification and solid objects in multiphase flows. In combination, our two methods allow rapid, multiscale characterization of the hydrodynamics and molecular evolution of polymers in realistic processing geometries.

36 MATERIALS SCIENCE↗
NASA NTRS2005

Investigation of Compressibility Effect for Aeropropulsive Shear Flows

Rocket Based Combined Cycle (RBCC) engines operate within a wide range of Mach numbers and altitudes. Fundamental fluid dynamic mechanisms involve complex choking, mass entrainment, stream mixing and wall interactions. The Propulsion Research Center at the University of Alabama in Huntsville is involved in an on- going experimental and numerical modeling study of non-axisymmetric ejector-based combined cycle propulsion systems. This paper attempts to address the modeling issues related to mixing, shear layer/wall interaction in a supersonic Strutjet/ejector flow field. Reynolds Averaged Navier-Stokes (RANS) solutions incorporating turbulence models are sought and compared to experimental measurements to characterize detailed flow dynamics. The effect of compressibility on fluids mixing and wall interactions were investigated using an existing CFD methodology. The compressibility correction to conventional incompressible two- equation models is found to be necessary for the supersonic mixing aspect of the ejector flows based on 2-D simulation results. 3-D strut-base flows involving flow separations were also investigated.

Balasubramanyam, M. S.↗
NASA NTRS2013

Sample Analysis at Mars Instrument Simulator

The Sample Analysis at Mars Instrument Simulator (SAMSIM) is a numerical model dedicated to plan and validate operations of the Sample Analysis at Mars (SAM) instrument on the surface of Mars. The SAM instrument suite, currently operating on the Mars Science Laboratory (MSL), is an analytical laboratory designed to investigate the chemical and isotopic composition of the atmosphere and volatiles extracted from solid samples. SAMSIM was developed using Matlab and Simulink libraries of MathWorks Inc. to provide MSL mission planners with accurate predictions of the instrument electrical, thermal, mechanical, and fluid responses to scripted commands. This tool is a first example of a multi-purpose, full-scale numerical modeling of a flight instrument with the purpose of supplementing or even eliminating entirely the need for a hardware engineer model during instrument development and operation. SAMSIM simulates the complex interactions that occur between the instrument Command and Data Handling unit (C&DH) and all subsystems during the execution of experiment sequences. A typical SAM experiment takes many hours to complete and involves hundreds of components. During the simulation, the electrical, mechanical, thermal, and gas dynamics states of each hardware component are accurately modeled and propagated within the simulation environment at faster than real time. This allows the simulation, in just a few minutes, of experiment sequences that takes many hours to execute on the real instrument. The SAMSIM model is divided into five distinct but interacting modules: software, mechanical, thermal, gas flow, and electrical modules. The software module simulates the instrument C&DH by executing a customized version of the instrument flight software in a Matlab environment. The inputs and outputs to this synthetic C&DH are mapped to virtual sensors and command lines that mimic in their structure and connectivity the layout of the instrument harnesses. This module executes, and thus validates, complex command scripts prior to their up-linking to the SAM instrument. As an output, this module generates synthetic data and message logs at a rate that is similar to the actual instrument.

Benna, Mehdi↗
NASA NTRS1986

Transport of absolute angular momentum in quasi-axisymmetric equatorial jet streams

It is well known that prograde equatorial jet stresses cannot occur in an axisymmetric inviscid fluid, owing to the constraints of local angular momentum conservation. For a viscous fluid, the constraints of mass conservation prevent the formation of any local maximum of absolute angular momentum (m) without a means of transferring m against its gradient (delta m) in the meridional plane. The circumstances under which m can be diffused up-gradient by normal molecular viscosity are derived, and illustrated with reference to numerical simulations of axisymmetric flows in a cylindrical annulus. Viscosity is shown to act so as to tend to expel m from the interior outwards from the rotation axis. Such an effect can produce local super-rotation even in a mechanically isolated fluid. The tendency of viscosity to result in the expulsion of m is shown to be analogous in certain respects to a vorticity-mixing hypothesis for the effects of non-axisymmetric eddies of the zonally-averaged flow. It is shown how the advective and diffusive transport of m by non-axisymmetric eddies can be represented by the Transformed Eulerian Mean meridional circulation and the Eliassen-Palm (EP) flux of Andrews and McIntyre respectively, in the zonal mean. Constraints on the form and direction of the EP flux in an advective/diffusive flow for such eddies are derived, by analogy with similar constraints on the diffusive flux of m due to viscosity.

Read, P. L.↗
NASA NTRS2015

Fluid Shifts

INTRODUCTION: Mechanisms responsible for the ocular structural and functional changes that characterize the visual impairment and intracranial pressure (ICP) syndrome (VIIP) are unclear, but hypothesized to be secondary to the cephalad fluid shift experienced in spaceflight. This study will relate the fluid distribution and compartmentalization associated with long-duration spaceflight with VIIP symptoms. We also seek to determine whether the magnitude of fluid shifts during spaceflight, as well as the VIIP-related effects of those shifts, can be predicted preflight with acute hemodynamic manipulations, and also if lower body negative pressure (LBNP) can reverse the VIIP effects. METHODS: Physiologic variables will be examined pre-, in- and post-flight in 10 International Space Station crewmembers including: fluid compartmentalization (D2O and NaBr dilution); interstitial tissue thickness (ultrasound); vascular dimensions and dynamics (ultrasound and MRI (including cerebrospinal fluid pulsatility)); ocular measures (optical coherence tomography, intraocular pressure, ultrasound); and ICP measures (tympanic membrane displacement, otoacoustic emissions). Pre- and post-flight measures will be assessed while upright, supine and during 15 deg head-down tilt (HDT). In-flight measures will occur early and late during 6 or 12 month missions. LBNP will be evaluated as a countermeasure during HDT and during spaceflight. RESULTS: The first two crewmembers are in the preflight testing phase. Preliminary results characterize the acute fluid shifts experienced from upright, to supine and HDT postures (increased stroke volume, jugular dimensions and measures of ICP) which are reversed with 25 millimeters Hg LBNP. DISCUSSION: Initial results indicate that acute cephalad fluid shifts may be related to VIIP symptoms, but also may be reversible by LBNP. The effect of a chronic fluid shift has yet to be evaluated. Learning Objectives: Current spaceflight VIIP research is described, including novel hardware and countermeasures.

Stenger, M. B.↗
DOE OSTI2021

Predicting Fluid Flow Regime, Permeability, and Diffusivity in Mudrocks from Multiscale Pore Characterisation

In geoenergy applications, mudrocks prevent fluids to leak from temporary (H 2 , CH 4 ) or permanent (CO 2 , radioactive waste) storage/disposal sites and serve as a source and reservoir for unconventional oil and gas. Understanding transport properties integrated with dominant fluid flow mechanisms in mudrocks is essential to better predict the performance of mudrocks within these applications. In this study, small-angle neutron scattering (SANS) experiments were conducted on 71 samples from 13 different sets of mudrocks across the globe to capture the pore structure of nearly the full pore size spectrum (2 nm–5 μm). We develop fractal models to predict transport properties (permeability and diffusivity) based on the SANS-derived pore size distributions. The results indicate that transport phenomena in mudrocks are intrinsically pore size-dependent. Depending on hydrostatic pore pressures, transition flow develops in micropores, slip flow in meso- and macropores, and continuum flow in larger macropores. Fluid flow regimes progress towards larger pore sizes during reservoir depletion or smaller pore sizes during fluid storage, so when pressure is decreased or increased, respectively. Capturing the heterogeneity of mudrocks by considering fractal dimension and tortuosity fractal dimension for defined pore size ranges, fractal models integrate apparent permeability with slip flow, Darcy permeability with continuum flow, and gas diffusivity with diffusion flow in the matrix. This new model of pore size-dependent transport and integrated transport properties using fractal models yields a systematic approach that can also inform multiscale multi-physics models to better understand fluid flow and transport phenomena in mudrocks on the reservoir and basin scale.

36 MATERIALS SCIENCE↗