Mathematical Inverse Problem of Magnetic Field from Exponentially Varying Conductive Ground
|
|
- Gwen Dalton
- 5 years ago
- Views:
Transcription
1 Applied Mathematical Sciences, Vol. 6, 212, no. 113, Mathematical Inverse Problem of Magnetic Field from Exponentially Varying Conductive Ground Warin Sripanya Faculty of Science and Technology, Nakhon Pathom Rajabhat University Nakhon Pathom 73, Thailand Suabsagun Yooyuanyong Department of Mathematics, Faculty of Science Silpakorn University, Nakhon Pathom 73, Thailand Abstract We derive an analytical solution of the steady state magnetic field due to a direct current source on a multilayered earth with a layer having exponentially varying conductivity. Our variation in conductivity is realistic and can be generalized to all cases of exponential profiles. The Hankel transform is introduced to our problem and analytical result is obtained. Our solution is achieved by solving a boundary value problem in the wave number domain and then transforming the solution back to the spatial domain. The effects of magnetic fields are plotted and compared to show the behavior in response to different ground structures at many depths while some parameters are approximately given. The curves of magnetic fields show some significance to the variation of conductivity. An inverse problem via the use of the Levenberg-Marquardt optimization technique is introduced for finding the conductivity parameters of the ground. The optimal result of our model is close to the true value with percentage errors of our two conductivity parameters less than.9% and 1.7%. Mathematics Subject Classification: 86A25 Keywords: inverse problem, magnetic field, direct current, exponentially varying conductivity Corresponding author.
2 564 W. Sripanya and S. Yooyuanyong 1 Introduction Many authors have investigated the nature of the resistivity response resulting from a direct current source on a heterogeneous ground whose electrical conductivity varies exponentially with depth. Paul and Banerjee [9] studied the problem of computing apparent resistivity for a half-space model. Similar results, but for a layered earth, were reported by Banerjee et al. [2], Kim and Lee [7], Stoyer and Wait [12]. Throughout these investigations, the electrical conductivity was assumed, for simplicity, to be exponentially dependent upon depth, denoted by σ (z) = α exp (βz), where α and β are the parameters that define the conductivity profile. Unfortunately, these works are not sufficiently general about the variation of conductivity to be used for many applications. This approach is realistic only for a very short range of depth. Fortuitously, Yooyuanyong [14] introduced a new approach for enhanced investigation of a varying conductive ground structure. The general form of the variation of conductivity was developed and defined by σ (z) = σ + (σ 1 σ ) exp ( βz), where β, σ and σ 1 are positive real number parameters that describe the conductivity profile. In this approach, the electrical conductivity varies exponentially and tends to a basement value σ as z tends to infinity. Chumchob [6] used this conductivity variation to study the problem of computing electromagnetic response for a multilayered earth structure. In this paper, the electrical exploration method based on the measurement of static magnetic fields associated with noninductive current flow (Sripanya and Yooyuanyong [11]) is applied to investigate the structure of earth. We develop an analytical solution of the steady state magnetic field for the problem of a horizontally stratified layered earth with a layer having exponentially varying conductivity presented by Yooyuanyong [14]. The Hankel transform is introduced to our problem and analytical result is obtained. The inversion process, using the Levenberg-Marquardt algorithm, is conducted to estimate the conductivity parameters of the ground. 2 Model and Basic Equations In our geometric model, a point source of direct current I is located at the interface between two half-spaces. The half-space above the interface (z < ) is the region of air with conductivity approximately equal to zero, whereas the half-space below the interface (z > ) is an n-layered horizontally stratified earth with depths to the layers h 1, h 2,..., h n 1 (the lowermost layer extending to infinity) measured from the ground surface, where n 2 is an integer. Each layer has conductivity as a function of depth, i.e., σ k (z) for layer 1 k n.
3 Mathematical inverse problem of magnetic field Magnetic Field Resulting from a DC Source The general steady state Maxwell s equations in the frequency domain (Chen and Oldenburg [4, 5]) can be used to determine the magnetic field for this problem, namely E =, H = σe, (1) where E is the vector electric field, H is the vector magnetic field and σ is the conductivity of the medium. Since the problem is axisymmetric and H has only the azimuthal component in cylindrical coordinates, for simplicity, we use H to represent the azimuthal component in the following derivations, and we now have 2 H z 2 + σ z ( ) 1 H σ z + 2 H r r H r 1 H =. (2) r2 The Hankel transform (Ali and Kalla [1]) is introduced and defined by and H (λ, z) = H (r, z) = λrh (r, z) J 1 (λr) dr (3) H (λ, z) J 1 (λr) dλ, (4) where J 1 is the Bessel function of the first kind of order one. Taking the transformation on both sides of equation (2), we obtain 2 H z 2 + σ z ( ) 1 H σ z λ2 H =. (5) 2.2 Boundary Conditions The magnetic field in each layer can be obtained by taking the inverse Hankel transform to the solution of equation (5), which satisfies the following boundary conditions: 1. The azimuthal component of the magnetic field needs to be continuous on each of the boundary planes in the earth, i.e., lim z h k H k = lim z h + k H k+1. (6) 2. The radial component of the electric field needs to be continuous on each of the boundary planes in the earth, i.e., lim z h k 1 H k σ k z = lim z h + k 1 H k+1. (7) σ k+1 z
4 5642 W. Sripanya and S. Yooyuanyong 3. The vertical component of the current density must be zero at the free surface except in an infinitesimal neighborhood around a current source (Kim and Lee [7]), i.e., 1 r r (rh 1) = I z= 2π δ (r) r. (8) 4. The azimuthal component of the magnetic field tends to zero as z tends to infinity, i.e., H n =. (9) lim z 3 Response of a Layer Having Exponentially Varying Conductivity For an exponentially varying conductive layer k, where 1 k n and n 2, the variation of conductivity is denoted by σ k (z) = c k + (a k c k ) e b kz, (1) where a k, b k and c k are positive real numbers. Hence, the equation for the magnetic field in each layer can be simplified by substituting equation (1) into (5) and we obtain 2 Hk + b k (σ k c k ) z 2 σ k The solution to the above equation is H k (λ, z) = ψ k ( Ak e ζ 1,k 2 F 1 (α 1,k, β 1,k ; γ 1,k ; ξ k ) H k z λ2 Hk =. (11) + B k e ζ 2,k 2 F 1 (α 2,k, β 2,k ; γ 2,k ; ξ k ) ), (12) where ϱ k α i,k = 3 2 λ + ( 1) i, β i,k = 3 b k 2b k 2 + λ + ( 1) i ϱ k, b k 2b k γ i,k = 1 + ( 1) i ϱ k, ϱ k = b 2k b + 4λ2, ξ k = c ke b kz, k c k a k ψ k = ( ( a k + c k e b k z 1 )) ( ) 2 z, ζi,k = ( 1) i ϱ k b k 2 (13) and 2 F 1 is the ordinary hypergeometric function. The unknown coefficients A k and B k are arbitrary constants, which can be determined by using the
5 Mathematical inverse problem of magnetic field 5643 propagator matrix technique applied to the above boundary conditions (see Sripanya and Yooyuanyong [11]). Thus, the magnetic field in layer k is H k (r, z) = ψ k ( Ak e ζ 1,k 2 F 1 (α 1,k, β 1,k ; γ 1,k ; ξ k ) + B k e ζ 2,k 2 F 1 (α 2,k, β 2,k ; γ 2,k ; ξ k ) ) J 1 (λr) dλ. (14) In the case of an exponentially varying conductive half-space, the magnetic field can be written as H (r, z) = I 2π 2F 1 (α, β; γ; ξ) ψ e ζ J 1 (λr) dλ, (15) 2F 1 (α, β; γ; ξ ) ψ where all the variables are given as above without subscript k while i = 1, and the subscript refers to z =. 4 Numerical Experiments In our numerical experiments, we calculate the magnetic fields due to a direct current source on two types of electrically conductive earth structures. Both of the example ground models are the heterogeneous conductive half-spaces. The first model has exponentially increasing conductivity with depth, whereas the second model has exponentially decreasing conductivity. The values of the model parameters are given in Table 1. Chave s algorithm [3] is used for numerically calculating the inverse Hankel transform of the magnetic field solutions. The special functions are computed by using the Numerical Recipes source codes (Press et al. [1]). The results from both models are plotted to show the behavior of magnetic fields against source-receiver spacing r at different depths z =,.5, 1,..., 1 metres as shown in Figure 1. Table 1: Model parameters used in our sample tests. Parameters Model a ( S m 1) b ( m 1) c ( S m 1)
6 5644 W. Sripanya and S. Yooyuanyong.16 Model 1.14 Magnetic Field Intensity (A/m) Source Receiver Spacing (m).16 Model 2.14 Magnetic Field Intensity (A/m) Source Receiver Spacing (m) Figure 1: Behavior of magnetic fields from synthesis models in our sample tests. 5 Inversion Process In our inverse model example, we simulate the reflection of magnetic radiation data from our forward model of practical interest. The magnetic fields are generated by the forward problem of the first model in our sample tests. Random errors up to 3% are superimposed on the scaled magnetic fields to simulate the set of real data. The iterative procedure using the Levenberg-Marquardt method (Press et al. [1]) is applied to estimate the model parameters of conductivity variation. The model parameter a is a conductivity of the earth s surface, which can be assumed to be known from the measurement. We start the iterative process to find the values of the conductivity parameters with initial guess values b =.1 m 1 and c =.1 S m 1. The inversion method leads to the optimal values of the parameters b and c in which their percentage errors are less than only 1.7% and.9%, respectively, after using 138 iterations. The graphs of the true and estimated conductivity models are plotted as shown in Figure 2.
7 Mathematical inverse problem of magnetic field Conductivity (S/m) True Model Estimated Model Initial Guess Model Depth (m) Figure 2: Graphs of conductivity σ against depth z for our inverse model example. 6 Discussions and Conclusions Analytical solution of the steady state magnetic field for exponentially conductive ground profile is derived. Our variation in conductivity is realistic and can be generalized to all cases of exponential profiles. The effects of magnetic fields are plotted and compared to show the behavior in response to different ground structures at many depths (see Figure 1). The magnetic curves of the model having exponentially increasing conductivity with depth (Model 1) are quite different from the model having exponentially decreasing conductivity (Model 2). The magnitude of magnetic field from the first model is much higher than the magnetic field from the second model at the same depth and source-receiver spacing. We observe that if the conductive ground has a high electrical conductivity, it will lead to the large sized magnitude of magnetic field. This means that the differences of curves for our models are depended on the variation of conductivity. An inverse problem via the use of an optimization technique is introduced for finding the conductivity parameters of the ground. The model of a simple case for the ground structure is used to investigate the conductivity profile. The iterative procedure using the Levenberg-Marquardt method is applied to estimate the model parameters of conductivity variation. The optimal result of our model converges to the true value with percentage errors of b and c less than only 1.7% and.9%, respectively, after using 138 iterations. The graphs of the true and estimated conductivity models are plotted as shown in Figure 2. We see that the graph of the estimated model is close to the true model. The inversion method leads to very good result and has high speed of convergence. This illustrates the advantage in using Levenberg-Marquardt method which gives the result much better than using another method of inversion (e.g.,
8 5646 W. Sripanya and S. Yooyuanyong Oldenburg [8], Vozoff and Jupp [13]). References [1] I. Ali and S. Kalla, A generalized Hankel transform and its use for solving certain partial differential equations, J. Austral. Math. Soc. Ser. B, 41 (1999), [2] B. Banerjee, B.J. Sengupta and B.P. Pal, Apparent resistivity of a multilayered earth with a layer having exponentiality varying conductivity, Geophys. Prosp., 28 (198), [3] A.D. Chave, Numerical integration of related Hankel transforms by quadrature and continued fraction expansion, Geophysics, 48 (1983), [4] J. Chen and D.W. Oldenburg, Magnetic and electrical fields of direct currents in a layered earth, Expl. Geophys., 35 (24), [5] J. Chen and D.W. Oldenburg, A new formula to compute apparent resistivities from marine magnetometric resistivity data, Geophysics, 71 (26), [6] N. Chumchob, Mathematical modelling of electromagnetic response on a multilayered earth with a layer having exponentially varying conductivity, Thai J. Math., 2 (24), [7] H.-S. Kim and K. Lee, Response of a multilayered earth with layers having exponentially varying resistivities, Geophysics, 61 (1996), [8] D.W. Oldenburg, The interpretation of direct current resistivity measurements, Geophysics, 43 (1978), [9] M.K. Paul and B. Banerjee, Electrical potentials due to a point source upon models of continuously varying conductivity, Pure Appl. Geophys., 8 (197), [1] W.H. Press, S.A. Teukolsky, W.T. Vetterling and B.P. Flannery, Numerical Recipes in Fortran 77: The Art of Scientific Computing, Fortran Numerical Recipes Vol. 1, 2nd ed., Cambridge Univ. Press, [11] W. Sripanya and S. Yooyuanyong, Mathematical inverse problem of magnetic field from heterogeneous earth structures, Far East J. Math. Sci. (to appear).
9 Mathematical inverse problem of magnetic field 5647 [12] C.H. Stoyer and J.R. Wait, Resistivity probing of an exponential earth with a homogeneous overburden, Geoexpl., 15 (1977), [13] K. Vozoff and D.L.B. Jupp, Joint inversion of geophysical data, Geophys. J. R. astr. Soc., 42 (1975), [14] S. Yooyuanyong, Electromagnetic response over a varying conductive ground, Songklanakarin J. Sci. Technol., 22 (2), Received: June, 212
PAijpam.eu MATHEMATICAL MODEL OF MAGNETOMETRIC RESISTIVITY SOUNDING FOR A CONDUCTIVE MEDIUM WITH A POSITIVELY SKEWED BULGE OVERBURDEN
International Journal of Pure and Applied Mathematics Volume 9 No. 2 214, 153-164 ISSN: 1311-88 printed version); ISSN: 1314-3395 on-line version) url: http://www.ijpam.eu doi: http://dx.doi.org/1.12732/ijpam.v9i2.5
More informationMathematical Modeling of Resistivity Probing of Ancient City at Prapathom Chadee, Nakhon Pathom, Thailand
Applied Mathematical Sciences, Vol. 7, 13, no. 98, 4847-4856 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/1.1988/ams.13.37353 Mathematical Modeling of Resistivity Probing of Ancient City at Prapathom
More informationMathematical 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 informationCalculation of Cylindrical Functions using Correction of Asymptotic Expansions
Universal Journal of Applied Mathematics & Computation (4), 5-46 www.papersciences.com Calculation of Cylindrical Functions using Correction of Asymptotic Epansions G.B. Muravskii Faculty of Civil and
More informationFast computation of Hankel Transform using orthonormal exponential approximation of complex kernel function
Fast computation of Hankel Transform using orthonormal exponential approximation of complex kernel function Pravin K Gupta 1, Sri Niwas 1 and Neeta Chaudhary 2 1 Department of Earth Sciences, Indian Institute
More informationinter.noise 2000 The 29th International Congress and Exhibition on Noise Control Engineering August 2000, Nice, FRANCE
Copyright SFA - InterNoise 2000 1 inter.noise 2000 The 29th International Congress and Exhibition on Noise Control Engineering 27-30 August 2000, Nice, FRANCE I-INCE Classification: 0.0 BACKSCATTERING
More informationMATHEMATICAL INVERSE PROBLEM OF ELECTRIC AND VERTICAL MAGNETIC FIELDS TO RESOLVE THE SALINE SOIL STRUCTURE
Journal of Mathematical Sciences: Advances and Applications Volume 6, Number,, Pages 57-85 MATHEMATICAL INVERSE PROBLEM OF ELECTRIC AND VERTICAL MAGNETIC FIELDS TO RESOLVE THE SALINE SOIL STRUCTURE Faculty
More informationA new formula to compute apparent resistivities from marine magnetometric resistivity data
GEOPHYSICS, VOL. 71, NO. 3 MAY-JUNE 006; P. G73 G81, 10 FIGS. 10.1190/1.194513 A new formula to compute apparent resistivities from marine magnetometric resistivity data Jiuping Chen 1 Douglas W. Oldenburg
More informationD034 On the Effects of Anisotropy in Marine CSEM
D034 On the Effects of Anisotropy in Marine CSEM L.O. Løseth* (Statoil ASA, formerly Dept. of Physics, NTNU), B. Ursin (Dept. of Appl. Geophys., NTNU) & L. Amundsen (Statoil ASA) SUMMARY In marine CSEM
More informationIntroduction to Condensed Matter Physics
Introduction to Condensed Matter Physics Diffraction I Basic Physics M.P. Vaughan Diffraction Electromagnetic waves Geometric wavefront The Principle of Linear Superposition Diffraction regimes Single
More informationUNWINDING FROM CONIC PACKAGES
DAAAM INTERNATIONAL SCIENTIFIC BOOK 2011 pp. 445-452 CHAPTER 36 UNWINDING FROM CONIC PACKAGES PRACEK, S.; SLUGA, F. & FRANKEN, G. Abstract: Mathematical modeling can be used to simulate the unwinding of
More informationAnomaly effects of arrays for 3d geoelectrical resistivity imaging using orthogonal or parallel 2d profiles
African Journal of Environmental Science and Technology Vol. 4(7), pp. 446-454, July 2010 Available online at http://www.academicjournals.org/ajest ISSN 1991-637X 2010 Academic Journals Full Length Research
More informationUniversity of Warwick institutional repository:
University of Warwick institutional repository: http://go.warwick.ac.uk/wrap This paper is made available online in accordance with publisher policies. Please scroll down to view the document itself. Please
More informationTHE SCATTERING FROM AN ELLIPTIC CYLINDER IRRADIATED BY AN ELECTROMAGNETIC WAVE WITH ARBITRARY DIRECTION AND POLARIZATION
Progress In Electromagnetics Research Letters, Vol. 5, 137 149, 2008 THE SCATTERING FROM AN ELLIPTIC CYLINDER IRRADIATED BY AN ELECTROMAGNETIC WAVE WITH ARBITRARY DIRECTION AND POLARIZATION Y.-L. Li, M.-J.
More informationNEAR FIELD REPRESENTATIONS OF THE ACOUSTIC GREEN S FUNCTION IN A SHALLOW OCEAN WITH FLUID-LIKE SEABED
Georgian Mathematical Journal Volume 4 2007, Number, 09 22 NEAR FIELD REPRESENTATIONS OF THE ACOUSTIC GREEN S FUNCTION IN A SHALLOW OCEAN WITH FLUID-LIKE SEAED ROERT GILERT AND MIAO-JUNG OU Abstract. In
More informationModified Fourier Spectral Method for Non-periodic CFD
Modified Fourier Spectral Method for Non-periodic CFD Huankun Fu Chaoqun Liu Technical Report 2011-14 http://www.uta.edu/math/preprint/ Modified Fourier Spectral Method for Non-periodic CFD Huankun Fu
More informationGreen s Tensors for the Electromagnetic Field in Cylindrical Coordinates
Chapter 25 Green s Tensors for the Electromagnetic Field in Cylindrical Coordinates The solutions of the vector Helmholtz equation in three dimensions can be expressed by a complete set of vector fields
More informationSongklanakarin Journal of Science and Technology SJST R1 KANYAMEE. Numerical methods for finding multiplicative inverses of a modulo N
Songklanakarin Journal of Science and Technology SJST-0-0.R KANYAMEE Numerical methods for finding multiplicative inverses of a modulo N Journal: Songklanakarin Journal of Science and Technology Manuscript
More informationOptimization 2. CS5240 Theoretical Foundations in Multimedia. Leow Wee Kheng
Optimization 2 CS5240 Theoretical Foundations in Multimedia Leow Wee Kheng Department of Computer Science School of Computing National University of Singapore Leow Wee Kheng (NUS) Optimization 2 1 / 38
More informationSolutions to Laplace s Equation in Cylindrical Coordinates and Numerical solutions. ρ + (1/ρ) 2 V
Solutions to Laplace s Equation in Cylindrical Coordinates and Numerical solutions Lecture 8 1 Introduction Solutions to Laplace s equation can be obtained using separation of variables in Cartesian and
More information3D MODELING OF EARTHQUAKE CYCLES OF THE XIANSHUIHE FAULT, SOUTHWESTERN CHINA
3D MODELING OF EARTHQUAKE CYCLES OF THE XIANSHUIHE FAULT, SOUTHWESTERN CHINA Li Xiaofan MEE09177 Supervisor: Bunichiro Shibazaki ABSTRACT We perform 3D modeling of earthquake generation of the Xianshuihe
More informationMagnetotelluric (MT) Method
Magnetotelluric (MT) Method Dr. Hendra Grandis Graduate Program in Applied Geophysics Faculty of Mining and Petroleum Engineering ITB Geophysical Methods Techniques applying physical laws (or theory) to
More informationA Study of the Thermal Properties of a One. Dimensional Lennard-Jones System
A Study of the Thermal Properties of a One Dimensional Lennard-Jones System Abstract In this study, the behavior of a one dimensional (1D) Lennard-Jones (LJ) system is simulated. As part of this research,
More information1 Solutions in cylindrical coordinates: Bessel functions
1 Solutions in cylindrical coordinates: Bessel functions 1.1 Bessel functions Bessel functions arise as solutions of potential problems in cylindrical coordinates. Laplace s equation in cylindrical coordinates
More informationWKB and Numerical Compound Matrix Methods for Solving the Problem of Everted Neo-Hookean Spherical Shell Morteza Sanjaranipour 1, Hamed Komeyli 2
Journal of mathematics and computer science (04), 77-90 WKB and Numerical Compound Matrix Methods for Solving the Problem of Everted Neo-Hookean Spherical Shell Morteza Sanjaranipour, Hamed Komeyli Faculty
More informationElectromagnetic Waves Retarded potentials 2. Energy and the Poynting vector 3. Wave equations for E and B 4. Plane EM waves in free space
Electromagnetic Waves 1 1. Retarded potentials 2. Energy and the Poynting vector 3. Wave equations for E and B 4. Plane EM waves in free space 1 Retarded Potentials For volume charge & current = 1 4πε
More informationNon-linear least-squares and chemical kinetics. An improved method to analyse monomer-excimer decay data
Journal of Mathematical Chemistry 21 (1997) 131 139 131 Non-linear least-squares and chemical kinetics. An improved method to analyse monomer-excimer decay data J.P.S. Farinha a, J.M.G. Martinho a and
More information1 Introduction. 2 Binary shift map
Introduction This notes are meant to provide a conceptual background for the numerical construction of random variables. They can be regarded as a mathematical complement to the article [] (please follow
More informationELECTROMAGNETIC SCATTERING FROM A CHIRAL- COATED NIHILITY CYLINDER
Progress In Electromagnetics Research Letters, Vol. 18, 41 5, 21 ELECTROMAGNETIC SCATTERING FROM A CHIRAL- COATED NIHILITY CYLINDER S. Ahmed and Q. A. Naqvi Department of Electronics Quaid-i-Azam University
More informationRadiation pattern in homogeneous and transversely isotropic attenuating media
Radiation pattern in homogeneous and transversely isotropic attenuating media Satish Sinha*, Sergey Abaseyev** and Evgeni Chesnokov** *Rajiv Gandhi Institute of Petroleum Technology, Rae Bareli, UP 229010
More informationTHEORETICAL ANALYSIS AND NUMERICAL CALCULATION OF LOOP-SOURCE TRANSIENT ELECTROMAGNETIC IMAGING
CHINESE JOURNAL OF GEOPHYSICS Vol.47, No.2, 24, pp: 379 386 THEORETICAL ANALYSIS AND NUMERICAL CALCULATION OF LOOP-SOURCE TRANSIENT ELECTROMAGNETIC IMAGING XUE Guo-Qiang LI Xiu SONG Jian-Ping GUO Wen-Bo
More informationNumerical Analysis of Plate Heat Exchanger Performance in Co-Current Fluid Flow Configuration
Numerical Analysis of Plate Heat Exchanger Performance in Co-Current Fluid Flow Configuration H. Dardour, S. Mazouz, and A. Bellagi Abstract For many industrial applications plate heat exchangers are demonstrating
More informationRobert W. Vestrum and R. James Brown
Group versus phase velocity in laboratory measurements on orthorhombic samples Robert W. Vestrum and R. James Brown ABSTRACT In laboratory measurements of traveltimes across anisotropic materials, there
More informationELECTRICAL CALCULATIONS SOFTWARE. R.Ubic 1996
Page 1 of 5 ELECTRICAL CALCULATIONS SOFTWARE R.Ubic 1996 These programmes are designed to calculate the dielectric constant of a microwave resonator from its geometry, the geometry of the test cavity,
More informationLOCAL DIFFERENTIAL QUADRATURE METHOD FOR ELLIPTIC EQUATIONS IN IRREGULAR DOMAINS
LOCAL DIFFERENTIAL QUADRATURE METHOD FOR ELLIPTIC EQUATIONS IN IRREGULAR DOMAINS L.H. Shen and D.L. Young * Department of Civil Engineering and Hydrotech Research Institute National Taiwan University Taipei,
More informationEFFICIENT 3D POISSON SOLVERS FOR SPACE-CHARGE SIMULATION
Abstract EFFICIENT 3D POISSON SOLVERS FOR SPACE-CHARGE SIMULATION Three-dimensional Poisson solver plays an important role in the self-consistent space-charge simulation. In this paper, we present several
More informationContents. 1 Basic Equations 1. Acknowledgment. 1.1 The Maxwell Equations Constitutive Relations 11
Preface Foreword Acknowledgment xvi xviii xix 1 Basic Equations 1 1.1 The Maxwell Equations 1 1.1.1 Boundary Conditions at Interfaces 4 1.1.2 Energy Conservation and Poynting s Theorem 9 1.2 Constitutive
More informationPropagation of Radio Frequency Waves Through Density Filaments
PSFC/JA-15-13 Propagation of Radio Frequency Waves Through Density Filaments A. K. Ram and K. Hizanidis a May 015 a National Technical University of Athens (part of HELLAS) School of Electrical and Computer
More informationComputational Modeling of Geoelectrical Soundings using PML-FDTD
120 Computational Modeling of Geoelectrical Soundings using PML-FDTD Luísa F. Ribeiro 1 and Marcela S. Novo 2 1,2 Federal University of Bahia UFBA, Electrical Engineering Department, Salvador BA, Brazil
More informationA Brief Introduction to Magnetotellurics and Controlled Source Electromagnetic Methods
A Brief Introduction to Magnetotellurics and Controlled Source Electromagnetic Methods Frank Morrison U.C. Berkeley With the help of: David Alumbaugh Erika Gasperikova Mike Hoversten Andrea Zirilli A few
More informationA Symbolic Numeric Environment for Analyzing Measurement Data in Multi-Model Settings (Extended Abstract)
A Symbolic Numeric Environment for Analyzing Measurement Data in Multi-Model Settings (Extended Abstract) Christoph Richard 1 and Andreas Weber 2? 1 Institut für Theoretische Physik, Universität Tübingen,
More informationGEOTECHNICAL SITE INVESTIGATION OF SUBSURFACES BY SURFACE WAVES CONSIDERING HIGHER MODES
13 th World Conference on Earthquake Engineering Vancouver, B.C., Canada August 1-6, 2004 Paper No. 2371 GEOECHNICAL SIE INVESIGAION OF SUBSURFACES BY SURFACE WAVES CONSIDERING HIGHER MODES Frank WUKE
More informationBessel function - Wikipedia, the free encyclopedia
Bessel function - Wikipedia, the free encyclopedia Bessel function Page 1 of 9 From Wikipedia, the free encyclopedia In mathematics, Bessel functions, first defined by the mathematician Daniel Bernoulli
More informationNumerical Analysis of Electromagnetic Fields
Pei-bai Zhou Numerical Analysis of Electromagnetic Fields With 157 Figures Springer-Verlag Berlin Heidelberg New York London Paris Tokyo Hong Kong Barcelona Budapest Contents Part 1 Universal Concepts
More informationA NOTE ON DIMENSION REDUCTION AND FINITE ENERGY LOCALIZED WAVE SOLUTIONS TO THE 3-D KLEIN-GORDON AND SCALAR WAVE EQUATIONS. PART II.
Progress In Electromagnetics Research, PIER 27, 357 365, 2 A NOTE ON DIMENSION REDUCTION AND FINITE ENERGY LOCALIZED WAVE SOLUTIONS TO THE 3-D KLEIN-GORDON AND SCALAR WAVE EQUATIONS. PART II. XWAVE-TYPE
More informationAn Evaluation of the Force Reduction Factor in the Force-Based Seismic Design
An Evaluation of the Force Reduction Factor in the Force-Based Seismic Design Gakuho Watanabe and Kazuhiko Kawashima Tokyo Institute of Technology, O-Okayama, Meguro, Tokyo, Japan, 5-55 ABSTRACT This paper
More informationOn the paper Role of potentials in the Aharonov-Bohm effect. A. M. Stewart
On the paper Role of potentials in the Aharonov-Bohm effect A. M. Stewart Emeritus Faculty, The Australian National University, Canberra, ACT 0200, Australia. http://www.anu.edu.au/emeritus/members/pages/a_stewart/
More informationELECTROMAGNETIC FIELD OF A HORIZONTAL IN- FINITELY LONG MAGNETIC LINE SOURCE OVER THE EARTH COATED WITH A DIELECTRIC LAYER
Progress In Electromagnetics Research Letters, Vol. 31, 55 64, 2012 ELECTROMAGNETIC FIELD OF A HORIZONTAL IN- FINITELY LONG MAGNETIC LINE SOURCE OVER THE EARTH COATED WITH A DIELECTRIC LAYER Y.-J. Zhi
More informationAn Alternative Exact Solution Method for the Reduced Wave Equation with a Variable Coefficient
Applied Mathematical Sciences, Vol. 1, 2007, no. 43, 2101-2109 An Alternative Exact Solution Method for the Reduced Wave Equation with a Variable Coefficient Adil Misir Department of Mathematics Faculty
More informationξ Da nc n = c Pe ξ = c o ( θ ) (95) c ξ = 0 (96) θ = 0 ; c= c i ξ ( ) (97) (98) Pe z = u L E z
5. THE AXIAL DISPERSION MODEL The axial dispersion model can readily be extended to apply to turbulent flow conditions. Under such conditions the velocity profile is almost flat and averaging with respect
More informationSpherical Wave Approach (SWA) for the electromagnetic modelling of 3D GPR scenarios. Nicola Tedeschi
Spherical Wave Approach (SWA) for the electromagnetic modelling of 3D GPR scenarios Nicola Tedeschi Outline Solutions of the vector wave equation Spherical wave approach Spherical reflected and transmitted
More informationOn spherical-wave scattering by a spherical scatterer and related near-field inverse problems
IMA Journal of Applied Mathematics (2001) 66, 539 549 On spherical-wave scattering by a spherical scatterer and related near-field inverse problems C. ATHANASIADIS Department of Mathematics, University
More informationQuasi-static Vertical Magnetic Field of a Large Horizontal Circular Loop Located at the Earth s Surface
Progress In Electromagnetics Research Letters, Vol. 6, 9 34, 16 Quasi-static Vertical Magnetic Field of a Large Horizontal Circular Loop Located at the Earth s Surface Mauro Parise * Abstract In this work,
More informationScattering of ECRF waves by edge density fluctuations and blobs
PSFC/JA-14-7 Scattering of ECRF waves by edge density fluctuations and blobs A. K. Ram and K. Hizanidis a June 2014 Plasma Science and Fusion Center, Massachusetts Institute of Technology Cambridge, MA
More informationRadiation by a dielectric wedge
Radiation by a dielectric wedge A D Rawlins Department of Mathematical Sciences, Brunel University, United Kingdom Joe Keller,Cambridge,2-3 March, 2017. We shall consider the problem of determining the
More informationDEVELOPMENT OF TOOL FOR THE IDENTIFICATION OF STIFFNESS AND DAMPING COEFFICIENTS OF JOURNAL BEARING
DEVELOPMENT OF TOOL FOR THE IDENTIFICATION OF STIFFNESS AND DAMPING COEFFICIENTS OF JOURNAL BEARING Sampath Kumar B 1 and G.Satish Babu 2 1 P.G student, Department of mechanical engineering, JNTUH College
More informationKirchhoff prestack depth migration in simple models of various anisotropy
Kirchhoff prestack depth migration in simple models of various anisotropy Václav Bucha Department of Geophysics, Faculty of Mathematics and Physics, Charles University, Ke Karlovu, 6 Praha, Czech Republic,
More informationElectromagnetic Scattering from an Anisotropic Uniaxial-coated Conducting Sphere
Progress In Electromagnetics Research Symposium 25, Hangzhou, China, August 22-26 43 Electromagnetic Scattering from an Anisotropic Uniaxial-coated Conducting Sphere You-Lin Geng 1,2, Xin-Bao Wu 3, and
More informationSlow Photons in Vacuum as Elementary Particles. Chander Mohan Singal
Ref ETOP98 Slow Photons in Vacuum as Elementary Particles Chander Mohan Singal Department of Physics, Indian Institute of Technology-Delhi, Hauz Khas, New Delhi-1116, INDIA E-Mail: drcmsingal@yahoocom
More information2009 Winton 1 Distributi ( ons 2) (2)
Distributions ib i (2) 2 IV. Triangular Distribution ib ti Known values The minimum (a) The mode (b - the most likely value of the pdf) The maimum (c) f() probability density function (area under the curve
More informationGreen s functions for planarly layered media
Green s functions for planarly layered media Massachusetts Institute of Technology 6.635 lecture notes Introduction: Green s functions The Green s functions is the solution of the wave equation for a point
More informationEstimating the Pure Characteristics Demand Model: A Computational Note
Estimating the Pure Characteristics Demand Model: A Computational Note Minjae Song School of Economics, Georgia Institute of Technology April, 2006 Abstract This paper provides details of the computational
More informationImpedance of a single-turn coil due to a double-layered sphere with varying properties
Impedance of a single-turn coil due to a double-layered sphere with varying properties A. A. Kolyshkin Rémi Vaillancourt CRM-293 June 994 Department of Applied Mathematics, Riga Technical University, Riga,
More informationPEMC PARABOLOIDAL REFLECTOR IN CHIRAL MEDIUM SUPPORTING POSITIVE PHASE VELOC- ITY AND NEGATIVE PHASE VELOCITY SIMULTANE- OUSLY
Progress In Electromagnetics Research Letters, Vol. 10, 77 86, 2009 PEMC PARABOLOIDAL REFLECTOR IN CHIRAL MEDIUM SUPPORTING POSITIVE PHASE VELOC- ITY AND NEGATIVE PHASE VELOCITY SIMULTANE- OUSLY T. Rahim
More informationAnalysis of Type-II Progressively Hybrid Censored Data
Analysis of Type-II Progressively Hybrid Censored Data Debasis Kundu & Avijit Joarder Abstract The mixture of Type-I and Type-II censoring schemes, called the hybrid censoring scheme is quite common in
More informationRELATION BETWEEN SMALL FUNCTIONS WITH DIFFERENTIAL POLYNOMIALS GENERATED BY MEROMORPHIC SOLUTIONS OF HIGHER ORDER LINEAR DIFFERENTIAL EQUATIONS
Kragujevac Journal of Mathematics Volume 38(1) (2014), Pages 147 161. RELATION BETWEEN SMALL FUNCTIONS WITH DIFFERENTIAL POLYNOMIALS GENERATED BY MEROMORPHIC SOLUTIONS OF HIGHER ORDER LINEAR DIFFERENTIAL
More informationPRELIMINARIES FOR HYPERGEOMETRIC EQUATION. We will only consider differential equations with regular singularities in this lectures.
PRELIMINARIES FOR HYPERGEOMETRIC EQUATION EDMUND Y.-M. CHIANG Abstract. We give a brief introduction to some preliminaries for Gauss hypergeometric equations. We will only consider differential equations
More informationSolutions of the Schrödinger equation for an attractive 1/r 6 potential
PHYSICAL REVIEW A VOLUME 58, NUMBER 3 SEPTEMBER 998 Solutions of the Schrödinger equation for an attractive /r 6 potential Bo Gao Department of Physics and Astronomy, University of Toledo, Toledo, Ohio
More informationPH36010: Numerical Methods - Evaluating the Lorenz Attractor using Runge-Kutta methods Abstract
PH36010: Numerical Methods - Evaluating the Lorenz Attractor using Runge-Kutta methods Mr. Benjamen P. Reed (110108461) IMPACS, Aberystwyth University January 31, 2014 Abstract A set of three coupled ordinary
More informationNumerical Analysis: Solving Systems of Linear Equations
Numerical Analysis: Solving Systems of Linear Equations Mirko Navara http://cmpfelkcvutcz/ navara/ Center for Machine Perception, Department of Cybernetics, FEE, CTU Karlovo náměstí, building G, office
More informationSOUND PROPAGATION CHARACTERISTICS IN ARCTIC OCEAN CALCULATED BY ELASTIC PE METHOD USING ROTATED PADE APPROXIMATION
SOUND PROPAGATION CHARACTERISTICS IN ARCTIC OCEAN CALCULATED BY ELASTIC PE METHOD USING ROTATED PADE APPROXIMATION PACS REFERENCE: 43.30.Bp Tsuchiya Takenobu; Endoh Nobuyuki; Anada Tetsuo Faculty of Engineering,
More information4 Power Series Solutions: Frobenius Method
4 Power Series Solutions: Frobenius Method Now the ODE adventure takes us to series solutions for ODEs, a technique A & W, that is often viable, valuable and informative. These can be readily applied Sec.
More informationDr. Parveen Lata Department of Basic and Applied Sciences, Punjabi University, Patiala, Punjab, India.
International Journal of Theoretical and Applied Mechanics. ISSN 973-685 Volume 12, Number 3 (217) pp. 435-443 Research India Publications http://www.ripublication.com Linearly Distributed Time Harmonic
More informationUSAGE OF NUMERICAL METHODS FOR ELECTROMAGNETIC SHIELDS OPTIMIZATION
October 4-6, 2007 - Chiinu, Rep.Moldova USAGE OF NUMERICAL METHODS FOR ELECTROMAGNETIC SHIELDS OPTIMIZATION Ionu- P. NICA, Valeriu Gh. DAVID, /tefan URSACHE Gh. Asachi Technical University Iai, Faculty
More informationEigenfunction Expansions of Source-excited Electromagnetic Fields on Open Cylindrical Guiding Structures in Unbounded Gyrotropic Media
128 Eigenfunction Expansions of Source-excited Electromagnetic Fields on Open Cylindrical Guiding Structures in Unbounded Gyrotropic Media A. V. Kudrin 1, E. Yu. Petrov 1, and T. M. Zaboronkova 2 1 University
More informationStrong Convergence of Two Iterative Algorithms for a Countable Family of Nonexpansive Mappings in Hilbert Spaces
International Mathematical Forum, 5, 2010, no. 44, 2165-2172 Strong Convergence of Two Iterative Algorithms for a Countable Family of Nonexpansive Mappings in Hilbert Spaces Jintana Joomwong Division of
More informationElectromagnetic scattering. Graduate Course Electrical Engineering (Communications) 1 st Semester, Sharif University of Technology
Electroagnetic scattering Graduate Course Electrical Engineering (Counications) 1 st Seester, 1388-1389 Sharif University of Technology Contents of lecture 5 Contents of lecture 5: Scattering fro a conductive
More informationModified Bessel functions : Iα, Kα
Modified Bessel functions : Iα, Kα The Bessel functions are valid even for complex arguments x, and an important special case is that of a purely imaginary argument. In this case, the solutions to the
More informationCOMBINED EFFECTS OF RADIATION AND JOULE HEATING WITH VISCOUS DISSIPATION ON MAGNETOHYDRODYNAMIC FREE CONVECTION FLOW AROUND A SPHERE
Suranaree J. Sci. Technol. Vol. 20 No. 4; October - December 2013 257 COMBINED EFFECTS OF RADIATION AND JOULE HEATING WITH VISCOUS DISSIPATION ON MAGNETOHYDRODYNAMIC FREE CONVECTION FLOW AROUND A SPHERE
More informationElectromagnetic Waves
Electromagnetic Waves Maxwell s equations predict the propagation of electromagnetic energy away from time-varying sources (current and charge) in the form of waves. Consider a linear, homogeneous, isotropic
More informationA novel difference schemes for analyzing the fractional Navier- Stokes equations
DOI: 0.55/auom-207-005 An. Şt. Univ. Ovidius Constanţa Vol. 25(),207, 95 206 A novel difference schemes for analyzing the fractional Navier- Stokes equations Khosro Sayevand, Dumitru Baleanu, Fatemeh Sahsavand
More informationSensors & Transducers 2015 by IFSA Publishing, S. L.
Sensors & Transducers 2015 by IFSA Publishing, S. L. http://www.sensorsportal.com Improved Performance in Soil Resistivity Measurement through Nonlinear Data Fitting 2* T. Islam and 1 Z. Chik 1 Department
More information8 STOCHASTIC SIMULATION
8 STOCHASTIC SIMULATIO 59 8 STOCHASTIC SIMULATIO Whereas in optimization we seek a set of parameters x to minimize a cost, or to maximize a reward function J( x), here we pose a related but different question.
More information1 The formation and analysis of optical waveguides
1 The formation and analysis of optical waveguides 1.1 Introduction to optical waveguides Optical waveguides are made from material structures that have a core region which has a higher index of refraction
More informationTime-varying electromagnetic field
Chapter 5 Time-varying electromagnetic field In this chapter we take up the time-varying electromagnetic field. By this study, we will be able to understand in a systematic way that actually the electric
More informationAvailable online at ScienceDirect. Procedia Engineering 144 (2016 )
Available online at www.sciencedirect.com ScienceDirect Procedia Engineering 144 (016 ) 170 177 1th International Conference on Vibration Problems, ICOVP 015 Propagation of Torsional Surface Wave in Anisotropic
More informationThe use of quad quad resistivity in helicopter electromagnetic mapping
GEOPHYSICS, VOL. 67, NO. 2 (MARCH-APRIL 2002); P. 459 467, 8 FIGS., 1 TABLE. 10.1190/1.1468605 The use of quad quad resistivity in helicopter electromagnetic mapping Haoping Huang and Douglas C. Fraser
More informationPropagation of Radio Frequency Waves Through Fluctuations in Plasmas
PSFC/JA-15- Propagation of Radio Frequency Waves Through Fluctuations in Plasmas A. K. Ram K. Hizanidis a and S. Valvis a a National Technical University of Athens (part of HELLAS) School of Electrical
More informationISSN: X (p); (e)
TORSIONA SURFACE WAVE IN SEF-REINFORCED AYER SANDWICHED BETWEEN TWO VISCO-EASIC HAF-SPACES UNDER THE INITIA STRESS Nidhi Dewangan, Sanjeev A. Sahu and Soniya Chaudhary S.G.G. Govt P.G. College, Kurud,
More informationarxiv:astro-ph/ v1 14 Jan 2002
Accretion onto a Black Hole in the Presence of Bremsstrahlung Cooling - Parameter Space Study Sivakumar G. Manickam and Sandip K. Chakrabarti arxiv:astro-ph/020203v 4 Jan 2002 S. N. Bose National Centre
More informationLecture 16: Relaxation methods
Lecture 16: Relaxation methods Clever technique which begins with a first guess of the trajectory across the entire interval Break the interval into M small steps: x 1 =0, x 2,..x M =L Form a grid of points,
More informationCHAPTER 4 ELECTROMAGNETIC WAVES IN CYLINDRICAL SYSTEMS
CHAPTER 4 ELECTROMAGNETIC WAVES IN CYLINDRICAL SYSTEMS The vector Helmholtz equations satisfied by the phasor) electric and magnetic fields are where. In low-loss media and for a high frequency, i.e.,
More informationLevenberg-Marquardt Method for Solving the Inverse Heat Transfer Problems
Journal of mathematics and computer science 13 (2014), 300-310 Levenberg-Marquardt Method for Solving the Inverse Heat Transfer Problems Nasibeh Asa Golsorkhi, Hojat Ahsani Tehrani Department of mathematics,
More information1. Comparison of stability analysis to previous work
. Comparison of stability analysis to previous work The stability problem (6.4) can be understood in the context of previous work. Benjamin (957) and Yih (963) have studied the stability of fluid flowing
More information9 Fluid Instabilities
9. Stability of a shear flow In many situations, gaseous flows can be subject to fluid instabilities in which small perturbations can rapidly flow, thereby tapping a source of free energy. An important
More informationON THE COEFFICIENTS OF AN ASYMPTOTIC EXPANSION RELATED TO SOMOS QUADRATIC RECURRENCE CONSTANT
Applicable Analysis and Discrete Mathematics available online at http://pefmath.etf.bg.ac.yu Appl. Anal. Discrete Math. x (xxxx), xxx xxx. doi:10.2298/aadmxxxxxxxx ON THE COEFFICIENTS OF AN ASYMPTOTIC
More informationLine adsorption in a mean-field density functional model
Chemistry Physical & Theoretical Chemistry fields Okayama University Year 2006 Line adsorption in a mean-field density functional model Kenichiro Koga Benjamin Widom Okayama University, koga@cc.okayama-u.ac.jp
More informationTRADITIONALLY, the singularity present in the method of
1200 IEEE TRANSACTIONS ON ANTENNAS AND PROPAGATION, VOL. 54, NO. 4, APRIL 2006 Evaluation Integration of the Thin Wire Kernel Donald R. Wilton, Fellow, IEEE, Nathan J. Champagne, Senior Member, IEEE Abstract
More informationSurface Waves and Free Oscillations. Surface Waves and Free Oscillations
Surface waves in in an an elastic half spaces: Rayleigh waves -Potentials - Free surface boundary conditions - Solutions propagating along the surface, decaying with depth - Lamb s problem Surface waves
More informationLimits at Infinity. Horizontal Asymptotes. Definition (Limits at Infinity) Horizontal Asymptotes
Limits at Infinity If a function f has a domain that is unbounded, that is, one of the endpoints of its domain is ±, we can determine the long term behavior of the function using a it at infinity. Definition
More information