Flow through Variable Permeability Porous Layers
|
|
- Denis Howard
- 5 years ago
- Views:
Transcription
1 Adv. Theor. Appl. Mech., Vol. 4, 011, no. 3, Flow through Vrible Permebility Porous Lyers M. H. Hmdn nd M. T. Kmel eprtment of Mthemticl Sciences University of New Brunswick P.O. Box 5050, Sint John, New Brunswick, EL 4L5, Cnd Abstrct Plne, Poiseuille flow through permebility-strtified porous lyer is nlyzed nd Brinkmn s eqution is solved for the given configurtion. A new dimensionless number is introduced for flow considered, nd n expression is obtined for the rtio of the driving pressure grdient through the vrible permebility medium to tht of the Nvier-Stokes flow in the sme configurtion. Keywords: Brinkmn s Eqution, Permebility Strtifiction, Poiseuille Flow 1 Introduction Fluid flow through porous medi is governed by vrious models depending on the porous microstructure, presence of mcroscopic boundries, type of flow considered, nd the presence or bsence of microscopic inertil ects, [7]. Populr mong the governing models is Brinkmn s eqution, [4], which incorportes viscous sher ects into rcy s lw to fcilitte the study of flow through porous medi possessing mcroscopic boundries. Vlidity, limittions nd uses of Brinkmn s eqution hve been discussed by vrious uthors (cf. [11, 1, 14, 16, 19, 0]). In the presence of solid, mcroscopic boundry, viscous forces ply prt in viscous seprtion of the flow, [3], where viscous sher ects re dominnt in thin lyer ner the boundry while the flow is rcy-like in the medium. In ddition, sher stress ects re importnt in the study of Nvier-Stokes flow through chnnel over porous lyer. This problem hs been nlyzed by lrge number of investigtors (cf. [5, 6, 10, 11, 13, 14, 16, 17] nd the references therein). Rudrih [16] concluded tht Brinkmn s eqution is more pproprite model when the porous lyer is of finite depth, while Prvzini et. l., [14],
2 136 M. H. Hmdn nd M. T. Kmel concluded tht when Brinkmn s eqution is used, three distinct flow regimes rise, depending on rcy number, nmely: free flow regime (for rcy number greter thn unity); Brinkmn regime (for rcy number less thn 6 unity nd greter thn10 ); nd rcy regime (for rcy number less thn 6 10 ). Their investigtion, [14], emphsized tht the Brinkmn regime is trnsition zone between the free nd the rcy flows. Nield, [1], rgued tht the use of Brinkmn s viscous sher term requires redefinition of the porosity ner solid boundry, while some uthors, [9, 17], emphsized the need for vrible permebility Brinkmn model. Typicl studies in this field of flow through nd over porous lyers consider the porous medium to hve constnt porosity nd permebility. However, two situtions rise in this connection: 1) In the presence of solid boundry, there is need to redefine porosity nd permebility due to the increse in permebility nd porosity ner solid wll, thus giving rise to chnneling ects. Chnneling ects hve n impct on sher stress nd het trnsfer ner the wll, (cf. [9, 10, 15]). ) In flow over porous lyers, there exists permebility discontinuity t the interfce between the flow regions. This hs n impct on momentum nd het trnsfer cross the interfce, (cf. [5, 9, 10, 15]). The bove situtions gve rise to the need for vrible permebility models to be used in connection with Brinkmn s eqution or in connection with models of flow through porous medi which employ viscous sher term (such s the rcy-forchheimer-brinkmn model, [, 7]). Fluid flow through porous medi of vrible porosity nd permebility, nd its pplictions, hs been discussed by number of investigtors (c.f. [15, 17, 18]). Permebility is usully modeled in number of wys depending on the porous structure nd the ppliction sought. For exmple, Rees nd Pop, [15], employed the following model in their nlysis of verticl free convection in porous medi: K y / L = K + ( K K ) e ] (1) w wherein L = length scle over which permebility vries; mbient medium; nd K w is permebility t the wll. K is permebility of Alloui et.l., [1], employed n exponentil model of the following form in their nlysis of convection in binry mixtures: cy K = e () where, when c is smll (c pproches zero) then the permebility behves like the liner function K = 1+cy.
3 Flow through vrible permebility porous lyers 137 In their study of vortex instbility of mixed convection flow, Hssnien et.l., [8], employed n exponentil model of the form: K y / γ = K [1 de ] (3) + where K is permebility t the edge of boundry lyer, d is constnt nd γ is prmeter. In modeling flow through free-spce chnnel bounded by porous lyer, Nield nd Kuznetsov, [13], voided the permebility discontinuity t the interfce between the regions by introducing porous lyer of vrible permebility sndwiched between the chnnel nd the porous lyer. They modeled the permebility in such wy tht the reciprocl of the permebility vries linerly cross the lyer, nmely K 1 = KL / y (4) where L is the overll porous medium thickness, is relted to defining the thickness of the trnsition lyer, K 1 is the trnsition lyer permebility, nd K is the constnt permebility of the underlying porous lyer. While permebility vrition in the bove models is intended to reflect permebility increse ner the wll or ner the interfce between two flow regimes, in some pplictions (such s in the design of synthetic nd rtificil porous medi) it my be desired to hve low permebility (nd more compction leding to lower porosity) in regions ner solid wlls, nd higher permebility in regions closer to the center of the flow domin. At the sme time, one needs to void permebility discontinuity. Idelly, the porous medium should hve different compction ner the mcroscopic boundry so tht porosity reduction is lso continuous function. In this work, we idelize the porous lyer by introducing permebility strtifiction in the medium, nd ssuming tht the permebility increses smoothly from zero on the mcroscopic boundry to imum t center of the porous domin. We illustrte this concept of permebility strtifiction by solving Brinkmn s eqution with vrible permebility under Poiseuille flow. The volumetric flow rte is obtined for the vrible permebility medi nd compred with the volumetric flow rte in the bsence of porous mtrix in order to obtin reltionship between the driving pressure grdients. A dimensionless number ssocited with the vrible permebility model rises in the formultion, nd is discussed in this work.
4 138 M. H. Hmdn nd M. T. Kmel We envision this permebility strtifiction to be of utility in the study of flow through rtificil porous medi, including the design of rtificil bone structures, nd in the nlysis of deformble medi (such s compressed consolidted medi). Anlysis Consider the stedy flow through porous structure tht is governed by the eqution of continuity nd Brinkmn s eqution of the following forms, respectively: v r = 0 (5) r r p + v v = 0 (6) k wherein: v r is the superficil verge velocity, k is the permebility, p is the pressure, is the bse-fluid viscosity, nd is the ective viscosity (tht is, the viscosity of the fluid in the porous medium). Now, consider the unidirectionl, plne Poiseuille flow in the configurtion shown in Fig. 1. A porous medium of vrible permebility is plced between the two solid wlls, nd fluid is driven through the infinite chnnel in the x-direction by constnt pressure grdient. y Solid Wll Porous Structure 0 x Solid Wll Fig. 1. imensionl Representtive Sketch On the solid wlls, no-slip velocity condition is imposed. We ssume tht the chnnel, of depth, is filled with porous mtrix cross which permebility is continuously vrying in such wy tht permebility is strtified in the verticl direction nd flls to zero on the solid boundries. The ide here is to cst the governing eqution in dimensionless form in order to rrive t permebility definition tht overcomes permebility discontinuity ner the wlls. For the prllel flow t hnd the norml velocity component vnishes nd equtions (5) nd (6) re replced by:
5 Flow through vrible permebility porous lyers 139 d u u = C dy k( y)...(7) where 1 dp C =. (8) dx If the permebility k(y) is constnt (sy k), then eqution (8) dmits the following velocity profile tht stisfies the no-slip condition t the solid wlls: C u( y) = ( 1+ cosh βy + { coth β cosechβ} sinh βy) (9) β where β =. (10) k Now, if the permebility is vrible function of y, the solution to eqution (7) is not redily vilble nd it depends on the form tht k(y) tkes. Let us ssume tht the permebility is qudrtic, prbolic function of y, with vlue of zero t the upper nd lower wlls, nd with imum vlue of permebility ttined long the horizontl centerline of the chnnel (y = /). Now, letting k ( y) = y + by + c (11) nd using the conditions k(0)=k()=0, we obtin k( y) = ( y y). (1) Since k(y) reches imum, sy k, t y = /, we see tht k = (13) 4 nd k(y) tkes the form 4k k( y) = ( y y). (14)
6 140 M. H. Hmdn nd M. T. Kmel Eqution (7) then tkes the following form with the help of (14); d u + dy 4 k ( y u = C. (15) y) y k( y) Now, introducing the dimensionless vribles: Y = ; K( Y ) =, we k obtin the following dimensionless permebility function: K( Y ) = 4( Y Y ). (16) The chnnel becomes dimensionless of unit depth, nd eqution (15) tkes the following form: d u dy + 4 k ( Y u = C. (17) Y ) In the bsence of generl solution to (17), we ssume u to be of the following qudrtic form, with imum velocity occurring t the point of 1 imum permebility, nmely ty = : u( Y ) = A( Y Y ). (18) Using eqution (18) in (17), we obtin: A = 4 8k Ck + / C = H (19) where H H is dimensionless number defined by: = =. (0) k ϑk nd ϑ =. The velocity profile (18) then tkes the form:
7 Flow through vrible permebility porous lyers 141 C u ( Y ) = 4 ( Y Y ). (1) 8 + H Since u ttins its imum, u, t 1 Y =, eqution (1) yields: C u =. 8 + H () efining U = u, the following dimensionless form of U is obtined from u equtions (1) nd (): U( Y ) = 4( Y Y ). (3) We note tht cross the chnnel, K(Y) hs imum vlue of 1, but k(y) hs imum vlue of k. Furthermore, the inlet velocity U(Y) hs imum vlue of 1, while the dimensionl inlet velocity u(y) hs imum vlue u which depends on the chnnel height,, the viscosity coicients nd, the pressure grdient, nd the imum permebility k. In other words, u is function of the driving pressure grdient nd the dimensionless number H. 3 Volumetric Flow Rte Using eqution (1), we obtin the following volumetric flow rte cross the chnnel: 1 C Q = udy = ( ) = u H 3 (4) By comprison, in the bsence of porous mtrix the plne Poiseuille Nvier- Stokes flow is governed by d 1 v dy dp = (5) dx which dmits the following solution tht stisfies no-slip conditions on the lower nd upper chnnel wlls:
8 14 M. H. Hmdn nd M. T. Kmel dp1 v = ( Y dx Y ) (6) dp where 1 is the (constnt) driving pressure grdient in the bsence of the porous dx mtrix. The imum velocity, v, is ttined t mid-chnnel, 1 Y =, nd tkes the form dp1 v = (7) 8 dx nd the volumetric flow rte, (6) s Q NS, for the Nvier-Stokes cse is obtined from 1 dp1 Q NS = vdy = = v 0 1 dx 3. (8) Equtions (4) nd (8) yield the following rtio between the volumetric flow rtes: Q u dp / dx = =. (9) QNS v [ 1+ H /8] dp1 / dx If the flow rtes re the sme, eqution (9) provides the following expression for the pressure grdient cross the porous lyer tht is needed to mintin the sme flow rte s tht of the Nvier-Stokes flow rte: dp dx [1 + H /8] dp1 =. (30) dx 4 Comments on the H Number The H number ppering in eqution (0) is dimensionless number whose vlue depends on the chnnel depth,, the imum permebility in given porous medium, nd the rtioϑ =. The chrcteristic velocity, u, is defined in terms of H in eqution (). etermintion of H is tied to the
9 Flow through vrible permebility porous lyers 143 knowledge of the rtioϑ =. If =, then H =. Furthermore, k H 0 for k 1. >> Nele nd Nder, [11], reported some vlues for the rtio. For the ske of illustrtion, we construct Tble 1 to illustrte the reltionship between u nd H for given permebility rtio, nd the ect of H on the reltionship between the Nvier-Stokes driving pressure grdient nd the corresponding pressure grdient tht drives the flow in the permebility-strtified porous chnnel. The expected decrese in the mgnitude of the driving pressure grdient through the porous medium with decresing ective viscosity is reflected in Tble 1. Mteril Fometl 16 Aloxite 0.01 H 16k 100 k u dp dp1 / dx dx Ck k 8k Ck k k Tble 1. Vlues of H nd relted quntities bsed on reported in [11]. 5 Conclusion In the cse of pressure driven, Poiseuille flow through permebility-strtified porous lyer, n expression for the velocity profile hs been properly introduced in this work so tht it reflects the dependence of the velocity on medium nd fluid properties. The profile is ssocited with choice of smooth permebility function tht overcomes permebility discontinuity s we move closer to the solid, bounding wlls. It is recommended tht the developed profile nd permebility function (or vritions thereof) be dopted in the nlysis of Brinkmn-type flow through rtificil nd industril porous mterils. A dimensionless number, H, hs been introduced in this nlysis. This number my prove to be useful in the nlysis of generl vrible permebility medi. It my lso be of some use in estimting imum velocity in given porous medium. Some estimted vlues
10 144 M. H. Hmdn nd M. T. Kmel of H re given in Tble 1 bsed on suggested vlues of the rtio of ective viscosity to fluid velocity. References [1] Z. Alloui, R. Benncer, nd P. Vsseur, Vrible permebility ect on convection in binry mixtures sturting porous lyer, Het Mss Trnsfer, 45(009), [] M.M. Awrtni nd M.H. Hmdn, Fully developed flow through porous chnnel bounded by flt pltes, Applied Mth. Comput., 169(005), [3] R.M. Brron nd M.H. Hmdn,Viscous seprtion of flow through porous medi into two dimensionl sink, Europen Journl of Mechnics/B: Fluids, 11#6(199), [4] H.C. Brinkmn, A Clcultion of the viscous force exerted by flowing fluid on dense swrm of prticles, Appl. Scientific Res., A1(1947), [5] M. Chndesris nd. Jmet, Boundry conditions t fluid porous interfce: n priori estimtion of the stress jump boundry conditions, Int. J. Het Mss Trnsf., 50(007), [6] R.A. Ford nd M.H. Hmdn, Coupled prllel flow through composite porous lyers, Appl. Mth. Comput., 97(1998), [7] M.H. Hmdn, M.T. Kmel nd H.I. Siyym, A permebility function for Brinkmn s eqution, In Mthemticl Methods, System Theory nd Control, 009, , WSEAS Publictions. [8] I.A. Hssnien, A.A. Slm, nd A.M. Eliw, Vrible permebility ect on vortex instbility of mixed convection flow in semi-infinite porous medium bounded by horizontl surfce, Applied Mth. Comput., 146(003), [9] M. Kviny, Principles of Het Trnsfer in Porous Medi, Second Edition, Spinger, Mechnicl Engineering Series, [10] A.V. Kuznetsov, Infuence of the stress jump condition t the porousmedium/cler-fluid interfce on flow t porous wll, Int. Commun. Het Mss Trnsfer, 4(1997), [11] G. Nele, nd W. Nder, Prcticl significnce of Brinkmn s extension of rcy s lw: coupled prllel flows within chnnel nd bounding porous medium, Cndin J. Chem. Eng., 5(1974),
11 Flow through vrible permebility porous lyers 145 [1].A. Nield, The Limittions of the Brinkmn-Forchheimer eqution in modeling flow in sturted porous medium nd t n interfce, Int. J. Het Fluid Flow. 1#3(1991), [13]. A. Nield nd A. V. Kuznetsov, The ect of trnsition lyer between fluid nd porous medium: sher flow in chnnel, Trnsport Porous Medi, 78(009), [14] M. Prvzini, V. Nssehi, R. J. Wkemn, nd M. H. R. Ghoreishy, Finite element modelling of flow through porous medium between two prllel pltes using the Brinkmn eqution, Trnsport in Porous Medi. 63(006), [15].A.S. Rees nd I. Pop, Verticl free convection in porous medium with vrible permebility ects, Interntionl Journl of Het nd Mss Trnsfer 43(000), [16] N. Rudrih, Flow pst porous lyers nd their stbility, Encyclopedi of Fluid Mechnics, Slurry Flow Technology, Gulf Publishing. Chpter 14, 1986, [17] M. Shroui nd M. Kviny, Slip nd no-slip velocity boundry conditions t interfce of porous, plin medi, Int. J. Het Mss Trnsfer. 35(199), [18] K. Vfi, Convective flow nd het trnsfer in vrible porosity medi, J. Fluid Mech., 147(1984), [19] F.W. Wiegel, Fluid Flow through Porous Mcromoleculr Systems, Lecture Notes in Physics, No.11, Springer-Verlg, Berlin, Heidelberg, New York, [0] W.O. Willims, Constitutive equtions for flow of n incompressible viscous fluid through porous medium, Qurterly Appl. Mth. October, (1978), Received: June, 010
Entropy ISSN
Entropy 006, 8[], 50-6 50 Entropy ISSN 099-4300 www.mdpi.org/entropy/ ENTROPY GENERATION IN PRESSURE GRADIENT ASSISTED COUETTE FLOW WITH DIFFERENT THERMAL BOUNDARY CONDITIONS Abdul Aziz Deprtment of Mechnicl
More informationOscillatory Mixed Convention in a Porous Medium
Journl of Applied Fluid Mechnics, Vol. 7, No., pp. -5,. Avilble online t www.jfmonline.net, ISSN 75-57, EISSN 75-65. Oscilltory Mixed Convention in Porous Medium M. Jn, S. L. Mji, S. Ds, R. N. Jn nd S.
More informationFlow in porous media
Red: Ch 2. nd 2.2 PART 4 Flow in porous medi Drcy s lw Imgine point (A) in column of wter (figure below); the point hs following chrcteristics: () elevtion z (2) pressure p (3) velocity v (4) density ρ
More informationFig. 1. Open-Loop and Closed-Loop Systems with Plant Variations
ME 3600 Control ystems Chrcteristics of Open-Loop nd Closed-Loop ystems Importnt Control ystem Chrcteristics o ensitivity of system response to prmetric vritions cn be reduced o rnsient nd stedy-stte responses
More informationConsequently, the temperature must be the same at each point in the cross section at x. Let:
HW 2 Comments: L1-3. Derive the het eqution for n inhomogeneous rod where the therml coefficients used in the derivtion of the het eqution for homogeneous rod now become functions of position x in the
More informationTravelling Profile Solutions For Nonlinear Degenerate Parabolic Equation And Contour Enhancement In Image Processing
Applied Mthemtics E-Notes 8(8) - c IN 67-5 Avilble free t mirror sites of http://www.mth.nthu.edu.tw/ men/ Trvelling Profile olutions For Nonliner Degenerte Prbolic Eqution And Contour Enhncement In Imge
More informationHYDROMAGNETIC FLOW OF CASSON FLUID THROUGH A VERTICAL DEFORMABLE POROUS STRATUM WITH VISCOUS DISSIPATION AND CHEMICAL REACTION
Interntionl Journl o Mechnicl Engineering nd Technology (IJMET) Volume 9, Issue 0, October 08, pp. 846 854, Article ID: IJMET_09_0_088 Avilble online t http://www.ieme.com/ijmet/issues.sp?jtype=ijmet&vtype=9&itype=0
More informationSummary: Method of Separation of Variables
Physics 246 Electricity nd Mgnetism I, Fll 26, Lecture 22 1 Summry: Method of Seprtion of Vribles 1. Seprtion of Vribles in Crtesin Coordintes 2. Fourier Series Suggested Reding: Griffiths: Chpter 3, Section
More informationStudies on Nuclear Fuel Rod Thermal Performance
Avilble online t www.sciencedirect.com Energy Procedi 1 (1) 1 17 Studies on Nucler Fuel od herml Performnce Eskndri, M.1; Bvndi, A ; Mihndoost, A3* 1 Deprtment of Physics, Islmic Azd University, Shirz
More informationpotentials A z, F z TE z Modes We use the e j z z =0 we can simply say that the x dependence of E y (1)
3e. Introduction Lecture 3e Rectngulr wveguide So fr in rectngulr coordintes we hve delt with plne wves propgting in simple nd inhomogeneous medi. The power density of plne wve extends over ll spce. Therefore
More informationA ROTATING DISC IN CONSTANT PURE SHEAR BY S. KUMAR AND C. V. JOGA RAO
A ROTATING DISC IN CONSTANT PURE SHEAR BY S. KUMAR AND C. V. JOGA RAO (Deprtment of Aeronuticl Engineering, Indin Institute of Science, Bnglore-3) Received April 25, 1954 SUMMARY The disc of constnt pure
More informationA Permeability Function for Brinkman s Equation
Proceedings of the th WSEAS Int. Conf. on MATHEMATICAL METHOS, COMPUTATIONAL TECHNIQUES AN INTELLIGENT SYSTEMS A Permeability Function for Brinman s Equation M.H. HAMAN, M.T. KAMEL, and H.I. SIYYAM 3,
More informationAn inverse steady state thermal stresses in a thin clamped circular plate with internal heat generation
Americn Journl of Engineering Reserch (AJER) e-issn : 2320-0847 p-issn : 2320-0936 Volume-02, Issue-10, pp-276-281 www.jer.org Reserch Pper Open Access An inverse stedy stte therml stresses in thin clmped
More informationMath 124A October 04, 2011
Mth 4A October 04, 0 Viktor Grigoryn 4 Vibrtions nd het flow In this lecture we will derive the wve nd het equtions from physicl principles. These re second order constnt coefficient liner PEs, which model
More informationJackson 2.26 Homework Problem Solution Dr. Christopher S. Baird University of Massachusetts Lowell
Jckson 2.26 Homework Problem Solution Dr. Christopher S. Bird University of Msschusetts Lowell PROBLEM: The two-dimensionl region, ρ, φ β, is bounded by conducting surfces t φ =, ρ =, nd φ = β held t zero
More information* * * is the adiabatic sound speed I. INTRODUCTION II. GENERAL EQUATIONS. 76, Hamburg D-22301, Germany.
Influence of grzing flow nd dissiption effects on the coustic boundry conditions t lined wll Yves Aurégn, Rudolf Strobinski, ) nd Vincent Pgneux Lbortoire d Acoustique de l Université du Mine (I. A. M.),
More informationNew Expansion and Infinite Series
Interntionl Mthemticl Forum, Vol. 9, 204, no. 22, 06-073 HIKARI Ltd, www.m-hikri.com http://dx.doi.org/0.2988/imf.204.4502 New Expnsion nd Infinite Series Diyun Zhng College of Computer Nnjing University
More informationME 309 Fluid Mechanics Fall 2006 Solutions to Exam3. (ME309_Fa2006_soln3 Solutions to Exam 3)
Fll 6 Solutions to Exm3 (ME39_F6_soln3 Solutions to Exm 3) Fll 6. ( pts totl) Unidirectionl Flow in Tringulr Duct (A Multiple-Choice Problem) We revisit n old friend, the duct with n equilterl-tringle
More informationChapter 6 Electrostatic Boundary Value Problems. Dr. Talal Skaik
Chpter 6 Electrosttic Boundry lue Problems Dr. Tll Skik 1 1 Introduction In previous chpters, E ws determined by coulombs lw or Guss lw when chrge distribution is known, or potentil is known throughout
More informationForces from Strings Under Tension A string under tension medites force: the mgnitude of the force from section of string is the tension T nd the direc
Physics 170 Summry of Results from Lecture Kinemticl Vribles The position vector ~r(t) cn be resolved into its Crtesin components: ~r(t) =x(t)^i + y(t)^j + z(t)^k. Rtes of Chnge Velocity ~v(t) = d~r(t)=
More informationTime Dependent Slip MHD Flow Due to Stretching or Shrinking Vertical Sheet with Thermal Radiation
J. Appl. Environ. Biol. Sci., 8(1)217-224, 2018 2018, TextRod Publiction ISSN: 2090-4274 Journl of Applied Environmentl nd Biologicl Sciences www.textrod.com Time Dependent Slip MHD Flow Due to Stretching
More informationVadose Zone Hydrology
Objectives Vdose Zone Hydrology 1. Review bsic concepts nd terminology of soil physics. 2. Understnd the role of wter-tble dynmics in GW-SW interction. Drcy s lw is useful in region A. Some knowledge of
More informationCredibility Hypothesis Testing of Fuzzy Triangular Distributions
666663 Journl of Uncertin Systems Vol.9, No., pp.6-74, 5 Online t: www.jus.org.uk Credibility Hypothesis Testing of Fuzzy Tringulr Distributions S. Smpth, B. Rmy Received April 3; Revised 4 April 4 Abstrct
More informationTHE INTERVAL LATTICE BOLTZMANN METHOD FOR TRANSIENT HEAT TRANSFER IN A SILICON THIN FILM
ROMAI J., v.9, no.2(2013), 173 179 THE INTERVAL LATTICE BOLTZMANN METHOD FOR TRANSIENT HEAT TRANSFER IN A SILICON THIN FILM Alicj Piseck-Belkhyt, Ann Korczk Institute of Computtionl Mechnics nd Engineering,
More informationConservation Law. Chapter Goal. 5.2 Theory
Chpter 5 Conservtion Lw 5.1 Gol Our long term gol is to understnd how mny mthemticl models re derived. We study how certin quntity chnges with time in given region (sptil domin). We first derive the very
More informationHT Module 2 Paper solution. Module 2. Q6.Discuss Electrical analogy of combined heat conduction and convection in a composite wall.
HT Module 2 Pper solution Qulity Solutions wwwqulitytutorilin Module 2 Q6Discuss Electricl nlogy of combined het conduction nd convection in composite wll M-16-Q1(c)-5m Ans: It is frequently convient to
More information13.4 Work done by Constant Forces
13.4 Work done by Constnt Forces We will begin our discussion of the concept of work by nlyzing the motion of n object in one dimension cted on by constnt forces. Let s consider the following exmple: push
More informationSolutions to Problems in Merzbacher, Quantum Mechanics, Third Edition. Chapter 7
Solutions to Problems in Merzbcher, Quntum Mechnics, Third Edition Homer Reid April 5, 200 Chpter 7 Before strting on these problems I found it useful to review how the WKB pproimtion works in the first
More informationEuler, Ioachimescu and the trapezium rule. G.J.O. Jameson (Math. Gazette 96 (2012), )
Euler, Iochimescu nd the trpezium rule G.J.O. Jmeson (Mth. Gzette 96 (0), 36 4) The following results were estblished in recent Gzette rticle [, Theorems, 3, 4]. Given > 0 nd 0 < s
More informationTANDEM QUEUE WITH THREE MULTISERVER UNITS AND BULK SERVICE WITH ACCESSIBLE AND NON ACCESSBLE BATCH IN UNIT III WITH VACATION
Indin Journl of Mthemtics nd Mthemticl Sciences Vol. 7, No., (June ) : 9-38 TANDEM QUEUE WITH THREE MULTISERVER UNITS AND BULK SERVICE WITH ACCESSIBLE AND NON ACCESSBLE BATCH IN UNIT III WITH VACATION
More informationPart I: Basic Concepts of Thermodynamics
Prt I: Bsic Concepts o Thermodynmics Lecture 4: Kinetic Theory o Gses Kinetic Theory or rel gses 4-1 Kinetic Theory or rel gses Recll tht or rel gses: (i The volume occupied by the molecules under ordinry
More informationChapter 4 Contravariance, Covariance, and Spacetime Diagrams
Chpter 4 Contrvrince, Covrince, nd Spcetime Digrms 4. The Components of Vector in Skewed Coordintes We hve seen in Chpter 3; figure 3.9, tht in order to show inertil motion tht is consistent with the Lorentz
More informationTHERMAL EXPANSION COEFFICIENT OF WATER FOR VOLUMETRIC CALIBRATION
XX IMEKO World Congress Metrology for Green Growth September 9,, Busn, Republic of Kore THERMAL EXPANSION COEFFICIENT OF WATER FOR OLUMETRIC CALIBRATION Nieves Medin Hed of Mss Division, CEM, Spin, mnmedin@mityc.es
More informationDeriving hydraulic conductivity function from soil column tests
Deriving hydrulic conductivity function from soil column tests Yvonne Lins, Mri Dtchev & Tom Schnz DFG Fo 444 - TP 4 Experimentelle und theoretische Untersuchungen teilgesättigter Reibungsmterilien Content
More informationKINEMATICS OF RIGID BODIES
KINEMTICS OF RIGID ODIES Introduction In rigid body kinemtics, e use the reltionships governing the displcement, velocity nd ccelertion, but must lso ccount for the rottionl motion of the body. Description
More informationPlates on elastic foundation
Pltes on elstic foundtion Circulr elstic plte, xil-symmetric lod, Winkler soil (fter Timoshenko & Woinowsky-Krieger (1959) - Chpter 8) Prepred by Enzo Mrtinelli Drft version ( April 016) Introduction Winkler
More informationCalculus of Variations
Clculus of Vritions Com S 477/577 Notes) Yn-Bin Ji Dec 4, 2017 1 Introduction A functionl ssigns rel number to ech function or curve) in some clss. One might sy tht functionl is function of nother function
More informationP 3 (x) = f(0) + f (0)x + f (0) 2. x 2 + f (0) . In the problem set, you are asked to show, in general, the n th order term is a n = f (n) (0)
1 Tylor polynomils In Section 3.5, we discussed how to pproximte function f(x) round point in terms of its first derivtive f (x) evluted t, tht is using the liner pproximtion f() + f ()(x ). We clled this
More informationApplications of Bernoulli s theorem. Lecture - 7
Applictions of Bernoulli s theorem Lecture - 7 Prcticl Applictions of Bernoulli s Theorem The Bernoulli eqution cn be pplied to gret mny situtions not just the pipe flow we hve been considering up to now.
More informationLECTURE 1. Introduction. 1. Rough Classiæcation of Partial Diæerential Equations
LECTURE 1 Introduction 1. Rough Clssiction of Prtil Dierentil Equtions A prtil dierentil eqution is eqution relting function of n vribles x 1 ;::: ;x n, its prtil derivtives, nd the coordintes x =èx 1
More informationDISTRIBUTION OF SUB AND SUPER HARMONIC SOLUTION OF MATHIEU EQUATION WITHIN STABLE ZONES
Fifth ASME Interntionl Conference on Multibody Systems, Nonliner Dynmics nd Control Symposium on Dynmics nd Control of Time-Vrying nd Time-Dely Systems nd Structures September 2-2, 05, Long Bech, Cliforni,
More informationReview of Calculus, cont d
Jim Lmbers MAT 460 Fll Semester 2009-10 Lecture 3 Notes These notes correspond to Section 1.1 in the text. Review of Clculus, cont d Riemnn Sums nd the Definite Integrl There re mny cses in which some
More informationModule 2: Rate Law & Stoichiomtery (Chapter 3, Fogler)
CHE 309: Chemicl Rection Engineering Lecture-8 Module 2: Rte Lw & Stoichiomtery (Chpter 3, Fogler) Topics to be covered in tody s lecture Thermodynmics nd Kinetics Rection rtes for reversible rections
More informationFirst midterm topics Second midterm topics End of quarter topics. Math 3B Review. Steve. 18 March 2009
Mth 3B Review Steve 18 Mrch 2009 About the finl Fridy Mrch 20, 3pm-6pm, Lkretz 110 No notes, no book, no clcultor Ten questions Five review questions (Chpters 6,7,8) Five new questions (Chpters 9,10) No
More informationCOMPARISON OF DIFFERENT PROCEDURES TO PREDICT UNSATURATED SOIL SHEAR STRENGTH. S.K. Vanapalli and D.G. Fredlund 1
COMPARISON OF DIFFERENT PROCEDURES TO PREDICT UNSATURATED SOIL SHEAR STRENGTH S.K. Vnplli nd D.G. Fredlund 1 Abstrct: Severl procedures hve been proposed in the recent yers to predict the sher strength
More informationFactors affecting the phonation threshold pressure and frequency
3SC Fctors ffecting the phontion threshold pressure nd frequency Zhoyn Zhng School of Medicine, University of Cliforni Los Angeles, CA, USA My, 9 57 th ASA Meeting, Portlnd, Oregon Acknowledgment: Reserch
More information4.4 Areas, Integrals and Antiderivatives
. res, integrls nd ntiderivtives 333. Ares, Integrls nd Antiderivtives This section explores properties of functions defined s res nd exmines some connections mong res, integrls nd ntiderivtives. In order
More informationSet up Invariable Axiom of Force Equilibrium and Solve Problems about Transformation of Force and Gravitational Mass
Applied Physics Reserch; Vol. 5, No. 1; 013 ISSN 1916-9639 E-ISSN 1916-9647 Published by Cndin Center of Science nd Eduction Set up Invrible Axiom of orce Equilibrium nd Solve Problems bout Trnsformtion
More informationProblems for HW X. C. Gwinn. November 30, 2009
Problems for HW X C. Gwinn November 30, 2009 These problems will not be grded. 1 HWX Problem 1 Suppose thn n object is composed of liner dielectric mteril, with constnt reltive permittivity ɛ r. The object
More informationINTRODUCTION. The three general approaches to the solution of kinetics problems are:
INTRODUCTION According to Newton s lw, prticle will ccelerte when it is subjected to unblnced forces. Kinetics is the study of the reltions between unblnced forces nd the resulting chnges in motion. The
More informationMAC-solutions of the nonexistent solutions of mathematical physics
Proceedings of the 4th WSEAS Interntionl Conference on Finite Differences - Finite Elements - Finite Volumes - Boundry Elements MAC-solutions of the nonexistent solutions of mthemticl physics IGO NEYGEBAUE
More informationPart II. Analysis of PDE
Prt II. Anlysis of PDE Lecture notes for MA342H P. Krgeorgis pete@mths.tcd.ie 1/38 Second-order liner equtions Consider the liner opertor L(u) which is defined by L(u) = 1 u xx + 2 u yy +b 1 u x +b 2 u
More informationChapter H1: Introduction, Heat Equation
Nme Due Dte: Problems re collected on Wednesdy. Mth 3150 Problems Hbermn Chpter H1 Submitted work. Plese submit one stpled pckge per problem set. Lbel ech problem with its corresponding problem number,
More informationA Brief Note on Quasi Static Thermal Stresses In A Thin Rectangular Plate With Internal Heat Generation
Americn Journl of Engineering Reserch (AJER) 13 Americn Journl of Engineering Reserch (AJER) e-issn : 3-847 p-issn : 3-936 Volume-, Issue-1, pp-388-393 www.jer.org Reserch Pper Open Access A Brief Note
More informationREGULARITY OF NONLOCAL MINIMAL CONES IN DIMENSION 2
EGULAITY OF NONLOCAL MINIMAL CONES IN DIMENSION 2 OVIDIU SAVIN AND ENICO VALDINOCI Abstrct. We show tht the only nonlocl s-miniml cones in 2 re the trivil ones for ll s 0, 1). As consequence we obtin tht
More informationEmden-Fowler Equation and Inverse Analysis of Simple Flows through Variable Permeability Porous Layers
4 Int. J. Open Problems Compt. Math., Vol. 0, No., June 07 ISSN 998-66; Copyright ICSRS Publication, 07 www.i-csrs.org Emden-Fowler Equation and Inverse Analysis of Simple Flows through Variable Permeability
More informationMultiple Integrals. Review of Single Integrals. Planar Area. Volume of Solid of Revolution
Multiple Integrls eview of Single Integrls eding Trim 7.1 eview Appliction of Integrls: Are 7. eview Appliction of Integrls: Volumes 7.3 eview Appliction of Integrls: Lengths of Curves Assignment web pge
More informationFirst Law of Thermodynamics. Control Mass (Closed System) Conservation of Mass. Conservation of Energy
First w of hermodynmics Reding Problems 3-3-7 3-0, 3-5, 3-05 5-5- 5-8, 5-5, 5-9, 5-37, 5-0, 5-, 5-63, 5-7, 5-8, 5-09 6-6-5 6-, 6-5, 6-60, 6-80, 6-9, 6-, 6-68, 6-73 Control Mss (Closed System) In this section
More informationOn the Linear Stability of Compound Capillary Jets
ILASS Americs, th Annul Conference on Liquid Atomiztion nd Spry Systems, Chicgo, IL, My 7 On the Liner Stbility of Compound Cpillry Jets Mksud (Mx) Ismilov, Stephen D Heister School of Aeronutics nd Astronutics,
More informationThe Moving Center of Mass of a Leaking Bob
The Moving Center of Mss of Leking Bob rxiv:1002.956v1 [physics.pop-ph] 21 Feb 2010 P. Arun Deprtment of Electronics, S.G.T.B. Khls College University of Delhi, Delhi 110 007, Indi. Februry 2, 2010 Abstrct
More informationState space systems analysis (continued) Stability. A. Definitions A system is said to be Asymptotically Stable (AS) when it satisfies
Stte spce systems nlysis (continued) Stbility A. Definitions A system is sid to be Asymptoticlly Stble (AS) when it stisfies ut () = 0, t > 0 lim xt () 0. t A system is AS if nd only if the impulse response
More information1 2-D Second Order Equations: Separation of Variables
Chpter 12 PDEs in Rectngles 1 2-D Second Order Equtions: Seprtion of Vribles 1. A second order liner prtil differentil eqution in two vribles x nd y is A 2 u x + B 2 u 2 x y + C 2 u y + D u 2 x + E u +
More informationCBE 291b - Computation And Optimization For Engineers
The University of Western Ontrio Fculty of Engineering Science Deprtment of Chemicl nd Biochemicl Engineering CBE 9b - Computtion And Optimiztion For Engineers Mtlb Project Introduction Prof. A. Jutn Jn
More informationIntro to Nuclear and Particle Physics (5110)
Intro to Nucler nd Prticle Physics (5110) Feb, 009 The Nucler Mss Spectrum The Liquid Drop Model //009 1 E(MeV) n n(n-1)/ E/[ n(n-1)/] (MeV/pir) 1 C 16 O 0 Ne 4 Mg 7.7 14.44 19.17 8.48 4 5 6 6 10 15.4.41
More informationUndergraduate Research
Undergrdute Reserch A Trigonometric Simpson s Rule By Ctherine Cusimno Kirby nd Sony Stnley Biogrphicl Sketch Ctherine Cusimno Kirby is the dughter of Donn nd Sm Cusimno. Originlly from Vestvi Hills, Albm,
More informationTransient Aspects of Heat Flux Bifurcation in Porous Media: An Exact Solution
Kun Yng School of Energy nd Power Engineering, Huzhong University of Science nd Technology, Wuhn 430074, PR Chin; Deprtment of Mechnicl Engineering, University of Cliforni, Riverside, Riverside, CA 95-045
More informationVadose Zone Hydrology
Objectives Vdose Zone Hydrology. Review bsic concepts nd terminology of soil physics. 2. Understnd the role of wter-tble dynmics in GW-SW interction. Wter storge in unsturted soil Minerl surfces hve uneven
More informationTheoretische Physik 2: Elektrodynamik (Prof. A.-S. Smith) Home assignment 4
WiSe 1 8.1.1 Prof. Dr. A.-S. Smith Dipl.-Phys. Ellen Fischermeier Dipl.-Phys. Mtthis Sb m Lehrstuhl für Theoretische Physik I Deprtment für Physik Friedrich-Alexnder-Universität Erlngen-Nürnberg Theoretische
More informationThe First Fundamental Theorem of Calculus. If f(x) is continuous on [a, b] and F (x) is any antiderivative. f(x) dx = F (b) F (a).
The Fundmentl Theorems of Clculus Mth 4, Section 0, Spring 009 We now know enough bout definite integrls to give precise formultions of the Fundmentl Theorems of Clculus. We will lso look t some bsic emples
More informationSCHEME OF WORK FOR IB MATHS STANDARD LEVEL
Snnrpsgymnsiet Lott Hydén Mthemtics, Stndrd Level Curriculum SCHEME OF WORK FOR IB MATHS STANDARD LEVEL Min resource: Mthemtics for the interntionl student, Mthemtics SL, Hese PART 1 Sequences nd Series
More informationIntroduction to Finite Element Method
Introduction to Finite Element Method Introductory Course on Multiphysics Modelling TOMASZ G. ZIELIŃSKI bluebox.ippt.pn.pl/ tzielins/ Tble of Contents 1 Introduction 1 1.1 Motivtion nd generl concepts.............
More informationEffects of dry density on soil water characteristic curve of clay
5th Interntionl Conference on Civil, Architecturl nd Hydrulic Engineering (ICCAHE 2016) Effects of dry density on soil wter chrcteristic curve of cly Hu Mengling, byo Hilin, cren Jinxi School of Architecture
More informationDefinition of Continuity: The function f(x) is continuous at x = a if f(a) exists and lim
Mth 9 Course Summry/Study Guide Fll, 2005 [1] Limits Definition of Limit: We sy tht L is the limit of f(x) s x pproches if f(x) gets closer nd closer to L s x gets closer nd closer to. We write lim f(x)
More informationu t = k 2 u x 2 (1) a n sin nπx sin 2 L e k(nπ/l) t f(x) = sin nπx f(x) sin nπx dx (6) 2 L f(x 0 ) sin nπx 0 2 L sin nπx 0 nπx
Chpter 9: Green s functions for time-independent problems Introductory emples One-dimensionl het eqution Consider the one-dimensionl het eqution with boundry conditions nd initil condition We lredy know
More informationThe Wave Equation I. MA 436 Kurt Bryan
1 Introduction The Wve Eqution I MA 436 Kurt Bryn Consider string stretching long the x xis, of indeterminte (or even infinite!) length. We wnt to derive n eqution which models the motion of the string
More informationBest Approximation. Chapter The General Case
Chpter 4 Best Approximtion 4.1 The Generl Cse In the previous chpter, we hve seen how n interpolting polynomil cn be used s n pproximtion to given function. We now wnt to find the best pproximtion to given
More informationTerminal Velocity and Raindrop Growth
Terminl Velocity nd Rindrop Growth Terminl velocity for rindrop represents blnce in which weight mss times grvity is equl to drg force. F 3 π3 ρ L g in which is drop rdius, g is grvittionl ccelertion,
More informationMath 113 Exam 2 Practice
Mth 3 Exm Prctice Februry 8, 03 Exm will cover 7.4, 7.5, 7.7, 7.8, 8.-3 nd 8.5. Plese note tht integrtion skills lerned in erlier sections will still be needed for the mteril in 7.5, 7.8 nd chpter 8. This
More informationContinuous Random Variables
STAT/MATH 395 A - PROBABILITY II UW Winter Qurter 217 Néhémy Lim Continuous Rndom Vribles Nottion. The indictor function of set S is rel-vlued function defined by : { 1 if x S 1 S (x) if x S Suppose tht
More informationNumerical Integration
Chpter 5 Numericl Integrtion Numericl integrtion is the study of how the numericl vlue of n integrl cn be found. Methods of function pproximtion discussed in Chpter??, i.e., function pproximtion vi the
More informationIN GAUSSIAN INTEGERS X 3 + Y 3 = Z 3 HAS ONLY TRIVIAL SOLUTIONS A NEW APPROACH
INTEGERS: ELECTRONIC JOURNAL OF COMBINATORIAL NUMBER THEORY 8 (2008), #A2 IN GAUSSIAN INTEGERS X + Y = Z HAS ONLY TRIVIAL SOLUTIONS A NEW APPROACH Elis Lmpkis Lmpropoulou (Term), Kiprissi, T.K: 24500,
More informationA REVIEW OF CALCULUS CONCEPTS FOR JDEP 384H. Thomas Shores Department of Mathematics University of Nebraska Spring 2007
A REVIEW OF CALCULUS CONCEPTS FOR JDEP 384H Thoms Shores Deprtment of Mthemtics University of Nebrsk Spring 2007 Contents Rtes of Chnge nd Derivtives 1 Dierentils 4 Are nd Integrls 5 Multivrite Clculus
More informationu( t) + K 2 ( ) = 1 t > 0 Analyzing Damped Oscillations Problem (Meador, example 2-18, pp 44-48): Determine the equation of the following graph.
nlyzing Dmped Oscilltions Prolem (Medor, exmple 2-18, pp 44-48): Determine the eqution of the following grph. The eqution is ssumed to e of the following form f ( t) = K 1 u( t) + K 2 e!"t sin (#t + $
More informationA Model of two mutually interacting Species with Mortality Rate for the Second Species
Avilble online t www.pelgireserchlibrry.com Advnces in Applied Science Reserch, 0, 3 ():757-764 ISS: 0976-860 CODE (USA): AASRFC A Model of two mutully intercting Species with Mortlity Rte for the Second
More informationEnergy Bands Energy Bands and Band Gap. Phys463.nb Phenomenon
Phys463.nb 49 7 Energy Bnds Ref: textbook, Chpter 7 Q: Why re there insultors nd conductors? Q: Wht will hppen when n electron moves in crystl? In the previous chpter, we discussed free electron gses,
More informationMA Handout 2: Notation and Background Concepts from Analysis
MA350059 Hndout 2: Nottion nd Bckground Concepts from Anlysis This hndout summrises some nottion we will use nd lso gives recp of some concepts from other units (MA20023: PDEs nd CM, MA20218: Anlysis 2A,
More informationSimulated Performance of Packed Bed Solar Energy Storage System having Storage Material Elements of Large Size - Part I
The Open Fuels & Energy Science Journl, 2008, 1, 91-96 91 Open Access Simulted Performnce of Pcked Bed Solr Energy Storge System hving Storge Mteril Elements of Lrge Size - Prt I Rnjit Singh *,1, R.P.
More informationKirchhoff and Mindlin Plates
Kirchhoff nd Mindlin Pltes A plte significntly longer in two directions compred with the third, nd it crries lod perpendiculr to tht plne. The theory for pltes cn be regrded s n extension of bem theory,
More informationMath 1B, lecture 4: Error bounds for numerical methods
Mth B, lecture 4: Error bounds for numericl methods Nthn Pflueger 4 September 0 Introduction The five numericl methods descried in the previous lecture ll operte by the sme principle: they pproximte the
More informationSolutions to Supplementary Problems
Solutions to Supplementry Problems Chpter 8 Solution 8.1 Step 1: Clculte the line of ction ( x ) of the totl weight ( W ).67 m W = 5 kn W 1 = 16 kn 3.5 m m W 3 = 144 kn Q 4m Figure 8.10 Tking moments bout
More informationScienceDirect. Thermal convection measurements inside aluminum foam and comparison to existing analytical solutions
Avilble online t www.sciencedirect.com ScienceDirect Procedi Mterils Science 4 (14 ) 341 346 8th Interntionl Conference on Porous Metls nd Metllic Foms, Metfom 13 Therml convection mesurements inside luminum
More informationMath Lecture 23
Mth 8 - Lecture 3 Dyln Zwick Fll 3 In our lst lecture we delt with solutions to the system: x = Ax where A is n n n mtrix with n distinct eigenvlues. As promised, tody we will del with the question of
More informationOrdinary differential equations
Ordinry differentil equtions Introduction to Synthetic Biology E Nvrro A Montgud P Fernndez de Cordob JF Urchueguí Overview Introduction-Modelling Bsic concepts to understnd n ODE. Description nd properties
More information1 Line Integrals in Plane.
MA213 thye Brief Notes on hpter 16. 1 Line Integrls in Plne. 1.1 Introduction. 1.1.1 urves. A piece of smooth curve is ssumed to be given by vector vlued position function P (t) (or r(t) ) s the prmeter
More informationChapter 1. Basic Concepts
Socrtes Dilecticl Process: The Þrst step is the seprtion of subject into its elements. After this, by deþning nd discovering more bout its prts, one better comprehends the entire subject Socrtes (469-399)
More informationQuestion 1: Figure 1: Schematic
Question : θ Figure : Schemtic Consider chnnel of height with rectngulr cross section s shown in the sketch. A hinged plnk of length L < nd t n ngle θ is locted t the center of the chnnel. You my ssume
More informationLecture 19: Continuous Least Squares Approximation
Lecture 19: Continuous Lest Squres Approximtion 33 Continuous lest squres pproximtion We begn 31 with the problem of pproximting some f C[, b] with polynomil p P n t the discrete points x, x 1,, x m for
More informationChapter 1. Chapter 1 1
Chpter Chpter : Signls nd Systems... 2. Introduction... 2.2 Signls... 3.2. Smpling... 4.2.2 Periodic Signls... 0.2.3 Discrete-Time Sinusoidl Signls... 2.2.4 Rel Exponentil Signls... 5.2.5 Complex Exponentil
More informationEquations of Motion. Figure 1.1.1: a differential element under the action of surface and body forces
Equtions of Motion In Prt I, lnce of forces nd moments cting on n component ws enforced in order to ensure tht the component ws in equilirium. Here, llownce is mde for stresses which vr continuousl throughout
More informationADVANCEMENT OF THE CLOSELY COUPLED PROBES POTENTIAL DROP TECHNIQUE FOR NDE OF SURFACE CRACKS
ADVANCEMENT OF THE CLOSELY COUPLED PROBES POTENTIAL DROP TECHNIQUE FOR NDE OF SURFACE CRACKS F. Tkeo 1 nd M. Sk 1 Hchinohe Ntionl College of Technology, Hchinohe, Jpn; Tohoku University, Sendi, Jpn Abstrct:
More information