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

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 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.

Generalized Newtonian fluid flow in porous media

Single-fluid-phase porous medium systems are typically modeled at an averaged length scale termed the macroscale, and Darcy's law is typically relied upon as an approximation of the momentum equation under Stokes flow conditions. Standard approaches for modeling macroscale single-fluid-phase flow of generalized Newtonian fluids (GNFs) extend the standard Newtonian model based upon Darcy's law using an effective viscosity and assuming that the intrinsic permeability is invariant with respect to fluid properties. This approach results in a need to perform an experiment for a non-Newtonian fluid, the introduction of effective parameters that are not tied to known microscale physics, and uncertainty regarding the dependencies of the fitted empirical parameter on system properties. We use the thermodynamically constrained averaging theory (TCAT) to examine the formulation and closure of a macroscale model for GNF flow that is consistent with microscale conservation principles and the second law of thermodynamics. A direct connection between microscale and macroscale quantities is used to formulate an expression for interphase momentum transfer for GNF flow in porous medium systems. Darcy's law is shown to approximate momentum transfer from the fluid phase to the solid phase. Momentum transfer is found to depend on the viscosity at the solid surface, which is only invariant for Newtonian flow. TCAT is used to derive a macroscale equation for the hydraulic resistance based on accessible fluid and solid properties. This hydraulic resistance may be used in the same way that hydraulic conductivity is typically used to model flow at the macroscale, and it includes parameters that can be calculated a priori, without the need to carry out microscale simulations, or experiments, for any GNF. The TCAT approach is validated for four model isotropic and anisotropic media and five Cross-model fluids. The traditional shift factor and effective viscosity are related to the newly derived TCAT model, shedding new light on this common empirical approach. The results from this work form a basis for the modeling of GNF flow in porous medium systems under Stokes flow, which is predictive given the rheological properties of the GNF and the resistance observed for Newtonian flow.

Microscale modeling of nondilute flow and transport in porous medium systems

Nondilute transport occurs routinely in porous medium systems. Experimental observations have revealed effects that seemingly depend upon density, viscosity, velocity, and chemical activity. Macroscale models based upon averaged behavior over many pores have been relied upon to describe such systems to date, which require parametrization of important physical phenomena in material coefficients. To advance fundamental understanding of these complex systems, we examine nondilute transport from a fundamental microscale, or pore-scale, continuum modeling perspective. We approximate the solution of a model based upon the variable-density Navier-Stokes equations and a nondilute species transport equation. Known dependencies of the densities, viscosities, chemical activity, and diffusion for a salt solution on chemical composition are included in the model. Microscale model solutions are averaged to the macroscale and compared with extant experimental observations. Investigation of the effects of various physical phenomena on the microscale velocity distribution and the observed macroscale dispersion are considered using dimensional analysis and constrained simulations. Simulation results are used to explain observed experimental results in light of underlying mechanisms. Conditions under which the various physicochemical effects investigated are important are revealed.

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.

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.

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.

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.

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.

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: 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.