Investigating the role of tortuosity in the Kozeny-Carman equation

Size: px
Start display at page:

Download "Investigating the role of tortuosity in the Kozeny-Carman equation"

Transcription

1 Investigating the role of tortuosity in the Kozeny-Carman equation Rebecca Allen, Shuyu Sun King Abdullah University of Science and Technology Sept 30, / 52

2 The objective of this work is to evaluate the role of tortuosity on fitting simulation data to the Kozeny-Carman equation. Outline: Review of Kozeny-Carman eqn. Review of different tortuosity definitions Obtaining permeability and tortuosity from pore-scale modeling Example geometries: in-line array, staggered-array, overlapping squares Fitting simulation data to Kozeny-Carman eqn. Conclusion 2 / 52

3 Kozeny-Carman equation derived from theory by treating porous media as comprised of parallel and uniform channels relates permeability to pore-structure properties: k = φ3 cs 2 = φ3 βτ 2 S 2 (1) Carman, 1937 & 1939: Shape of cross-section Kozeny constant, c Figure: Tortuosity: τ = λ/l Circle 2 Ellipse (major/minor = 2) 2.13 Ellipse (major/minor = 10) 2.45 Rectangle (width/height = 1) 1.78 Rectangle (width/height = 2) 1.94 Rectangle (width/height = 10) 2.65 Rectangle (width/height = ) 3 In granular beds, how can we compute tortuosity? Is it a function of porosity? 3 / 52

4 Past work on measuring tortuosity Many authors have theoretically or empirically derived tortuosity as a function of porosity. Figure: Porosity vs tortuosity trends from literature: generally, τ 1, and τ 1 as φ 1. 4 / 52

5 Study Samples Considered Tortuosity vs Porosity Fit Tortuosity Type Maxwell, 1881 array of spheres (3D), dilute τ = 1 + suspension, non-conducting electrical (1 φ) 2 conductivity Rayleigh, 1892 array of cylinders (2D) τ = 2 φ diffusion? Mackie & Meares, 2 φ diffusion of electrolytes in membrane τ = ( 1955 φ )2 diffusion Weissberg, 1963 bed of uniform spheres (applicable to overlapping, non-uniform spheres) τ = lnφ diffusion Kim et al, 1987 isotropic system, 0 < φ < 0.5 τ = φ 0.4 diffusion Koponen et al, 1996 Koponen et al, 1997 Matyka et al, 2008 Duda et al, 2011 Pisani, 2011 Liu & Kitanidis, D random overlapping mono-sized squares, 0.5 < φ < 1 2D random overlapping mono-sized squares, 0.4 < φ < 1 2D random overlapping mono-sized squares 2D freely overlapping squares random, partial overlapping shapes isotropic grain (spherical), staggered 0.25 < φ < 0.5 τ = (1 φ) 1 φ τ = 1 + a (φ φc ) m, a = 0.65, m = 0.19 τ 1 R φ S τ 1 (1 φ) γ, γ = 1/2 τ = 1 1 α(1 φ), α=shape factor τ = φ 1 m , m = 1.28 hydraulic hydraulic hydraulic hydraulic diffusion electrical conductivity * as referenced in Ochoa-Tapia et al 1994 ** as referenced in Shen & Chen 2007, and Boudreau / 52

6 However recent work states that any definition of tortuosity (i.e., hydraulic and diffusive) is not a function of porosity but rather a function of the pore geometry only and is a tensorial property (Valdes-Parada et al. 2011, Liu & Kitanidis 2013). 6 / 52

7 Different forms of tortuosity The tortuous nature of fluid flow through porous media: The tortuous nature of a solute diffusing through porous media: τ h = path traveled by fluid domain unit length τ d = path traveled by diffusing solute domain unit length 7 / 52

8 Different forms of tortuosity The tortuous nature of fluid flow through porous media: The tortuous nature of a solute diffusing through porous media: τ h = path traveled by fluid domain unit length τ d = path traveled by diffusing solute domain unit length The Kozeny-Carman equation uses hydraulic tortuosity: k = φ3 cτ h 2 S 2 = φ3 βs 2 Effective diffusivity is the binary diffusion scaled by diffusive tortuosity D eff = D τ d = Dτ d 8 / 52

9 Different forms of tortuosity The tortuous nature of fluid flow through porous media: The tortuous nature of a solute diffusing through porous media: τ h = path traveled by fluid domain unit length τ d = path traveled by diffusing solute domain unit length The Kozeny-Carman equation uses hydraulic tortuosity: k = φ3 cτ h 2 S 2 = φ3 βs 2 Past Work on measuring τ h and/or fitting simulated data to KC-eqn: overlapping squares (Koponen et al 1996/7, Matyka et al 2008, Duda et al 2011) 3D body-centered-cubic unit cells (Ebrahimi Khabbazi et al 2013) Effective diffusivity is the binary diffusion scaled by diffusive tortuosity D eff = D τ d = Dτ d Past Work on measuring τ d or D eff : in cellular geometry (Ochoa et al. 1987) in porous wick (Beyhaghi & Pillai 2011) in packed beds and unconsolidated porous media (Kim et al. 1987, Quintard 1993) 9 / 52

10 Different forms of tortuosity The tortuous nature of fluid flow through porous media: The tortuous nature of a solute diffusing through porous media: τ h = path traveled by fluid domain unit length τ d = path traveled by diffusing solute domain unit length The Kozeny-Carman equation uses hydraulic tortuosity: k = φ3 cτ h 2 S 2 = φ3 βs 2 Effective diffusivity is the binary diffusion scaled by diffusive tortuosity D eff = D τ d = Dτ d τ h τ d 10 / 52

11 Work flow of this research represent porous media in domain, using simple geometries do pore-scale modeling: solve Stokes flow obtain pressure and velocity fields compute permeability tensor and hydraulic tortuosity solve Closure Variable problem obtain diffusive tortuosity tensor compare the hydraulic and diffusive tortuosity results fit simulation results (permeability, hydraulic tortuosity) to Kozeny-Carman equation 11 / 52

12 Solve Stokes flow obtain velocities, permeability Wang et al Solve Stokes flow in pore space of media obtain p and v. Pressure and flow fields in 0.05mm x 0.05mm REV, µ 2 v p + ρg = 0 (2) v = 0 (3) 2. Assuming Darcy s equation is valid (u = K µ ( p + ρg)), compute K from u = φ v ( p = 0 in REV). (a) Flux in X-direction (b) Flux in Y-direction and the resulting K in m 2 is» kxx k xy = k yx k yy» 6.03x x x x / 52

13 Solve Stokes flow obtain velocities, hydraulic tortuosity The hydraulic tortuosity can be computed from fluid velocity fields, vx 2 + vy 2 vx 2 + vy 2 τ hx =, τ hy = v x v y (4) In other words, vx (i, j) 2 + v y (i, j) 2 vx (i, j) 2 + v y (i, j) 2 τ hx = i,j v x (i, j) i,j, τ hy = i,j v y (i, j) i,j * Koponen et al 1997, Duda et al / 52

14 Solve closure problem obtain diffusive tortuosity diffusive transport at pore-scale: c t = (D m c) (5) use theory of volume averaging to define c f local spatial deviation c = c c f c is a linear function of c f, so is given by c = b c f (6) b is a vector field that maps c f onto c closed form of local volume average transport eqn., where φ c f t = (φd eff c f ) (7) ( D eff := D m I + 1 ) n bda V f A fs (8) by convention, D eff /D m = 1/τ 14 / 52

15 Solve closure problem obtain diffusive tortuosity The closure problem is given by: 2 b = 0 (9) n fs b = n fs at A fs (10) b(r + I i ) = b(r) i = 1, 2 (11) b f = 0 (12) Figure: Closure Problem The tortuosity (factor) tensor is τ = I + 1 V f A fs n fs bda (13) * Valdes-Parada et al. 2011, Beyhaghi & Pillai 2011, Quintard 1993, Kim et al. 1987, Ochoa et al / 52

16 Solve closure problem obtain diffusive tortuosity In 2D, the tensor components are τ = " τ xx τ yx # " τ xy V τ = f RA n x b x da fs 1 yy V f 1 V f RA fs n x b y da R A fs n y b x da V f RA fs n y b y da # (14) The Laplace equation 2 b = 0 is used to solve the closure variable b, where b = b x i + b y j. In 2D, these spatial derivatives are 2 b x = x 2 b y = x bx x by x «+ y «+ y bx y by y «= 0 (15) «= 0 (16) In this work, we solve these PDEs using finite difference. * Valdes-Parada et al. 2011, Beyhaghi & Pillai 2011, Quintard 1993, Kim et al. 1987, Ochoa et al / 52

17 Example 1: In-line array of uniform shapes (a) φ = 0.32 (b) φ = 0.54 (c) φ = 0.71 (d) φ = 0.87 (e) φ = / 52

18 In-line array of circles; φ = 0.71 Hydraulic tortuosity: Diffusive tortuosity: (a) Flow in x-dir. (b) Flow in y-dir. Figure: Magnitude of pressure gradients (i.e., velocities) with hydraulic flow lines 18 / 52

19 In-line array of circles; φ = 0.71 Hydraulic tortuosity: Diffusive tortuosity: (a) Flow in x-dir. (b) Flow in y-dir. Figure: Magnitude of pressure gradients (i.e., velocities) with hydraulic flow lines τ hx = vx 2 + vy 2 / v x = 1.02 τ hy = vx 2 + vy 2 / v y = / 52

20 In-line array of circles; φ = 0.71 Hydraulic tortuosity: Diffusive tortuosity: (a) Flow in x-dir. (b) Flow in y-dir. (a) b x (b) b y Figure: Magnitude of pressure gradients (i.e., velocities) with hydraulic flow lines Figure: Closure variable fields used to obtain surface intergrals τ hx = vx 2 + vy 2 / v x = 1.02 τ hy = vx 2 + vy 2 / v y = / 52

21 In-line array of circles; φ = 0.71 Hydraulic tortuosity: Diffusive tortuosity: (a) Flow in x-dir. (b) Flow in y-dir. (a) b x (b) b y Figure: Magnitude of pressure gradients (i.e., velocities) with hydraulic flow lines Figure: Closure variable fields used to obtain surface intergrals τ hx = vx 2 + vy 2 / v x = 1.02 τ hy = vx 2 + vy 2 / v y = 1.02 ) 1 τ d x = (1 + 1Vf n x b x da = 1.29 A fs ) 1 τ d y = (1 + 1Vf n y b y da = 1.29 A fs 21 / 52

22 In-line array of circles; φ = 0.71 Hydraulic tortuosity: Diffusive tortuosity: (a) Flow in x-dir. (b) Flow in y-dir. q (a) c x x 2 + c x y 2 (b) q c y x 2 c + y 2 y Figure: Magnitude of pressure gradients (i.e., velocities) with hydraulic flow lines Figure: Magnitude of concentration gradients with diffusive flow lines τ hx = τ hy = τ d x = τ d y = / 52

23 In-line array of circles and squares; 0 < φ < 1 (a) Circle (b) Square Figure: Porosity vs. tortuosity: diffusive tortuosity trend follows analytical solution from Rayleigh 1892 until inadequate mesh refinement, while hydraulic tortuosity is independent of porosity for this geometry. 23 / 52

24 Example 2: Staggered-array of uniform shapes (a) φ = 0.36 (b) φ = 0.64 (c) φ = 0.84 (d) φ = / 52

25 Staggered-array of squares; φ = Hydraulic tortuosity: Diffusive tortuosity: (a) Flow in x-dir. (b) Flow in y-dir. Figure: Magnitude of pressure gradients (i.e., velocities) with hydraulic flow lines τ hx = 1.37 τ hy = 1.06 (a) X-Prob (b) Y-Prob Figure: Magnitude of concentration gradients with diffusive flow lines τ d xx = 1.37 τ d yy = / 52

26 Staggered-array of squares; 0.35 < φ < 1 (a) τ xx (b) τ yy Figure: Tortuosity vs. porosity trends for staggered-array of squares: hydraulic flow is predominantly parallel to driving force, unless an obstacle prevents a linear flow path. 26 / 52

27 Example 3: Randomly distributed, freely overlapping squares (a) φ = 0.53 (b) φ = 0.61 (c) φ = 0.70 (d) φ = 0.80 (e) φ = 0.90 Figure: Various pore-structure geometries. Notice isolated fluid sites are filled in, causing larger, non-uniform shapes. Square length = 0.01 x domain length (Koza et al 2000, Duda et al 2011). 27 / 52

28 Overlapping Squares; φ = 0.7 Hydraulic tortuosity: Diffusive tortuosity: q (a) p x x 2 + p x 2 y (b) Hydr. lines q (a) c x x 2 + c x 2 y (b) Diff. lines Figure: Magnitude of pressure gradients (i.e., velocities) and hydraulic flow lines for X-problem Figure: Magnitude of concentration gradients and diffusive flow lines for X-problem τ hx = τ d x = / 52

29 Overlapping Squares; 0.45 < φ < 1 (a) Plot of full simulated data (b) Close-up plot of simulated data Figure: Porosity vs tortuosity for pore-structures of overlapping squares: diffusive tortuosity >> than hydraulic tortuosity at low porosities 29 / 52

30 Overlapping Squares: Anisotropy of τ h and τ d (a) φ = (b) φ = Figure: Anisotropic ratio for hydraulic and diffusive tortuosity: the degree of anisotropy for hydraulic flow and diffusive flow 1 for higher porosities, while the diffusive flow is anisotropic at lower porosities. Figure: Pore-structures of lower porosity exhibit more isolated fluid site that are filled it, thus creating non-uniform distribution of solids in overlapping squares configuration. This impacts degree of anisotropy. 30 / 52

31 Computing Kozeny constant for Ex. 3: overlapping squares 31 / 52

32 Overlapping Squares: Kozeny constant (a) Permeability k (b) Hydraulic tortuosity τ (c) Spec. surface S Figure: Simulated data The Kozeny constant c, or shape factor β, can be computed by where S = A f s V tot k = φ3 cs 2 = φ3 βτ 2 S 2 (17) (i.e., V tot is the fluid and solid space combined). 32 / 52

33 Overlapping Squares: Kozeny constant (a) Kozeny constant, c (b) Shape factor, β Figure: Parameters that fit simulated data to Kozeny-Carman equation. Kozeny constant and shape factor do not change within a porosity interval of 0.65 < φ < / 52

34 Conclusions Past work has proposed τ = τ(φ), however recent work has suggested this is not true, rather τ = τ(pore geometry) only 34 / 52

35 Conclusions Past work has proposed τ = τ(φ), however recent work has suggested this is not true, rather τ = τ(pore geometry) only Hydraulic tortuosity can be computed from fluid velocity field, while diffusive tortuosity can be computed from closure variable problem 35 / 52

36 Conclusions Past work has proposed τ = τ(φ), however recent work has suggested this is not true, rather τ = τ(pore geometry) only Hydraulic tortuosity can be computed from fluid velocity field, while diffusive tortuosity can be computed from closure variable problem Hydraulic tortuosity diffusive tortuosity 36 / 52

37 Conclusions Past work has proposed τ = τ(φ), however recent work has suggested this is not true, rather τ = τ(pore geometry) only Hydraulic tortuosity can be computed from fluid velocity field, while diffusive tortuosity can be computed from closure variable problem Hydraulic tortuosity diffusive tortuosity The diffusive tortuosity is able to capture the anisotropic nature of pore geometry more readily that hydraulic tortuosity 37 / 52

38 Conclusions Past work has proposed τ = τ(φ), however recent work has suggested this is not true, rather τ = τ(pore geometry) only Hydraulic tortuosity can be computed from fluid velocity field, while diffusive tortuosity can be computed from closure variable problem Hydraulic tortuosity diffusive tortuosity The diffusive tortuosity is able to capture the anisotropic nature of pore geometry more readily that hydraulic tortuosity The Kozeny-Carman eqn. proposes a relationship between permeability, porosity, specific surface area, and hydraulic tortuosity 38 / 52

39 Questions? 39 / 52

40 Appendix 40 / 52

41 Obtaining full permeability tensor The compute the full permeability tensor, consider the pore structure: Figure: Pore structure example, φ = 0.82 We assume this pore structure is a representative elementary volume (REV), which implies the properties measured from this sample represents the properties of the porous media the sample came from. The rest of the domain is made up of repeating structures of this figure. As such, we use periodic boundary conditions on all external sides. 41 / 52

42 Slow and steady flow is modeled through the pore space by Stokes equation, 0 = ρg i + p x i + µ 2 v i x j x j (18) where the solution for pressure p and velocity v i are obtained for a specified flow driving force, g i. The spatial average pore-space velocity ˆv i can be computed for each flow scenario. (a) Driving force in X- (b) Driving force in Y- direction: (g x, g y ) = direction: (g x, g y ) = (1, 0). (ˆv x, ˆv y ) = (?,?) (0, 1). (ˆv x, ˆv y ) = (?,?) Figure: Pressure and velocity solutions after solving Stokes flow through 42 / 52

43 Then, we use Darcy s equation to describe flow at the macro-scale, u i = k ( ) ij P + ρg i (19) µ x i where u i is the Darcy s velocity of the pore structure (u i = ˆv i /φ), and where k ij and P/ x i are with respect to the whole structure (i.e., the REV). In 2D, Darcy s equation is [ ux u y ] = 1 [ kxx k xy µ k yx k yy ] ([ P x P y ] + ρ [ gx and since we assume periodic boundaries on all external sides of the pore structure, the pressure gradients are P x REV = 0, thus equation 20 becomes [ ] ux = ρ µ u y g y ]) (20) P y REV = 0 (21) [ ] [ ] kxx k xy gx k yx k yy g y (22) 43 / 52

44 When (g x, g y ) = (1, 0), equation 22 is ] [ˆvx /φ = ˆv y /φ [ kxx k yx ] (23) and when (g x, g y ) = (0, 1), equation 22 is ] [ˆvx /φ = ˆv y /φ [ kxy k yy ] (24) where ρ/µ = 1, and u i = ˆv i /φ. In our pore structure example, the resulting permeability tensor is [ ] [ ] kxx k xy = (25) k yx k yy To conclude, the full permeability tensor can be computed by using the spatial average pore-space velocity ˆv i obtained by solving Stokes flow for two different driving force scenarios. 44 / 52

45 Solve closure problem obtain diffusive tortuosity diffusive transport at pore-scale: c t = (D m c) (26) use theory of volume averaging to define c f local spatial deviation c = c c f transport of local volume average, c f : φ c f t )) = (φd m ( c f + 1Vf n cda A fs c is a linear function of c f, so is given by b is a vector field that maps c f onto c (27) c = b c f (28) * Valdes-Parada et al. 2011, Beyhaghi & Pillai 2011, Quintard 1993, Kim et al. 1987, Ochoa et al / 52

46 take the surface integral of c = b c f substitute it into the local volume average transport eqn. to get its closed form, where φ c f t = (φd eff c f ) (29) ( D eff := D m I + 1 ) n bda V f A fs (30) by convention, D eff /D m = 1/τ, thus the tortuosity (or tortuosity factor) is computed by ) 1 τ = τ = (I + 1Vf n bda A fs (31) 46 / 52

47 In-line array of circles; φ = 0.71 Hydraulic tortuosity: Diffusive tortuosity: (a) Flow in x-dir. (b) Flow in y-dir. Figure: Magnitude of pressure gradients (i.e., velocities) with hydraulic flow lines (a) b x/ x (b) b y / y Figure: Closure variable gradients used to obtain volume integrals τ hx = vx 2 + vy 2 / v x = 1.02 τ hy = vx 2 + vy 2 / v y = 1.02 τ d x = (1 + 1Vf τ d y = (1 + 1Vf V f V f ) 1 b x x dv b y y dv ) 1 47 / 52

48 In-line array of circles; φ = 0.71 Hydraulic tortuosity: Diffusive tortuosity: (a) Flow in x-dir. (b) Flow in y-dir. (a) c x = b x + xs (b) c y = b y + ys Figure: Magnitude of pressure gradients (i.e., velocities) with hydraulic flow lines Figure: Concentration fields τ hx = vx 2 + vy 2 / v x = 1.02 τ hy = vx 2 + vy 2 / v y = 1.02 ( 1 τ d x = V f ( 1 τ d y = V f V f V f ) 1 c x x dv c y y dv ) 1 48 / 52

49 Figure: Diffusive tortuosity, τ yy 49 / 52 Staggered-array of squares; φ = q c (a) b x (b) c x = b x + xs (c) x x diffusive flow lines 2 + c x 2 y with Figure: Diffusive tortuosity, τ xx q c (a) b y (b) c y = b y + ys (c) y x diffusive flow lines 2 c + y 2 y with

50 Overlapping Squares; φ = 0.7 q (a) b x (b) c x = b x + xs (c) c x x 2 + c x 2 y (d) Diffu. lines Figure: Scalar fields (closure variable b x, concentration c x ), magnitude of concentration gradients, and diffusive flow lines for X-problem [ τ xx τ yx τ xy τ yy ] = I + 1Vf τ xx = 1 τ d x = V f V f V f bdv f = 1 V f b x x dv f = 1 V f V f V f cdv f c x x dv f 50 / 52

51 Overlapping Squares; φ = 0.7 q (a) b x (b) c x = b x + xs (c) c x x 2 + c x 2 y (d) Diffu. lines Figure: Scalar fields (closure variable b x, concentration c x ), magnitude of concentration gradients, and diffusive flow lines for X-problem V f i,j 1 V f i,j b x x δv f 1 b y V f x δv f i,j b x y δv f V f i,j b y y δv f = 1 V f 1 V f i,j i,j c x x δv f c x y δv f 1 c y V f x δv f i,j 1 c y V f y δv f i,j 51 / 52

52 Overlapping Squares: Anisotropy of k and τ h Figure: Comparison between the degree of anisotropy of permeability and hydraulic tortuosity: permeability is more anisotropic than hydraulic tortuosity at lower porosities 52 / 52

THE ANALYSIS OF THE RELATIONS BETWEEN POROSITY AND TORTUOSITY IN GRANULAR BEDS

THE ANALYSIS OF THE RELATIONS BETWEEN POROSITY AND TORTUOSITY IN GRANULAR BEDS Technical Sciences, 2017, 20(1), 75 85 THE ANALYSIS OF THE RELATIONS BETWEEN POROSITY AND TORTUOSITY IN GRANULAR BEDS Department of Mechanics and Basics of Machine Construction University of Warmia and

More information

Two-dimensional flow in a porous medium with general anisotropy

Two-dimensional flow in a porous medium with general anisotropy Two-dimensional flow in a porous medium with general anisotropy P.A. Tyvand & A.R.F. Storhaug Norwegian University of Life Sciences 143 Ås Norway peder.tyvand@umb.no 1 Darcy s law for flow in an isotropic

More information

Hydraulic properties of porous media

Hydraulic properties of porous media PART 5 Hydraulic properties of porous media Porosity Definition: Void space: n V void /V total total porosity e V void /V solid Primary porosity - between grains Secondary porosity - fracture or solution

More information

1 Modeling Immiscible Fluid Flow in Porous Media

1 Modeling Immiscible Fluid Flow in Porous Media Excerpts from the Habilitation Thesis of Peter Bastian. For references as well as the full text, see http://cox.iwr.uni-heidelberg.de/people/peter/pdf/bastian_habilitationthesis.pdf. Used with permission.

More information

VI. Porous Media. Lecture 34: Transport in Porous Media

VI. Porous Media. Lecture 34: Transport in Porous Media VI. Porous Media Lecture 34: Transport in Porous Media 4/29/20 (corrected 5/4/2 MZB) Notes by MIT Student. Conduction In the previous lecture, we considered conduction of electricity (or heat conduction

More information

How are calculus, alchemy and forging coins related?

How are calculus, alchemy and forging coins related? BMOLE 452-689 Transport Chapter 8. Transport in Porous Media Text Book: Transport Phenomena in Biological Systems Authors: Truskey, Yuan, Katz Focus on what is presented in class and problems Dr. Corey

More information

Transfer Equations: An Attempt to Pose an Optimization Problem. Project for CE291 Henry Kagey

Transfer Equations: An Attempt to Pose an Optimization Problem. Project for CE291 Henry Kagey Transfer Equations: An Attempt to Pose an Optimization Problem Project for CE291 Henry Kagey Background System Solar Disinfection of Greywater The goal of this study is to define the mass transfer in a

More information

KOZENY-CARMAN EQUATION REVISITED. Jack Dvorkin Abstract

KOZENY-CARMAN EQUATION REVISITED. Jack Dvorkin Abstract KOZENY-CARMAN EQUATION REVISITED Jack Dvorkin -- 009 Abstract The Kozeny-Carman equation is often presented as permeability versus porosity, grain size, and tortuosity. When it is used to estimate permeability

More information

An Immersed Boundary Method for Computing Anisotropic Permeability of Structured Porous Media

An Immersed Boundary Method for Computing Anisotropic Permeability of Structured Porous Media An Immersed Boundary Method for Computing Anisotropic Permeability of Structured Porous Media D.J. Lopez Penha a, B.J. Geurts a, S. Stolz a,b, M. Nordlund b a Dept. of Applied Mathematics, University of

More information

q v = - K h = kg/ν units of velocity Darcy's Law: K = kρg/µ HYDRAULIC CONDUCTIVITY, K Proportionality constant in Darcy's Law

q v = - K h = kg/ν units of velocity Darcy's Law: K = kρg/µ HYDRAULIC CONDUCTIVITY, K Proportionality constant in Darcy's Law Darcy's Law: q v - K h HYDRAULIC CONDUCTIVITY, K m/s K kρg/µ kg/ν units of velocity Proportionality constant in Darcy's Law Property of both fluid and medium see D&S, p. 62 HYDRAULIC POTENTIAL (Φ): Φ g

More information

Permeability of Dual-Structured Porous Media

Permeability of Dual-Structured Porous Media 56 The Open Transport Phenomena Journal, 2011, 3, 56-61 Permeability of Dual-Structured Porous Media Open Access Ehsan Khajeh * and Daan M. Maijer Department of Materials Engineering, The University of

More information

Fluid Particle Interactions Basics

Fluid Particle Interactions Basics Fluid Particle Interactions Basics Outline u 1. Single sphere: Stokes and turbulent drag 2. Many spheres: Stokes flow 3. Many spheres: intermediate Re 4. Many spheres: terminal velocity? 1. Single sphere:

More information

Mobility of Power-law and Carreau Fluids through Fibrous Media

Mobility of Power-law and Carreau Fluids through Fibrous Media Mobility of Power-law and Carreau Fluids through Fibrous Media Setareh Shahsavari, Gareth H. McKinley Department of Mechanical Engineering, Massachusetts Institute of Technology September 3, 05 Abstract

More information

Getting started: CFD notation

Getting started: CFD notation PDE of p-th order Getting started: CFD notation f ( u,x, t, u x 1,..., u x n, u, 2 u x 1 x 2,..., p u p ) = 0 scalar unknowns u = u(x, t), x R n, t R, n = 1,2,3 vector unknowns v = v(x, t), v R m, m =

More information

Number of pages in the question paper : 05 Number of questions in the question paper : 48 Modeling Transport Phenomena of Micro-particles Note: Follow the notations used in the lectures. Symbols have their

More information

arxiv: v1 [physics.flu-dyn] 26 Mar 2012

arxiv: v1 [physics.flu-dyn] 26 Mar 2012 How to Calculate ortuosity Easily? Maciej Matyka and Zbigniew Koza arxiv:1203.5646v1 [physics.flu-dyn] 26 Mar 2012 Institute of heoretical Physics, University of Wrocław, pl. M. Borna 9, 50-204 Wrocław,

More information

Permeability and fluid transport

Permeability and fluid transport Permeability and fluid transport Thermal transport: Fluid transport: q = " k # $p with specific discharge (filter velocity) q [m s 1 ] pressure gradient p [N m 3 ] dynamic viscosity η [N s m 2 ] (intrinsic)

More information

Studies on flow through and around a porous permeable sphere: II. Heat Transfer

Studies on flow through and around a porous permeable sphere: II. Heat Transfer Studies on flow through and around a porous permeable sphere: II. Heat Transfer A. K. Jain and S. Basu 1 Department of Chemical Engineering Indian Institute of Technology Delhi New Delhi 110016, India

More information

Random packings of spiky particles: Geometry and transport properties

Random packings of spiky particles: Geometry and transport properties Random packings of spiky particles: Geometry and transport properties I. Malinouskaya a,b, V.V. Mourzenko b, J.-F. Thovert b, P.M. Adler a a. Institut P, Département Fluides, Thermique, Combustion, ENSMA

More information

ON THE VALIDITY OF THE CARMAN-KOZENY EQUATION IN RANDOM FIBROUS MEDIA

ON THE VALIDITY OF THE CARMAN-KOZENY EQUATION IN RANDOM FIBROUS MEDIA II International Conference on Particle-based Methods Fundamentals and Applications PARTICLES 011 E. Oñate and D.R.J. Owen (Eds) ON THE VALIDITY OF THE CARMAN-KOZENY EQUATION IN RANDOM FIBROUS MEDIA K.

More information

Introduction to Aspects of Multiscale Modeling as Applied to Porous Media

Introduction to Aspects of Multiscale Modeling as Applied to Porous Media Introduction to Aspects of Multiscale Modeling as Applied to Porous Media Part III Todd Arbogast Department of Mathematics and Center for Subsurface Modeling, Institute for Computational Engineering and

More information

Number of pages in the question paper : 06 Number of questions in the question paper : 48 Modeling Transport Phenomena of Micro-particles Note: Follow the notations used in the lectures. Symbols have their

More information

NUMERICAL MODELING OF FLOW THROUGH DOMAINS WITH SIMPLE VEGETATION-LIKE OBSTACLES

NUMERICAL MODELING OF FLOW THROUGH DOMAINS WITH SIMPLE VEGETATION-LIKE OBSTACLES XIX International Conference on Water Resources CMWR 2012 University of Illinois at Urbana-Champaign June 17-22,2012 NUMERICAL MODELING OF FLOW THROUGH DOMAINS WITH SIMPLE VEGETATION-LIKE OBSTACLES Steven

More information

Fluid Flow Fluid Flow and Permeability

Fluid Flow Fluid Flow and Permeability and Permeability 215 Viscosity describes the shear stresses that develop in a flowing fluid. V z Stationary Fluid Velocity Profile x Shear stress in the fluid is proportional to the fluid velocity gradient.

More information

In all of the following equations, is the coefficient of permeability in the x direction, and is the hydraulic head.

In all of the following equations, is the coefficient of permeability in the x direction, and is the hydraulic head. Groundwater Seepage 1 Groundwater Seepage Simplified Steady State Fluid Flow The finite element method can be used to model both steady state and transient groundwater flow, and it has been used to incorporate

More information

Three-Dimensional Simulation of Mixing Flow in a Porous Medium with Heat and Mass Transfer in a Moisture Recovery System

Three-Dimensional Simulation of Mixing Flow in a Porous Medium with Heat and Mass Transfer in a Moisture Recovery System 12 th Fluid Dynamics Conference, Babol Noshirvani University of Technology, 28-30 April 2009 Three-Dimensional Simulation of Mixing Flow in a Porous Medium with Heat and Mass Transfer in a Moisture Recovery

More information

Evaporation-driven transport and precipitation of salt in porous media: A multi-domain approach

Evaporation-driven transport and precipitation of salt in porous media: A multi-domain approach Evaporation-driven transport and precipitation of salt in porous media: A multi-domain approach Vishal Jambhekar Karen Schmid, Rainer Helmig Department of Hydromechanics and Hydrosystemmodeling EGU General

More information

Chapter 1. Viscosity and the stress (momentum flux) tensor

Chapter 1. Viscosity and the stress (momentum flux) tensor Chapter 1. Viscosity and the stress (momentum flux) tensor Viscosity and the Mechanisms of Momentum Transport 1.1 Newton s law of viscosity ( molecular momentum transport) 1.2 Generalization of Newton

More information

Geometry modeling of open-cell foams for efficient fluid flow and heat transfer computations using modified Kelvin cells

Geometry modeling of open-cell foams for efficient fluid flow and heat transfer computations using modified Kelvin cells E. Werzner, M.A.A. Mendes, J. Storm, M. Laurinat, S. Ray, D. Trimis Chair of Gas and Heat Technology Institute of Thermal Engineering Gustav-Zeuner-Str. 7 09596 Freiberg Germany Geometry modeling of open-cell

More information

Mathematical Tripos Part IA Lent Term Example Sheet 1. Calculate its tangent vector dr/du at each point and hence find its total length.

Mathematical Tripos Part IA Lent Term Example Sheet 1. Calculate its tangent vector dr/du at each point and hence find its total length. Mathematical Tripos Part IA Lent Term 205 ector Calculus Prof B C Allanach Example Sheet Sketch the curve in the plane given parametrically by r(u) = ( x(u), y(u) ) = ( a cos 3 u, a sin 3 u ) with 0 u

More information

Tortuosity-porosity relation in porous media flow

Tortuosity-porosity relation in porous media flow Tortuosity-porosity relation in porous media flow Maciej Matyka, 1 Arzhang Khalili, 2,3 and Zbigniew Koza 1 1 Institute of Theoretical Physics, University of Wrocław, pl. M. Borna 9, 50-204 Wrocław, Poland

More information

THE ROCK PHYSICS HANDBOOK

THE ROCK PHYSICS HANDBOOK THE ROCK PHYSICS HANDBOOK TOOLS FOR SEISMIC ANALYSIS IN POROUS MEDIA Gary Mavko Tapan Mukerji Jack Dvorkin Stanford University Stanford University Stanford University CAMBRIDGE UNIVERSITY PRESS CONTENTS

More information

arxiv: v1 [physics.flu-dyn] 21 Jan 2009

arxiv: v1 [physics.flu-dyn] 21 Jan 2009 Finite-size anisotropy in statistically uniform porous media Zbigniew Koza and Maciej Matyka Institute of Theoretical Physics,University of Wroc law, pl. M. Borna 9, 5-4 Wroc law, Poland Arzhang Khalili

More information

Anisotropic permeabilities evolution of reservoir rocks under pressure:

Anisotropic permeabilities evolution of reservoir rocks under pressure: Extended reserves Clean refining Fuel-efficient vehicles Diversified fuels Controlled CO 2 Anisotropic permeabilities evolution : New experimental and numerical approaches (1) Dautriat J. 1-2*, Gland N.

More information

PORE-SCALE PHASE FIELD MODEL OF TWO-PHASE FLOW IN POROUS MEDIUM

PORE-SCALE PHASE FIELD MODEL OF TWO-PHASE FLOW IN POROUS MEDIUM Excerpt from the Proceedings of the COMSOL Conference 2010 Paris PORE-SCALE PHASE FIELD MODEL OF TWO-PHASE FLOW IN POROUS MEDIUM Igor Bogdanov 1*, Sylvain Jardel 1, Anis Turki 1, Arjan Kamp 1 1 Open &

More information

Effects of Viscous Dissipation on Unsteady Free Convection in a Fluid past a Vertical Plate Immersed in a Porous Medium

Effects of Viscous Dissipation on Unsteady Free Convection in a Fluid past a Vertical Plate Immersed in a Porous Medium Transport in Porous Media (2006) 64: 1 14 Springer 2006 DOI 10.1007/s11242-005-1126-6 Effects of Viscous Dissipation on Unsteady Free Convection in a Fluid past a Vertical Plate Immersed in a Porous Medium

More information

Available online at ScienceDirect. Energy Procedia 78 (2015 ) th International Building Physics Conference, IBPC 2015

Available online at   ScienceDirect. Energy Procedia 78 (2015 ) th International Building Physics Conference, IBPC 2015 Available online at www.sciencedirect.com ScienceDirect Energy Procedia 78 (2015 ) 537 542 6th International Building Physics Conference, IBPC 2015 Characterization of fibrous insulating materials in their

More information

Darcy s law in 3-D. K * xx K * yy K * zz

Darcy s law in 3-D. K * xx K * yy K * zz PART 7 Equations of flow Darcy s law in 3-D Specific discarge (vector) is calculated by multiplying te ydraulic conductivity (second-order tensor) by te ydraulic gradient (vector). We obtain a general

More information

A Statistical Correlation Between Permeability, Porosity, Tortuosity and Conductance

A Statistical Correlation Between Permeability, Porosity, Tortuosity and Conductance Transp Porous Med https://doi.org/10.1007/s11242-017-0983-0 A Statistical Correlation Between Permeability, Porosity, Tortuosity and Conductance S. M. Rezaei Niya 1 A. P. S. Selvadurai 1 Received: 7 March

More information

Experimental Study of Darcy and Non-Darcy Flow in Porous Media

Experimental Study of Darcy and Non-Darcy Flow in Porous Media International Journal of Engineering and Technology Volume No., December, 0 Experimental Study of Darcy and Non-Darcy Flow in Porous Media Oyinkepreye D. Orodu, Favour A. Makinde, Kale B. Orodu Department

More information

Advanced Heat and Mass Transfer by Amir Faghri, Yuwen Zhang, and John R. Howell

Advanced Heat and Mass Transfer by Amir Faghri, Yuwen Zhang, and John R. Howell Heat Transfer Heat transfer rate by conduction is related to the temperature gradient by Fourier s law. For the one-dimensional heat transfer problem in Fig. 1.8, in which temperature varies in the y-

More information

Streamline calculations. Lecture note 2

Streamline calculations. Lecture note 2 Streamline calculations. Lecture note 2 February 26, 2007 1 Recapitulation from previous lecture Definition of a streamline x(τ) = s(τ), dx(τ) dτ = v(x,t), x(0) = x 0 (1) Divergence free, irrotational

More information

UNCLASSI FILED L10160 DEFENSE DOCUMENTATION CENTER FOR SCIENTIFIC AND TECHNICAL INFORMATION CAMERON STAMTIOI, ALEXANDRIA, VIRGINIA UNCLASSIFIED

UNCLASSI FILED L10160 DEFENSE DOCUMENTATION CENTER FOR SCIENTIFIC AND TECHNICAL INFORMATION CAMERON STAMTIOI, ALEXANDRIA, VIRGINIA UNCLASSIFIED UNCLASSI FILED AD L10160 DEFENSE DOCUMENTATION CENTER FOR SCIENTIFIC AND TECHNICAL INFORMATION CAMERON STAMTIOI, ALEXANDRIA, VIRGINIA UNCLASSIFIED NOTICE: When government or other drawings, specifications

More information

ENGI Partial Differentiation Page y f x

ENGI Partial Differentiation Page y f x ENGI 344 4 Partial Differentiation Page 4-0 4. Partial Differentiation For functions of one variable, be found unambiguously by differentiation: y f x, the rate of change of the dependent variable can

More information

Modeling as a tool for understanding the MEA. Henrik Ekström Utö Summer School, June 22 nd 2010

Modeling as a tool for understanding the MEA. Henrik Ekström Utö Summer School, June 22 nd 2010 Modeling as a tool for understanding the MEA Henrik Ekström Utö Summer School, June 22 nd 2010 COMSOL Multiphysics and Electrochemistry Modeling The software is based on the finite element method A number

More information

M98-P2 (formerly C98-P1) Non-Newtonian Fluid Flow through Fabrics Matthew W. Dunn Philadelphia University

M98-P2 (formerly C98-P1) Non-Newtonian Fluid Flow through Fabrics Matthew W. Dunn Philadelphia University 1 Non-Newtonian Fluid Flow through Fabrics Matthew W. Dunn Philadelphia University http://spike.philacol.edu/perm/ Goal Statement The overall objectives of this program are to Model fabric porosity based

More information

Darcy's Law. Laboratory 2 HWR 531/431

Darcy's Law. Laboratory 2 HWR 531/431 Darcy's Law Laboratory HWR 531/431-1 Introduction In 1856, Henry Darcy, a French hydraulic engineer, published a report in which he described a series of experiments he had performed in an attempt to quantify

More information

MATH COURSE NOTES - CLASS MEETING # Introduction to PDEs, Fall 2011 Professor: Jared Speck

MATH COURSE NOTES - CLASS MEETING # Introduction to PDEs, Fall 2011 Professor: Jared Speck MATH 8.52 COURSE NOTES - CLASS MEETING # 6 8.52 Introduction to PDEs, Fall 20 Professor: Jared Speck Class Meeting # 6: Laplace s and Poisson s Equations We will now study the Laplace and Poisson equations

More information

Effect of fibre shape on transverse thermal conductivity of unidirectional composites

Effect of fibre shape on transverse thermal conductivity of unidirectional composites Sādhanā Vol. 4, Part 2, April 25, pp. 53 53. c Indian Academy of Sciences Effect of fibre shape on transverse thermal conductivity of unidirectional composites B RAGHAVA RAO,, V RAMACHANDRA RAJU 2 and

More information

MATH COURSE NOTES - CLASS MEETING # Introduction to PDEs, Spring 2018 Professor: Jared Speck

MATH COURSE NOTES - CLASS MEETING # Introduction to PDEs, Spring 2018 Professor: Jared Speck MATH 8.52 COURSE NOTES - CLASS MEETING # 6 8.52 Introduction to PDEs, Spring 208 Professor: Jared Speck Class Meeting # 6: Laplace s and Poisson s Equations We will now study the Laplace and Poisson equations

More information

This section develops numerically and analytically the geometric optimisation of

This section develops numerically and analytically the geometric optimisation of 7 CHAPTER 7: MATHEMATICAL OPTIMISATION OF LAMINAR-FORCED CONVECTION HEAT TRANSFER THROUGH A VASCULARISED SOLID WITH COOLING CHANNELS 5 7.1. INTRODUCTION This section develops numerically and analytically

More information

MODELLING MULTIPHASE FLOWS OF DISCRETE PARTICLES IN VISCOELASTIC FLUIDS

MODELLING MULTIPHASE FLOWS OF DISCRETE PARTICLES IN VISCOELASTIC FLUIDS MODELLING MULTIPHASE FLOWS OF DISCRETE PARTICLES IN VISCOELASTIC FLUIDS R. RIBEIRO 1, C. FERNANDES 1, S.A. FAROUGHI 2, G.H. McKINLEY 2 AND J.M. NÓBREGA 1 1 Institute for Polymers and Composites/i3N, University

More information

DISPERSION IN POROUS MEDIA FOR MULTICOMPONENT SYSTEMS. Introduction. Theory

DISPERSION IN POROUS MEDIA FOR MULTICOMPONENT SYSTEMS. Introduction. Theory DISPERSION IN POROUS MEDIA FOR MULTICOMPONENT SYSTEMS Quintard M., IMFT, Toulouse-France Bletzacker, L., IMFT, Toulouse-France Chenu D., IMFT, Toulouse-France Whitaker S., UCD, Davis Introduction In this

More information

Dissolution and precipitation during flow in porous media

Dissolution and precipitation during flow in porous media 1/25 Class project for GEOS 692: Transport processes and physical properties of rocks Dissolution and precipitation during flow in porous media Gry Andrup-Henriksen Fall 2006 1 2/25 Outline Introduction

More information

Modeling seismic wave propagation during fluid injection in a fractured network: Effects of pore fluid pressure on time-lapse seismic signatures

Modeling seismic wave propagation during fluid injection in a fractured network: Effects of pore fluid pressure on time-lapse seismic signatures Modeling seismic wave propagation during fluid injection in a fractured network: Effects of pore fluid pressure on time-lapse seismic signatures ENRU LIU, SERAFEIM VLASTOS, and XIANG-YANG LI, Edinburgh

More information

Diffusion and Adsorption in porous media. Ali Ahmadpour Chemical Eng. Dept. Ferdowsi University of Mashhad

Diffusion and Adsorption in porous media. Ali Ahmadpour Chemical Eng. Dept. Ferdowsi University of Mashhad Diffusion and Adsorption in porous media Ali Ahmadpour Chemical Eng. Dept. Ferdowsi University of Mashhad Contents Introduction Devices used to Measure Diffusion in Porous Solids Modes of transport in

More information

I. Borsi. EMS SCHOOL ON INDUSTRIAL MATHEMATICS Bedlewo, October 11 18, 2010

I. Borsi. EMS SCHOOL ON INDUSTRIAL MATHEMATICS Bedlewo, October 11 18, 2010 : an : an (Joint work with A. Fasano) Dipartimento di Matematica U. Dini, Università di Firenze (Italy) borsi@math.unifi.it http://web.math.unifi.it/users/borsi porous EMS SCHOOL ON INDUSTRIAL MATHEMATICS

More information

Table of Contents. Preface... xiii

Table of Contents. Preface... xiii Preface... xiii PART I. ELEMENTS IN FLUID MECHANICS... 1 Chapter 1. Local Equations of Fluid Mechanics... 3 1.1. Forces, stress tensor, and pressure... 4 1.2. Navier Stokes equations in Cartesian coordinates...

More information

Stokes Flow in a Slowly Varying Two-Dimensional Periodic Pore

Stokes Flow in a Slowly Varying Two-Dimensional Periodic Pore Transport in Porous Media 6: 89 98, 997. 89 c 997 Kluwer Academic Publishers. Printed in the Netherlands. Stokes Flow in a Slowly Varying Two-Dimensional Periodic Pore PETER K. KITANIDIS? and BRUCE B.

More information

Fluid flow over and through a regular bundle of rigid fibres

Fluid flow over and through a regular bundle of rigid fibres Under consideration for publication in J. Fluid Mech. 1 Fluid flow over and through a regular bundle of rigid fibres Giuseppe A. Zampogna and Alessandro Bottaro DICCA, Scuola Politecnica, Università di

More information

The velocity distribution in a random porous medium

The velocity distribution in a random porous medium Engineering Conferences International ECI Digital Archives 5th International Conference on Porous Media and Their Applications in Science, Engineering and Industry Refereed Proceedings Summer 6-24-2014

More information

water L v i Chapter 4 Saturation

water L v i Chapter 4 Saturation 4. Resistivity The presence of hydrocarbons is identified by the electrical resistance of the formation. These electrical properties of rocks depend on the pore geometry and fluid distribution. That is,

More information

PRACTICE PROBLEMS. Please let me know if you find any mistakes in the text so that i can fix them. 1. Mixed partial derivatives.

PRACTICE PROBLEMS. Please let me know if you find any mistakes in the text so that i can fix them. 1. Mixed partial derivatives. PRACTICE PROBLEMS Please let me know if you find any mistakes in the text so that i can fix them. 1.1. Let Show that f is C 1 and yet How is that possible? 1. Mixed partial derivatives f(x, y) = {xy x

More information

(a) The points (3, 1, 2) and ( 1, 3, 4) are the endpoints of a diameter of a sphere.

(a) The points (3, 1, 2) and ( 1, 3, 4) are the endpoints of a diameter of a sphere. MATH 4 FINAL EXAM REVIEW QUESTIONS Problem. a) The points,, ) and,, 4) are the endpoints of a diameter of a sphere. i) Determine the center and radius of the sphere. ii) Find an equation for the sphere.

More information

Introduction of Partial Differential Equations and Boundary Value Problems

Introduction of Partial Differential Equations and Boundary Value Problems Introduction of Partial Differential Equations and Boundary Value Problems 2009 Outline Definition Classification Where PDEs come from? Well-posed problem, solutions Initial Conditions and Boundary Conditions

More information

SCA : A STRUCTURAL MODEL TO PREDICT TRANSPORT PROPERTIES OF GRANULAR POROUS MEDIA Guy Chauveteau, IFP, Yuchun Kuang IFP and Marc Fleury, IFP

SCA : A STRUCTURAL MODEL TO PREDICT TRANSPORT PROPERTIES OF GRANULAR POROUS MEDIA Guy Chauveteau, IFP, Yuchun Kuang IFP and Marc Fleury, IFP SCA2003-53: A STRUCTURAL MODEL TO PREDICT TRANSPORT PROPERTIES OF GRANULAR POROUS MEDIA Guy Chauveteau, IFP, Yuchun Kuang IFP and Marc Fleury, IFP This paper was prepared for presentation at the International

More information

2. Modeling of shrinkage during first drying period

2. Modeling of shrinkage during first drying period 2. Modeling of shrinkage during first drying period In this chapter we propose and develop a mathematical model of to describe nonuniform shrinkage of porous medium during drying starting with several

More information

On the origin of Darcy s law 1

On the origin of Darcy s law 1 Chapter 1 On the origin of Darcy s law 1 Cyprien Soulaine csoulain@stanford.edu When one thinks about porous media, the immediate concepts that come to mind are porosity, permeability and Darcy s law.

More information

Flow Modeling & Transport in Fiber Mats used as Reinforcement in Polymer Composites

Flow Modeling & Transport in Fiber Mats used as Reinforcement in Polymer Composites Flow Modeling & Transport in Fiber Mats used as Reinforcement in Polymer Composites Krishna M. Pillai Associate Professor Department of Mechanical Engineering A few facts about Milwaukee and Wisconsin:

More information

A Study on Numerical Solution to the Incompressible Navier-Stokes Equation

A Study on Numerical Solution to the Incompressible Navier-Stokes Equation A Study on Numerical Solution to the Incompressible Navier-Stokes Equation Zipeng Zhao May 2014 1 Introduction 1.1 Motivation One of the most important applications of finite differences lies in the field

More information

1.72, Groundwater Hydrology Prof. Charles Harvey Lecture Packet #4: Continuity and Flow Nets

1.72, Groundwater Hydrology Prof. Charles Harvey Lecture Packet #4: Continuity and Flow Nets 1.7, Groundwater Hydrology Prof. Charles Harvey Lecture Packet #4: Continuity and Flow Nets Equation of Continuity Our equations of hydrogeology are a combination of o Conservation of mass o Some empirical

More information

Non-Newtonian Flows in Porous Media: upscaling problems

Non-Newtonian Flows in Porous Media: upscaling problems 6/26/8 4th Cargèse Summer School, 208 Non-Newtonian Flows in Porous Media: upscaling problems https://www.dropbox.com/s/mcgg0ifpogsznv2/non_newtonian_v00.pdf?dl=0 Davit Y., Zami-Pierre F.,2, de Loubens

More information

Porosity-Permeability Relations in Granular, Fibrous and Tubular Geometries

Porosity-Permeability Relations in Granular, Fibrous and Tubular Geometries Porosity-Permeability Relations in Granular, Fibrous and Tubular Geometries November 21 st Feng Xiao, Xiaolong Yin Colorado School of Mines American Physical Society, 64 th Annual DFD Meeting, Baltimore,

More information

Chapter 5. The Differential Forms of the Fundamental Laws

Chapter 5. The Differential Forms of the Fundamental Laws Chapter 5 The Differential Forms of the Fundamental Laws 1 5.1 Introduction Two primary methods in deriving the differential forms of fundamental laws: Gauss s Theorem: Allows area integrals of the equations

More information

DRAFT. Analytical Model for the Effective Thermal Conductivity Matrix of Anisotropic Composites

DRAFT. Analytical Model for the Effective Thermal Conductivity Matrix of Anisotropic Composites DRAFT Proceedings of IMECE 25 ASME 21 International Mechanical Engineering Congress & Exposition November 12-18 21, Vancouver, British Columbia Canada IMECE21-466 Analytical Model for the Effective Thermal

More information

PERMEABILITY AND THE STRUCTURE OF POROSITY IN PARTICULATE MATERIALS

PERMEABILITY AND THE STRUCTURE OF POROSITY IN PARTICULATE MATERIALS PERMEABILITY AND THE STRUCTURE OF POROSITY IN PARTICULATE MATERIALS A thesis submitted for fulfilment of the requirements for the award of the degree of Doctor of Philosophy from The University of Newcastle

More information

Permeability calculation of sphere-packed porous media using dissipative particle dynamics

Permeability calculation of sphere-packed porous media using dissipative particle dynamics Desalination and Water Treatment www.deswater.com 34 (011) 77 83 October 1944-3994/1944-3986 011 Desalination Publications. All rights reserved doi: 10/5004/dwt.011.80 Permeability calculation of sphere-packed

More information

Study on Estimation of Hydraulic Conductivity of Porous Media Using Drag Force Model Jashandeep Kaur, M. A. Alam

Study on Estimation of Hydraulic Conductivity of Porous Media Using Drag Force Model Jashandeep Kaur, M. A. Alam 26 IJSRSET Volume 2 Issue 3 Print ISSN : 2395-99 Online ISSN : 2394-499 Themed Section: Engineering and Technology Study on Estimation of Hydraulic Conductivity of Porous Media Using Drag Force Model Jashandeep

More information

An Application of Path-Percolation Theory and Lattice-Boltzmann Model on Mass Transfer in Inhomogeneous Porous Media

An Application of Path-Percolation Theory and Lattice-Boltzmann Model on Mass Transfer in Inhomogeneous Porous Media University of Tennessee, Knoxville Trace: Tennessee Research and Creative Exchange Doctoral Dissertations Graduate School 8-2015 An Application of Path-Percolation Theory and Lattice-Boltzmann Model on

More information

GORE: Modeling Filter Compression, Fluid Generation and Parallelizing the Potts Model

GORE: Modeling Filter Compression, Fluid Generation and Parallelizing the Potts Model GORE: Modeling Filter Compression, Fluid Generation and Parallelizing the Potts Model Mathematical Problems in Industry, June 12 17, 2016 June 17, 2016 (Mathematical Problems in Industry, June 12 17, 2016)

More information

NON-DARCY POROUS MEDIA FLOW IN NO-SLIP AND SLIP REGIMES

NON-DARCY POROUS MEDIA FLOW IN NO-SLIP AND SLIP REGIMES THERMAL SCIENCE, Year 2012, Vol. 16, No. 1, pp. 167-176 167 NON-DARCY POROUS MEDIA FLOW IN NO-SLIP AND SLIP REGIMES by Antonio F. MIGUEL Geophysics Centre of Evora & Department of Physics, University of

More information

On pore fluid pressure and effective solid stress in the mixture theory of porous media

On pore fluid pressure and effective solid stress in the mixture theory of porous media On pore fluid pressure and effective solid stress in the mixture theory of porous media I-Shih Liu Abstract In this paper we briefly review a typical example of a mixture of elastic materials, in particular,

More information

Modelling mass and heat transfer in a porous structure

Modelling mass and heat transfer in a porous structure Modelling mass and heat transfer in a porous structure Microclimate in a hygiene product Master s Thesis within the Innovative and Sustainable Chemical Engineering programme EMMA PERSSON Department of

More information

Determination of Locally Varying Directions through Mass Moment of Inertia Tensor

Determination of Locally Varying Directions through Mass Moment of Inertia Tensor Determination of Locally Varying Directions through Mass Moment of Inertia Tensor R. M. Hassanpour and C.V. Deutsch Centre for Computational Geostatistics Department of Civil and Environmental Engineering

More information

Site Characterization & Hydrogeophysics

Site Characterization & Hydrogeophysics Site Characterization & Hydrogeophysics (Source: Matthew Becker, California State University) Site Characterization Definition: quantitative description of the hydraulic, geologic, and chemical properties

More information

Lecture 10. (2) Functions of two variables. Partial derivatives. Dan Nichols February 27, 2018

Lecture 10. (2) Functions of two variables. Partial derivatives. Dan Nichols February 27, 2018 Lecture 10 Partial derivatives Dan Nichols nichols@math.umass.edu MATH 233, Spring 2018 University of Massachusetts February 27, 2018 Last time: functions of two variables f(x, y) x and y are the independent

More information

Flow and Transport. c(s, t)s ds,

Flow and Transport. c(s, t)s ds, Flow and Transport 1. The Transport Equation We shall describe the transport of a dissolved chemical by water that is traveling with uniform velocity ν through a long thin tube G with uniform cross section

More information

Contents. MATH 32B-2 (18W) (L) G. Liu / (TA) A. Zhou Calculus of Several Variables. 1 Multiple Integrals 3. 2 Vector Fields 9

Contents. MATH 32B-2 (18W) (L) G. Liu / (TA) A. Zhou Calculus of Several Variables. 1 Multiple Integrals 3. 2 Vector Fields 9 MATH 32B-2 (8W) (L) G. Liu / (TA) A. Zhou Calculus of Several Variables Contents Multiple Integrals 3 2 Vector Fields 9 3 Line and Surface Integrals 5 4 The Classical Integral Theorems 9 MATH 32B-2 (8W)

More information

MATLAB Solution of Flow and Heat Transfer through a Porous Cooling Channel and the Conjugate Heat Transfer in the Surrounding Wall

MATLAB Solution of Flow and Heat Transfer through a Porous Cooling Channel and the Conjugate Heat Transfer in the Surrounding Wall MATLAB Solution of Flow and Heat Transfer through a Porous Cooling Channel and the Conjugate Heat Transfer in the Surrounding Wall James Cherry, Mehmet Sözen Grand Valley State University, cherryj1@gmail.com,

More information

Uniformity of the Universe

Uniformity of the Universe Outline Universe is homogenous and isotropic Spacetime metrics Friedmann-Walker-Robertson metric Number of numbers needed to specify a physical quantity. Energy-momentum tensor Energy-momentum tensor of

More information

Inverse problems in tumor growth modelling

Inverse problems in tumor growth modelling MC2 Damiano Lombardi Angelo Iollo Thierry Colin Olivier Saut Marseille, 18-09-09 Outline 1) Introduction: phenomenological model for tumor growth models 2) Models based on mixture theory: a simplified

More information

LBM mesoscale modelling of porous media

LBM mesoscale modelling of porous media Advanced Computational Methods in Heat Transfer X 59 LBM mesoscale modelling of porous media A. C. M. Sousa 1,2 & A. Nabovati 1 1 Department of Mechanical Engineering, University of New Brunswick, Fredericton,

More information

MATH20411 PDEs and Vector Calculus B

MATH20411 PDEs and Vector Calculus B MATH2411 PDEs and Vector Calculus B Dr Stefan Güttel Acknowledgement The lecture notes and other course materials are based on notes provided by Dr Catherine Powell. SECTION 1: Introctory Material MATH2411

More information

Simulating Fluid-Fluid Interfacial Area

Simulating Fluid-Fluid Interfacial Area Simulating Fluid-Fluid Interfacial Area revealed by a pore-network model V. Joekar-Niasar S. M. Hassanizadeh Utrecht University, The Netherlands July 22, 2009 Outline 1 What s a Porous medium 2 Intro to

More information

Erratum to: Anisotropic Permeability of Trabecular Bone and Its Relationship to Fabric and Architecture: A Computational Study

Erratum to: Anisotropic Permeability of Trabecular Bone and Its Relationship to Fabric and Architecture: A Computational Study Annals of Biomedical Engineering, Vol. 45, No. 6, June 2017 (Ó 2017) pp. 1555 1559 DOI: 10.1007/s10439-017-1835-3 Erratum Erratum to: Anisotropic Permeability of Trabecular Bone and Its Relationship to

More information

Four Verification Cases for PORODRY

Four Verification Cases for PORODRY Four Verification Cases for PORODRY We designed four different verification cases to validate different aspects of our numerical solution and its code implementation, and these are summarized in Table

More information

Appendix: Nomenclature

Appendix: Nomenclature Appendix: Nomenclature 209 Appendix: Nomenclature Chapter 2: Roman letters A amplitude of surface wave m A j (complex) amplitude of tidal component [various] a right ascension rad or a magnitude of tidal

More information

MULTIVARIABLE CALCULUS

MULTIVARIABLE CALCULUS MULTIVARIABLE CALCULUS JOHN QUIGG Contents 13.1 Three-Dimensional Coordinate Systems 2 13.2 Vectors 3 13.3 The Dot Product 5 13.4 The Cross Product 6 13.5 Equations of Lines and Planes 7 13.6 Cylinders

More information

Calculation of Permeability Tensors for Unstructured Grid Blocks

Calculation of Permeability Tensors for Unstructured Grid Blocks Calculation of Permeability Tensors for Unstructured Grid Blocs R. M. Hassanpour, O. Leuangthong and C.V. Deutsch Centre for Computational Geostatistics Department of Civil and Environmental Engineering

More information

Maxwell s Equations:

Maxwell s Equations: Course Instructor Dr. Raymond C. Rumpf Office: A-337 Phone: (915) 747-6958 E-Mail: rcrumpf@utep.edu Maxwell s Equations: Physical Interpretation EE-3321 Electromagnetic Field Theory Outline Maxwell s Equations

More information