Thermodynamically Constrained Averaging Theory Approach for Modeling Flow and Transport Phenomena in Porous Medium Systems: 4. Species Transport Fundamentals

Modeling flow of Carreau fluids in porous media

Carreau fluids occur routinely in porous medium systems for a range of applications, and the dependence of the viscosity for such fluids on the rate of strain tensor poses challenges to modeling at an averaged macroscale. Traditional approaches for macroscale modeling such flows have relied upon experimental observations of flows for generalized Newtonian fluids (GNFs) and a phenomenological approach referred to herein as the shift factor. A recently developed approach based upon averaging conservation and thermodynamic equations from the microscale for Cross model GNFs is extended to the case of Carreau fluids and shown to predict the flow through both isotropic and anisotropic media accurately without the need for GNF-flow experiments. The model is formulated in terms of rheological properties, a standard Newtonian resistance tensor, and a length-scale tensor, which does require estimation. An approach based upon measures of the morphology and topology of the pore space is developed to approximate this length-scale tensor. Thus, this work provides the missing components needed to predict Carreau GNF macroscale flow with only rheological information for the fluid and analysis of the pore morphology and topology independent of any fluid flow experiments. Accuracy of predictions based upon this approach is quantified, and extension to other GNFs is straightforward.

A Hierarchy of Probability, Fluid and Generalized Densities for the Eulerian Velocivolumetric Description of Fluid Flow, for New Families of Conservation Laws

The Reynolds transport theorem occupies a central place in continuum mechanics, providing a generalized integral conservation equation for the transport of any conserved quantity within a fluid or material volume, which can be connected to its corresponding differential equation. Recently, a more generalized framework was presented for this theorem, enabling parametric transformations between positions on a manifold or in any generalized coordinate space, exploiting the underlying continuous multivariate (Lie) symmetries of a vector or tensor field associated with a conserved quantity. We explore the implications of this framework for fluid flow systems, based on an Eulerian velocivolumetric (position-velocity) description of fluid flow. The analysis invokes a hierarchy of five probability density functions, which by convolution are used to define five fluid densities and generalized densities relevant to this description. We derive 11 formulations of the generalized Reynolds transport theorem for different choices of the coordinate space, parameter space and density, only the first of which is commonly known. These are used to generate a table of integral and differential conservation laws applicable to each formulation, for eight important conserved quantities (fluid mass, species mass, linear momentum, angular momentum, energy, charge, entropy and probability). The findings substantially expand the set of conservation laws for the analysis of fluid flow and dynamical systems.

A continuum mechanical framework for modeling tumor growth and treatment in two- and three-phase systems

The growth and treatment of tumors is an important problem to society that involves the manifestation of cellular phenomena at length scales on the order of centimeters. Continuum mechanical approaches are being increasingly used to model tumors at the largest length scales of concern. The issue of how to best connect such descriptions to smaller-scale descriptions remains open. We formulate a framework to derive macroscale models of tumor behavior using the thermodynamically constrained averaging theory (TCAT), which provides a firm connection with the microscale and constraints on permissible forms of closure relations. We build on developments in the porous medium mechanics literature to formulate fundamental entropy inequality expressions for a general class of three-phase, compositional models at the macroscale. We use the general framework derived to formulate two classes of models, a two-phase model and a three-phase model. The general TCAT framework derived forms the basis for a wide range of potential models of varying sophistication, which can be derived, approximated, and applied to understand not only tumor growth but also the effectiveness of various treatment modalities.

In Silico Prediction of Food Properties: A Multiscale Perspective

Several open software packages have popularized modeling and simulation strategies at the food product scale. Food processing and key digestion steps can be described in 3D using the principles of continuum mechanics. However, compared to other branches of engineering, the necessary transport, mechanical, chemical, and thermodynamic properties have been insufficiently tabulated and documented. Natural variability, accented by food evolution during processing and deconstruction, requires considering composition and structure-dependent properties. This review presents practical approaches where the premises for modeling and simulation start at a so-called “microscopic” scale where constituents or phase properties are known. The concept of microscopic or ground scale is shown to be very flexible from atoms to cellular structures. Zooming in on spatial details tends to increase the overall cost of simulations and the integration over food regions or time scales. The independence of scales facilitates the reuse of calculations and makes multiscale modeling capable of meeting food manufacturing needs. On one hand, new image-modeling strategies without equations or meshes are emerging. On the other hand, complex notions such as compositional effects, multiphase organization, and non-equilibrium thermodynamics are naturally incorporated in models without linearization or simplifications. Multiscale method’s applicability to hierarchically predict food properties is discussed with comprehensive examples relevant to food science, engineering and packaging. Entropy-driven properties such as transport and sorption are emphasized to illustrate how microscopic details bring new degrees of freedom to explore food-specific concepts such as safety, bioavailability, shelf-life and food formulation. Routes for performing spatial and temporal homogenization with and without chemical details are developed. Creating a community sharing computational codes, force fields, and generic food structures is the next step and should be encouraged. This paper provides a framework for the transfer of results from other fields and the development of methods specific to the food domain.

Toward a New Generation of Two-Fluid Flow Models Based on the Thermodynamically-Constrained Averaging Theory

Traditional models of two-fluid flow through porous media at the macroscale have existed for nearly a century. These phenomenological models are not firmly connected to the microscale; thermodynamic constraints are not enforced; empirical closure relations are well known to be hysteretic; fluid pressures are typically assumed to be in a local equilibrium state with fluid saturations; and important quantities such as interfacial and curvilinear geometric extents, tensions, and curvatures, known to be important from microscale studies, do not explicitly appear in traditional macroscale models. Despite these shortcomings, the traditional model for two-fluid flow in porous media has been extensively studied to develop efficient numerical approximation methods, experimental and surrogate measure parameterization approaches, and convenient pre- and post-processing environments; and they have been applied in a large number of applications from a variety of fields. The thermodynamically constrained averaging theory (TCAT) was developed to overcome the limitations associated with traditional approaches, and we consider here issues associated with the closure of this new generation of models. It has been shown that a hysteretic-free state equation exists based upon integral geometry that relates changes in volume fractions, capillary pressure, interfacial areas, and the Euler characteristic. We show an analysis of how this state equation can be parameterized with a relatively small amount of data. We also formulate a state equation for resistance coefficients that we show to be hysteretic free, unlike traditional relative permeability models. Lastly, we comment on the open issues remaining for this new generation of models.

A monolithic multiphase porous medium framework for (a-)vascular tumor growth

We present a dynamic vascular tumor model combining a multiphase porous medium framework for avascular tumor growth in a consistent Arbitrary Lagrangian Eulerian formulation and a novel approach to incorporate angiogenesis. The multiphase model is based on Thermodynamically Constrained Averaging Theory and comprises the extracellular matrix as a porous solid phase and three fluid phases: (living and necrotic) tumor cells, host cells and the interstitial fluid. Angiogenesis is modeled by treating the neovasculature as a proper additional phase with volume fraction or blood vessel density. This allows us to define consistent inter-phase exchange terms between the neovasculature and the interstitial fluid. As a consequence, transcapillary leakage and lymphatic drainage can be modeled. By including these important processes we are able to reproduce the increased interstitial pressure in tumors which is a crucial factor in drug delivery and, thus, therapeutic outcome. Different coupling schemes to solve the resulting five-phase problem are realized and compared with respect to robustness and computational efficiency. We find that a fully monolithic approach is superior to both the standard partitioned and a hybrid monolithic-partitioned scheme for a wide range of parameters. The flexible implementation of the novel model makes further extensions (e.g., inclusion of additional phases and species) straightforward.

Thermodynamically Constrained Averaging Theory: Principles, Model Hierarchies, and Deviation Kinetic Energy Extensions

The thermodynamically constrained averaging theory (TCAT) is a comprehensive theory used to formulate hierarchies of multiphase, multiscale models that are closed based upon the second law of thermodynamics. The rate of entropy production is posed in terms of the product of fluxes and forces of dissipative processes. The attractive features of TCAT include consistency across disparate length scales; thermodynamic consistency across scales; the inclusion of interfaces and common curves as well as phases; the development of kinematic equations to provide closure relations for geometric extent measures; and a structured approach to model building. The elements of the TCAT approach are shown; the ways in which each of these attractive features emerge from the TCAT approach are illustrated; and a review of the hierarchies of models that have been formulated is provided. Because the TCAT approach is mathematically involved, we illustrate how this approach can be applied by leveraging existing components of the theory that can be applied to a wide range of applications. This can result in a substantial reduction in formulation effort compared to a complete derivation while yielding identical results. Lastly, we note the previous neglect of the deviation kinetic energy, which is not important in slow porous media flows, formulate the required equations to extend the theory, and comment on applications for which the new components would be especially useful. This work should serve to make TCAT more accessible for applications, thereby enabling higher fidelity models for applications such as turbulent multiphase flows.

Translational models of tumor angiogenesis: A nexus of in silico and in vitro models

Emerging evidence shows that endothelial cells are not only the building blocks of vascular networks that enable oxygen and nutrient delivery throughout a tissue but also serve as a rich resource of angiocrine factors. Endothelial cells play key roles in determining cancer progression and response to anti-cancer drugs. Furthermore, the endothelium-specific deposition of extracellular matrix is a key modulator of the availability of angiocrine factors to both stromal and cancer cells. Considering tumor vascular network as a decisive factor in cancer pathogenesis and treatment response, these networks need to be an inseparable component of cancer models. Both computational and in vitro experimental models have been extensively developed to model tumor-endothelium interactions. While informative, they have been developed in different communities and do not yet represent a comprehensive platform. In this review, we overview the necessity of incorporating vascular networks for both in vitro and in silico cancer models and discuss recent progresses and challenges of in vitro experimental microfluidic cancer vasculature-on-chip systems and their in silico counterparts. We further highlight how these two approaches can merge together with the aim of presenting a predictive combinatorial platform for studying cancer pathogenesis and testing the efficacy of single or multi-drug therapeutics for cancer treatment.

An overview of multiphase cartilage mechanical modelling and its role in understanding function and pathology

There is a long history of mathematical and computational modelling with the objective of understanding the mechanisms governing cartilage׳s remarkable mechanical performance. Nonetheless, despite sophisticated modelling development, simulations of cartilage have consistently lagged behind structural knowledge and thus the relationship between structure and function in cartilage is not fully understood. However, in the most recent generation of studies, there is an emerging confluence between our structural knowledge and the structure represented in cartilage modelling. This raises the prospect of further refinement in our understanding of cartilage function and also the initiation of an engineering-level understanding for how structural degradation and ageing relates to cartilage dysfunction and pathology, as well as informing the potential design of prospective interventions. Aimed at researchers entering the field of cartilage modelling, we thus review the basic principles of cartilage models, discussing the underlying physics and assumptions in relatively simple settings, whilst presenting the derivation of relatively parsimonious multiphase cartilage models consistent with our discussions. We proceed to consider modern developments that start aligning the structure captured in the models with observed complexities. This emphasises the challenges associated with constitutive relations, boundary conditions, parameter estimation and validation in cartilage modelling programmes. Consequently, we further detail how both experimental interrogations and modelling developments can be utilised to investigate and reduce such difficulties before summarising how cartilage modelling initiatives may improve our understanding of cartilage ageing, pathology and intervention.

Predicting the growth of glioblastoma multiforme spheroids using a multiphase porous media model

Tumor spheroids constitute an effective in vitro tool to investigate the avascular stage of tumor growth. These three-dimensional cell aggregates reproduce the nutrient and proliferation gradients found in the early stages of cancer and can be grown with a strict control of their environmental conditions. In the last years, new experimental techniques have been developed to determine the effect of mechanical stress on the growth of tumor spheroids. These studies report a reduction in cell proliferation as a function of increasingly applied stress on the surface of the spheroids. This work presents a specialization for tumor spheroid growth of a previous more general multiphase model. The equations of the model are derived in the framework of porous media theory, and constitutive relations for the mass transfer terms and the stress are formulated on the basis of experimental observations. A set of experiments is performed, investigating the growth of U-87MG spheroids both freely growing in the culture medium and subjected to an external mechanical pressure induced by a Dextran solution. The growth curves of the model are compared to the experimental data, with good agreement for both the experimental settings. A new mathematical law regulating the inhibitory effect of mechanical compression on cancer cell proliferation is presented at the end of the paper. This new law is validated against experimental data and provides better results compared to other expressions in the literature.

Statistical mechanics of unsaturated porous media

We explore a mean-field theory of fluid imbibition and drainage through permeable porous solids. In the limit of vanishing inertial and viscous forces, the theory predicts the hysteretic "retention curves" relating the capillary pressure applied across a connected domain to its degree of saturation in wetting fluid in terms of known surface energies and void space geometry. To avoid complicated calculations, we adopt the simplest statistical mechanics, in which a pore interacts with its neighbors through narrow openings called "necks," while being either full or empty of wetting fluid. We show how the main retention curves can be calculated from the statistical distribution of two dimensionless parameters λ and α measuring the specific areas of, respectively, neck cross section and wettable pore surface relative to pore volume. The theory attributes hysteresis of these curves to collective first-order phase transitions. We illustrate predictions with a porous domain consisting of a random packing of spheres, show that hysteresis strength grows with λ and weakens as the distribution of α broadens, and reproduce the behavior of Haines jumps observed in recent experiments on an ordered pore network.

An analysis platform for multiscale hydrogeologic modeling with emphasis on hybrid multiscale methods

One of the most significant challenges faced by hydrogeologic modelers is the disparity between the spatial and temporal scales at which fundamental flow, transport, and reaction processes can best be understood and quantified (e.g., microscopic to pore scales and seconds to days) and at which practical model predictions are needed (e.g., plume to aquifer scales and years to centuries). While the multiscale nature of hydrogeologic problems is widely recognized, technological limitations in computation and characterization restrict most practical modeling efforts to fairly coarse representations of heterogeneous properties and processes. For some modern problems, the necessary level of simplification is such that model parameters may lose physical meaning and model predictive ability is questionable for any conditions other than those to which the model was calibrated. Recently, there has been broad interest across a wide range of scientific and engineering disciplines in simulation approaches that more rigorously account for the multiscale nature of systems of interest. In this article, we review a number of such approaches and propose a classification scheme for defining different types of multiscale simulation methods and those classes of problems to which they are most applicable. Our classification scheme is presented in terms of a flowchart (Multiscale Analysis Platform), and defines several different motifs of multiscale simulation. Within each motif, the member methods are reviewed and example applications are discussed. We focus attention on hybrid multiscale methods, in which two or more models with different physics described at fundamentally different scales are directly coupled within a single simulation. Very recently these methods have begun to be applied to groundwater flow and transport simulations, and we discuss these applications in the context of our classification scheme. As computational and characterization capabilities continue to improve, we envision that hybrid multiscale modeling will become more common and also a viable alternative to conventional single-scale models in the near future.

A tumor growth model with deformable ECM

Existing tumor growth models based on fluid analogy for the cells do not generally include the extracellular matrix (ECM), or if present, take it as rigid. The three-fluid model originally proposed by the authors and comprising tumor cells (TC), host cells (HC), interstitial fluid (IF) and an ECM, considered up to now only a rigid ECM in the applications. This limitation is here relaxed and the deformability of the ECM is investigated in detail. The ECM is modeled as a porous solid matrix with Green-elastic and elasto-visco-plastic material behavior within a large strain approach. Jauman and Truesdell objective stress measures are adopted together with the deformation rate tensor. Numerical results are first compared with those of a reference experiment of a multicellular tumor spheroid (MTS) growing in vitro, then three different tumor cases are studied: growth of an MTS in a decellularized ECM, growth of a spheroid in the presence of host cells and growth of a melanoma. The influence of the stiffness of the ECM is evidenced and comparison with the case of a rigid ECM is made. The processes in a deformable ECM are more rapid than in a rigid ECM and the obtained growth pattern differs. The reasons for this are due to the changes in porosity induced by the tumor growth. These changes are inhibited in a rigid ECM. This enhanced computational model emphasizes the importance of properly characterizing the biomechanical behavior of the malignant mass in all its components to correctly predict its temporal and spatial pattern evolution.

A two-phase model of plantar tissue: a step toward prediction of diabetic foot ulceration

A new computational model, based on the thermodynamically constrained averaging theory, has been recently proposed to predict tumor initiation and proliferation. A similar mathematical approach is proposed here as an aid in diabetic ulcer prevention. The common aspects at the continuum level are the macroscopic balance equations governing the flow of the fluid phase, diffusion of chemical species, tissue mechanics, and some of the constitutive equations. The soft plantar tissue is modeled as a two-phase system: a solid phase consisting of the tissue cells and their extracellular matrix, and a fluid one (interstitial fluid and dissolved chemical species). The solid phase may become necrotic depending on the stress level and on the oxygen availability in the tissue. Actually, in diabetic patients, peripheral vascular disease impacts tissue necrosis; this is considered in the model via the introduction of an effective diffusion coefficient that governs transport of nutrients within the microvasculature. The governing equations of the mathematical model are discretized in space by the finite element method and in time domain using the θ-Wilson Method. While the full mathematical model is developed in this paper, the example is limited to the simulation of several gait cycles of a healthy foot.

USNCTAM perspectives on mechanics in medicine

Over decades, the theoretical and applied mechanics community has developed sophisticated approaches for analysing the behaviour of complex engineering systems. Most of these approaches have targeted systems in the transportation, materials, defence and energy industries. Applying and further developing engineering approaches for understanding, predicting and modulating the response of complicated biomedical processes not only holds great promise in meeting societal needs, but also poses serious challenges. This report, prepared for the US National Committee on Theoretical and Applied Mechanics, aims to identify the most pressing challenges in biological sciences and medicine that can be tackled within the broad field of mechanics. This echoes and complements a number of national and international initiatives aiming at fostering interdisciplinary biomedical research. This report also comments on cultural/educational challenges. Specifically, this report focuses on three major thrusts in which we believe mechanics has and will continue to have a substantial impact. (i) Rationally engineering injectable nano/microdevices for imaging and therapy of disease. Within this context, we discuss nanoparticle carrier design, vascular transport and adhesion, endocytosis and tumour growth in response to therapy, as well as uncertainty quantification techniques to better connect models and experiments. (ii) Design of biomedical devices, including point-of-care diagnostic systems, model organ and multi-organ microdevices, and pulsatile ventricular assistant devices. (iii) Mechanics of cellular processes, including mechanosensing and mechanotransduction, improved characterization of cellular constitutive behaviour, and microfluidic systems for single-cell studies.

Averaging Theory for Description of Environmental Problems: What Have We Learned?

Advances in Water Resources has been a prime archival source for implementation of averaging theories in changing the scale at which processes of importance in environmental modeling are described. Thus in celebration of the 35th year of this journal, it seems appropriate to assess what has been learned about these theories and about their utility in describing systems of interest. We review advances in understanding and use of averaging theories to describe porous medium flow and transport at the macroscale, an averaged scale that models spatial variability, and at the megascale, an integral scale that only considers time variation of system properties. We detail physical insights gained from the development and application of averaging theory for flow through porous medium systems and for the behavior of solids at the macroscale. We show the relationship between standard models that are typically applied and more rigorous models that are derived using modern averaging theory. We discuss how the results derived from averaging theory that are available can be built upon and applied broadly within the community. We highlight opportunities and needs that exist for collaborations among theorists, numerical analysts, and experimentalists to advance the new classes of models that have been derived. Lastly, we comment on averaging developments for rivers, estuaries, and watersheds.

A multiphase model for three-dimensional tumor growth

Several mathematical formulations have analyzed the time-dependent behaviour of a tumor mass. However, most of these propose simplifications that compromise the physical soundness of the model. Here, multiphase porous media mechanics is extended to model tumor evolution, using governing equations obtained via the Thermodynamically Constrained Averaging Theory (TCAT). A tumor mass is treated as a multiphase medium composed of an extracellular matrix (ECM); tumor cells (TC), which may become necrotic depending on the nutrient concentration and tumor phase pressure; healthy cells (HC); and an interstitial fluid (IF) for the transport of nutrients. The equations are solved by a Finite Element method to predict the growth rate of the tumor mass as a function of the initial tumor-to-healthy cell density ratio, nutrient concentration, mechanical strain, cell adhesion and geometry. Results are shown for three cases of practical biological interest such as multicellular tumor spheroids (MTS) and tumor cords. First, the model is validated by experimental data for time-dependent growth of an MTS in a culture medium. The tumor growth pattern follows a biphasic behaviour: initially, the rapidly growing tumor cells tend to saturate the volume available without any significant increase in overall tumor size; then, a classical Gompertzian pattern is observed for the MTS radius variation with time. A core with necrotic cells appears for tumor sizes larger than 150 μm, surrounded by a shell of viable tumor cells whose thickness stays almost constant with time. A formula to estimate the size of the necrotic core is proposed. In the second case, the MTS is confined within a healthy tissue. The growth rate is reduced, as compared to the first case - mostly due to the relative adhesion of the tumor and healthy cells to the ECM, and the less favourable transport of nutrients. In particular, for tumor cells adhering less avidly to the ECM, the healthy tissue is progressively displaced as the malignant mass grows, whereas tumor cell infiltration is predicted for the opposite condition. Interestingly, the infiltration potential of the tumor mass is mostly driven by the relative cell adhesion to the ECM. In the third case, a tumor cord model is analyzed where the malignant cells grow around microvessels in a 3D geometry. It is shown that tumor cells tend to migrate among adjacent vessels seeking new oxygen and nutrient. This model can predict and optimize the efficacy of anticancer therapeutic strategies. It can be further developed to answer questions on tumor biophysics, related to the effects of ECM stiffness and cell adhesion on tumor cell proliferation.

Thermodynamically Constrained Averaging Theory Approach for Modeling Flow and Transport Phenomena in Porous Medium Systems: 8. Interface and Common Curve Dynamics

This work is the eighth in a series that develops the fundamental aspects of the thermodynamically constrained averaging theory (TCAT) that allows for a systematic increase in the scale at which multiphase transport phenomena is modeled in porous medium systems. In these systems, the explicit locations of interfaces between phases and common curves, where three or more interfaces meet, are not considered at scales above the microscale. Rather, the densities of these quantities arise as areas per volume or length per volume. Modeling of the dynamics of these measures is an important challenge for robust models of flow and transport phenomena in porous medium systems, as the extent of these regions can have important implications for mass, momentum, and energy transport between and among phases, and formulation of a capillary pressure relation with minimal hysteresis. These densities do not exist at the microscale, where the interfaces and common curves correspond to particular locations. Therefore, it is necessary for a well-developed macroscale theory to provide evolution equations that describe the dynamics of interface and common curve densities. Here we point out the challenges and pitfalls in producing such evolution equations, develop a set of such equations based on averaging theorems, and identify the terms that require particular attention in experimental and computational efforts to parameterize the equations. We use the evolution equations developed to specify a closed two-fluid-phase flow model.

Thermodynamically Constrained Averaging Theory Approach for Modeling Flow and Transport Phenomena in Porous Medium Systems: 7. Single-Phase Megascale Flow Models

This work is the seventh in a series that introduces and employs the thermodynamically constrained averaging theory (TCAT) for modeling flow and transport in multiscale porous medium systems. This paper expands the previous analyses in the series by developing models at a scale where spatial variations within the system are not considered. Thus the time variation of variables averaged over the entire system is modeled in relation to fluxes at the boundary of the system. This implementation of TCAT makes use of conservation equations for mass, momentum, and energy as well as an entropy balance. Additionally, classical irreversible thermodynamics is assumed to hold at the microscale and is averaged to the megascale, or system scale. The fact that the local equilibrium assumption does not apply at the megascale points to the importance of obtaining closure relations that account for the large-scale manifestation of small-scale variations. Example applications built on this foundation are suggested to stimulate future work.

Thermodynamically Constrained Averaging Theory Approach for Modeling Flow and Transport Phenomena in Porous Medium Systems: 5. Single-Fluid-Phase Transport

This work is the fifth in a series of papers on the thermodynamically constrained averaging theory (TCAT) approach for modeling flow and transport phenomena in multiscale porous medium systems. The general TCAT framework and the mathematical foundation presented in previous works are used to develop models that describe species transport and single-fluid-phase flow through a porous medium system in varying physical regimes. Classical irreversible thermodynamics formulations for species in fluids, solids, and interfaces are developed. Two different approaches are presented, one that makes use of a momentum equation for each entity along with constitutive relations for species diffusion and dispersion, and a second approach that makes use of a momentum equation for each species in an entity. The alternative models are developed by relying upon different approaches to constrain an entropy inequality using mass, momentum, and energy conservation equations. The resultant constrained entropy inequality is simplified and used to guide the development of closed models. Specific instances of dilute and non-dilute systems are examined and compared to alternative formulation approaches.

Thermodynamically Constrained Averaging Theory Approach for Modeling Flow and Transport Phenomena in Porous Medium Systems: 5. Single-Fluid-Phase Transport

This work is the fifth in a series of papers on the thermodynamically constrained averaging theory (TCAT) approach for modeling flow and transport phenomena in multiscale porous medium systems. The general TCAT framework and the mathematical foundation presented in previous works are used to develop models that describe species transport and single-fluid-phase flow through a porous medium system in varying physical regimes. Classical irreversible thermodynamics formulations for species in fluids, solids, and interfaces are developed. Two different approaches are presented, one that makes use of a momentum equation for each entity along with constitutive relations for species diffusion and dispersion, and a second approach that makes use of a momentum equation for each species in an entity. The alternative models are developed by relying upon different approaches to constrain an entropy inequality using mass, momentum, and energy conservation equations. The resultant constrained entropy inequality is simplified and used to guide the development of closed models. Specific instances of dilute and non-dilute systems are examined and compared to alternative formulation approaches.

Thermodynamically Constrained Averaging Theory Approach for Modeling Flow and Transport Phenomena in Porous Medium Systems: 5. Single-Fluid-Phase Transport

This work is the fifth in a series of papers on the thermodynamically constrained averaging theory (TCAT) approach for modeling flow and transport phenomena in multiscale porous medium systems. The general TCAT framework and the mathematical foundation presented in previous works are used to develop models that describe species transport and single-fluid-phase flow through a porous medium system in varying physical regimes. Classical irreversible thermodynamics formulations for species in fluids, solids, and interfaces are developed. Two different approaches are presented, one that makes use of a momentum equation for each entity along with constitutive relations for species diffusion and dispersion, and a second approach that makes use of a momentum equation for each species in an entity. The alternative models are developed by relying upon different approaches to constrain an entropy inequality using mass, momentum, and energy conservation equations. The resultant constrained entropy inequality is simplified and used to guide the development of closed models. Specific instances of dilute and non-dilute systems are examined and compared to alternative formulation approaches.

Thermodynamically Constrained Averaging Theory Approach for Modeling Flow and Transport Phenomena in Porous Medium Systems: 4. Species Transport Fundamentals

This work is the fourth in a series of papers on the thermodynamically constrained averaging theory (TCAT) approach for modeling flow and transport phenomena in multiscale porous medium systems. The general TCAT framework and the mathematical foundation presented in previous works are built upon by formulating macroscale models for conservation of mass, momentum, and energy, and the balance of entropy for a species in a phase volume, interface, and common curve. In addition, classical irreversible thermodynamic relations for species in entities are averaged from the microscale to the macroscale. Finally, we comment on alternative approaches that can be used to connect species and entity conservation equations to a constrained system entropy inequality, which is a key component of the TCAT approach. The formulations detailed in this work can be built upon to develop models for species transport and reactions in a variety of multiphase systems.

Thermodynamically Constrained Averaging Theory Approach for Modeling Flow and Transport Phenomena in Porous Medium Systems: 4. Species Transport Fundamentals

This work is the fourth in a series of papers on the thermodynamically constrained averaging theory (TCAT) approach for modeling flow and transport phenomena in multiscale porous medium systems. The general TCAT framework and the mathematical foundation presented in previous works are built upon by formulating macroscale models for conservation of mass, momentum, and energy, and the balance of entropy for a species in a phase volume, interface, and common curve. In addition, classical irreversible thermodynamic relations for species in entities are averaged from the microscale to the macroscale. Finally, we comment on alternative approaches that can be used to connect species and entity conservation equations to a constrained system entropy inequality, which is a key component of the TCAT approach. The formulations detailed in this work can be built upon to develop models for species transport and reactions in a variety of multiphase systems.

A multiscale theoretical model for diffusive mass transfer in cellular biological media

An integrated methodology is developed for the theoretical analysis of solute transport and reaction in cellular biological media, such as tissues, microbial flocs, and biofilms. First, the method of local spatial averaging with a weight function is used to establish the equation which describes solute conservation at the cellular biological medium scale, starting with a continuum-based formulation of solute transport at finer spatial scales. Second, an effective-medium model is developed for the self-consistent calculation of the local diffusion coefficient in the cellular biological medium, including the effects of the structural heterogeneity of the extra-cellular space and the reversible adsorption to extra-cellular polymers. The final expression for the local effective diffusion coefficient is: D(Abeta)=lambda(beta)D(Aupsilon), where D(Aupsilon) is the diffusion coefficient in water, and lambda(beta) is a function of the composition and fundamental geometric and physicochemical system properties, including the size of solute molecules, the size of extra-cellular polymer fibers, and the mass permeability of the cell membrane. Furthermore, the analysis sheds some light on the function of the extra-cellular hydrogel as a diffusive barrier to solute molecules approaching the cell membrane, and its implications on the transport of chemotherapeutic agents within a cellular biological medium. Finally, the model predicts the qualitative trend as well as the quantitative variability of a large number of published experimental data on the diffusion coefficient of oxygen in cell-entrapping gels, microbial flocs, biofilms, and mammalian tissues.

Pore-scale investigation of viscous coupling effects for two-phase flow in porous media

Recent studies have revealed that viscous coupling effects in immiscible two-phase flow, caused by momentum transfer between the two fluid phases, can be important in porous medium systems. In this work, we use a three-dimensional parallel processing version of a two-fluid-phase lattice Boltzmann (LB) model to investigate this phenomenon. A multiple-relaxation-time (MRT) approximation of the LB equations is used in the simulator, which leads to a viscosity-independent velocity field. We validate our model by verifying the velocity profile for two-phase flow through a channel with a square cross section. We then simulate co-current flow through a sphere-pack porous medium and obtain correlations of the relative permeabilities as a function of capillary number, wettability, and the fluid viscosities. The results are qualitatively consistent with experimental observations. In addition, we calculate the generalized permeability coefficients and show that the coupling coefficients are significant and the matrix is nonsymmetric. We also find a strong correlation between the relative permeability and interfacial area between fluids, indicating that both the common extension of Darcy's Law and the generalized formulation accounting for viscous coupling effects do not provide adequate insight into two-phase flow processes in porous media. This work lends additional support for the hypothesis that interfacial area is a key variable for multiphase flow in porous medium systems.