Two-dimensional flow in a porous medium with general anisotropy
|
|
- Nathaniel Young
- 6 years ago
- Views:
Transcription
1 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
2 Darcy s law for flow in an isotropic porous medium µ 0= p v K p is pressure μ is dynamic viscosity K is permeability Henry Darcy (1956) formulated this linear law v is velocity where a driving force is balanced by a resistance. (vector) Analogous linear constitutive laws are: Ohm s law of electric conduction Fourier s law of heat conduction Fick s law of mass diffusion Essentially these are first-order Taylor expansions
3 Darcy s law for flow in an isotropic porous medium: Newton s. law µ 0= p v K mass x acceleration = pressure force + resistance force (taken per volume, with acceleration term neglected) Basic experiment for obtaining Darcy s law 3
4 Basic 1D experiment for obtaining Darcy s law Tube (a)-(b) filled with a porous medium. A p p A p Q = Kfall Δp=p = K Pressure µ L µ a x -pb. Tube length L. Cross-section area A. Volume flux Q. a b A pa pb A p Q= K = K µ L µ x 4 This is the observed relationship between pressure fall and volume flux. The permeability K is introduced as a proportionality constant.
5 5 Darcy s law for flow in an isotropic porous medium µ 0= p v K p is pressure μ is dynamic viscosity K is Follows by generalizing to 3Dpermeability v is velocity the empirical relationship (vector) A pa pb A p Q= K = K µ L µ x by introducing as velocity the specific flux vx=q/a. Alternatively one may work with the average pore velocity nq/a (n is porosity).
6 Darcy s law for 1D flow in an anisotropic porous medium The 1D version of Darcy s law A pa pb A p Q= K = K µ L µ x is rewritten for an anisotropic medium as Qx K xx p vx = = A µ x 6 where Kxx is the permeability for flow in the x direction as driven by a gradient in the x direction
7 7 Darcy s law for 1D flow in an anisotropic porous medium Darcy s law for an anisotropic porous medium with 1D pressure gradient p/ x=δp/lx Qx K xx p vx = = Ax µ x Qy K yx p vy = = Ay µ x K zx p Qz vz = = Az µ x vz vy x p/ x vx Kzx Kyx Kxx Kyx, Kzx are the permeabilities for flows in the y and z directions as driven by a gradient in the x direction. Given by Kyx/Kxx=vy/vx, Kzx/Kxx=vz/vx
8 Darcy s law for 3D flow in an anisotropic porous medium Darcy s law for an anisotropic porous medium with 3D pressure gradient ( p/ x, p/ y, p/ z) µ ( vx vy K11 v z ) = K 1 K 31 or, in index form 8 K1 K K 3 K13 p / x K 3 p / y K 33 p / z p µ vi = K ij xj Kij (i=1,,3) (j=1,,3) are the permeabilities for flow in the x1=x, x=y and x3=z directions as driven by gradients in the xj (j=1,,3) directions
9 3D anisotropic porous media We consider a homogeneous porous medium with general 3D anisotropy The anisotropy can be caused by fibres with preferred directions The anisotropy can be caused by periodic layering These cases may have transverse isotropy: One permeability along fibres or across layer. permeability in perpendicular direction the Another Our analysis will assume general anisotropy, with arbitrary 9 directions of principal axes
10 Flow in a 3D anisotropic porous medium Coordinate axes x,y,z. Principal axes of permeability X, Y, Z, with permeabilities KI, KII, KIII. The two coordinate systems are linked by the three Euler angles α, β, γ KIII KII KI 10
11 11 Flow in a 3D anisotropic porous medium Principal axes of permeability are X, Y, Z. The permeability matrix is KI K II 0 KIII 0 0 K III KII KI in the X, Y, Z system. The matrix is transformed to the x,y,z system
12 1 D flow in a 3D anisotropic porous medium The permeability matrix in the x,y,z system K11 K 1 K 31 K1 K K 3 K13 K 3 K 33 KIII KII KI is given by the three principal permeabilites and the Euler angles (formulas omitted) Now we restrict our attention to general D flow ( / z=0)
13 D flow in a 3D anisotropic porous medium We introduce the pressure gradient that drives a D flow (Gx,Gy) = - ( p/ x, p/ y) Note that this flow has a passive component in the z direction, even though there is no pressure gradient in the z direction Darcy s law for the 3D porous medium is written as v = K G = G K Here K is the 3D permeability 13 tensor (matrix), which is symmetric. The velocity vector v has three components (vx,vy,vz) The idea is now to introduce an effective D description for the flow in the (x,y) plane. K1 and K are defined as effective D principal permeabilities, with axes rotated an angle φ compared with the coordinate axes x and y φ y K1 K φ x
14 D flow in a 3D anisotropic porous medium An effective D description for the flow in the (x,y) plane takes K1 and K as effective D principal permeabilities, with axes rotated an angle φ Darcy s law for D flow in a 3D porous medium with (Gx,Gy) = ( p/ x,- p/ y) v = K G = G K The velocity vector v has three components (vx,vy,vz), but we are only interested in its D projection (vx,vy), which is given by the D formula φ K1 cos ϕ + K sin ϕ = (1 / )( K1 K ) sin ϕ 14 K eff K1 K (v x, v y ) = G K eff y φ (1 / )( K1 K ) sin ϕ K1 sin ϕ + K cos ϕ x
15 D flow in a 3D anisotropic porous medium The equations governing K1, K and φ are shown here. They can easily be solved numerically. K1 cos ϕ + K sin ϕ = K I a11 + K II a1 + K III a31 K1 sin ϕ + K cos ϕ = K I a1 + K II a + K III a3 (1 / )( K1 K ) sin ϕ = K I a11a1 + K II a1a + K III a31a3 φ The coefficients consist of products of of the rotation matrix 15 a11 a1 a1 a a 31 a3 yelements K a13 a3 = a33 cos γ cos α cos β sin α sin γ sin γ cos α cos β sin α cos γ sin β sin α K1 φ cos γ sin α + cos β cos α sin γ sin γ sin α + cos β cos α cos γ sin β cos α sin γ sin β cos γ sin β cos β x
16 Numerical results for KI=1, KII=, KIII=3 α β γ ϕ K1 K π /6 0 π /6 π /3 1 π / 6 π / 6 π / π / 6 π / 3 π / KIII φ KII y K1 K φ 16 KI x
17 17 Numerical results for KI=3, KII=, KIII=1 α β π /6 0 π /6 π /8 π /6 π /4 π / 6 3π / 8 γ ϕ K1 K π / 4 5π / 1 3 π / π / π / Increasing β from zero, we see how the effective permeabilities K1 and K are influenced more and more from KIII. The values of these effective permeabilities are between the extremal 3D values
18 Original motivation for this work We investigated drainage flow from a D anisotropic porous medium into a ditch. General D permeability matrix y K xx K xy K yx K yy x Our question was: Will our D drainage theory for a 18 uniformly draining ditch be valid for a class of porous media with 3D anisotropy? The answer is yes. It is valid for general 3D anisotropy, defining the effective D permeability matrix presented here
19 Summary and conclusions A D flow in a general 3D anisotropic porous medium does not involve all six independent components of the full permeability tensor (matrix). It involves five, and reduces them to three. Three independent components of an effective D permeability tensor are sufficient to describe the D flow in the plane. These three are expressed by the principal values K1, K and the tilt angle φ of the principal axes. One disadvantage with this approach is that it overlooks the passive flow perpendicular to the x,y plane of the D flow Another disadvantage is that we cannot reconstruct the full 19 permeability tensor. Even if we took the passive flow components into account, we only have 5 equations to determine the 6 independent tensor components. The single tensor component that we cannot find with D flow is Kzz
20 APPENDIX: Mathematica program for calculating the effective permeability matrix K eff K1 cos ϕ + K sin ϕ = (1 / )( K1 K ) sin ϕ 0 (1 / )( K1 K ) sin ϕ K1 sin ϕ + K cos ϕ
1871. twadaa t, 30 cta. pat Haa;fe,ttaw Spiritism. From Uis luport of tie vision, and in U e n i e h t i a d i W A C h r f i
V < > X Q x X > >! 5> V3 23 3 - - - : -- { - -- (!! - - - -! :- 4 -- : -- -5--4 X -
More informationA DARK GREY P O N T, with a Switch Tail, and a small Star on the Forehead. Any
Y Y Y X X «/ YY Y Y ««Y x ) & \ & & } # Y \#$& / Y Y X» \\ / X X X x & Y Y X «q «z \x» = q Y # % \ & [ & Z \ & { + % ) / / «q zy» / & / / / & x x X / % % ) Y x X Y $ Z % Y Y x x } / % «] «] # z» & Y X»
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 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 informationoenofc : COXT&IBCTOEU. AU skaacst sftwer thsa4 aafcekr will be ehat«s«ai Bi. C. W. JUBSSOS. PERFECT THBOUGH SDFFEBISG. our
x V - --- < x x 35 V? 3?/ -V 3 - ) - - [ Z8 - & Z - - - - - x 0-35 - 3 75 3 33 09 33 5 \ - - 300 0 ( -? 9 { - - - -- - < - V 3 < < - - Z 7 - z 3 - [ } & _ 3 < 3 ( 5 7< ( % --- /? - / 4-4 - & - % 4 V 2
More informationTwo Posts to Fill On School Board
Y Y 9 86 4 4 qz 86 x : ( ) z 7 854 Y x 4 z z x x 4 87 88 Y 5 x q x 8 Y 8 x x : 6 ; : 5 x ; 4 ( z ; ( ) ) x ; z 94 ; x 3 3 3 5 94 ; ; ; ; 3 x : 5 89 q ; ; x ; x ; ; x : ; ; ; ; ; ; 87 47% : () : / : 83
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 informationCounty Council Named for Kent
\ Y Y 8 9 69 6» > 69 ««] 6 : 8 «V z 9 8 x 9 8 8 8?? 9 V q» :: q;; 8 x () «; 8 x ( z x 9 7 ; x >«\ 8 8 ; 7 z x [ q z «z : > ; ; ; ( 76 x ; x z «7 8 z ; 89 9 z > q _ x 9 : ; 6? ; ( 9 [ ) 89 _ ;»» «; x V
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 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 informationA. H. Hall, 33, 35 &37, Lendoi
7 X x > - z Z - ----»»x - % x x» [> Q - ) < % - - 7»- -Q 9 Q # 5 - z -> Q x > z»- ~» - x " < z Q q»» > X»? Q ~ - - % % < - < - - 7 - x -X - -- 6 97 9
More information3D and Planar Constitutive Relations
3D and Planar Constitutive Relations A School on Mechanics of Fibre Reinforced Polymer Composites Knowledge Incubation for TEQIP Indian Institute of Technology Kanpur PM Mohite Department of Aerospace
More informationProblem Set 2 Due Tuesday, September 27, ; p : 0. (b) Construct a representation using five d orbitals that sit on the origin as a basis: 1
Problem Set 2 Due Tuesday, September 27, 211 Problems from Carter: Chapter 2: 2a-d,g,h,j 2.6, 2.9; Chapter 3: 1a-d,f,g 3.3, 3.6, 3.7 Additional problems: (1) Consider the D 4 point group and use a coordinate
More informationLOWELL WEEKLY JOURNAL
Y -» $ 5 Y 7 Y Y -Y- Q x Q» 75»»/ q } # ]»\ - - $ { Q» / X x»»- 3 q $ 9 ) Y q - 5 5 3 3 3 7 Q q - - Q _»»/Q Y - 9 - - - )- [ X 7» -» - )»? / /? Q Y»» # X Q» - -?» Q ) Q \ Q - - - 3? 7» -? #»»» 7 - / Q
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 informationLOWELL WEEKLY JOURNAL
Y G y G Y 87 y Y 8 Y - $ X ; ; y y q 8 y $8 $ $ $ G 8 q < 8 6 4 y 8 7 4 8 8 < < y 6 $ q - - y G y G - Y y y 8 y y y Y Y 7-7- G - y y y ) y - y y y y - - y - y 87 7-7- G G < G y G y y 6 X y G y y y 87 G
More informationLOWELL JOURNAL. MUST APOLOGIZE. such communication with the shore as Is m i Boimhle, noewwary and proper for the comfort
- 7 7 Z 8 q ) V x - X > q - < Y Y X V - z - - - - V - V - q \ - q q < -- V - - - x - - V q > x - x q - x q - x - - - 7 -» - - - - 6 q x - > - - x - - - x- - - q q - V - x - - ( Y q Y7 - >»> - x Y - ] [
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 informationi r-s THE MEMPHIS, TENN., SATURDAY. DEGfMBER
N k Q2 90 k ( < 5 q v k 3X3 0 2 3 Q :: Y? X k 3 : \ N 2 6 3 N > v N z( > > :}9 [ ( k v >63 < vq 9 > k k x k k v 6> v k XN Y k >> k < v Y X X X NN Y 2083 00 N > N Y Y N 0 \ 9>95 z {Q ]k3 Q k x k k z x X
More informationII&Ij <Md Tmlaiiiiiit, aad once in Ihe y a w Teataa m i, the vmb thatalmta oot Uiaapirit world. into as abode or wotld by them- CooTBOtioa
382 4 7 q X
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 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 informationChapter 2: Heat Conduction Equation
-1 General Relation for Fourier s Law of Heat Conduction - Heat Conduction Equation -3 Boundary Conditions and Initial Conditions -1 General Relation for Fourier s Law of Heat Conduction (1) The rate of
More informationIntroduction to Seismology Spring 2008
MIT OpenCourseWare http://ocw.mit.edu 1.510 Introduction to Seismology Spring 008 For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms. 1.510 Introduction to
More informationLOWELL WEEKLY JOURNAL
: Y J G V $ 5 V V G Y 2 25 Y 2» 5 X # VG q q q 6 6 X J 6 $3 ( 6 2 6 2 6 25 3 2 6 Y q 2 25: JJ JJ < X Q V J J Y J Q V (» Y V X Y? G # V Y J J J G J»Y ) J J / J Y Y X ({ G #? J Y ~» 9? ) < ( J VY Y J G (
More informationIn this section, mathematical description of the motion of fluid elements moving in a flow field is
Jun. 05, 015 Chapter 6. Differential Analysis of Fluid Flow 6.1 Fluid Element Kinematics In this section, mathematical description of the motion of fluid elements moving in a flow field is given. A small
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 informationGG655/CEE623 Groundwater Modeling. Aly I. El-Kadi
GG655/CEE63 Groundwater Modeling Model Theory Water Flow Aly I. El-Kadi Hydrogeology 1 Saline water in oceans = 97.% Ice caps and glaciers =.14% Groundwater = 0.61% Surface water = 0.009% Soil moisture
More informationRotational Motion. Chapter 4. P. J. Grandinetti. Sep. 1, Chem P. J. Grandinetti (Chem. 4300) Rotational Motion Sep.
Rotational Motion Chapter 4 P. J. Grandinetti Chem. 4300 Sep. 1, 2017 P. J. Grandinetti (Chem. 4300) Rotational Motion Sep. 1, 2017 1 / 76 Angular Momentum The angular momentum of a particle with respect
More information1. FUNDAMENTAL CONCEPTS AND MATH REVIEW
1. FUNDAMENTAL CONCEPTS AND MATH REVIEW 1.1. Introduction Here we provide for your reading pleasure a review of some of the math concepts used in part of this course. Most of this falls under the heading
More informationCoupled free-flow and porous media flow: a numerical and experimental investigation
Coupled free-flow and porous media flow: a numerical and experimental investigation Master s Thesis Pavan Cornelissen 3863514 Supervisors: Kilian Weishaupt, MSc prof. dr. ir. Rainer Helmig prof. dr. ir.
More informationMOMENTUM TRANSPORT Velocity Distributions in Turbulent Flow
TRANSPORT PHENOMENA MOMENTUM TRANSPORT Velocity Distributions in Turbulent Flow Introduction to Turbulent Flow 1. Comparisons of laminar and turbulent flows 2. Time-smoothed equations of change for incompressible
More informationMath Exam 2, October 14, 2008
Math 96 - Exam 2, October 4, 28 Name: Problem (5 points Find all solutions to the following system of linear equations, check your work: x + x 2 x 3 2x 2 2x 3 2 x x 2 + x 3 2 Solution Let s perform Gaussian
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 informationAn introduction to Solid State NMR and its Interactions
An introduction to Solid State NMR and its Interactions From tensor to NMR spectra CECAM Tutorial September 9 Calculation of Solid-State NMR Parameters Using the GIPAW Method Thibault Charpentier - CEA
More information' Liberty and Umou Ono and Inseparablo "
3 5? #< q 8 2 / / ) 9 ) 2 ) > < _ / ] > ) 2 ) ) 5 > x > [ < > < ) > _ ] ]? <
More informationProblem 1 (Exercise 2.2, Monograph)
MS&E 314/CME 336 Assignment 2 Conic Linear Programming January 3, 215 Prof. Yinyu Ye 6 Pages ASSIGNMENT 2 SOLUTIONS Problem 1 (Exercise 2.2, Monograph) We prove the part ii) of Theorem 2.1 (Farkas Lemma
More informationChem8028(1314) - Spin Dynamics: Spin Interactions
Chem8028(1314) - Spin Dynamics: Spin Interactions Malcolm Levitt see also IK m106 1 Nuclear spin interactions (diamagnetic materials) 2 Chemical Shift 3 Direct dipole-dipole coupling 4 J-coupling 5 Nuclear
More informationCoalbed Methane Properties
Coalbed Methane Properties Subtopics: Permeability-Pressure Relationship Coal Compressibility Matrix Shrinkage Seidle and Huitt Palmer and Mansoori Shi and Durucan Constant Exponent Permeability Incline
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 informationOWELL WEEKLY JOURNAL
Y \»< - } Y Y Y & #»»» q ] q»»»>) & - - - } ) x ( - { Y» & ( x - (» & )< - Y X - & Q Q» 3 - x Q Y 6 \Y > Y Y X 3 3-9 33 x - - / - -»- --
More informationLOWELL. MICHIGAN, OCTOBER morning for Owen J. Howard, M last Friday in Blodpett hospital.
G GG Y G 9 Y- Y 77 8 Q / x -! -} 77 - - # - - - - 0 G? x? x - - V - x - -? : : - q -8 : : - 8 - q x V - - - )?- X - - 87 X - ::! x - - -- - - x -- - - - )0 0 0 7 - - 0 q - V -
More information2. Second-order Linear Ordinary Differential Equations
Advanced Engineering Mathematics 2. Second-order Linear ODEs 1 2. Second-order Linear Ordinary Differential Equations 2.1 Homogeneous linear ODEs 2.2 Homogeneous linear ODEs with constant coefficients
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 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 informationSimplified Analytical Model of a Six-Degree-of-Freedom Large-Gap Magnetic Suspension System
NASA Technical Memorandum 112868 Simplified Analytical Model of a Six-Degree-of-Freedom Large-Gap Magnetic Suspension System Nelson J. Groom Langley Research Center, Hampton, Virginia June 1997 National
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 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 informationTransmission Lines and E. M. Waves Prof. R. K. Shevgaonkar Department of Electrical Engineering Indian Institute of Technology, Bombay
Transmission Lines and E. M. Waves Prof. R. K. Shevgaonkar Department of Electrical Engineering Indian Institute of Technology, Bombay Lecture 18 Basic Laws of Electromagnetics We saw in the earlier lecture
More informationLecture 8 Analyzing the diffusion weighted signal. Room CSB 272 this week! Please install AFNI
Lecture 8 Analyzing the diffusion weighted signal Room CSB 272 this week! Please install AFNI http://afni.nimh.nih.gov/afni/ Next lecture, DTI For this lecture, think in terms of a single voxel We re still
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 informationg(t) = f(x 1 (t),..., x n (t)).
Reading: [Simon] p. 313-333, 833-836. 0.1 The Chain Rule Partial derivatives describe how a function changes in directions parallel to the coordinate axes. Now we shall demonstrate how the partial derivatives
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 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 informationSummary of Fourier Optics
Summary of Fourier Optics Diffraction of the paraxial wave is described by Fresnel diffraction integral, u(x, y, z) = j λz dx 0 dy 0 u 0 (x 0, y 0 )e j(k/2z)[(x x 0) 2 +(y y 0 ) 2 )], Fraunhofer diffraction
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 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 informationCSE 167: Introduction to Computer Graphics Lecture #2: Linear Algebra Primer
CSE 167: Introduction to Computer Graphics Lecture #2: Linear Algebra Primer Jürgen P. Schulze, Ph.D. University of California, San Diego Fall Quarter 2016 Announcements Monday October 3: Discussion Assignment
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 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 informationPolynomials. In many problems, it is useful to write polynomials as products. For example, when solving equations: Example:
Polynomials Monomials: 10, 5x, 3x 2, x 3, 4x 2 y 6, or 5xyz 2. A monomial is a product of quantities some of which are unknown. Polynomials: 10 + 5x 3x 2 + x 3, or 4x 2 y 6 + 5xyz 2. A polynomial is a
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 informationProblem Set 2 Due Thursday, October 1, & & & & # % (b) Construct a representation using five d orbitals that sit on the origin as a basis:
Problem Set 2 Due Thursday, October 1, 29 Problems from Cotton: Chapter 4: 4.6, 4.7; Chapter 6: 6.2, 6.4, 6.5 Additional problems: (1) Consider the D 3h point group and use a coordinate system wherein
More informationJUST THE MATHS UNIT NUMBER ORDINARY DIFFERENTIAL EQUATIONS 3 (First order equations (C)) A.J.Hobson
JUST THE MATHS UNIT NUMBER 15.3 ORDINARY DIFFERENTIAL EQUATIONS 3 (First order equations (C)) by A.J.Hobson 15.3.1 Linear equations 15.3.2 Bernouilli s equation 15.3.3 Exercises 15.3.4 Answers to exercises
More informationCSE 167: Introduction to Computer Graphics Lecture #2: Linear Algebra Primer
CSE 167: Introduction to Computer Graphics Lecture #2: Linear Algebra Primer Jürgen P. Schulze, Ph.D. University of California, San Diego Spring Quarter 2016 Announcements Project 1 due next Friday at
More informationThe Sphere OPTIONAL - I Vectors and three dimensional Geometry THE SPHERE
36 THE SPHERE You must have played or seen students playing football, basketball or table tennis. Football, basketball, table tennis ball are all examples of geometrical figures which we call "spheres"
More informationModule 2: Governing Equations and Hypersonic Relations
Module 2: Governing Equations and Hypersonic Relations Lecture -2: Mass Conservation Equation 2.1 The Differential Equation for mass conservation: Let consider an infinitely small elemental control volume
More informationExercise: concepts from chapter 5
Reading: Fundamentals of Structural Geology, Ch 5 1) Study the oöids depicted in Figure 1a and 1b. Figure 1a Figure 1b Figure 1. Nearly undeformed (1a) and significantly deformed (1b) oöids with spherulitic
More informationClosed-Form Solution Of Absolute Orientation Using Unit Quaternions
Closed-Form Solution Of Absolute Orientation Using Unit Berthold K. P. Horn Department of Computer and Information Sciences November 11, 2004 Outline 1 Introduction 2 3 The Problem Given: two sets of corresponding
More informationSIMULATION AND INTERPRETATION OF BOREHOLE GEOPHYSICAL MEASUREMENTS USING hp FINTE ELEMENTS
SIMULATION AND INTERPRETATION OF BOREHOLE GEOPHYSICAL MEASUREMENTS USING hp FINTE ELEMENTS hp-fem team: D. Pardo, M. J. Nam, L. Demkowicz, C. Torres-Verdín, V. M. Calo, M. Paszynski, and P. J. Matuszyk
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 Preliminary Math Concept of Stress Stress Components Equilibrium
More informationAM 106/206: Applied Algebra Madhu Sudan 1. Lecture Notes 12
AM 106/206: Applied Algebra Madhu Sudan 1 Lecture Notes 12 October 19 2016 Reading: Gallian Chs. 27 & 28 Most of the applications of group theory to the physical sciences are through the study of the symmetry
More informationHEAT TRANSFER MATERIALS. in COMPOSITE SEIICHI NOMURA A. HAJI-SHEIKH. The University of Texas at Arlington
HEAT TRANSFER in COMPOSITE MATERIALS SEIICHI NOMURA A. HAJI-SHEIKH The University of Texas at Arlington Heat Transfer in Composite Materials D EStech Publications, Inc. 439 North Duke Street Lancaster,
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 informationSimulation study of density-driven natural convection mechanism in isotropic and anisotropic brine aquifers using a black oil reservoir simulator
Available online at www.sciencedirect.com Energy Procedia 37 (23 ) 5562 5569 GHGT- Simulation study of density-driven natural convection mechanism in isotropic and anisotropic brine aquifers using a black
More informationGG303 Lecture 6 8/27/09 1 SCALARS, VECTORS, AND TENSORS
GG303 Lecture 6 8/27/09 1 SCALARS, VECTORS, AND TENSORS I Main Topics A Why deal with tensors? B Order of scalars, vectors, and tensors C Linear transformation of scalars and vectors (and tensors) II Why
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 informationWater in Soil Sections in Craig
Water in Soil Sections 2.1-2.6 in Craig Outlines Introduction Darcy s Law Volume of water flowing per unit time Measuring K in laboratory Seepage Theory Flow Net Introduction All soils are permeable materials,
More informationStrain page 1. Strain. Paul Bons Mineralogie & Geodynamik, Eberhard Karls Universität Tübingen
page 1 Paul Bons Mineralogie & Geodynamik, Eberhard Karls Universität Tübingen Figure 1. Population of Asaphus trilobites before and after a homogeneous deformation event. The amount of strain is visualised
More informationThis report shows the capabilities of Tdyn for modelling the fluid flow through porous media.
Flow through porous media 1 Flow through porous media Introduction This report shows the capabilities of Tdyn for modelling the fluid flow through porous media. Modelling approach Modeling of flows through
More informationMA6351-TRANSFORMS AND PARTIAL DIFFERENTIAL EQUATIONS. Question Bank. Department of Mathematics FATIMA MICHAEL COLLEGE OF ENGINEERING & TECHNOLOGY
MA6351-TRANSFORMS AND PARTIAL DIFFERENTIAL EQUATIONS Question Bank Department of Mathematics FATIMA MICHAEL COLLEGE OF ENGINEERING & TECHNOLOGY MADURAI 625 020, Tamilnadu, India 1. Define PDE and Order
More informationSpin Interactions. Giuseppe Pileio 24/10/2006
Spin Interactions Giuseppe Pileio 24/10/2006 Magnetic moment µ = " I ˆ µ = " h I(I +1) " = g# h Spin interactions overview Zeeman Interaction Zeeman interaction Interaction with the static magnetic field
More informationAbility of Darcy s Law for Extension in Two- Phase Flow for Sedimentary Medium in Capillary Non-equilibrium Situations
Research Article imedpub Journals http://www.imedpub.com Resources, Recycling and Waste Management Ability of Darcy s Law for Extension in Two- Phase Flow for Sedimentary Medium in Capillary Non-equilibrium
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 informationALGEBRAIC GEOMETRY HOMEWORK 3
ALGEBRAIC GEOMETRY HOMEWORK 3 (1) Consider the curve Y 2 = X 2 (X + 1). (a) Sketch the curve. (b) Determine the singular point P on C. (c) For all lines through P, determine the intersection multiplicity
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 informationREFERENCE: CROFT & DAVISON CHAPTER 20 BLOCKS 1-3
IV ORDINARY DIFFERENTIAL EQUATIONS REFERENCE: CROFT & DAVISON CHAPTER 0 BLOCKS 1-3 INTRODUCTION AND TERMINOLOGY INTRODUCTION A differential equation (d.e.) e) is an equation involving an unknown function
More information3 Constitutive Relations: Macroscopic Properties of Matter
EECS 53 Lecture 3 c Kamal Sarabandi Fall 21 All rights reserved 3 Constitutive Relations: Macroscopic Properties of Matter As shown previously, out of the four Maxwell s equations only the Faraday s and
More informationA short review of continuum mechanics
A short review of continuum mechanics Professor Anette M. Karlsson, Department of Mechanical ngineering, UD September, 006 This is a short and arbitrary review of continuum mechanics. Most of this material
More informationNumerical Modelling in Geosciences. Lecture 6 Deformation
Numerical Modelling in Geosciences Lecture 6 Deformation Tensor Second-rank tensor stress ), strain ), strain rate ) Invariants quantities independent of the coordinate system): - First invariant trace:!!
More information6. 3D Kinematics DE2-EA 2.1: M4DE. Dr Connor Myant
DE2-EA 2.1: M4DE Dr Connor Myant 6. 3D Kinematics Comments and corrections to connor.myant@imperial.ac.uk Lecture resources may be found on Blackboard and at http://connormyant.com Contents Three-Dimensional
More informationReservoir Simulator Practical
Reservoir Simulator Practical Course Notes 2012 Philipp Lang IZR Room 403 Tel 3004 philipp.lang@unileoben.ac.at for further information please refer to the accompanying document Info Sheet & Course Logistics
More informationCOMPOSITE PLATE THEORIES
CHAPTER2 COMPOSITE PLATE THEORIES 2.1 GENERAL Analysis of composite plates is usually done based on one of the following the ries. 1. Equivalent single-layer theories a. Classical laminate theory b. Shear
More informationCHAPTER 2 THERMAL EFFECTS IN STOKES SECOND PROBLEM FOR UNSTEADY MICROPOLAR FLUID FLOW THROUGH A POROUS
CHAPTER THERMAL EFFECTS IN STOKES SECOND PROBLEM FOR UNSTEADY MICROPOLAR FLUID FLOW THROUGH A POROUS MEDIUM. Introduction The theory of micropolar fluids introduced by Eringen [34,35], deals with a class
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 informationChapter 3 Permeability
3.2 Darcy s Law In 1856, Darcy investigated the flow of water through sand filters for water purification. His experimental apparatus is shown in Figure 3.11. By empirical observation Figure 3.11 Schematic
More informationReflection of SV- Waves from the Free Surface of a. Magneto-Thermoelastic Isotropic Elastic. Half-Space under Initial Stress
Mathematica Aeterna, Vol. 4, 4, no. 8, 877-93 Reflection of SV- Waves from the Free Surface of a Magneto-Thermoelastic Isotropic Elastic Half-Space under Initial Stress Rajneesh Kakar Faculty of Engineering
More informationNotes 4: Differential Form of the Conservation Equations
Low Speed Aerodynamics Notes 4: Differential Form of the Conservation Equations Deriving Conservation Equations From the Laws of Physics Physical Laws Fluids, being matter, must obey the laws of Physics.
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 information