DEVELOPMENT OF A MULTI-FEED P-T WELLBORE MODEL FOR GEOTHERMAL WELLS
|
|
- Silvester Morrison
- 5 years ago
- Views:
Transcription
1 PROCEEDINGS, Thirty-First Workshop on Geothermal Reservoir Engineering Stanford University, Stanford, California, January 3-February 1, 6 SGP-TR-179 DEVELOPMENT OF MULTI-FEED P-T WELLBORE MODEL FOR GEOTHERML WELLS Murat Cinar, Mustafa Onur and bdurrahman Satman Petroleum and Natural Gas Engineering Department Istanbul Tehnial University Maslak Istanbul TURKEY s: inarmura@itu.edu.tr, onur@itu.edu.tr, mdsatman@itu.edu.tr BSTRCT wellbore pressure temperature (p-t) model for geothermal wells with multiple feed zones is presented. The flow rate, temperature and pressure may vary depending on the feed zones. The fluid onsidered is water ontaining low amount of CO. Espeially in Turkey, geothermal waters ontain onsiderable amounts of CO dissolved in water. The presene of CO leads to two phase flow at higher pressure onditions than the vapor pressure of water. The wellbore model onsiders all these onditions. The model assumes that the flow is in steady state and Henry s law is valid for CO s partial pressure alulations. The main inputs of the model are reservoir properties and the depths of the feed points. The main outputs are pressure and temperature profiles. If measured profiles of temperature and pressure are available then the model an be used to estimate the mass fration of CO. Dynami temperature profiles measured at some wells in fyon Omer-Geek geothermal field in Turkey showed multi-feed wellbore behavior. The multi-feed p-t wellbore model developed in this study was applied to analyze the p-t profiles of one of the wells, the well F-1. There are two feed points in the wellbore, one of whih is ooler. The field data from the F-1 well was analyzed with the model and a good math between the alulated and measured data was obtained. lso the CO ontent was estimated from the math. This paper disusses the analysis proedure as well as the results obtained from the modeling approah. INTRODUCTION Prediting flowing pressure and temperature profiles in geothermal wells is a very important issue to exploit geothermal resoure in a rational way. High temperatures and flow rates are ommonly enountered in geothermal wells. These onditions inrease the errors in the measurements. In addition, the presene of CO, frequently found in geothermal fluids auses sale deposition of alium arbonate in the wellbore where the fluid is flashing. In order to take the neessary preautions, the depth of flashing has to be known (Satman and lkan, 1989). Compliated geologial strutures of geothermal fields make wellbore flow a muh more ompliated problem. For example, feeds with different temperatures at different depths may be enountered. In the following work, a pressure-temperature (p-t) model for geothermal wells with multiple feed points is presented. Based on the assumptions made, the model disussed is valid for shallow wells of water dominated systems ontaining low amount of dissolved CO. DESCRIPTION OF THE MODEL Many of the geothermal systems are water dominated systems produing a two-phase mixture of liquid and gas. In Turkey the majority of the geothermal systems is water-dominated systems inluding CO dissolved in water. These systems produe a two phase mixture, vapor phase of whih is formed mostly by CO, at the wellbore. The alulation of pressure and temperature profiles for those systems has to onsider the following effets: 1. Flashing point of CO and mass transfer between phases.. Potential and fritional pressure losses for both single and two phase onditions 3. Heat loss to the surrounding formations. SINGLE PHSE FLOW In this setion, the equations for pressure and temperature drop for flow of single-phase water are presented.
2 Pressure Profile Calulations General pressure drop equation for fluid flow is given by Equation 1. dp dp dp dp = + + dz dz dz dz F E (1) Here the first term on the right hand side represents pressure drop due to frition, seond term represents pressure drop due to elevation hanges, and the last term represents pressure drop due to aeleration. Pressure loss due to aeleration is muh smaller than pressure loss due to elevation hange and frition and so an be ignored. Therefore Equation 1 beomes, dp dp dp = + dz dz dz Here, F dp f = dz F g d υ ρ dp ρg sinθ = dz E g ssuming the well is vertial (sinθ =1), dp fυ ρ ρg = dz g d g E () (3) (4) (5) By using Equation 5, pressure drop for single phase water flow an be alulated. Temperature Profile Calulations For prediting temperature profile in the wellbore in single phase onditions, Ramey s lassial method is used (Ramey, 1961). Ramey s method for alulating temperature distribution in the wellbore is based on the following two simplifying assumptions: 1. Heat flows radially away from the wellbore. Heat flow through various thermal resistanes in the immediate viinity of the wellbore is rapid ompared to heat flow in the formation and an be represented by steady state solution. Ramey s approah was investigated by many researhers and validity of the model was examined (Hagoort, 4; lves et al., 199; Wu and Pruess, 199). In Hagoort s reent work, it is observed that Ramey s solution is an exellent approximation exept for an early transient period. For geothermal wells the equation for temperature distribution along the well is given by, y ( ) ( 1 ) ( ) T = T ay + a e + T T e (6) bh / y/ bh Here is the diffusion depth and for geothermal wells is defined as, () t () = wf t π k (7) In Equation 7, f(t) is the dimensionless time funtion representing the transient heat transfer to the formation. In literature there are different orrelations to alulate this term (Ramey, 1961; Hagoort,4; Hasan and Kabir, 1994). The orrelation proposed here is found by urve fitting the data given by Ramey et al. (see Ramey Table 6.1, 1961). The dimensionless time funtion an be defined as, () = ( D ) x Log( t ) x Log( t ) Log f t x Log t D D 4 for 1 td 1 (8) where, αt td = r Log( f(t) ) ( ) ' garwal - Ramey Hagoort Proposed (Eq.8) (9) - 4 Log(t D ) Figure. 1. Comparison of different orrelations in the literature for f(t).
3 The dimensionless time funtion may be found from Figure 1. By using Equations 6, 7, 8 and 9 temperature distribution along the well for the flow of single-phase water an be estimated. MULTI FEED WELLBORE MODEL Figure.. Representation of the multiple feed (Tokita and Itoi, 4). Representation of a multi-feed ase is shown in Figure. In the model it is assumed that the flow from the deepest feed zone ours first and forms a pressure gradient in the well. That pressure gradient stimulates the other shallower feed zones to flow onseutively. Here flow rate will inrease after eah feed zone, and the temperature of the fluid flowing to the wellbore may be different from the fluid flowing in the well. ssuming that thermal equilibrium is established instantaneously, new temperature at whih fluids are mixed at the feed point an be alulated by the following formula. Qtotal = Qt + Q f (1) T mix = q ρt + q ρ T t t t f f f ( qf + qt) ρ f ( Tmix) (11) Here subsript f stands for the feed properties, and t stands for the properties of the fluid in the tubing. TWO PHSE FLOW The fluid flowing in the well forms a pressure gradient. Pressure dereases while the fluid rises in the well. s a result of pressure drop the dissolved CO in the fluid leaves the water and forms a seond phase. Developing the model it is assumed that the gas phase is only omposed of CO. Furthermore, Henry s law is assumed to be valid for partial pressure alulations of CO. Under those assumptions, the following method is used for alulating the two-phase pressure and temperature profiles. Pressure Profile Calulations Two phase alulations are ompliated by the fat that two phases exist simultaneously. The interfae between phases an be formed in many different patterns whih is alled flow pattern. To make aurate alulations, flow pattern and the riteria for flow pattern transitions have to be determined orretly. In the model it is assumed that flow is in the bubble flow regime. The harateristi of this flow pattern is the flowing of gas in the liquid as disrete bubbles. In two-phase flow, phases flow with different veloities. The gas phase, having lower density relatively, flows faster than the liquid phase. s a result the faster flowing gas phase holds up the liquid phase. This phenomenon is alled liquid hold up and the differene in the veloity is alled slip (Hasan and Kabir, ). Using homogenous flow equation for pressure drop alulations in bubble flow redues omplexity without introduing any signifiant inauray (Hasan and Kabir, ). dp f g dv = dz g d g g dz mυmρm ρm υmρm m (1) Fluid aeleration during bubble flow is very small so the last term in Equation 1 is negleted. Hene the equation beomes, dp fmυmρm ρmg = dz g d g (13) The density of the mixture an be alulated by the following equations, ( 1 ) ρ = ρ f + ρ f (14) f m L L g L L vsg = 1 C υ + m υ ( L g) g ρ ρ σ υ = 1.53 ρ L 1/ 4 (15) (16) The flashing point pressure for a given temperature an be alulated by the help of Henry s law that is given by,
4 pi = KH ni (17) For CO, 18 p = x K 44 p = p mpa CO CO HW CO CO where, ( ) ( ) ( ) ln x1 xt x1 xt T = T C KHW = + x xt The flashing point pressure an be estimated by, ( 1 ) fp CO CO s (18) (19) p = p + n p () For the vapor pressure alulations of water ntoine equation is used (Ohe, 1976). Temperature Profile Calulations Temperature profile alulations are governed by the equation proposed by lves (lves et al., 199) whih redues to Ramey's method under appropriate assumptions and is given by, z ( ) ( 1 ) ( ) T = T az + a e + T T e bh 1 g + ηρ g e Jg gjp / z/ bh z/ ( 1 ) (6) For the single-phase onditions the above equation redues to Equation 6. For flow rates higher than 5 lb m /se, ρηg term an be taken as (Sagar et al., 1991). Hene the equation redues to the following one: L/ ( bh ) ( 1 ) L / g L/ ( T Tbh ) e ( 1 e ) T = T al + a e + gj P (7) By using Equation 7 temperature distribution along the well for two phase onditions an be estimated. P-T MODEL Two governing equations (temperature and pressure equations) are oupled through the dependene of the fluid density on temperature for single-phase and on both temperature and pressure for two-phase. z/ ( bh θ) θ( 1 ) T = T azsin + asin e + Here, 1 ρ J z / z/ ( T Tbh ) e + Φ( 1 e ) P dp dz (1) The wellbore depth profile is broken up into a number of segments and the pressure and temperature equations applied suessfully to eah segment. The first segment is hosen as the bottom of the well and from that referene point the alulations are performed up to the wellhead for eah segment. dp g d dp ρ p Sin ρυ υ Φ= ρη θ dz g g dz dz 1 T z η = yg yl p z + T p () (3) Pressure loss due to elevation hange is the most important omponent of pressure drop equation of vertial flow and pressure losses due to aeleration and frition are muh more smaller ompared to pressure loss due to elevation hange so we an take, dp g Sin dz g ρ θ (4) ρυ dυ g dz (5) By assuming the well is vertial Equation 1 is redued to the following form, MODEL VLIDTION Dynami temperature profiles measured at some wells in the fyon Omer-Geek geothermal field in Turkey showed multi-feed wellbore behavior. The multi-feed p-t wellbore model developed in this study was applied to analyze the p-t profiles of one of the wells, the well F-1. There are two feed points in the wellbore, one of whih is ooler. The F-1 well is 695 ft deep. ording to the measurements the reservoir water ontains.4 % CO by mass. The field data from the F-1 well was analyzed using the model developed in this study and a good math between the alulated and measured data was obtained. By mathing the field data with the alulated ones, the properties of the feed zones (temperature and flow rate) and CO ontent were estimated. The temperatures of feed zones are found by the analysis of a stati p-t test onduted in the well.
5 Figure 3 shows the stati and dynami temperature profiles. In the stati temperature profile there is a sharp transition between 5 F and 13.8 F whih is an indiation of the upper and ooler feed point. Here the temperature of the seond feed is found to be 13.8 F. To alulate the flow rate of the seond feed Equation 11 is used. Hene the only unknown in Equation 11 is the flow rate of the feed and is estimated to be.53 ft 3 /se. Temperature, F Calulated values Measured values The CO ontent was found to be.4 % from the math of alulated and dynami p-t test data. lso the flashing point was found to be approximately 38 ft. The results are presented in Figures 4 and 5. Depth, ft 4 CONCLUSIONS In this paper, a wellbore p-t model for geothermal wells with multiple feed points is presented. The model has the ability to alulate pressure and temperature profiles along the well for both single and two phase onditions in the presene of CO. Model an be used to alulate CO ontent by using field data. In addition a orrelation for prediting dimensionless time funtion is proposed. good math between alulated and measured data is obtained. 6 Figure. 4. Comparison of alulated and measured temperature data. Pressure, psia 1 3 Temperature, F Calulated values Measured values Depth, ft 4 Depth, ft Stati Temperature Profile Dynami Temperature Profile 8 Figure. 5. Comparison of alulated and measured pressure data. Figure. 3. Stati and dynami temperature profiles. NOMECLTURE a = geothermal gradient, F / ft = diffusion depth, ft p = fluid heat apaity (onstant pressure), Btu / (lb m - F) d = diameter, ft f = frition fator f g = void fration f L = liquid hold up
6 f(t) = dimensionless time funtion g = aeleration of gravity, ft / se g = unit onversion fator,3.17(lb m - ft) / lb f / se J = mehanial equivalent of heat, 778 ft-lb fore / BTU k = thermal ondutivity of formation, 33.6 Btu / (ft-day- F) K H = Henry s oeffiient. n = mole fration p = pressure, psia p fp = flashing point pressure p s = vapor pressure q = flow rate, ft 3 / se Q = Heat, Btu / lb m r = outside radius of asing, ft t = time, days t D = Fourier time, dimensionless T = temperature, F T bh = bottom hole temperature T = temperature of fluid in tubing w = mass flow rate, lb m / hr x = mass fration. y = no slip hold up. z = Length, ft. z = gas ompressibility fator α = thermal diffusivity of formation,.96 ft /gün ρ = density, lb m / ft 3 υ = veloity of the fluid, ft / se υ sg = superfiial in-situ veloity, ft / se υ = terminal rise veloity, ft / se θ = inlination angle with horizontal Φ = dimensionless orretion parameter η = Joule Thomson oeffiient σ = surfae tension, lb m / se SUBSCRIPTS m = mixture g = gaseous L = liquid W = water Hasan,.R., Kabir, C.S. (1994), spets of Wellbore Heat Transfer During Two-Phase Flow, SPE Prodution & Failities, ugust, Hasan,.R., Kabir, C.S. (), Fluid Flow and Heat Transfer in Wellbores, Soiety of Petroleum Engineers, Texas. Ohe, S. (1976), Computer ided Data Book of Vapor Pressure, Data Book Publishing Co., Tokyo. Ramey, H.J., Jr. (1981) Reservoir Engineering ssessment of Geothermal Systems, Stanford University, Stanford, C. Ramey, H.J., Jr. (1961) Wellbore Heat Transmission, IME Transations, Sagar, R.K., Dotty, D.R., and Shmidt, Z. (1991) Prediting Temperatures in a Flowing Oil Well, SPE Prodution Engineering, 1991, Satman,., lkan, H. (1989) Modeling of Wellbore Flow and Calite Deposition for Geothermal Wells in the Presene of CO, SPE Paper No: Tokita, H., Itoi, R., (4) Development of the MULFEWS Multi-Feed Wellbore Simulator, 9th Workshop on Geothermal Reservoir Engineering, Stanford University, Stanford, California, 6-8 January. Yu-Shu Wu and Pruess, K. (199) n nalytial Solution for Wellbore Heat Transmission in Layered Formations, SPE Reservoir Engineering, November, SUPERSCRIPTS = average REFERENCES lves, I.N., lhanati, F.J., Shoham, O. (199), Unified Model for Prediting Flowing Temperature Distribution in Wellbores and Pipelines, SPE Prodution Engineering, November, Hagoort, J. (4), Ramey s Wellbore Heat Transmission Revisited, SPE Journal, Deember,
Millennium Relativity Acceleration Composition. The Relativistic Relationship between Acceleration and Uniform Motion
Millennium Relativity Aeleration Composition he Relativisti Relationship between Aeleration and niform Motion Copyright 003 Joseph A. Rybzyk Abstrat he relativisti priniples developed throughout the six
More informationProcess engineers are often faced with the task of
Fluids and Solids Handling Eliminate Iteration from Flow Problems John D. Barry Middough, In. This artile introdues a novel approah to solving flow and pipe-sizing problems based on two new dimensionless
More informationHomework Set 4. gas B open end
Homework Set 4 (1). A steady-state Arnold ell is used to determine the diffusivity of toluene (speies A) in air (speies B) at 298 K and 1 atm. If the diffusivity is DAB = 0.0844 m 2 /s = 8.44 x 10-6 m
More informationNUMERICAL SIMULATION OF ATOMIZATION WITH ADAPTIVE JET REFINEMENT
Paper ID ILASS8--7 ILASS 28 Sep. 8-, 28, Como Lake, Italy A44 NUMERICAL SIMULATION OF ATOMIZATION WITH ADAPTIVE JET REFINEMENT Anne Bagué, Daniel Fuster, Stéphane Popinet + & Stéphane Zaleski Université
More informationSOME FUNDAMENTAL ASPECTS OF COMPRESSIBLE FLOW
SOE FUNDAENAL ASECS OF CORESSIBLE FLOW ah number gas veloity mah number, speed of sound a a R < : subsoni : transoni > : supersoni >> : hypersoni art three : ah Number 7 Isentropi flow in a streamtube
More informationCALCULATION AND USE OF STEAM/WATER RELATIVE PERMEABILITIES IN GEOTHERMAL RESERVOIRS KIYOSHI SHINOHAR4. June 1978 A WSTER OF SCIENCE REPORT
CACUATION AND USE OF STEAM/WATER REATIVE PERMEABIITIES IN GEOTHERMA RESERVOIRS i A WSTER OF SCIENCE REPORT BY KIYOSHI SHINOHAR4 i I June 1978 ABSTRACT A new method to alulate the steam/water relative permeabilities
More informationDeveloping Excel Macros for Solving Heat Diffusion Problems
Session 50 Developing Exel Maros for Solving Heat Diffusion Problems N. N. Sarker and M. A. Ketkar Department of Engineering Tehnology Prairie View A&M University Prairie View, TX 77446 Abstrat This paper
More informationAnalysis of discretization in the direct simulation Monte Carlo
PHYSICS OF FLUIDS VOLUME 1, UMBER 1 OCTOBER Analysis of disretization in the diret simulation Monte Carlo iolas G. Hadjionstantinou a) Department of Mehanial Engineering, Massahusetts Institute of Tehnology,
More informationWave Propagation through Random Media
Chapter 3. Wave Propagation through Random Media 3. Charateristis of Wave Behavior Sound propagation through random media is the entral part of this investigation. This hapter presents a frame of referene
More informationNatural Convection Experiment Measurements from a Vertical Surface
OBJECTIVE Natural Convetion Experiment Measurements from a Vertial Surfae 1. To demonstrate te basi priniples of natural onvetion eat transfer inluding determination of te onvetive eat transfer oeffiient.
More informationCavity flow with surface tension past a flat plate
Proeedings of the 7 th International Symposium on Cavitation CAV9 Paper No. ## August 7-, 9, Ann Arbor, Mihigan, USA Cavity flow with surfae tension past a flat plate Yuriy Savhenko Institute of Hydromehanis
More informationBINARY RANKINE CYCLE OPTIMIZATION Golub, M., Koscak-Kolin, S., Kurevija, T.
BINARY RANKINE CYCLE OPTIMIZATION Golub, M., Kosak-Kolin, S., Kurevija, T. Faulty of Mining, Geology and Petroleum Engineering Department of Petroleum Engineering Pierottijeva 6, Zagreb 0 000, Croatia
More informationA Cubic Equation of State for Reservoir Fluids
6th IASME/WSEAS International Conferene on HEAT TRANSFER, THERMAL ENGINEERING and ENVIRONMENT (HTE'08) Rhodes, Greee, August 0-, 008 A Cubi Equation of State for Reservoir Fluids AMIR AHMAD SHIRAZIMANESH,
More informationIn this problem, we are given the following quantities: We want to find: Equations and basic calculations:
.1 It takes. million tons of oal per year to a 1000-W power plant that operates at a apaity fator of 70%. If the heating value of the oal is 1,000 /lb, alulate the plant s effiieny and the heat rate. In
More informationReview for Exam #2. Specific Heat, Thermal Conductivity, and Thermal Diffusivity. Conduction
Review for Exam # Speifi Heat, Thermal Condutivity, and Thermal Diffusivity Speifi heat ( p ) A measure of how muh energy is required to raise the temperature of an objet Thermal ondutivity (k) A measure
More informationEffect of Different Types of Promoters on Bed Expansion in a Gas-Solid Fluidized Bed with Varying Distributor Open Areas
Journal of Chemial Engineering of Japan, Vol. 35, No. 7, pp. 681 686, 2002 Short Communiation Effet of Different Types of Promoters on Bed Expansion in a Gas-Solid Fluidized Bed with Varying Distributor
More information( ) ( ) Volumetric Properties of Pure Fluids, part 4. The generic cubic equation of state:
CE304, Spring 2004 Leture 6 Volumetri roperties of ure Fluids, part 4 The generi ubi equation of state: There are many possible equations of state (and many have been proposed) that have the same general
More informationCOMBINED PROBE FOR MACH NUMBER, TEMPERATURE AND INCIDENCE INDICATION
4 TH INTERNATIONAL CONGRESS OF THE AERONAUTICAL SCIENCES COMBINED PROBE FOR MACH NUMBER, TEMPERATURE AND INCIDENCE INDICATION Jiri Nozika*, Josef Adame*, Daniel Hanus** *Department of Fluid Dynamis and
More informationDYNAMICS OF CAPILLARY FLOW AND TRANSPORT PROPERTIES IN CARBONATE SEDIMENTARY FORMATION BY TIME-CONTROLLED POROSIMETRY
SCA4-44 /5 DYNAMICS OF CAPILLARY FLOW AND TRANSPORT PROPERTIES IN CARBONATE SEDIMENTARY FORMATION BY TIME-CONTROLLED POROSIMETRY A. Cerepi, Institut EGID-Bordeaux 3, Université Mihel de Montaigne,, allée
More informationPart G-4: Sample Exams
Part G-4: Sample Exams 1 Cairo University M.S.: Eletronis Cooling Faulty of Engineering Final Exam (Sample 1) Mehanial Power Engineering Dept. Time allowed 2 Hours Solve as muh as you an. 1. A heat sink
More informationCRITICAL EXPONENTS TAKING INTO ACCOUNT DYNAMIC SCALING FOR ADSORPTION ON SMALL-SIZE ONE-DIMENSIONAL CLUSTERS
Russian Physis Journal, Vol. 48, No. 8, 5 CRITICAL EXPONENTS TAKING INTO ACCOUNT DYNAMIC SCALING FOR ADSORPTION ON SMALL-SIZE ONE-DIMENSIONAL CLUSTERS A. N. Taskin, V. N. Udodov, and A. I. Potekaev UDC
More informationLOAD-RATIO DEPENDENCE ON FATIGUE LIFE OF COMPOSITES
LOAD-RATIO DEPENDENCE ON FATIGUE LIFE OF COMPOSITES Joakim Shön 1 and Anders F. Blom 1, 1 Strutures Department, The Aeronautial Researh Institute of Sweden Box 1101, SE-161 11 Bromma, Sweden Department
More informationA simple expression for radial distribution functions of pure fluids and mixtures
A simple expression for radial distribution funtions of pure fluids and mixtures Enrio Matteoli a) Istituto di Chimia Quantistia ed Energetia Moleolare, CNR, Via Risorgimento, 35, 56126 Pisa, Italy G.
More informationGreen s function for the wave equation
Green s funtion for the wave equation Non-relativisti ase January 2019 1 The wave equations In the Lorentz Gauge, the wave equations for the potentials are (Notes 1 eqns 43 and 44): 1 2 A 2 2 2 A = µ 0
More information2 How far? Equilibrium Answers
How far? Equilibrium Answers ratie: pages 37 39 1 Answer is D. Only a hange in temperature harges the value of the equilibrium onstant. Answer is D. [B] /[A] so [B] [A] or [B] [A] 1/ 3 Answer is B. Amounts
More informationGeneral Equilibrium. What happens to cause a reaction to come to equilibrium?
General Equilibrium Chemial Equilibrium Most hemial reations that are enountered are reversible. In other words, they go fairly easily in either the forward or reverse diretions. The thing to remember
More informationExperiment 03: Work and Energy
MASSACHUSETTS INSTITUTE OF TECHNOLOGY Physis Department Physis 8.01 Purpose of the Experiment: Experiment 03: Work and Energy In this experiment you allow a art to roll down an inlined ramp and run into
More informationDIGITAL DISTANCE RELAYING SCHEME FOR PARALLEL TRANSMISSION LINES DURING INTER-CIRCUIT FAULTS
CHAPTER 4 DIGITAL DISTANCE RELAYING SCHEME FOR PARALLEL TRANSMISSION LINES DURING INTER-CIRCUIT FAULTS 4.1 INTRODUCTION Around the world, environmental and ost onsiousness are foring utilities to install
More informationTheory. Coupled Rooms
Theory of Coupled Rooms For: nternal only Report No.: R/50/TCR Prepared by:. N. taey B.., MO Otober 00 .00 Objet.. The objet of this doument is present the theory alulations to estimate the reverberant
More informationELECTROMAGNETIC NORMAL MODES AND DISPERSION FORCES.
ELECTROMAGNETIC NORMAL MODES AND DISPERSION FORCES. All systems with interation of some type have normal modes. One may desribe them as solutions in absene of soures; they are exitations of the system
More informationRC DEEP BEAMS ANALYSIS CONSIDERING LOCALIZATION IN COMPRESSION
RC DEEP BEAMS ANAYSIS CONSIDERING OCAIZATION IN COMPRESSION Manakan ERTSAMATTIYAKU* 1, Torsak ERTSRISAKURAT* 1, Tomohiro MIKI* 1 and Junihiro NIWA* ABSTRACT: It has been found that RC deep beams usually
More informationThe Laws of Acceleration
The Laws of Aeleration The Relationships between Time, Veloity, and Rate of Aeleration Copyright 2001 Joseph A. Rybzyk Abstrat Presented is a theory in fundamental theoretial physis that establishes the
More informationCombined Electric and Magnetic Dipoles for Mesoband Radiation, Part 2
Sensor and Simulation Notes Note 53 3 May 8 Combined Eletri and Magneti Dipoles for Mesoband Radiation, Part Carl E. Baum University of New Mexio Department of Eletrial and Computer Engineering Albuquerque
More informationThe Hanging Chain. John McCuan. January 19, 2006
The Hanging Chain John MCuan January 19, 2006 1 Introdution We onsider a hain of length L attahed to two points (a, u a and (b, u b in the plane. It is assumed that the hain hangs in the plane under a
More informationUPPER-TRUNCATED POWER LAW DISTRIBUTIONS
Fratals, Vol. 9, No. (00) 09 World Sientifi Publishing Company UPPER-TRUNCATED POWER LAW DISTRIBUTIONS STEPHEN M. BURROUGHS and SARAH F. TEBBENS College of Marine Siene, University of South Florida, St.
More informationExperimental Investigation and FE Analysis of Fiber Woven Layered Composites under Dynamic Loading
2th International LS-DYNA Users Conferene Constitutive Modeling(2) xperimental Investigation and F Analysis of Fiber Woven Layered Composites under Dynami Loading Pavel A. Mossakovsky, Fedor K. Antonov,
More informationWhere as discussed previously we interpret solutions to this partial differential equation in the weak sense: b
Consider the pure initial value problem for a homogeneous system of onservation laws with no soure terms in one spae dimension: Where as disussed previously we interpret solutions to this partial differential
More informationA NORMALIZED EQUATION OF AXIALLY LOADED PILES IN ELASTO-PLASTIC SOIL
Journal of Geongineering, Vol. Yi-Chuan 4, No. 1, Chou pp. 1-7, and April Yun-Mei 009 Hsiung: A Normalized quation of Axially Loaded Piles in lasto-plasti Soil 1 A NORMALIZD QUATION OF AXIALLY LOADD PILS
More informationInter-fibre contacts in random fibrous materials: experimental verification of theoretical dependence on porosity and fibre width
J Mater Si (2006) 41:8377 8381 DOI 10.1007/s10853-006-0889-7 LETTER Inter-fibre ontats in random fibrous materials: experimental verifiation of theoretial dependene on porosity and fibre width W. J. Bathelor
More information23.1 Tuning controllers, in the large view Quoting from Section 16.7:
Lesson 23. Tuning a real ontroller - modeling, proess identifiation, fine tuning 23.0 Context We have learned to view proesses as dynami systems, taking are to identify their input, intermediate, and output
More information755. A survey of Hopf bifurcation analysis in nonlinear railway wheelset dynamics
755. A survey of Hopf bifuration analysis in nonlinear railway wheelset dynamis Hamid M. Sedighi, Kourosh H. Shirazi Department of Mehanial Engineering, Shahid Chamran University, Ahvaz, Iran E-mail: hmsedighi@gmail.om,
More informationAdvanced Computational Fluid Dynamics AA215A Lecture 4
Advaned Computational Fluid Dynamis AA5A Leture 4 Antony Jameson Winter Quarter,, Stanford, CA Abstrat Leture 4 overs analysis of the equations of gas dynamis Contents Analysis of the equations of gas
More informationFrequency Domain Analysis of Concrete Gravity Dam-Reservoir Systems by Wavenumber Approach
Frequeny Domain Analysis of Conrete Gravity Dam-Reservoir Systems by Wavenumber Approah V. Lotfi & A. Samii Department of Civil and Environmental Engineering, Amirkabir University of Tehnology, Tehran,
More informationHeat exchangers: Heat exchanger types:
Heat exhangers: he proess of heat exhange between two fluids that are at different temperatures and separated by a solid wall ours in many engineering appliations. he devie used to implement this exhange
More information13.Prandtl-Meyer Expansion Flow
3.Prandtl-eyer Expansion Flow This hapter will treat flow over a expansive orner, i.e., one that turns the flow outward. But before we onsider expansion flow, we will return to onsider the details of the
More informationChapter 14. The Concept of Equilibrium and the Equilibrium Constant. We have for the most part depicted reactions as going one way.
Chapter 14 The Conept of Equilibrium and the Equilibrium Constant In hapter 1 we dealt with Physial Equilibrium Physial Changes HO 2 (l) HO 2 (g) In hapter 14 we will learn about Chemial Equilibrium. We
More informationSubject: Introduction to Component Matching and Off-Design Operation % % ( (1) R T % (
16.50 Leture 0 Subjet: Introdution to Component Mathing and Off-Design Operation At this point it is well to reflet on whih of the many parameters we have introdued (like M, τ, τ t, ϑ t, f, et.) are free
More informationWavetech, LLC. Ultrafast Pulses and GVD. John O Hara Created: Dec. 6, 2013
Ultrafast Pulses and GVD John O Hara Created: De. 6, 3 Introdution This doument overs the basi onepts of group veloity dispersion (GVD) and ultrafast pulse propagation in an optial fiber. Neessarily, it
More informationWhat s New in ChemSep TM 6.8
What s New in ChemSep TM 6.8 January 2011 (Updated Marh 2011) Harry Kooijman and Ross Taylor In this doument we identify and desribe the most important new features in ChemSep. 1. New: GUI an diretly load
More information10.2 The Occurrence of Critical Flow; Controls
10. The Ourrene of Critial Flow; Controls In addition to the type of problem in whih both q and E are initially presribed; there is a problem whih is of pratial interest: Given a value of q, what fators
More informationUTC. Engineering 329. Proportional Controller Design. Speed System. John Beverly. Green Team. John Beverly Keith Skiles John Barker.
UTC Engineering 329 Proportional Controller Design for Speed System By John Beverly Green Team John Beverly Keith Skiles John Barker 24 Mar 2006 Introdution This experiment is intended test the variable
More informationBeams on Elastic Foundation
Professor Terje Haukaas University of British Columbia, Vanouver www.inrisk.ub.a Beams on Elasti Foundation Beams on elasti foundation, suh as that in Figure 1, appear in building foundations, floating
More informationChapter 3 Lecture 7. Drag polar 2. Topics. Chapter-3
hapter 3 eture 7 Drag polar Topis 3..3 Summary of lift oeffiient, drag oeffiient, pithing moment oeffiient, entre of pressure and aerodynami entre of an airfoil 3..4 Examples of pressure oeffiient distributions
More informationMass Transfer (Stoffaustausch) Fall 2012
Mass Transfer (Stoffaustaush) Fall Examination 9. Januar Name: Legi-Nr.: Edition Diffusion by E. L. Cussler: none nd rd Test Duration: minutes The following materials are not permitted at your table and
More informationFORCE DISTRIBUTION OF REINFORCED CONCRETE COUPLING BEAMS WITH DIAGONAL REINFORCEMENT
FORCE DISTRIBUTION OF REINFORCED CONCRETE COULING BEAMS WITH DIAGONAL REINFORCEMENT Yenny Nurhasanah Jurusan Teknik Sipil, Fakultas Teknik, Universitas Muhammadiyah Surakarta Jl. A. Yani Tromol os 1 abelan
More informationThree-dimensional Meso-scopic Analyses of Mortar and Concrete Model by Rigid Body Spring Model
Three-dimensional Meso-sopi Analyses of Mortar and Conrete Model by Rigid Body Spring Model K. Nagai, Y. Sato & T. Ueda Hokkaido University, Sapporo, Hokkaido, JAPAN ABSTRACT: Conrete is a heterogeneity
More informationSTUDY OF INHERENT FREQUENCY OF HELMHOLTZ RESONATOR
005 WJTA Amerian Waterjet Conferene August -3, 005! Houston, Texas Paper 6B-4 STUDY OF INHERENT FREQUENCY OF HELMHOLT RESONATOR Gong Weili An Liqian Cui Longlian Xie Guixin Shool of Mehanis, Arhiteture
More informationFig Review of Granta-gravel
0 Conlusion 0. Sope We have introdued the new ritial state onept among older onepts of lassial soil mehanis, but it would be wrong to leave any impression at the end of this book that the new onept merely
More informationDr G. I. Ogilvie Lent Term 2005
Aretion Diss Mathematial Tripos, Part III Dr G. I. Ogilvie Lent Term 2005 1.4. Visous evolution of an aretion dis 1.4.1. Introdution The evolution of an aretion dis is regulated by two onservation laws:
More informationCharacterizing Pure and Undefined Petroleum Components
International Journal of Engineering & ehnology IJE-IJENS Vol:10 No:0 8 Charaterizing Pure and Undefined Petroleum Components Hassan S. Naji King Abdulaziz University, Jeddah, Saudi Arabia Website: http://hnaji.au.edu.sa
More informationFailure Assessment Diagram Analysis of Creep Crack Initiation in 316H Stainless Steel
Failure Assessment Diagram Analysis of Creep Crak Initiation in 316H Stainless Steel C. M. Davies *, N. P. O Dowd, D. W. Dean, K. M. Nikbin, R. A. Ainsworth Department of Mehanial Engineering, Imperial
More informationSEDIMENT TRANSPORT CALCULATION CONSIDERING COHESIVE EFFECTS AND ITS APPLICATION TO WAVE-INDUCED TOPOGRAPHIC CHANGE
Proeedings of the 7 th International Conferene on Asian and Paifi Coasts (APAC 03) Bali, Indonesia, September 4-6, 03 SEDIMENT TRANSPORT CALCULATION CONSIDERING COHESIVE EFFECTS AND ITS APPLICATION TO
More information2. Mass transfer takes place in the two contacting phases as in extraction and absorption.
PRT 11- CONVECTIVE MSS TRNSFER 2.1 Introdution 2.2 Convetive Mass Transfer oeffiient 2.3 Signifiant parameters in onvetive mass transfer 2.4 The appliation of dimensional analysis to Mass Transfer 2.4.1
More informationTHE EQUATION CONSIDERING CONCRETE STRENGTH AND STIRRUPS FOR DIAGONAL COMPRESSIVE CAPACITY OF RC BEAM
- Tehnial Paper - THE EQUATION CONSIDERING CONCRETE STRENGTH AND STIRRUPS FOR DIAGONAL COMPRESSIE CAPACITY OF RC BEAM Patarapol TANTIPIDOK *, Koji MATSUMOTO *, Ken WATANABE *3 and Junihiro NIWA *4 ABSTRACT
More informationModeling of Threading Dislocation Density Reduction in Heteroepitaxial Layers
A. E. Romanov et al.: Threading Disloation Density Redution in Layers (II) 33 phys. stat. sol. (b) 99, 33 (997) Subjet lassifiation: 6.72.C; 68.55.Ln; S5.; S5.2; S7.; S7.2 Modeling of Threading Disloation
More informationSlenderness Effects for Concrete Columns in Sway Frame - Moment Magnification Method
Slenderness Effets for Conrete Columns in Sway Frame - Moment Magnifiation Method Slender Conrete Column Design in Sway Frame Buildings Evaluate slenderness effet for olumns in a sway frame multistory
More informationDetermination of the reaction order
5/7/07 A quote of the wee (or amel of the wee): Apply yourself. Get all the eduation you an, but then... do something. Don't just stand there, mae it happen. Lee Iaoa Physial Chemistry GTM/5 reation order
More informationIN-PLANE VIBRATIONS OF CURVED BEAMS WITH VARIABLE CROSS-SECTIONS CARRYING ADDITIONAL MASS
11 th International Conferene on Vibration Problems Z. Dimitrovová et al. (eds.) Lisbon, Portugal, 9-1 September 013 IN-PLANE VIBRATIONS OF CURVED BEAMS WITH VARIABLE CROSS-SECTIONS CARRYING ADDITIONAL
More informationFour-dimensional equation of motion for viscous compressible substance with regard to the acceleration field, pressure field and dissipation field
Four-dimensional equation of motion for visous ompressible substane with regard to the aeleration field, pressure field and dissipation field Sergey G. Fedosin PO box 6488, Sviazeva str. -79, Perm, Russia
More informationKINETICS OF IRON OXIDE DIRECT REDUCTION BY COAL E.R. ABRIL 1
KINETICS OF IRON OXIDE DIRECT REDUCTION BY COAL E.R. ABRIL 1 CIMM- Av.Velez Sarsfield 1561 C.P.5000 Córdoba, Argentina. aabril@intiemor.gov.ar Abstrat - A new interpretation to the kinetis of iron oxide
More informationShear-Friction Strength of RC Walls with 550 MPa Bars
Proeedings of the Tenth Paifi Conferene on Earthquake Engineering Building an Earthquake-Resilient Paifi 6-8 November 215, Sydney, Australia Shear-Frition Strength of RC Walls with 55 MPa Bars Jang-woon
More informationPASSIVE CONTROL OF BUILDING FRAMES BY MEANS OF LIQUID DAMPERS SEALED BY VISCOELASTIC MATERIAL
ASSIVE CONTROL OF BUILDING FRAMES BY MEANS OF LIQUID DAMERS SEALED BY VISCOELASTIC MATERIAL 98 Shigeya KAWAMATA, Naoki FUNAKI And Yoshiro ITOH 3 SUMMARY In the purpose of installing in passive vibration
More informationBEAMS: SHEARING STRESS
LECTURE Third Edition BEAMS: SHEARNG STRESS A. J. Clark Shool of Engineering Department of Civil and Environmental Engineering 14 Chapter 6.1 6.4 b Dr. brahim A. Assakkaf SPRNG 200 ENES 220 Mehanis of
More informationA model for measurement of the states in a coupled-dot qubit
A model for measurement of the states in a oupled-dot qubit H B Sun and H M Wiseman Centre for Quantum Computer Tehnology Centre for Quantum Dynamis Griffith University Brisbane 4 QLD Australia E-mail:
More informationTest of General Relativity Theory by Investigating the Conservation of Energy in a Relativistic Free Fall in the Uniform Gravitational Field
Test of General Relatiity Theory by Inestigating the Conseration of Energy in a Relatiisti Free Fall in the Uniform Graitational Field By Jarosla Hyneek 1 Abstrat: This paper inestigates the General Relatiity
More informationSlenderness Effects for Concrete Columns in Sway Frame - Moment Magnification Method
Slenderness Effets for Conrete Columns in Sway Frame - Moment Magnifiation Method Slender Conrete Column Design in Sway Frame Buildings Evaluate slenderness effet for olumns in a sway frame multistory
More informationNumerical simulation of a one-dimensional shock tube problem at supercritical fluid conditions
International Journal of Physial Sienes Vol. 3 (1), pp. 314-30, Deember, 008 Available online at http://www.aademijournals.org/ijps ISSN 199-1950 008 Aademi Journals Full ength esearh Paper Numerial simulation
More informationTHE SPANN VIBROACOUSTIC METHOD Revision A
THE SPNN VIBROCOUSTIC METHOD Revision By Tom Irvine Deember 15, 01 Email: tom@vibrationdata.om Figure 1. vionis Installation and Testing Introdution vionis omponents in airraft and launh vehiles may be
More information3 Tidal systems modelling: ASMITA model
3 Tidal systems modelling: ASMITA model 3.1 Introdution For many pratial appliations, simulation and predition of oastal behaviour (morphologial development of shorefae, beahes and dunes) at a ertain level
More information+Ze. n = N/V = 6.02 x x (Z Z c ) m /A, (1.1) Avogadro s number
In 1897, J. J. Thomson disovered eletrons. In 1905, Einstein interpreted the photoeletri effet In 1911 - Rutherford proved that atoms are omposed of a point-like positively harged, massive nuleus surrounded
More informationNUMERICAL SIMULATION OF SEA ICE GROWTH AND DESALINATION
17th International Symposium on Ie Saint Petersburg, Russia, 1-5 June 004 International Assoiation of Hydrauli Engineering and Researh Session I NUMERICAL SIMULATION OF SEA ICE GROWTH AND DESALINATION
More informationCALCULATION OF NONLINEAR TUNE SHIFT USING BEAM POSITION MEASUREMENT RESULTS
International Journal of Modern Physis A Vol. 24, No. 5 (2009) 974 986 World Sientifi Publishing Company CALCULATION OF NONLINEAR TUNE SHIFT USING BEAM POSITION MEASUREMENT RESULTS PAVEL SNOPOK, MARTIN
More informationEvaluation of effect of blade internal modes on sensitivity of Advanced LIGO
Evaluation of effet of blade internal modes on sensitivity of Advaned LIGO T0074-00-R Norna A Robertson 5 th Otober 00. Introdution The urrent model used to estimate the isolation ahieved by the quadruple
More informationNUMERICAL SIMULATION OF SEA ICE GROWTH AND DESALINATION
NUMERICAL SIMULATION OF SEA ICE GROWTH AND DESALINATION Chris Petrih 1, Pat J. Langhorne 1, Zhifa Sun 1 ABSTRACT We present a numerial model based on the finite volume method to simulate sea ie growth
More informationELECTROMAGNETIC WAVES
ELECTROMAGNETIC WAVES Now we will study eletromagneti waves in vauum or inside a medium, a dieletri. (A metalli system an also be represented as a dieletri but is more ompliated due to damping or attenuation
More informationINTERNATIONAL JOURNAL OF CIVIL AND STRUCTURAL ENGINEERING Volume 2, No 4, 2012
INTERNATIONAL JOURNAL OF CIVIL AND STRUCTURAL ENGINEERING Volume, No 4, 01 Copyright 010 All rights reserved Integrated Publishing servies Researh artile ISSN 0976 4399 Strutural Modelling of Stability
More informationu x u t Internal Waves
Internal Waves We now examine internal waves for the ase in whih there are two distint layers and in whih the lower layer is at rest. This is an approximation of the ase in whih the upper layer is muh
More informationCollinear Equilibrium Points in the Relativistic R3BP when the Bigger Primary is a Triaxial Rigid Body Nakone Bello 1,a and Aminu Abubakar Hussain 2,b
International Frontier Siene Letters Submitted: 6-- ISSN: 9-8, Vol., pp -6 Aepted: -- doi:.8/www.sipress.om/ifsl.. Online: --8 SiPress Ltd., Switzerland Collinear Equilibrium Points in the Relativisti
More informationEvaluation of a Dual-Load Nondestructive Testing System To Better Discriminate Near-Surface Layer Moduli
52 TRANSPORTATION RESEARCH RECORD 1355 Evaluation of a Dual-Load Nondestrutive Testing System To Better Disriminate Near-Surfae Layer Moduli REYNALDO ROQUE, PEDRO ROMERO, AND BYRON E. RUTH Theoretial analyses
More informationMass Transfer 2. Diffusion in Dilute Solutions
Mass Transfer. iffusion in ilute Solutions. iffusion aross thin films and membranes. iffusion into a semi-infinite slab (strength of weld, tooth deay).3 Eamples.4 ilute diffusion and onvetion Graham (85)
More informationECE-320 Linear Control Systems. Winter 2013, Exam 1. No calculators or computers allowed, you may leave your answers as fractions.
ECE-320 Linear Control Systems Winter 2013, Exam 1 No alulators or omputers allowed, you may leave your answers as frations. All problems are worth 3 points unless noted otherwise. Total /100 1 Problems
More informationSpeed-feedback Direct-drive Control of a Low-speed Transverse Flux-type Motor with Large Number of Poles for Ship Propulsion
Speed-feedbak Diret-drive Control of a Low-speed Transverse Flux-type Motor with Large Number of Poles for Ship Propulsion Y. Yamamoto, T. Nakamura 2, Y. Takada, T. Koseki, Y. Aoyama 3, and Y. Iwaji 3
More informationA Spatiotemporal Approach to Passive Sound Source Localization
A Spatiotemporal Approah Passive Sound Soure Loalization Pasi Pertilä, Mikko Parviainen, Teemu Korhonen and Ari Visa Institute of Signal Proessing Tampere University of Tehnology, P.O.Box 553, FIN-330,
More informationEFFECTS OF COUPLE STRESSES ON PURE SQUEEZE EHL MOTION OF CIRCULAR CONTACTS
-Tehnial Note- EFFECTS OF COUPLE STRESSES ON PURE SQUEEZE EHL MOTION OF CIRCULAR CONTACTS H.-M. Chu * W.-L. Li ** Department of Mehanial Engineering Yung-Ta Institute of Tehnology & Commere Ping-Tung,
More informationCity, University of London Institutional Repository
City Researh Online City, University of London Institutional Repository Citation: Labib, M., Moslehy, Y. & Ayoub, A. (07). Softening Coeffiient of Reinfored Conrete Elements Subjeted to Three-Dimensional
More informationEFFECTIVE STRESS LAW FOR THE PERMEABILITY OF CLAY-RICH SANDSTONES
SCA22-5 1/6 EFFECTIVE STRESS LAW FOR THE PERMEABILITY OF CLAY-RICH SANDSTONES Widad Al-Wardy and Robert W. Zimmerman Department of Earth Siene and Engineering Imperial College of Siene, Tehnology and Mediine
More informationCONVECTION AT A.MODEL ICE EDGE
JACK CALMAN CONVECTION AT A.MODEL ICE EDGE The flow pattern near the edge of a melting ie blok is modeled by heating a metal edge in a saltstratified fluid. An unexpetedly strong, horizontal boundary urrent
More informationSimple Considerations on the Cosmological Redshift
Apeiron, Vol. 5, No. 3, July 8 35 Simple Considerations on the Cosmologial Redshift José Franiso Garía Juliá C/ Dr. Maro Mereniano, 65, 5. 465 Valenia (Spain) E-mail: jose.garia@dival.es Generally, the
More informationDiffusion and compensating subsidence as limiting cases of a single flux parameterization. David Randall
! Revised September 014 9:49 AM! 1 Diffusion and ompensating subsidene as limiting ases of a single flux parameterization David Randall The tendeny of the average of a onservative variable h due to a vertial
More informationPOROUS CARBON PARTICLE COMBUSTION IN AIR
MCS 7 Chia Laguna, Cagliari, Sardinia, taly, 11-15, 11 POOUS CABON PATCLE COMBUSTON N A V. M. Gremyahkin grema@ipmnet.ru nstitute for Problems in Mehanis, AS, Mosow, ussia Abstrat Theoretial investigation
More information