How are calculus, alchemy and forging coins related?
|
|
- Matilda Sharp
- 5 years ago
- Views:
Transcription
1 BMOLE 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 J. Bishop Assistant Professor of Biomedical Engineering Principal Investigator of the Pharmacoengineering Laboratory: pharmacoengineering.com Dwight Look College of Engineering Texas A&M University Emerging Technologies Building Room 5016 College Station, TX How are calculus, alchemy and forging coins related? 1
2 Porosity, Tortuosity, and Available Volume Fraction Chapter 8, Section 2 2
3 Specific Surface = s = Total interface area Total volume Units = 1 length Porosity = ε = Void volume Total volume ε is dimensionless Void volume: total volume of void space in a porous medium Interface: border between solid and void spaces If total volume based on interstitial space: ε is larger than 0.9 If cells and vessels need to be considered: ε varies between 0.06 and 0.30 Note: In some tumors, ε is as high as 0.6 3
4 Porosity (ε) Porosity: a measure of the average void volume fraction in a specific region of porous medium Does not provide information on how pores are connected or number of pores available for water and solute transport ε = ε i + ε p + ε n Note: isolated pores are not accessible to external solvents and solutes: They can sometimes be considered part of the solid phase. In this case, void volume is the total volume of penetrable pores 4
5 Tortuosity (T) Path length between A and B in a porous medium is the distance between the points through connected pores L min is the shorted path length L is the straight-line distance between A and B T is always greater than or equal to unity T = ( L min L )2 5
6 Available Volume Fraction (K AV ) Not all penetrable pores are accessible to solutes: Dependent on molecular properties of the solutes Pore smaller than solute molecule or is surrounded by too small pores Available Volume: portion of accessible volume that can be occupied by the solute K AV = Available volume Total volume 6
7 K AV is molecule dependent and always smaller than porosity. This can be caused by 3 scenarios: 1)Centers of the solute molecules cannot reach the solid surface in the void space Difference between total void volume and the available volume can be estimated as the product of the area of the surface and distance ( ) between solute and surface 2)Some of the void space is smaller than the solute molecules 3)Inaccessibility of large penetrable pores surrounded by pores smaller than the solutes Generally, K AV decreases with size of solutes 7
8 Partition Coefficient (Φ): ratio of available volume to void volume Measure of solute partitioning at equilibrium between external solutions and void space in porous media Φ = K AV ε Radius of Gyration (R g ): size of flexible molecules R g 2 = 1 N i r ic 2 = 1 2N 2 i j r ij 2 r ic is the average distance between segment i and the center of mass of the polymer chain r ij is the average distance between segments i and j in the polymer chain N is the total number of segments in the polymer chain 8
9 When flexible polymer can be modelled as a linear chain of rigid rods (segments) and the length of a segment (l) is much less than the pore radius (l/a<<1): λ s is the ratio of R g to the radius of the spherical pore 1 Φ sphere = 6 π 2 n 2 exp[ n2 π 2 λ 2 s ] n=1 Φ plate = 8 π 2 n=0 λ p is the ratio of R g to the half-distance between plates 1 (2n + 1) 2 exp[ (2n + 1)2 π 2 λ 2 4 p ] Φ cylinder = 4 1 n=1 α n 2 exp[ α n 2 λ c 2 ] λ c is the ratio of R g to the radius of the cylindrical pore and α n are the roots of J 0 (α n )=0 J 0 is the Bessel function of the first kind of order zero Bessel Function = x d dx dy x dx + m2 x 2 X 2 y = 0 9
10 Exclusion Volume Some porous media are fiber matrices: space inside and near surface of fibers not available to solutes Exclusion Volume: size of the space r f and r s are the radii of the fiber and solute L is the length of the fiber Exclusion volume = π(r f + r s ) 2 L If minimum distance between fibers is larger than 2(r f + r s ), such as when fiber density is very low or when fibers are parallel: Exclusion volume fraction = π(r f + r s ) 2 LN V θ is the volume fraction of fibers K AV = 1 exclusion volume fraction = 1 θ = θ( r s r f + 1) 2 r s r f
11 Ogston Equation If the minimum distance between fibers is smaller than 2(r f + r s ): K AV = exp[ θ(1 + r s r f ) 2 ] Assumes: Fibers are rigid rods with random orientation Length of fibers is finite but much larger than r s There is no fiber overlap Concentration of solute is low (steric interaction between nearest neighbors of solute and fiber) Fiber radius much less than interfiber distance 11
12 When θ(1 + r s r f ) 2 is much less than unity, Ogston Equation reduces to previous equation. Porosity can be derived from the Ogston Equation by letting r s =0 When θ is much less than unity: ε = exp[ θ] ε = 1 θ Partition coefficient of solutes in liquid phase of fiber-matrix material: Φ = exp[ θ(1 + ( r s rf )) 2 ] 1 θ 12
13 8.3.2: Brinkman Equation 8.3.3: Squeeze Flow
14 Review: Darcy s Law o Describes fluid flow in porous media o Invalid for Non-Newtonian Fluids o We talked about δ, l, and L o L is the distance over which macroscopic changes of physical quantities have to be considered o V = K p o K is the hydraulic conductivity (a constant) o p is the gradient of hydrostatic pressure o V = K p = ϕ B ϕ L o ϕ B is volumetric flow in source o ϕ L is volumetric flow in sink o Generally ϕ B = ϕ L = 0 o Darcy s Law Neglects Friction within the fluid
15 Brinkman Equation ok = Kμ o k is the specific hydraulic permeability (usually units of nm 2 ) o Darcy s law is used when k is low: when k is much smaller than the square of L o When k is not low: Brinkman Equation obrinkman Equation: o μ 2 V 1 V p=0 K o Darcy s Law is a special case of this equation when the first term = 0 (V = K p) Example 8.7!! (The math is annoying and long and I m sorry) (It s not my fault) (Don t hate me)
16 Ex 8.7
17 Squeeeze Flow o Squeeze flow is the fluid flow caused by the relative movement of solid boundaries towards each other o Tissue deformation causes change in volume fraction of the interstitial space o Leads to fluid flow
18 Squeeze Flow ovh is the velocity of cell membrane movement or is the radius of the plate othese equations will be more useful and will make more sense in later discussions (hopefully) othe derivation of them is in ovelocity profiles: oz and r are cylindrical coordinates ovr is velocity in the r direction ovz is velocity in the z direction
19 Darcy s Law
20 Darcy s Law Flow rate is proportional to pressure gradient Valid in many porous media NOT valid for: Non-Newtonian fluids Newtonian liquids at high velocities Gasses at very low and high velocities
21 Darcy s Law Describing fluid flow in porous media Two Ways: Numerically solve governing equations for fluid flow in individual pores if structure is known Assume porous medium is a uniform material (Continuum Approach)
22 Three Length Scales 1. δ = average size of pores 2. L = distance which macroscopic changes of physical quantities must be considered (ie. Fluid velocity and pressure) 3. l = a small volume with dimension l Continuum Approach Requires: δ << l 0 << L
23 Darcy s Law l 0 = Representative Elementary Volume (REV l 0 << L Total volume of REV can give the averaging over a volume value v f = fluid velocity, velocity of each fluid particle AVERAGED over volume of the fluid phase v = velocity of each fluid particle averaged over REV V = εv f
24 Darcy s Law Law of Mass Conservation No Fluid Production = Source Fluid Consumption = Sink Mass Balance with velocities v = φ B - φ L Ev f = φ B - φ L Determined by Starling s Law φ B = rate of volumetric flow in sources φ L = rate of volumetric flow in sinks
25 Darcy s Law Momentum Balance in Porous Media v = -K p p = gradient of hydrostatic pressure K = hydraulic conductivity constant p = average quantity within the fluid phase in the REV
26 Darcy s Law Substituting the Equations to form: ( K p) = φ B - φ L Steady State: 2 * p = 0
27 8.4 Solute Transport in Porous Media
28 8.4.1 General Considerations Solute Transport in Porous Media
29 8.4.1 General Considerations There are 4 general problems using the continuum approach(darcy s Law) for analyzing transport of solutes through porous media 1)D eff 2) Solute velocity 3) Dispersion 4) Boundary Conditions
30 1) D eff Diffusion of solutes is characterized by the effective diffusion coefficient: D eff D eff in porous media < D eff in solutions Why?...
31 1) D eff cont. Factor effecting Diffusion in Porous Media: Connectedness of Pores Layer 100 D eff Layer Layer # As we near layer 100, available space for diffusion and D eff decrease exponentially
32 2) Solute Velocity Convective velocities Solvent velocities Solutes are hindered by porous structures Ex: filtering coffee grounds f = v s v f = solute velocity solvent velocity f is retardation coefficient
33 2) Solvent Velocity cont. σ = 1 f : reflection coefficient Characterizes the hindrance of convective transport across a membrane f is dependent on: Fluid velocity Solute size Pore Structure Flux of Convective transport across tissue: N s = v s C = fv f C Where v s is solute velocity and C is local concentration of solutes
34 3) Dispersion of Solutes
35 4) Boundary Conditions Concentrations of solutes can be discontinuous at interfases between solutions and porous media Thus, for two regions, 1 and 2 N 1 = N 2 N = flux and C 1 K AA1 = C 2 K AA2 K AA = area fraction available at inerfase for solute transport
36 When all 4 are considered Governing Equation for Transport of neutral molecule through porous media: C t + fv fc = D eff 2 C + φ B φ L + Q Flux of Convective Transport Isotropic, uniform and dispersion coefficient is negligible Solute vel. Reaction If we do consider dispersion coefficient: Deff = Deff + Disp. Coeff.
37 Governing Equation cont. C t + fv fc = D eff 2 C + φ B φ L + Q C, φ B, φ L, Q: avg. quantity unit tissue volume f, v f, D eff : avg. quantity unit volume of fluid phase
38 8.4.2 Effective Diffusion Coefficient in Hydrogels 3 Factors effecting D eff of uncharged solutes in hydrogels: 1) Diffusion coefficient of solutes in WATER (D 0 ) 2) Hydrodynamic Interactions between solute and surrounding solvent molecules, F 3) Tortuosity of diffusion pathways due to the steric exclusion of solutes in the matrix, S D eff = D 0 FS So, how do we solve for F and S?
39 Hydrodynamic Interactions, F F is a ratio: F = friction coeff. of solute in porous media friction coeff. of solute inwater Friction coeff. in water = 6πμr s μ = viscocity r s = radius of molecule F measures the enhancement of drag on solute molecule due to presence of polymeric fibers in water
40 Hydrodynamic Interactions, F cont. Two approaches to determining F: 1) Effectivemedium or Brinkman-medium 2) 3-D space with cylindrical fibers model
41 1) Effective-Medium or Brinkman-Medium Assumptions: Hydrogel is a uniform medium Spherical solute molecule Constant velocity Movement governed by Brinkman Equation F = 1 + r s k r s k 2 1 where φ B = φ L = 0
42 2) 3-D space with cylindrical fibers model Assumptions: Hydrogel is modeled as 3-D spaces filled with water and randomly placed cylindrical fibers Movement of spherical particles is determined by Stokes-Einstein Equation F α, θ = exp( a 1 θ a 2) pg. 426 done after normalization of 6πμr s
43 2) 3-D space with cylindrical fibers model cont. F α, θ = exp( a 1 θ a 2) Depends on: α = θ = fiber radii solute radii fiber volume hydrogel volume Where: a 1 = α α 2 a 2 = α α 2
44 Tortuosity Factor, S Depends on f a f a = (1 + 1 α )2 f a is the excluded volume fraction of solute in hydrogel if: Low fiber density Fibers arranged in parallel manner If f a < 0.7 S α, θ = exp 0.84f a 1.09
45 Effective Diffusion Coefficient: In Liquid-Filled Pore and Biological Tissue Chapter 8: Section
46 Effective Diffusion Coefficient in a Liquid- Filled Pore Depends on diffusion coefficient D 0 of solutes in water, hydrodynamic interactions between solute and solvent molecules, and steric exclusion of solutes near the walls of pores Assume entrance effect is negligible: v z = 2v m (1 r R )2 v z is the axial fluid velocity v m is the mean velocity in the pore r is the radial coordinate R is the radii of the cylinder 46
47 Assume the volume fraction of spherical solutes is less than 2λ 3 λ = a R a is the radii of the solute R is the radii of the cylinder Solute-Solute interactions become negligible When λ ->0, solute-pore interactions are negligible N z = D 0 C z + Cv z C z 0 r R a C = ቊ 0 R a < r R N z is the flux D 0 is the diffusion coefficient C is the solute concentration z is the axial coordinate 47
48 When λ does not approach 0: N s = D 0 K dc dz + GCv z K is the enhanced friction coefficient G is the lag coefficient K and G are functions of λ and r/r Flux averaged over the entire cross-sectional area: ഥN s = 2 R 2 න 0 R N s rdr Integrating the top flux equation: ഥN s = HD 0 dc dz + WCv m 48
49 H and W are called hydrodynamic resistance coefficients H = 2 R 2 න 0 R a 1 K rdr W = 4 R 2 න 0 R a G 1 r R 2 rdr D eff = HD 0 Centerline approximation: assume all spheres are distributed on the centerline position in the pore For λ < 0.4: K 1 λ, 0 = λ λ λ 5 G λ, 0 = λ λ 3 49
50 H λ = Φ λ λ λ 5 W λ = Φ 2 Φ λ λ 3 Partition coefficient of solute in the pore: Φ = 1 λ 2 For λ > 0.4: 50
51 Diffusion of spherical molecules between parallel plates: ഥN s = HD 0 dc dz + WCv m ഥN s = 1 h න 0 H = 1 h න 0 h N s dy h a 1 K dy h is the half-width of the slit W = 3 h a 2h න G 1 y 0 h 2 dy K and G depend on λ and y/h D eff = HD 0 H λ = Φ λ λ λ λ 5 + Oλ 6 W λ = Φ 2 3 Φ λ2 + Oλ 3 Φ = 1 λ 51
52 Effective Diffusion Coefficient in Biological Tissues D eff = b 1 M r b 2 M r is the molecular weight of the solutes b 1 and b 2 are functions of charge and shape of solutes, and structures of tissues The effects of tissue structures on D eff increases with the size of solutes 52
53 Fluid Transport in Poroelastic Materials Section 8.5 By: Hannah Humbert 53
54 Biological Tissues: Deformable Deformation can be linear or nonlinear "Extreme softness of brain matter in simple shear" 54
55 Poroelastic Response σ yy σ xy Porosity, ε y Pore Fluid Pressure, p x σ xx σ xx Total Stress = Effective Stress + Pore Pressure σ = τ + pi σ yy 55
56 Understanding the linear function of Stress If, I = then, σ = τ xx τ xy τ xz τ yx τ yy τ yz τ zx τ zy τ zz + p p p = 2μ G E + μ λ ei + pi Must be describing stresses that are NOT pore pressure related 56
57 Understanding the linear Function of Stress σ = 2μ G E + μ λ ei + pi We know this is pore pressure, and this is the only force happening INTERNALLY, right? We must look at the Lamé constants Then, these two terms must be describing the forces EXTERNALLY. How many external forces are there? Shear and Normal. Which term is which? 57
58 LamÉ Constants Named after French Mathematician, Gabriel Lamé Two constants which relate stress to strain in isotropic, elastic material Depend on the material and its temperature μ λ : first parameter (related to bulk modulus) μ G : second parameter (related to shear modulus) μ λ = K 2 3 μ G μ G = τ γ = F A L L 58
59 Understanding the linear Function of Stress σ = 2μ G E + μ λ ei + pi Describes EXTERNAL SHEAR FORCES Describes INTERNAL FORCES Describes EXTERNAL BULK FORCES Now we need to understand E and e 59
60 Understanding the linear Function of Stress E = strain tensor (vector gradient) E = 1 2 [ u + u T ] u = u x x, u y y, u z z ( u) T = u x x u y y u z z E = u x x u x x + u y y u x x + u z z u y y + u x x 2 u y y u y y + u z z u z z + u x x u z z + u y y 2 u z z e = volume dilation (scalar divergence) e = Tr E = u e = u = u x x + u y y + u z z 60
61 Understanding the linear Function of Stress So now we understand this: A tensor that s in units of stress σ = 2μ G E + μ λ ei + pi A scalar quantity that s multiplied by the identify matrix to make a tensor in units of pressure Now we can reduce mass conservation equations to get an equation we can use to solve problems like these! A scalar quantity that s multiplied by the identify matrix to make a tensor in units of stress 61
62 Deriving mass and Momentum equations mass conservation in the fluid phase is this: (ρ f ε) t + ρ f εv f = ρ f (φ B φ L ) mass conservation in the solid phase is this: [ρ s 1 ε ] t + ρ s (1 ε) u t = 0 62
63 Deriving mass and Momentum equations Summing the two equations together: (ρ f ε) + ρ t f εv f ρ f φ B φ L = [ρ s 1 ε ] t (ρ f ε) + ρ t f εv f ρ s 1 ε u t = ρ s 1 ε t + ρ s (1 ε) u t + ρ f φ B φ L εv f 1 ε u t = φ B φ L 63
64 Biot Law According to a paper written in 1984 (that I do not have access too): ε v f u t = K p σ = 0 Zienkiewicz, O. C., and T. Shiomi. "Dynamic behaviour of saturated porous media; the generalized Biot formulation and its numerical solution." International journal for numerical and analytical methods in geomechanics 8.1 (1984):
65 Using BIOTS law to substitute σ = μ G 2 u + (μ G + μ λ ) e + p = 0 (2μ G + μ λ ) 2 e = 2 p 65
66 If volume dilation and hydraulic conductivity are homogenous ε v f e t = K 2 p e t = K 2 p + φ B φ L e t = K(2μ G + μ λ ) 2 e + φ B φ L Coefficient of consolidation 66
Chapter 9: Differential Analysis of Fluid Flow
of Fluid Flow Objectives 1. Understand how the differential equations of mass and momentum conservation are derived. 2. Calculate the stream function and pressure field, and plot streamlines for a known
More informationChapter 9: Differential Analysis
9-1 Introduction 9-2 Conservation of Mass 9-3 The Stream Function 9-4 Conservation of Linear Momentum 9-5 Navier Stokes Equation 9-6 Differential Analysis Problems Recall 9-1 Introduction (1) Chap 5: Control
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 informationV (r,t) = i ˆ u( x, y,z,t) + ˆ j v( x, y,z,t) + k ˆ w( x, y, z,t)
IV. DIFFERENTIAL RELATIONS FOR A FLUID PARTICLE This chapter presents the development and application of the basic differential equations of fluid motion. Simplifications in the general equations and common
More information3D Elasticity Theory
3D lasticity Theory Many structural analysis problems are analysed using the theory of elasticity in which Hooke s law is used to enforce proportionality between stress and strain at any deformation level.
More informationNumerical Heat and Mass Transfer
Master Degree in Mechanical Engineering Numerical Heat and Mass Transfer 15-Convective Heat Transfer Fausto Arpino f.arpino@unicas.it Introduction In conduction problems the convection entered the analysis
More informationNANO 243/CENG 207 Course Use Only
L8. Drug Dispersion and Diffusion in Biological Systems (2) April 26, 2018 2. Diffusion in Water (3) Diffusion of small molecule in water drug molecules/solvents are identical to solvent Assume molecules
More informationChapter 1 Fluid Characteristics
Chapter 1 Fluid Characteristics 1.1 Introduction 1.1.1 Phases Solid increasing increasing spacing and intermolecular liquid latitude of cohesive Fluid gas (vapor) molecular force plasma motion 1.1.2 Fluidity
More informationCENG 501 Examination Problem: Estimation of Viscosity with a Falling - Cylinder Viscometer
CENG 501 Examination Problem: Estimation of Viscosity with a Falling - Cylinder Viscometer You are assigned to design a fallingcylinder viscometer to measure the viscosity of Newtonian liquids. A schematic
More informationLecture 8: Tissue Mechanics
Computational Biology Group (CoBi), D-BSSE, ETHZ Lecture 8: Tissue Mechanics Prof Dagmar Iber, PhD DPhil MSc Computational Biology 2015/16 7. Mai 2016 2 / 57 Contents 1 Introduction to Elastic Materials
More informationBiotransport: Principles
Robert J. Roselli Kenneth R. Diller Biotransport: Principles and Applications 4 i Springer Contents Part I Fundamentals of How People Learn (HPL) 1 Introduction to HPL Methodology 3 1.1 Introduction 3
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 informationMechanics PhD Preliminary Spring 2017
Mechanics PhD Preliminary Spring 2017 1. (10 points) Consider a body Ω that is assembled by gluing together two separate bodies along a flat interface. The normal vector to the interface is given by n
More informationFundamentals of Fluid Dynamics: Elementary Viscous Flow
Fundamentals of Fluid Dynamics: Elementary Viscous Flow Introductory Course on Multiphysics Modelling TOMASZ G. ZIELIŃSKI bluebox.ippt.pan.pl/ tzielins/ Institute of Fundamental Technological Research
More informationOutline. Definition and mechanism Theory of diffusion Molecular diffusion in gases Molecular diffusion in liquid Mass transfer
Diffusion 051333 Unit operation in gro-industry III Department of Biotechnology, Faculty of gro-industry Kasetsart University Lecturer: Kittipong Rattanaporn 1 Outline Definition and mechanism Theory of
More informationBasic Equations of Elasticity
A Basic Equations of Elasticity A.1 STRESS The state of stress at any point in a loaded bo is defined completely in terms of the nine components of stress: σ xx,σ yy,σ zz,σ xy,σ yx,σ yz,σ zy,σ zx,andσ
More informationMacroscopic theory Rock as 'elastic continuum'
Elasticity and Seismic Waves Macroscopic theory Rock as 'elastic continuum' Elastic body is deformed in response to stress Two types of deformation: Change in volume and shape Equations of motion Wave
More informationMechanics of Earthquakes and Faulting
Mechanics of Earthquakes and Faulting www.geosc.psu.edu/courses/geosc508 Standard Solids and Fracture Fluids: Mechanical, Chemical Effects Effective Stress Dilatancy Hardening and Stability Mead, 1925
More information12. Stresses and Strains
12. Stresses and Strains Finite Element Method Differential Equation Weak Formulation Approximating Functions Weighted Residuals FEM - Formulation Classification of Problems Scalar Vector 1-D T(x) u(x)
More informationGame Physics. Game and Media Technology Master Program - Utrecht University. Dr. Nicolas Pronost
Game and Media Technology Master Program - Utrecht University Dr. Nicolas Pronost Soft body physics Soft bodies In reality, objects are not purely rigid for some it is a good approximation but if you hit
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 informationSoft Bodies. Good approximation for hard ones. approximation breaks when objects break, or deform. Generalization: soft (deformable) bodies
Soft-Body Physics Soft Bodies Realistic objects are not purely rigid. Good approximation for hard ones. approximation breaks when objects break, or deform. Generalization: soft (deformable) bodies Deformed
More informationNDT&E Methods: UT. VJ Technologies CAVITY INSPECTION. Nondestructive Testing & Evaluation TPU Lecture Course 2015/16.
CAVITY INSPECTION NDT&E Methods: UT VJ Technologies NDT&E Methods: UT 6. NDT&E: Introduction to Methods 6.1. Ultrasonic Testing: Basics of Elasto-Dynamics 6.2. Principles of Measurement 6.3. The Pulse-Echo
More informationMultiple Integrals and Vector Calculus (Oxford Physics) Synopsis and Problem Sets; Hilary 2015
Multiple Integrals and Vector Calculus (Oxford Physics) Ramin Golestanian Synopsis and Problem Sets; Hilary 215 The outline of the material, which will be covered in 14 lectures, is as follows: 1. Introduction
More informationENGR Heat Transfer II
ENGR 7901 - Heat Transfer II Convective Heat Transfer 1 Introduction In this portion of the course we will examine convection heat transfer principles. We are now interested in how to predict the value
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 informationElements of Rock Mechanics
Elements of Rock Mechanics Stress and strain Creep Constitutive equation Hooke's law Empirical relations Effects of porosity and fluids Anelasticity and viscoelasticity Reading: Shearer, 3 Stress Consider
More informationChapter 2 Governing Equations
Chapter Governing Equations Abstract In this chapter fundamental governing equations for propagation of a harmonic disturbance on the surface of an elastic half-space is presented. The elastic media is
More informationLecture 8. Stress Strain in Multi-dimension
Lecture 8. Stress Strain in Multi-dimension Module. General Field Equations General Field Equations [] Equilibrium Equations in Elastic bodies xx x y z yx zx f x 0, etc [2] Kinematics xx u x x,etc. [3]
More informationChapter 1. Continuum mechanics review. 1.1 Definitions and nomenclature
Chapter 1 Continuum mechanics review We will assume some familiarity with continuum mechanics as discussed in the context of an introductory geodynamics course; a good reference for such problems is Turcotte
More informationPEAT SEISMOLOGY Lecture 2: Continuum mechanics
PEAT8002 - SEISMOLOGY Lecture 2: Continuum mechanics Nick Rawlinson Research School of Earth Sciences Australian National University Strain Strain is the formal description of the change in shape of a
More informationSamantha Ramirez, MSE. Stress. The intensity of the internal force acting on a specific plane (area) passing through a point. F 2
Samantha Ramirez, MSE Stress The intensity of the internal force acting on a specific plane (area) passing through a point. Δ ΔA Δ z Δ 1 2 ΔA Δ x Δ y ΔA is an infinitesimal size area with a uniform force
More informationRheology of Soft Materials. Rheology
Τ Thomas G. Mason Department of Chemistry and Biochemistry Department of Physics and Astronomy California NanoSystems Institute Τ γ 26 by Thomas G. Mason All rights reserved. γ (t) τ (t) γ τ Δt 2π t γ
More information3 2 6 Solve the initial value problem u ( t) 3. a- If A has eigenvalues λ =, λ = 1 and corresponding eigenvectors 1
Math Problem a- If A has eigenvalues λ =, λ = 1 and corresponding eigenvectors 1 3 6 Solve the initial value problem u ( t) = Au( t) with u (0) =. 3 1 u 1 =, u 1 3 = b- True or false and why 1. if A is
More informationFLUID MECHANICS PROF. DR. METİN GÜNER COMPILER
FLUID MECHANICS PROF. DR. METİN GÜNER COMPILER ANKARA UNIVERSITY FACULTY OF AGRICULTURE DEPARTMENT OF AGRICULTURAL MACHINERY AND TECHNOLOGIES ENGINEERING 1 5. FLOW IN PIPES 5.1.3. Pressure and Shear Stress
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 informationChapter 2 CONTINUUM MECHANICS PROBLEMS
Chapter 2 CONTINUUM MECHANICS PROBLEMS The concept of treating solids and fluids as though they are continuous media, rather thancomposedofdiscretemolecules, is one that is widely used in most branches
More informationTensor Visualization. CSC 7443: Scientific Information Visualization
Tensor Visualization Tensor data A tensor is a multivariate quantity Scalar is a tensor of rank zero s = s(x,y,z) Vector is a tensor of rank one v = (v x,v y,v z ) For a symmetric tensor of rank 2, its
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 informationShell Balances in Fluid Mechanics
Shell Balances in Fluid Mechanics R. Shankar Subramanian Department of Chemical and Biomolecular Engineering Clarkson University When fluid flow occurs in a single direction everywhere in a system, shell
More informationAE/ME 339. K. M. Isaac Professor of Aerospace Engineering. 12/21/01 topic7_ns_equations 1
AE/ME 339 Professor of Aerospace Engineering 12/21/01 topic7_ns_equations 1 Continuity equation Governing equation summary Non-conservation form D Dt. V 0.(2.29) Conservation form ( V ) 0...(2.33) t 12/21/01
More informationChemical Reaction Engineering Prof. Jayant Modak Department of Chemical Engineering Indian Institute of Science, Bangalore
Chemical Reaction Engineering Prof. Jayant Modak Department of Chemical Engineering Indian Institute of Science, Bangalore Lecture No. # 26 Problem solving : Heterogeneous reactions Friends, in last few
More informationElements of Continuum Elasticity. David M. Parks Mechanics and Materials II February 25, 2004
Elements of Continuum Elasticity David M. Parks Mechanics and Materials II 2.002 February 25, 2004 Solid Mechanics in 3 Dimensions: stress/equilibrium, strain/displacement, and intro to linear elastic
More information20. Rheology & Linear Elasticity
I Main Topics A Rheology: Macroscopic deformation behavior B Linear elasticity for homogeneous isotropic materials 10/29/18 GG303 1 Viscous (fluid) Behavior http://manoa.hawaii.edu/graduate/content/slide-lava
More informationContinuum Mechanics. Continuum Mechanics and Constitutive Equations
Continuum Mechanics Continuum Mechanics and Constitutive Equations Continuum mechanics pertains to the description of mechanical behavior of materials under the assumption that the material is a uniform
More informationRHEOLOGY Principles, Measurements, and Applications. Christopher W. Macosko
RHEOLOGY Principles, Measurements, and Applications I -56081-5'79~5 1994 VCH Publishers. Inc. New York Part I. CONSTITUTIVE RELATIONS 1 1 l Elastic Solid 5 1.1 Introduction 5 1.2 The Stress Tensor 8 1.2.1
More informationdynamics of f luids in porous media
dynamics of f luids in porous media Jacob Bear Department of Civil Engineering Technion Israel Institute of Technology, Haifa DOVER PUBLICATIONS, INC. New York Contents Preface xvii CHAPTER 1 Introduction
More informationMechanics of Earthquakes and Faulting
Mechanics of Earthquakes and Faulting www.geosc.psu.edu/courses/geosc508 Overview Milestones in continuum mechanics Concepts of modulus and stiffness. Stress-strain relations Elasticity Surface and body
More informationPart I.
Part I bblee@unimp . Introduction to Mass Transfer and Diffusion 2. Molecular Diffusion in Gasses 3. Molecular Diffusion in Liquids Part I 4. Molecular Diffusion in Biological Solutions and Gels 5. Molecular
More informationStress, Strain, Mohr s Circle
Stress, Strain, Mohr s Circle The fundamental quantities in solid mechanics are stresses and strains. In accordance with the continuum mechanics assumption, the molecular structure of materials is neglected
More informationMultiple Integrals and Vector Calculus: Synopsis
Multiple Integrals and Vector Calculus: Synopsis Hilary Term 28: 14 lectures. Steve Rawlings. 1. Vectors - recap of basic principles. Things which are (and are not) vectors. Differentiation and integration
More informationA First Course on Kinetics and Reaction Engineering Unit D and 3-D Tubular Reactor Models
Unit 34. 2-D and 3-D Tubular Reactor Models Overview Unit 34 describes two- and three-dimensional models for tubular reactors. One limitation of the ideal PFR model is that the temperature and composition
More informationFinite Element Method in Geotechnical Engineering
Finite Element Method in Geotechnical Engineering Short Course on + Dynamics Boulder, Colorado January 5-8, 2004 Stein Sture Professor of Civil Engineering University of Colorado at Boulder Contents Steps
More informationRotational & Rigid-Body Mechanics. Lectures 3+4
Rotational & Rigid-Body Mechanics Lectures 3+4 Rotational Motion So far: point objects moving through a trajectory. Next: moving actual dimensional objects and rotating them. 2 Circular Motion - Definitions
More informationSupporting Information. Technique for real-time measurements of endothelial permeability in a
Supporting Information Technique for real-time measurements of endothelial permeability in a microfluidic membrane chip using laser-induced fluorescence detection Edmond W.K. Young a,b,, Michael W.L. Watson
More informationFluid Mechanics II Viscosity and shear stresses
Fluid Mechanics II Viscosity and shear stresses Shear stresses in a Newtonian fluid A fluid at rest can not resist shearing forces. Under the action of such forces it deforms continuously, however small
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 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 informationPolymers Reactions and Polymers Production (3 rd cycle)
EQ, Q, DEQuim, DQuim nd semester 017/018, IST-UL Science and Technology of Polymers ( nd cycle) Polymers Reactions and Polymers Production (3 rd cycle) Lecture 5 Viscosity easurements of the viscosity
More informationRock Rheology GEOL 5700 Physics and Chemistry of the Solid Earth
Rock Rheology GEOL 5700 Physics and Chemistry of the Solid Earth References: Turcotte and Schubert, Geodynamics, Sections 2.1,-2.4, 2.7, 3.1-3.8, 6.1, 6.2, 6.8, 7.1-7.4. Jaeger and Cook, Fundamentals of
More informationInvestigating the role of tortuosity in the Kozeny-Carman equation
Investigating the role of tortuosity in the Kozeny-Carman equation Rebecca Allen, Shuyu Sun King Abdullah University of Science and Technology rebecca.allen@kaust.edu.sa, shuyu.sun@kaust.edu.sa Sept 30,
More informationChapter 4: Fluid Kinematics
Overview Fluid kinematics deals with the motion of fluids without considering the forces and moments which create the motion. Items discussed in this Chapter. Material derivative and its relationship to
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 informationPHASE TRANSITIONS IN SOFT MATTER SYSTEMS
OUTLINE: Topic D. PHASE TRANSITIONS IN SOFT MATTER SYSTEMS Definition of a phase Classification of phase transitions Thermodynamics of mixing (gases, polymers, etc.) Mean-field approaches in the spirit
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 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 informationDynamics of Glaciers
Dynamics of Glaciers McCarthy Summer School 01 Andy Aschwanden Arctic Region Supercomputing Center University of Alaska Fairbanks, USA June 01 Note: This script is largely based on the Physics of Glaciers
More informationBoundary Conditions in Fluid Mechanics
Boundary Conditions in Fluid Mechanics R. Shankar Subramanian Department of Chemical and Biomolecular Engineering Clarkson University The governing equations for the velocity and pressure fields are partial
More informationComputer Applications in Engineering and Construction Programming Assignment #9 Principle Stresses and Flow Nets in Geotechnical Design
CVEN 302-501 Computer Applications in Engineering and Construction Programming Assignment #9 Principle Stresses and Flow Nets in Geotechnical Design Date distributed : 12/2/2015 Date due : 12/9/2015 at
More informationUNIVERSITY OF SASKATCHEWAN ME MECHANICS OF MATERIALS I FINAL EXAM DECEMBER 13, 2008 Professor A. Dolovich
UNIVERSITY OF SASKATCHEWAN ME 313.3 MECHANICS OF MATERIALS I FINAL EXAM DECEMBER 13, 2008 Professor A. Dolovich A CLOSED BOOK EXAMINATION TIME: 3 HOURS For Marker s Use Only LAST NAME (printed): FIRST
More informationColloidal Suspension Rheology Chapter 1 Study Questions
Colloidal Suspension Rheology Chapter 1 Study Questions 1. What forces act on a single colloidal particle suspended in a flowing fluid? Discuss the dependence of these forces on particle radius. 2. What
More informationChapter 2: Basic Governing Equations
-1 Reynolds Transport Theorem (RTT) - Continuity Equation -3 The Linear Momentum Equation -4 The First Law of Thermodynamics -5 General Equation in Conservative Form -6 General Equation in Non-Conservative
More informationLectures on Applied Reactor Technology and Nuclear Power Safety. Lecture No 6
Lectures on Nuclear Power Safety Lecture No 6 Title: Introduction to Thermal-Hydraulic Analysis of Nuclear Reactor Cores Department of Energy Technology KTH Spring 2005 Slide No 1 Outline of the Lecture
More informationAGITATION AND AERATION
AGITATION AND AERATION Although in many aerobic cultures, gas sparging provides the method for both mixing and aeration - it is important that these two aspects of fermenter design be considered separately.
More informationContinuum mechanism: Stress and strain
Continuum mechanics deals with the relation between forces (stress, σ) and deformation (strain, ε), or deformation rate (strain rate, ε). Solid materials, rigid, usually deform elastically, that is the
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 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 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 informationIntroduction to Marine Hydrodynamics
1896 1920 1987 2006 Introduction to Marine Hydrodynamics (NA235) Department of Naval Architecture and Ocean Engineering School of Naval Architecture, Ocean & Civil Engineering First Assignment The first
More informationIntroduction to Fluid Mechanics
Introduction to Fluid Mechanics Tien-Tsan Shieh April 16, 2009 What is a Fluid? The key distinction between a fluid and a solid lies in the mode of resistance to change of shape. The fluid, unlike the
More informationELASTICITY (MDM 10203)
ELASTICITY () Lecture Module 3: Fundamental Stress and Strain University Tun Hussein Onn Malaysia Normal Stress inconstant stress distribution σ= dp da P = da A dimensional Area of σ and A σ A 3 dimensional
More informationM E 320 Professor John M. Cimbala Lecture 10
M E 320 Professor John M. Cimbala Lecture 10 Today, we will: Finish our example problem rates of motion and deformation of fluid particles Discuss the Reynolds Transport Theorem (RTT) Show how the RTT
More informationStress and Strains in Soil and Rock. Hsin-yu Shan Dept. of Civil Engineering National Chiao Tung University
Stress and Strains in Soil and Rock Hsin-yu Shan Dept. of Civil Engineering National Chiao Tung University Stress and Strain ε 1 1 2 ε 2 ε Dimension 1 2 0 ε ε ε 0 1 2 ε 1 1 2 ε 2 ε Plane Strain = 0 1 2
More informationFORMULA SHEET. General formulas:
FORMULA SHEET You may use this formula sheet during the Advanced Transport Phenomena course and it should contain all formulas you need during this course. Note that the weeks are numbered from 1.1 to
More informationCHAPTER 1 Fluids and their Properties
FLUID MECHANICS Gaza CHAPTER 1 Fluids and their Properties Dr. Khalil Mahmoud ALASTAL Objectives of this Chapter: Define the nature of a fluid. Show where fluid mechanics concepts are common with those
More informationChapter 3. Load and Stress Analysis
Chapter 3 Load and Stress Analysis 2 Shear Force and Bending Moments in Beams Internal shear force V & bending moment M must ensure equilibrium Fig. 3 2 Sign Conventions for Bending and Shear Fig. 3 3
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 informationExplaining and modelling the rheology of polymeric fluids with the kinetic theory
Explaining and modelling the rheology of polymeric fluids with the kinetic theory Dmitry Shogin University of Stavanger The National IOR Centre of Norway IOR Norway 2016 Workshop April 25, 2016 Overview
More informationBasic Fluid Mechanics
Basic Fluid Mechanics Chapter 6A: Internal Incompressible Viscous Flow 4/16/2018 C6A: Internal Incompressible Viscous Flow 1 6.1 Introduction For the present chapter we will limit our study to incompressible
More informationSEMM Mechanics PhD Preliminary Exam Spring Consider a two-dimensional rigid motion, whose displacement field is given by
SEMM Mechanics PhD Preliminary Exam Spring 2014 1. Consider a two-dimensional rigid motion, whose displacement field is given by u(x) = [cos(β)x 1 + sin(β)x 2 X 1 ]e 1 + [ sin(β)x 1 + cos(β)x 2 X 2 ]e
More informationAE/ME 339. Computational Fluid Dynamics (CFD) K. M. Isaac. Momentum equation. Computational Fluid Dynamics (AE/ME 339) MAEEM Dept.
AE/ME 339 Computational Fluid Dynamics (CFD) 9//005 Topic7_NS_ F0 1 Momentum equation 9//005 Topic7_NS_ F0 1 Consider the moving fluid element model shown in Figure.b Basis is Newton s nd Law which says
More informationMECH 5312 Solid Mechanics II. Dr. Calvin M. Stewart Department of Mechanical Engineering The University of Texas at El Paso
MECH 5312 Solid Mechanics II Dr. Calvin M. Stewart Department of Mechanical Engineering The University of Texas at El Paso Table of Contents Thermodynamics Derivation Hooke s Law: Anisotropic Elasticity
More informationEKC314: TRANSPORT PHENOMENA Core Course for B.Eng.(Chemical Engineering) Semester II (2008/2009)
EKC314: TRANSPORT PHENOMENA Core Course for B.Eng.(Chemical Engineering) Semester II (2008/2009) Dr. Mohamad Hekarl Uzir-chhekarl@eng.usm.my School of Chemical Engineering Engineering Campus, Universiti
More informationConvective Mass Transfer
Convective Mass Transfer Definition of convective mass transfer: The transport of material between a boundary surface and a moving fluid or between two immiscible moving fluids separated by a mobile interface
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 informationOne dimensional steady state diffusion, with and without source. Effective transfer coefficients
One dimensional steady state diffusion, with and without source. Effective transfer coefficients 2 mars 207 For steady state situations t = 0) and if convection is not present or negligible the transport
More information5 The Oldroyd-B fluid
5 The Oldroyd-B fluid Last time we started from a microscopic dumbbell with a linear entropic spring, and derived the Oldroyd-B equations: A u = u ρ + u u = σ 2 pi + η u + u 3 + u A A u u A = τ Note that
More informationStrain Transformation equations
Strain Transformation equations R. Chandramouli Associate Dean-Research SASTRA University, Thanjavur-613 401 Joint Initiative of IITs and IISc Funded by MHRD Page 1 of 8 Table of Contents 1. Stress transformation
More informationFoundations of. Colloid Science SECOND EDITION. Robert J. Hunter. School of Chemistry University of Sydney OXPORD UNIVERSITY PRESS
Foundations of Colloid Science SECOND EDITION Robert J. Hunter School of Chemistry University of Sydney OXPORD UNIVERSITY PRESS CONTENTS 1 NATURE OF COLLOIDAL DISPERSIONS 1.1 Introduction 1 1.2 Technological
More informationPLAXIS. Scientific Manual
PLAXIS Scientific Manual 2016 Build 8122 TABLE OF CONTENTS TABLE OF CONTENTS 1 Introduction 5 2 Deformation theory 7 2.1 Basic equations of continuum deformation 7 2.2 Finite element discretisation 8 2.3
More information