Investigating the role of tortuosity in the Kozeny-Carman equation
|
|
- Moses Stevens
- 6 years ago
- Views:
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
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 informationTwo-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 informationHydraulic 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 information1 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 informationVI. 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 informationHow 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 informationTransfer 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 informationKOZENY-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 informationAn 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 informationq 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 informationPermeability 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 informationFluid 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 informationMobility 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 informationGetting 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 informationNumber 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 informationarxiv: 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 informationPermeability 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 informationStudies 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 informationRandom 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 informationON 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 informationIntroduction 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 informationNumber 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 informationNUMERICAL 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 informationFluid 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 informationIn 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 informationThree-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 informationEvaporation-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 informationChapter 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 informationGeometry 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 informationMathematical 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 informationTortuosity-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 informationTHE 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 informationarxiv: 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 informationAnisotropic 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 informationPORE-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 informationEffects 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 informationAvailable 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 informationDarcy 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 informationA 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 informationExperimental 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 informationAdvanced 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 informationStreamline 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 informationUNCLASSI 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 informationENGI 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 informationModeling 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 informationM98-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 informationDarcy'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 informationMATH 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 informationEffect 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 informationMATH 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 informationThis 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 informationMODELLING 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 informationDISPERSION 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 informationDissolution 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 informationModeling 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 informationDiffusion 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 informationI. 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 informationTable 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 informationStokes 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 informationFluid 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 informationThe 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 informationwater 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 informationPRACTICE 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.
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 informationIntroduction 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 informationSCA : 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 information2. 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 informationOn 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 informationFlow 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 informationA 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 information1.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 informationNon-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 informationPorosity-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 informationChapter 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 informationDRAFT. 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 informationPERMEABILITY 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 informationPermeability 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 informationStudy 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 informationAn 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 informationGORE: 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 informationNON-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 informationOn 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 informationModelling 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 informationDetermination 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 informationSite Characterization & Hydrogeophysics
Site Characterization & Hydrogeophysics (Source: Matthew Becker, California State University) Site Characterization Definition: quantitative description of the hydraulic, geologic, and chemical properties
More informationLecture 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 informationFlow 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 informationContents. 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 informationMATLAB 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 informationUniformity 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 informationInverse 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 informationLBM 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 informationMATH20411 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 informationSimulating 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 informationErratum 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 informationFour 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 informationAppendix: 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 informationMULTIVARIABLE 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 informationCalculation 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 informationMaxwell 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