Effect of no-flow boundaries on interference testing. in fractured reservoirs

Size: px
Start display at page:

Download "Effect of no-flow boundaries on interference testing. in fractured reservoirs"

Transcription

1 Effect of no-flo boundaies on intefeence testing in factued esevois T.Aa. Jelmet 1 1 epatement of petoleum engineeing and applied geophysics,, Noegian Univesity of Science and Tecnology, NTNU. Tondheim, Noay Coesponding autho: tomj@ntnu.no ABSTRACT: The existence of a no-flo bounday is impotant fo esevoi management. This is because fluid cannot flo acoss a sealing bounday. In many cases no-flo boundaies do not sho up on seismic. Sub-seismic faults may be detected by ell testing. This study deals ith mathematical modeling of intefeence tests in factued esevois. The objective is to help locate sub-seismic faults to impove esevoi management. We concentate on factued netoks of factal type. An intefeence test involves at least to ells, an active ell and an obsevation ell. The latte is idle. The esulting mathematical model includes the single ell esponses as a special case. An equation to pedict the pessue behavio of intefeence and single ell tests in the pesence of a no-flo bounday has been deived. To exteme cases ee identified. These depend on the distance beteen the active ell and the fault hich could be infinitely lage (case 1 o zeo (case 2. All othe cases fall in beteen these. It as found that the influence of a no-flo bounday shos up as a tansition beteen the to limiting cuves. Initially the pessue esponse ill follo the cuve of case 1. Afte the tansition the pessue signatue ill mege ith cuve to. A lage distance leads to late tansition. KEYWORS : Well testing, factued esevoi, factal esevoi. 1. Intoduction We concentate on factue netoks of factal behavio. The factued esevoi is ithout seconday poosity. Such esevois ae chaacteized by factue netoks of diffeent achitectue. Well test intepetation may povide impotant input to the geological model of the esevoi. This is because pessue tansient data eflect esevoi dynamics. In case of an intefeence test, the entie volume in beteen ells is sampled. The taditional intepetation methodology is type cuve matching. A type cuve is a gaph of the solution to the dimensionless flo equation in a log-log coodinate system. If the ight model has been obtained, then the ell test esponse and the type cuve ae chaacteized by the same shape. In fact, if one is put on top the othe, the to cuves ill ovelap. This popety is the basis of classical type cuve analysis, Ealoughe (1977. Chang and Yotsos (1990 poposed a one-phase model of fluid flo in factal netoks. Thei solution is fo single ell tests. Beie (1990 applied the Chang and Yotsos factal model to a non-factued but disodeed inhomogeneous esevoi. Apillan et al. (1993

2 studied the pessue behavio of intefeence tests. Thei methodology as based on a similaity solution and did not accommodate ellboe stoage, skin o possible faults. 2. Methodology We assume an infinite esevoi and constant poduction in the active ell. The Laplace space solution of the govening equation, see Chang and Yotsos (1990, as obtained by classical theoy of diffeential equations. This in tun as tansfomed back to time domain by a numeical technique, Stehfest (1970. Much has been itten about the behavio of factals. The Society of Petoleum Enginees (SPE data base has been seached fo a solution to the poblem of an intefeence test ith ellboe stoage and skin at the active ell. We ee unable to locate one. Hence e decided to deive a solution athe than seach additional databases. The effect of a no-flo bounday as geneated by the method of images. Sandal et al. (1978 found the pessue in an obsevation ell due to a poduction ell ith ellboe stoage and skin fo a homogeneous esevoi. The same methodology is applicable to a factal esevoi. The esulting equation fo the Laplace tansfom of pessue behavio in the obsevation ell is: p ( ( 1 δ ν ( ξ Φ Φ K = s C s K S sk sk ( ν ( ξ + ν 1( ξ + ν 1( ξ (1 and in the poduction ell: p Kυ ( ξ + S skυ 1 ( ξ ( ξ + ( ξ + ( ξ = s C s K S sk sk ( υ υ 1 υ 1 (2 υ = 1 δ (3 ( θ ξ = 2/ + 2 s (4 ( θ Φ= 2 + /2 (5 δ = ( θ + 2 (6

3 hee: C is ellboe stoage constant is factal dimension d is Eucledian dimension S is skin K is modified Bessel function of second kind, of any eal ode ν θ is connectivity index is the distance beteen ells s is Laplace vaiable The dimensionless vaiables ae k k ( t = t = t (7 ϕμ ϕ μ θ c t ( c t obs = (8 d 2 kg pf = ( pi pf (9 ( d 2 d 2 d+ θ kg k G p( = ( pi p( = ( pi p( (10 Folloing Apillan (1990, e define ( t k TF = = t (11 c obs 2+ θ 2 ϕ( obs μ t obs ( ( p k G PF = = p p (12 d 2 obs obs d+ θ ( i obs k is the pemeability aound the active ell and k( is the pemeability at a distance fom the active ell. Well test intepetation may be caied out by type cuve matching as explained by Ealoughe (1977. d=1, 2 and 3 coesponds to flo in a Catesian, adial and spheical coodinate system espectively. The constant G assumes diffeent values depending on the Eucledian dimension: d = 1 G= A, d = 2 G= 2π h and d = 3 G= 4π. The distance beteen the active ell and the obsevation ell is: = Δ x +Δ y ( obs obs obs

4 No-flo boundaies ee simulated by the method of images. Suppose a no-flo fault is located at a distance Δy F aay fom the active ell. The fault is of infinite extension and uns paallel to the x-axis. The physical bounday may be eplaced by an image ell at the same distance, Δy F, behind the fault. Intefeence beteen the eal and fictitious (image ell ill ceate a no-flo bounday at the position of the fault. The total esponse of the to ells may be calculated by supeposition. This is because of the lineaity of the govening equations. The distance beteen the obsevation ell and image ell is: ( 2 = Δ x + Δ y ±Δ y (14 2 img obs img obs hee Δ yimg = 2Δ yf. The sign depends on hich ell, active o obsevation ell, has the shotest distance to the fault. It is positive hen the active ell is moe close to the fault than the obsevation ell and vice vesa. The pessue dop in the poduction ell in the pesence of a fault is: ( ( p t = p, t + p (, t (15 obs obs img 3. Results The effect of the distance to the bounday on the pessue behavio is shon in Fig. 1 and on the logaithmic pessue deivative in Fig. 2. A dimensionless distance of 2000 coesponds to 200 m fo a petoleum esevoi. This is because the adius of the ellboe is aound 10 cm. Fo a shalloe esevoi the ellboe adius could be geate. We find that the effect of the ellboe stoage constant and skin is negligible hen the dimensionless distance beteen the active and obsevation ell, obs, is As expected the effect of the fault shos up ealie ith deceasing distance, img. Table 1: Input data: Effect of the distance to a fault Run obs img d θ γ S C

5 PF=p / γ (2+Theta TF=t / obs Fig. 1. The effect of distance fom a fault on the Pessue Function PF' TF Fig. 2. The effect of distance fom a fault on the deivative of the Pessue Function 4. iscussion The factal model implies the assumption of an inhomogeneous esevoi. The ock popeties depend on the spatial coodinate in the fom of poe la expessions. In most cases the ock popeties ill be deceasing functions of the distance aay fom the active ell. This means that the pemeability and poosity at the obsevation ell ill be diffeent as seen fom the active ell and the image ell. This is because the to ells ae located at diffeent distances aay fom the obsevation ell. Then one could question the validity of the supeposition by simple addition of pd-functions. Appaently pd- and td- functions (fo the to ells ae scaled diffeent ays since they ae based on ock popeties evaluated at diffeent distances aay fom the eal and the imaginay image ell. The pimay dimensionless vaiables,

6 hoeve, ae scaled accoding to the ock popeties at the active ell,, see equations (7 and (10, left hand side. Hence diect addition of pd-functions is possible. A negative skin facto may be accounted fo by use the equivalent ellboe adius, e. Fo sylindical geomety, d=2, the elationship beteen equivalent and actual ellboe adius is given by: e = e S Suppose the skin facto of a stimulated ell is -3. Then the effective adius is 20 times lage than the adius dilled. Then the dimensionless distances to the obsevation ell and the image ell ae educed by the same facto. Unde such conditions e find that that the ellboe stoage facto and skin may have an influence on the pessue signatue. The ellboe stoage constant is moe impotant than the skin fo an intefeence test. 5. Conclusions and efeences A numeical model to simulate the effect of faults in a factal esevoi has been implemented. The numeical model depends on the Laplace space solution to the govening equation. The solution as conveted back to time domain by a numeical technique. Thee ae to limiting cuves fo a esevoi ith a fault. The loe one is chaacteized by no fault at all o equivalently a fault at an infinite distance. The uppe one is chaacteized by a fault in the immediate neighbohood of the active ell, Δ y F = 0. When a fault is placed anyhee in beteen the limiting positions, the esponse is chaacteized by a tansition fom the loe to the uppe cuve. The stat of the tansition peiod is delayed ith inceasing distance, Δ y. Hence thee is a possibility to estimate the distance to a sealing fault. F The ellboe stoage facto and skin of the active ell may have an effect on the pessue signatue in the obsevation ell in case the dimensionless distance beteen the to ells is in the ode of 10. This situation is likely to occu if the active ell has been stimulated. REFERENCES Beie, R. A., Pessue Tansient field data shoing factal esevoi stuctue. Pape SPE 21553, CIM/SPE Intenational technical meeting, Calgay, Albeta June Chang, J. & Yotsos, Y.C., Pessue tansient behavio of factal esevois. SPE Fomation Evaluation, (Mach: Ealoghe, R. C., 1977 Advances in ell testing. Monogaph Volume 5 of the Heny L. ohety seies, Society of Petoleum Enginees, Ne Yok, p. 107 Sandal, H., Hone, R. N., Ramey, H. J. & Williamson, J. W., Intefeence testing ith ellboe stoage and skin at the poducing ell. Pape SPE 7454, 83 d SPE Annual Meeting and Exhibition of the Society of Petoleum Enginees, Houston, Oct 1-4 Stehfest, H. Algoithm 368: Numeical invesion of Laplace tansfoms. -5, Communications of the ACM. Jan. (13, 47-49

An Investigation of Wellbore Storage and Skin Effect in Unsteady Liquid Flow: II. Finite Difference Treatment

An Investigation of Wellbore Storage and Skin Effect in Unsteady Liquid Flow: II. Finite Difference Treatment An Investigation of Wellboe Stoage and Sin Effect in Unsteady Liquid Flo: II. Finite Diffeence Teatment ROBERT A. WATTENBARGER* H. J. RAMEY, JR. MEMBERS A/ME MOBI L RESEARCH & DEVELOPMENT CORP. DALLAS,

More information

ANALYSIS OF PRESSURE VARIATION OF FLUID IN AN INFINITE ACTING RESERVOIR

ANALYSIS OF PRESSURE VARIATION OF FLUID IN AN INFINITE ACTING RESERVOIR Nigeian Jounal of Technology (NIJOTECH) Vol. 36, No. 1, Januay 2017, pp. 80 86 Copyight Faculty of Engineeing, Univesity of Nigeia, Nsukka, Pint ISSN: 0331-8443, Electonic ISSN: 2467-8821 www.nijotech.com

More information

T. Raja Rani. Military Technological College, Muscat, Oman. Abstract

T. Raja Rani. Military Technological College, Muscat, Oman. Abstract ISSN: 78-8 Vol. Issue, Octobe - Fee Convection ove a Vaying all Vetical Cylinde embedded in a Poous medium ith effect of Radiation, Vaiable Fluid Popeties and Statification. T. Raja Rani Militay Technological

More information

MAGNETIC FIELD AROUND TWO SEPARATED MAGNETIZING COILS

MAGNETIC FIELD AROUND TWO SEPARATED MAGNETIZING COILS The 8 th Intenational Confeence of the Slovenian Society fo Non-Destuctive Testing»pplication of Contempoay Non-Destuctive Testing in Engineeing«Septembe 1-3, 5, Potoož, Slovenia, pp. 17-1 MGNETIC FIELD

More information

PREDICTION OF THERMAL BREAKTHROUGH FROM TRACER TESTS. G. Michael Shook

PREDICTION OF THERMAL BREAKTHROUGH FROM TRACER TESTS. G. Michael Shook PROCEEDINGS, Tenty-Fouth Wokshop on Geothemal Resevoi Engineeing Stanfod Univesity, Stanfod, Califonia, Januay 25-27, 1999 SGP-TR-162 PREDICTION OF THERMAL BREAKTHROUGH FROM TRACER TESTS G. Michael Shook

More information

An Exact Solution of Navier Stokes Equation

An Exact Solution of Navier Stokes Equation An Exact Solution of Navie Stokes Equation A. Salih Depatment of Aeospace Engineeing Indian Institute of Space Science and Technology, Thiuvananthapuam, Keala, India. July 20 The pincipal difficulty in

More information

A Comment on Increasing Returns and Spatial. Unemployment Disparities

A Comment on Increasing Returns and Spatial. Unemployment Disparities The Society fo conomic Studies The nivesity of Kitakyushu Woking Pape Seies No.06-5 (accepted in Mach, 07) A Comment on Inceasing Retuns and Spatial nemployment Dispaities Jumpei Tanaka ** The nivesity

More information

Right-handed screw dislocation in an isotropic solid

Right-handed screw dislocation in an isotropic solid Dislocation Mechanics Elastic Popeties of Isolated Dislocations Ou study of dislocations to this point has focused on thei geomety and thei ole in accommodating plastic defomation though thei motion. We

More information

Liquid gas interface under hydrostatic pressure

Liquid gas interface under hydrostatic pressure Advances in Fluid Mechanics IX 5 Liquid gas inteface unde hydostatic pessue A. Gajewski Bialystok Univesity of Technology, Faculty of Civil Engineeing and Envionmental Engineeing, Depatment of Heat Engineeing,

More information

Solution and Type Curve Analysis of Fluid Flow Model for Fractal Reservoir

Solution and Type Curve Analysis of Fluid Flow Model for Fractal Reservoir Wold Jounal of Mechanics, 0,, 09-6 doi:0.436/wjm.0.507 Published Online Octobe 0 (htt://www.scirp.og/jounal/wjm) Solution and Tye Cuve Analysis of Fluid Flow Model fo Factal Resevoi Yulong Zhao, Liehui

More information

Fresnel Diffraction. monchromatic light source

Fresnel Diffraction. monchromatic light source Fesnel Diffaction Equipment Helium-Neon lase (632.8 nm) on 2 axis tanslation stage, Concave lens (focal length 3.80 cm) mounted on slide holde, iis mounted on slide holde, m optical bench, micoscope slide

More information

Econ 201: Problem Set 2 Answers

Econ 201: Problem Set 2 Answers Econ 0: Poblem Set Anses Instucto: Alexande Sollaci T.A.: Ryan Hughes Winte 08 Question (a) The fixed cost is F C = 4 and the total vaiable costs ae T CV (y) = 4y. (b) To anse this question, let x = (x,...,

More information

Notes on McCall s Model of Job Search. Timothy J. Kehoe March if job offer has been accepted. b if searching

Notes on McCall s Model of Job Search. Timothy J. Kehoe March if job offer has been accepted. b if searching Notes on McCall s Model of Job Seach Timothy J Kehoe Mach Fv ( ) pob( v), [, ] Choice: accept age offe o eceive b and seach again next peiod An unemployed oke solves hee max E t t y t y t if job offe has

More information

Keywords: Non-Darcy flow, Wellbore induced stresses, Analytical solution, Numerical models, Sensitivity analysis

Keywords: Non-Darcy flow, Wellbore induced stresses, Analytical solution, Numerical models, Sensitivity analysis The effect of Non-Dacy Flo on Induced Stesses Aound a boe in an anisotopic in-situ stess field Amin Tohidi 1, Ahmad Fahimifa 1,*, Vamegh asouli 1 Depatment of civil and envionmental Engineeing, Amikabi

More information

Contact impedance of grounded and capacitive electrodes

Contact impedance of grounded and capacitive electrodes Abstact Contact impedance of gounded and capacitive electodes Andeas Hödt Institut fü Geophysik und extateestische Physik, TU Baunschweig The contact impedance of electodes detemines how much cuent can

More information

NOTE. Some New Bounds for Cover-Free Families

NOTE. Some New Bounds for Cover-Free Families Jounal of Combinatoial Theoy, Seies A 90, 224234 (2000) doi:10.1006jcta.1999.3036, available online at http:.idealibay.com on NOTE Some Ne Bounds fo Cove-Fee Families D. R. Stinson 1 and R. Wei Depatment

More information

Mathematical Model of Magnetometric Resistivity. Sounding for a Conductive Host. with a Bulge Overburden

Mathematical Model of Magnetometric Resistivity. Sounding for a Conductive Host. with a Bulge Overburden Applied Mathematical Sciences, Vol. 7, 13, no. 7, 335-348 Mathematical Model of Magnetometic Resistivity Sounding fo a Conductive Host with a Bulge Ovebuden Teeasak Chaladgan Depatment of Mathematics Faculty

More information

PRESSURE DRAWDOWN EQUATIONS FOR MULTIPLE-WELL SYSTEMS IN CIRCULAR-CYLINDRICAL RESERVOIRS

PRESSURE DRAWDOWN EQUATIONS FOR MULTIPLE-WELL SYSTEMS IN CIRCULAR-CYLINDRICAL RESERVOIRS VOL. 8, NO. 7, JULY 013 ISSN 1819-6608 ARPN Jounal of Engineeing and Applied Sciences 006-013 Asian Reseac Publising Netwok (ARPN). All igts eseved. www.apnjounals.com PRESSURE RAWOWN EQUATIONS FOR MULTIPLE-WELL

More information

Pulse Neutron Neutron (PNN) tool logging for porosity Some theoretical aspects

Pulse Neutron Neutron (PNN) tool logging for porosity Some theoretical aspects Pulse Neuton Neuton (PNN) tool logging fo poosity Some theoetical aspects Intoduction Pehaps the most citicism of Pulse Neuton Neuon (PNN) logging methods has been chage that PNN is to sensitive to the

More information

STUDY ON 2-D SHOCK WAVE PRESSURE MODEL IN MICRO SCALE LASER SHOCK PEENING

STUDY ON 2-D SHOCK WAVE PRESSURE MODEL IN MICRO SCALE LASER SHOCK PEENING Study Rev. Adv. on -D Mate. shock Sci. wave 33 (13) pessue 111-118 model in mico scale lase shock peening 111 STUDY ON -D SHOCK WAVE PRESSURE MODEL IN MICRO SCALE LASER SHOCK PEENING Y.J. Fan 1, J.Z. Zhou,

More information

OSCILLATIONS AND GRAVITATION

OSCILLATIONS AND GRAVITATION 1. SIMPLE HARMONIC MOTION Simple hamonic motion is any motion that is equivalent to a single component of unifom cicula motion. In this situation the velocity is always geatest in the middle of the motion,

More information

Axisymmetric Stokes Flow past a Swarm of Porous Cylindrical Shells

Axisymmetric Stokes Flow past a Swarm of Porous Cylindrical Shells Jounal of Applied Fluid Mechanics Vol. 9 No. pp. 957-963 06. Available online at www.jafmonline.net ISSN 735-357 EISSN 735-365. Axisymmetic Stokes Flow past a Swam of Poous Cylindical Shells S. Deo and

More information

Compactly Supported Radial Basis Functions

Compactly Supported Radial Basis Functions Chapte 4 Compactly Suppoted Radial Basis Functions As we saw ealie, compactly suppoted functions Φ that ae tuly stictly conditionally positive definite of ode m > do not exist The compact suppot automatically

More information

20th Century Atomic Theory - Hydrogen Atom

20th Century Atomic Theory - Hydrogen Atom 0th Centuy Atomic Theoy - Hydogen Atom Ruthefod s scatteing expeiments (Section.5, pp. 53-55) in 1910 led to a nuclea model of the atom whee all the positive chage and most of the mass wee concentated

More information

The condition for maximum intensity by the transmitted light in a plane parallel air film is. For an air film, μ = 1. (2-1)

The condition for maximum intensity by the transmitted light in a plane parallel air film is. For an air film, μ = 1. (2-1) hapte Two Faby--Peot ntefeomete A Faby-Peot intefeomete consists of two plane paallel glass plates A and B, sepaated by a distance d. The inne sufaces of these plates ae optically plane and thinly silveed

More information

Solution to Problem First, the firm minimizes the cost of the inputs: min wl + rk + sf

Solution to Problem First, the firm minimizes the cost of the inputs: min wl + rk + sf Econ 0A Poblem Set 4 Solutions ue in class on Tu 4 Novembe. No late Poblem Sets accepted, so! This Poblem set tests the knoledge that ou accumulated mainl in lectues 5 to 9. Some of the mateial ill onl

More information

TheWaveandHelmholtzEquations

TheWaveandHelmholtzEquations TheWaveandHelmholtzEquations Ramani Duaiswami The Univesity of Mayland, College Pak Febuay 3, 2006 Abstact CMSC828D notes (adapted fom mateial witten with Nail Gumeov). Wok in pogess 1 Acoustic Waves 1.1

More information

Analysis of simple branching trees with TI-92

Analysis of simple branching trees with TI-92 Analysis of simple banching tees with TI-9 Dušan Pagon, Univesity of Maibo, Slovenia Abstact. In the complex plane we stat at the cente of the coodinate system with a vetical segment of the length one

More information

On the Sun s Electric-Field

On the Sun s Electric-Field On the Sun s Electic-Field D. E. Scott, Ph.D. (EE) Intoduction Most investigatos who ae sympathetic to the Electic Sun Model have come to agee that the Sun is a body that acts much like a esisto with a

More information

A dual-reciprocity boundary element method for axisymmetric thermoelastodynamic deformations in functionally graded solids

A dual-reciprocity boundary element method for axisymmetric thermoelastodynamic deformations in functionally graded solids APCOM & ISCM 11-14 th Decembe, 013, Singapoe A dual-ecipocity bounday element method fo axisymmetic themoelastodynamic defomations in functionally gaded solids *W. T. Ang and B. I. Yun Division of Engineeing

More information

15 Solving the Laplace equation by Fourier method

15 Solving the Laplace equation by Fourier method 5 Solving the Laplace equation by Fouie method I aleady intoduced two o thee dimensional heat equation, when I deived it, ecall that it taes the fom u t = α 2 u + F, (5.) whee u: [0, ) D R, D R is the

More information

Entropy and Free Energy: Predicting the direction of spontaneous change The approach to Chemical equilibrium

Entropy and Free Energy: Predicting the direction of spontaneous change The approach to Chemical equilibrium Lectue 8-9 Entopy and Fee Enegy: Pedicting the diection of spontaneous change The appoach to Chemical equilibium Absolute entopy and the thid law of themodynamics To define the entopy of a compound in

More information

PES 3950/PHYS 6950: Homework Assignment 6

PES 3950/PHYS 6950: Homework Assignment 6 PES 3950/PHYS 6950: Homewok Assignment 6 Handed out: Monday Apil 7 Due in: Wednesday May 6, at the stat of class at 3:05 pm shap Show all woking and easoning to eceive full points. Question 1 [5 points]

More information

Graphs of Sine and Cosine Functions

Graphs of Sine and Cosine Functions Gaphs of Sine and Cosine Functions In pevious sections, we defined the tigonometic o cicula functions in tems of the movement of a point aound the cicumfeence of a unit cicle, o the angle fomed by the

More information

Fluid flow in curved geometries: Mathematical Modeling and Applications

Fluid flow in curved geometries: Mathematical Modeling and Applications Fluid flow in cuved geometies: Mathematical Modeling and Applications D. Muhammad Sajid Theoetical Plasma Physics Division PINSTECH, P.O. Niloe, PAEC, Islamabad Mach 01-06, 010 Islamabad, Paistan Pesentation

More information

Introduction to Arrays

Introduction to Arrays Intoduction to Aays Page 1 Intoduction to Aays The antennas we have studied so fa have vey low diectivity / gain. While this is good fo boadcast applications (whee we want unifom coveage), thee ae cases

More information

Kunming, , R.P. China. Kunming, , R.P. China. *Corresponding author: Jianing He

Kunming, , R.P. China. Kunming, , R.P. China. *Corresponding author: Jianing He Applied Mechanics and Mateials Online: 2014-04-28 ISSN: 1662-7482, Vol. 540, pp 92-95 doi:10.4028/www.scientific.net/amm.540.92 2014 Tans Tech Publications, Switzeland Reseach on Involute Gea Undecutting

More information

, and the curve BC is symmetrical. Find also the horizontal force in x-direction on one side of the body. h C

, and the curve BC is symmetrical. Find also the horizontal force in x-direction on one side of the body. h C Umeå Univesitet, Fysik 1 Vitaly Bychkov Pov i teknisk fysik, Fluid Dynamics (Stömningsläa), 2013-05-31, kl 9.00-15.00 jälpmedel: Students may use any book including the textbook Lectues on Fluid Dynamics.

More information

Analytical Solutions for Confined Aquifers with non constant Pumping using Computer Algebra

Analytical Solutions for Confined Aquifers with non constant Pumping using Computer Algebra Poceedings of the 006 IASME/SEAS Int. Conf. on ate Resouces, Hydaulics & Hydology, Chalkida, Geece, May -3, 006 (pp7-) Analytical Solutions fo Confined Aquifes with non constant Pumping using Compute Algeba

More information

THE INFLUENCE OF THE MAGNETIC NON-LINEARITY ON THE MAGNETOSTATIC SHIELDS DESIGN

THE INFLUENCE OF THE MAGNETIC NON-LINEARITY ON THE MAGNETOSTATIC SHIELDS DESIGN THE INFLUENCE OF THE MAGNETIC NON-LINEARITY ON THE MAGNETOSTATIC SHIELDS DESIGN LIVIU NEAMŢ 1, ALINA NEAMŢ, MIRCEA HORGOŞ 1 Key wods: Magnetostatic shields, Magnetic non-lineaity, Finite element method.

More information

8 Separation of Variables in Other Coordinate Systems

8 Separation of Variables in Other Coordinate Systems 8 Sepaation of Vaiables in Othe Coodinate Systems Fo the method of sepaation of vaiables to succeed you need to be able to expess the poblem at hand in a coodinate system in which the physical boundaies

More information

Steady State and Transient Performance Analysis of Three Phase Induction Machine using MATLAB Simulations

Steady State and Transient Performance Analysis of Three Phase Induction Machine using MATLAB Simulations Intenational Jounal of Recent Tends in Engineeing, Vol, No., May 9 Steady State and Tansient Pefomance Analysis of Thee Phase Induction Machine using MATAB Simulations Pof. Himanshu K. Patel Assistant

More information

INFLUENCE OF GROUND INHOMOGENEITY ON WIND INDUCED GROUND VIBRATIONS. Abstract

INFLUENCE OF GROUND INHOMOGENEITY ON WIND INDUCED GROUND VIBRATIONS. Abstract INFLUENCE OF GROUND INHOMOGENEITY ON WIND INDUCED GROUND VIBRATIONS Mohammad Mohammadi, National Cente fo Physical Acoustics, Univesity of Mississippi, MS Caig J. Hicey, National Cente fo Physical Acoustics,

More information

STABILITY AND PARAMETER SENSITIVITY ANALYSES OF AN INDUCTION MOTOR

STABILITY AND PARAMETER SENSITIVITY ANALYSES OF AN INDUCTION MOTOR HUNGARIAN JOURNAL OF INDUSTRY AND CHEMISTRY VESZPRÉM Vol. 42(2) pp. 109 113 (2014) STABILITY AND PARAMETER SENSITIVITY ANALYSES OF AN INDUCTION MOTOR ATTILA FODOR 1, ROLAND BÁLINT 1, ATTILA MAGYAR 1, AND

More information

Physics 235 Chapter 5. Chapter 5 Gravitation

Physics 235 Chapter 5. Chapter 5 Gravitation Chapte 5 Gavitation In this Chapte we will eview the popeties of the gavitational foce. The gavitational foce has been discussed in geat detail in you intoductoy physics couses, and we will pimaily focus

More information

Chem 453/544 Fall /08/03. Exam #1 Solutions

Chem 453/544 Fall /08/03. Exam #1 Solutions Chem 453/544 Fall 3 /8/3 Exam # Solutions. ( points) Use the genealized compessibility diagam povided on the last page to estimate ove what ange of pessues A at oom tempeatue confoms to the ideal gas law

More information

A Most Useful Device of Studying Electrode Processes: The Rotating Disk Electrode

A Most Useful Device of Studying Electrode Processes: The Rotating Disk Electrode A Most Useful Device of Studying Electode Pocesses: The Rotating Disk Electode the theoetical basis Soma Vesztegom Laboatoy of Electochemisty & Electoanalytical Chemisty Eötvös Loánd Univesity of Budapest

More information

NATURAL CONVECTION HEAT TRANSFER WITHIN VERTICALLY ECCENTRIC DOMED SKYLIGHTS CAVITIES

NATURAL CONVECTION HEAT TRANSFER WITHIN VERTICALLY ECCENTRIC DOMED SKYLIGHTS CAVITIES Poceedings: Building Simulation 007 NATURAL CONVECTION HEAT TRANSFER WITHIN VERTICALLY ECCENTRIC DOMED SKYLIGHTS CAVITIES A. Satipi, A. Laouadi, D. Naylo 3, R. Dhib 4 Depatment of Mechanical and Industial

More information

Section 8.2 Polar Coordinates

Section 8.2 Polar Coordinates Section 8. Pola Coodinates 467 Section 8. Pola Coodinates The coodinate system we ae most familia with is called the Catesian coodinate system, a ectangula plane divided into fou quadants by the hoizontal

More information

DonnishJournals

DonnishJournals DonnishJounals 041-1189 Donnish Jounal of Educational Reseach and Reviews. Vol 1(1) pp. 01-017 Novembe, 014. http:///dje Copyight 014 Donnish Jounals Oiginal Reseach Pape Vecto Analysis Using MAXIMA Savaş

More information

Central Coverage Bayes Prediction Intervals for the Generalized Pareto Distribution

Central Coverage Bayes Prediction Intervals for the Generalized Pareto Distribution Statistics Reseach Lettes Vol. Iss., Novembe Cental Coveage Bayes Pediction Intevals fo the Genealized Paeto Distibution Gyan Pakash Depatment of Community Medicine S. N. Medical College, Aga, U. P., India

More information

On the radial derivative of the delta distribution

On the radial derivative of the delta distribution On the adial deivative of the delta distibution Fed Backx, Fank Sommen & Jasson Vindas Depatment of Mathematical Analysis, Faculty of Engineeing and Achitectue, Ghent Univesity Depatment of Mathematics,

More information

ELASTIC ANALYSIS OF CIRCULAR SANDWICH PLATES WITH FGM FACE-SHEETS

ELASTIC ANALYSIS OF CIRCULAR SANDWICH PLATES WITH FGM FACE-SHEETS THE 9 TH INTERNATIONAL CONFERENCE ON COMPOSITE MATERIALS ELASTIC ANALYSIS OF CIRCULAR SANDWICH PLATES WITH FGM FACE-SHEETS R. Sbulati *, S. R. Atashipou Depatment of Civil, Chemical and Envionmental Engineeing,

More information

Math 124B February 02, 2012

Math 124B February 02, 2012 Math 24B Febuay 02, 202 Vikto Gigoyan 8 Laplace s equation: popeties We have aleady encounteed Laplace s equation in the context of stationay heat conduction and wave phenomena. Recall that in two spatial

More information

A method for solving dynamic problems for cylindrical domains

A method for solving dynamic problems for cylindrical domains Tansactions of NAS of Azebaijan, Issue Mechanics, 35 (7), 68-75 (016). Seies of Physical-Technical and Mathematical Sciences. A method fo solving dynamic poblems fo cylindical domains N.B. Rassoulova G.R.

More information

A NEW VARIABLE STIFFNESS SPRING USING A PRESTRESSED MECHANISM

A NEW VARIABLE STIFFNESS SPRING USING A PRESTRESSED MECHANISM Poceedings of the ASME 2010 Intenational Design Engineeing Technical Confeences & Computes and Infomation in Engineeing Confeence IDETC/CIE 2010 August 15-18, 2010, Monteal, Quebec, Canada DETC2010-28496

More information

Why Professor Richard Feynman was upset solving the Laplace equation for spherical waves? Anzor A. Khelashvili a)

Why Professor Richard Feynman was upset solving the Laplace equation for spherical waves? Anzor A. Khelashvili a) Why Pofesso Richad Feynman was upset solving the Laplace equation fo spheical waves? Anzo A. Khelashvili a) Institute of High Enegy Physics, Iv. Javakhishvili Tbilisi State Univesity, Univesity St. 9,

More information

The Schwartzchild Geometry

The Schwartzchild Geometry UNIVERSITY OF ROCHESTER The Schwatzchild Geomety Byon Osteweil Decembe 21, 2018 1 INTRODUCTION In ou study of geneal elativity, we ae inteested in the geomety of cuved spacetime in cetain special cases

More information

I. CONSTRUCTION OF THE GREEN S FUNCTION

I. CONSTRUCTION OF THE GREEN S FUNCTION I. CONSTRUCTION OF THE GREEN S FUNCTION The Helmohltz equation in 4 dimensions is 4 + k G 4 x, x = δ 4 x x. In this equation, G is the Geen s function and 4 efes to the dimensionality. In the vey end,

More information

General Railgun Function

General Railgun Function Geneal ailgun Function An electomagnetic ail gun uses a lage Loentz foce to fie a pojectile. The classic configuation uses two conducting ails with amatue that fits between and closes the cicuit between

More information

KOEBE DOMAINS FOR THE CLASSES OF FUNCTIONS WITH RANGES INCLUDED IN GIVEN SETS

KOEBE DOMAINS FOR THE CLASSES OF FUNCTIONS WITH RANGES INCLUDED IN GIVEN SETS Jounal of Applied Analysis Vol. 14, No. 1 2008), pp. 43 52 KOEBE DOMAINS FOR THE CLASSES OF FUNCTIONS WITH RANGES INCLUDED IN GIVEN SETS L. KOCZAN and P. ZAPRAWA Received Mach 12, 2007 and, in evised fom,

More information

FORMULA FOR SPIRAL CLOUD-RAIN BANDS OF A TROPICAL CYCLONE. Boris S. Yurchak* University of Maryland, Baltimore County, Greenbelt, Maryland =, (2)

FORMULA FOR SPIRAL CLOUD-RAIN BANDS OF A TROPICAL CYCLONE. Boris S. Yurchak* University of Maryland, Baltimore County, Greenbelt, Maryland =, (2) PF.3 FORMULA FOR SPIRAL CLOUD-RAIN BANDS OF A TROPICAL CYCLONE Bois S. Yuchak* Univesity of Mayland, Baltimoe County, Geenbelt, Mayland. INTRODUCTION Spial cloud-ain bands (SCRB) ae some of the most distinguishing

More information

Rigid Body Dynamics 2. CSE169: Computer Animation Instructor: Steve Rotenberg UCSD, Winter 2018

Rigid Body Dynamics 2. CSE169: Computer Animation Instructor: Steve Rotenberg UCSD, Winter 2018 Rigid Body Dynamics 2 CSE169: Compute Animation nstucto: Steve Rotenbeg UCSD, Winte 2018 Coss Poduct & Hat Opeato Deivative of a Rotating Vecto Let s say that vecto is otating aound the oigin, maintaining

More information

Stress Intensity Factor

Stress Intensity Factor S 47 Factue Mechanics http://imechanicaog/node/7448 Zhigang Suo Stess Intensity Facto We have modeled a body by using the linea elastic theoy We have modeled a cack in the body by a flat plane, and the

More information

Velocimetry Techniques and Instrumentation

Velocimetry Techniques and Instrumentation AeE 344 Lectue Notes Lectue # 05: elocimety Techniques and Instumentation D. Hui Hu Depatment of Aeospace Engineeing Iowa State Univesity Ames, Iowa 500, U.S.A Methods to Measue Local Flow elocity - Mechanical

More information

4. Electrodynamic fields

4. Electrodynamic fields 4. Electodynamic fields D. Rakhesh Singh Kshetimayum 1 4.1 Intoduction Electodynamics Faaday s law Maxwell s equations Wave equations Lenz s law Integal fom Diffeential fom Phaso fom Bounday conditions

More information

To Feel a Force Chapter 7 Static equilibrium - torque and friction

To Feel a Force Chapter 7 Static equilibrium - torque and friction To eel a oce Chapte 7 Chapte 7: Static fiction, toque and static equilibium A. Review of foce vectos Between the eath and a small mass, gavitational foces of equal magnitude and opposite diection act on

More information

Math 2263 Solutions for Spring 2003 Final Exam

Math 2263 Solutions for Spring 2003 Final Exam Math 6 Solutions fo Sping Final Exam ) A staightfowad appoach to finding the tangent plane to a suface at a point ( x, y, z ) would be to expess the cuve as an explicit function z = f ( x, y ), calculate

More information

7.2. Coulomb s Law. The Electric Force

7.2. Coulomb s Law. The Electric Force Coulomb s aw Recall that chaged objects attact some objects and epel othes at a distance, without making any contact with those objects Electic foce,, o the foce acting between two chaged objects, is somewhat

More information

Chapter Introduction to Finite Element Methods

Chapter Introduction to Finite Element Methods Chapte 1.4 Intoduction to Finite Element Methods Afte eading this chapte, you should e ale to: 1. Undestand the asics of finite element methods using a one-dimensional polem. In the last fifty yeas, the

More information

Force between two parallel current wires and Newton s. third law

Force between two parallel current wires and Newton s. third law Foce between two paallel cuent wies and Newton s thid law Yannan Yang (Shanghai Jinjuan Infomation Science and Technology Co., Ltd.) Abstact: In this pape, the essence of the inteaction between two paallel

More information

arxiv: v2 [quant-ph] 30 May 2013

arxiv: v2 [quant-ph] 30 May 2013 Entanglement in the Gove s Seach Algoithm Shantanav Chakaboty, Subhashish Banejee, Satyabata Adhikai, and Atul Kuma Indian Institute of Technology Jodhpu, Jodhpu-34011, India axiv:13054454v [quant-ph]

More information

arxiv: v1 [physics.pop-ph] 3 Jun 2013

arxiv: v1 [physics.pop-ph] 3 Jun 2013 A note on the electostatic enegy of two point chages axiv:1306.0401v1 [physics.pop-ph] 3 Jun 013 A C Tot Instituto de Física Univesidade Fedeal do io de Janeio Caixa Postal 68.58; CEP 1941-97 io de Janeio,

More information

2 Governing Equations

2 Governing Equations 2 Govening Equations This chapte develops the govening equations of motion fo a homogeneous isotopic elastic solid, using the linea thee-dimensional theoy of elasticity in cylindical coodinates. At fist,

More information

1) Consider an object of a parabolic shape with rotational symmetry z

1) Consider an object of a parabolic shape with rotational symmetry z Umeå Univesitet, Fysik 1 Vitaly Bychkov Pov i teknisk fysik, Fluid Mechanics (Stömningsläa), 01-06-01, kl 9.00-15.00 jälpmedel: Students may use any book including the tetbook Lectues on Fluid Dynamics.

More information

Application of homotopy perturbation method to the Navier-Stokes equations in cylindrical coordinates

Application of homotopy perturbation method to the Navier-Stokes equations in cylindrical coordinates Computational Ecology and Softwae 5 5(): 9-5 Aticle Application of homotopy petubation method to the Navie-Stokes equations in cylindical coodinates H. A. Wahab Anwa Jamal Saia Bhatti Muhammad Naeem Muhammad

More information

Surveillance Points in High Dimensional Spaces

Surveillance Points in High Dimensional Spaces Société de Calcul Mathématique SA Tools fo decision help since 995 Suveillance Points in High Dimensional Spaces by Benad Beauzamy Januay 06 Abstact Let us conside any compute softwae, elying upon a lage

More information

MODULE 5 ADVANCED MECHANICS GRAVITATIONAL FIELD: MOTION OF PLANETS AND SATELLITES VISUAL PHYSICS ONLINE

MODULE 5 ADVANCED MECHANICS GRAVITATIONAL FIELD: MOTION OF PLANETS AND SATELLITES VISUAL PHYSICS ONLINE VISUAL PHYSICS ONLIN MODUL 5 ADVANCD MCHANICS GRAVITATIONAL FILD: MOTION OF PLANTS AND SATLLITS SATLLITS: Obital motion of object of mass m about a massive object of mass M (m

More information

( ) [ ] [ ] [ ] δf φ = F φ+δφ F. xdx.

( ) [ ] [ ] [ ] δf φ = F φ+δφ F. xdx. 9. LAGRANGIAN OF THE ELECTROMAGNETIC FIELD In the pevious section the Lagangian and Hamiltonian of an ensemble of point paticles was developed. This appoach is based on a qt. This discete fomulation can

More information

Effect of drag on the performance for an efficient wind turbine blade design

Effect of drag on the performance for an efficient wind turbine blade design Available online at www.sciencediect.com Enegy Pocedia 18 (01 ) 404 415 Abstact Effect of dag on the pefomance fo an efficient wind tubine blade design D. Eng. Ali H. Almukhta Univesity of Technology Email-

More information

A Relativistic Electron in a Coulomb Potential

A Relativistic Electron in a Coulomb Potential A Relativistic Electon in a Coulomb Potential Alfed Whitehead Physics 518, Fall 009 The Poblem Solve the Diac Equation fo an electon in a Coulomb potential. Identify the conseved quantum numbes. Specify

More information

A Backward Identification Problem for an Axis-Symmetric Fractional Diffusion Equation

A Backward Identification Problem for an Axis-Symmetric Fractional Diffusion Equation Mathematical Modelling and Analysis Publishe: Taylo&Fancis and VGTU Volume 22 Numbe 3, May 27, 3 32 http://www.tandfonline.com/tmma https://doi.og/.3846/3926292.27.39329 ISSN: 392-6292 c Vilnius Gediminas

More information

LET a random variable x follows the two - parameter

LET a random variable x follows the two - parameter INTERNATIONAL JOURNAL OF MATHEMATICS AND SCIENTIFIC COMPUTING ISSN: 2231-5330, VOL. 5, NO. 1, 2015 19 Shinkage Bayesian Appoach in Item - Failue Gamma Data In Pesence of Pio Point Guess Value Gyan Pakash

More information

From Gravitational Collapse to Black Holes

From Gravitational Collapse to Black Holes Fom Gavitational Collapse to Black Holes T. Nguyen PHY 391 Independent Study Tem Pape Pof. S.G. Rajeev Univesity of Rocheste Decembe 0, 018 1 Intoduction The pupose of this independent study is to familiaize

More information

DIRECT AND INVERSE PROBLEMS FOR THE EVERSION OF A SPHERICAL SHELL

DIRECT AND INVERSE PROBLEMS FOR THE EVERSION OF A SPHERICAL SHELL THE PUBLISHING HOUSE POCEEDINGS OF THE OMANIAN ACADEMY, Seies OF THE OMANIAN ACADEMY Volume, Numbe /, pp. 33 DIECT AND INVESE POBLEMS FO THE EVESION OF A SPHEICAL SHELL aluca I. OŞC Coina S. DAPACA Depatment

More information

Solution of a Spherically Symmetric Static Problem of General Relativity for an Elastic Solid Sphere

Solution of a Spherically Symmetric Static Problem of General Relativity for an Elastic Solid Sphere Applied Physics eseach; Vol. 9, No. 6; 7 ISSN 96-969 E-ISSN 96-9647 Published by Canadian Cente of Science and Education Solution of a Spheically Symmetic Static Poblem of Geneal elativity fo an Elastic

More information

Gauss s Law Simulation Activities

Gauss s Law Simulation Activities Gauss s Law Simulation Activities Name: Backgound: The electic field aound a point chage is found by: = kq/ 2 If thee ae multiple chages, the net field at any point is the vecto sum of the fields. Fo a

More information

Duality between Statical and Kinematical Engineering Systems

Duality between Statical and Kinematical Engineering Systems Pape 00, Civil-Comp Ltd., Stiling, Scotland Poceedings of the Sixth Intenational Confeence on Computational Stuctues Technology, B.H.V. Topping and Z. Bittna (Editos), Civil-Comp Pess, Stiling, Scotland.

More information

We are IntechOpen, the world s leading publisher of Open Access books Built by scientists, for scientists. International authors and editors

We are IntechOpen, the world s leading publisher of Open Access books Built by scientists, for scientists. International authors and editors We ae IntechOpen, the old s leading publishe of Open Access books Built by scientists, fo scientists 3,7 8,.7 M Open access books available Intenational authos and editos Donloads Ou authos ae among the

More information

Introduction to Nuclear Forces

Introduction to Nuclear Forces Intoduction to Nuclea Foces One of the main poblems of nuclea physics is to find out the natue of nuclea foces. Nuclea foces diffe fom all othe known types of foces. They cannot be of electical oigin since

More information

B. Spherical Wave Propagation

B. Spherical Wave Propagation 11/8/007 Spheical Wave Popagation notes 1/1 B. Spheical Wave Popagation Evey antenna launches a spheical wave, thus its powe density educes as a function of 1, whee is the distance fom the antenna. We

More information

Radial Inflow Experiment:GFD III

Radial Inflow Experiment:GFD III Radial Inflow Expeiment:GFD III John Mashall Febuay 6, 003 Abstact We otate a cylinde about its vetical axis: the cylinde has a cicula dain hole in the cente of its bottom. Wate entes at a constant ate

More information

Asymptotic Model of the Flow in Metallic and Elastomeric Progressive Cavity Pump

Asymptotic Model of the Flow in Metallic and Elastomeric Progressive Cavity Pump Poceedings of COBEM 009 Copyight c 009 by ABCM 0th Intenational Congess of Mechanical Engineeing Novembe 15-0, 009, Gamado, RS, Bazil Asymptotic Model of the Flow in Metallic and Elastomeic Pogessive Cavity

More information

A PERTURBATION EXPANSION APPROACH TO SOLVING THE ELECTROMAGNETIC INDUCTION PROBLEM IN THREE DIMENSIONS

A PERTURBATION EXPANSION APPROACH TO SOLVING THE ELECTROMAGNETIC INDUCTION PROBLEM IN THREE DIMENSIONS A PERTURBATION EXPANSION APPROACH TO SOLVING THE ELECTROMAGNETIC INDUCTION PROBLEM IN THREE DIMENSIONS By Nancy H. Natek Submitted in Patial Fulfillment of the Requiements fo the Maste of Science in Geophysics

More information

1 Similarity Analysis

1 Similarity Analysis ME43A/538A/538B Axisymmetic Tubulent Jet 9 Novembe 28 Similaity Analysis. Intoduction Conside the sketch of an axisymmetic, tubulent jet in Figue. Assume that measuements of the downsteam aveage axial

More information

Lecture 8 - Gauss s Law

Lecture 8 - Gauss s Law Lectue 8 - Gauss s Law A Puzzle... Example Calculate the potential enegy, pe ion, fo an infinite 1D ionic cystal with sepaation a; that is, a ow of equally spaced chages of magnitude e and altenating sign.

More information

Improvement in Accuracy for Design of Multidielectric Layers Microstrip Patch Antenna

Improvement in Accuracy for Design of Multidielectric Layers Microstrip Patch Antenna 498 Impovement in Accuacy fo Design of Multidielectic Layes Micostip Patch Antenna Sami Dev Gupta*, Anvesh Gag and Anuag P. Saan Jaypee Institute of Infomation Technology Univesity Noida, Utta Padesh,

More information

Chemo-mechanical coupling and damage enhanced dissolution at intergranular contact

Chemo-mechanical coupling and damage enhanced dissolution at intergranular contact Numeical Models in Geomechanics-NUMOG IX-Pande & Pietuszczak (eds) 4 Taylo & Fancis Goup, London, ISBN 9 589 636 X Chemo-mechanical coupling and damage enhanced dissolution at integanula contact T. Hueckel

More information

Chapter 5. Uniform Circular Motion. a c =v 2 /r

Chapter 5. Uniform Circular Motion. a c =v 2 /r Chapte 5 Unifom Cicula Motion a c =v 2 / Unifom cicula motion: Motion in a cicula path with constant speed s v 1) Speed and peiod Peiod, T: time fo one evolution Speed is elated to peiod: Path fo one evolution:

More information

Casimir-Polder potential for parallel metallic plates in background of a conical defect

Casimir-Polder potential for parallel metallic plates in background of a conical defect Amenian Jounal of Physics, 7, vol., issue, pp. 8-9 Casimi-Polde potential fo paallel metallic plates in backgound of a conical defect A.Kh. Gigoyan Institute of Applied Poblems in Physics NAS RA, 5 Nesessian

More information