Two Kinds of Parabolic Equation algorithms in the Computational Electromagnetics
|
|
- Joleen White
- 5 years ago
- Views:
Transcription
1 Avaiabe onine at Procedia Engineering 9 (0) Internationa Workshop on Information and Eectronics Engineering (IWIEE) Two Kinds of Paraboic Equation agorithms in the Computationa Eectromagnetics Kong Meng a,chen Mingsheng a,hu Yujuan b,wu Xianiang a a Department of physics and Eectronic Engineering,,Hefei Norma University,Hefei,3060,China. b Department of Pubic Computer teaching,hefei Norma University,Hefei,3060, China. Abstract Sip-step Fourier Transform (SSFT) and Finite Difference (FD) agorithms are introduced to sove Paraboic Equation (PE). SSFT is appied to anayze the probem of radio propagation in vacuum. Based on Spit-step Padé approximation, a high order FD agorithm is introduced to compute the Radar Cross Section (RCS). According to the resuts, the differences of the two kinds of PE method have been discussed. 0 Pubished by Esevier Ltd. Seection and/or peer-review under responsibiity of Harbin University of Science and Technoogy. Open access under CC BY-NC-ND icense. Keywords: Paraboic Equation, Sip-step Fourier Transform, Finite Difference, Radio Propagation, Radar Cross Section. Introduction As an approximated form of the wave equation, the Paraboic Equation (PE) method modes energy propagating in a cone centered on the paraxia direction. Recenty, many researchers have appied the paraboic equation to sove the probems of radio propagation and eectromagnetic scattering []. It was reaized that PE method may aso be adopted as an efficient numerica too of eectromagnetic fied cacuation, and it may bridging the gap between rigorous numerica methods ike the method of moments (MOM) or finite difference time domain (FDTD) [], and asymptotica methods based on ray-tracing or physica optics (PO). Using PE method can avoid the imitation of Corresponding author. Chen Mingsheng Te.: ; E-mai address: cmsh@ahu.edu.cn Pubished by Esevier Ltd. doi:0.06/j.proeng Open access under CC BY-NC-ND icense.
2 46 Kong Meng et a. / Procedia Engineering 9 (0) CPU time and memory required by these rigorous methods. It aso can reduce the error caused by the high-frequency approximation. Recenty, the usua toos for soving the PE are spit-step Fourier transform (SSFT) and finite difference (FD) agorithms. SSFT using the FFT technique is a frequency-domain agorithm, when it is used to sove the PE, the step size of SSFT is amost free from the restrictions and can be seected as a arge discritized eement. Therefore, SSFT agorithm of PE is suitabe for arge-scae wave propagation probems. SSFT agorithm is compicated when deaing with irreguar boundary, so it is not easy to be used in the cacuation of eectromagnetic scattering probems with compex structures. FD agorithm is carried out through the mesh of the computationa domain, and accompishes the computation according to the eectromagnetic fied on the adjacent grid points. Since precise mesh is used, FD agorithm can directy cacuate the fied on arbitrary boundaries, therefore it is more convenient to dea with irreguar boundaries, but the discritized step in FD agorithm is restricted by the waveength. Therefore one must take very sma step, which wi resut arge-scae matrix operations. Soving the high-frequency, arge-scae wave propagation probems based on PE, FD agorithm is sow and the computer memory consuming is aso high. FD agorithm is mainy used to cacuate the target RCS and forecast the eectromagnetic distribution of urban and other sma-scae wave propagation. [3]. SSFT agorithm of Paraboic Equation As we know, froward paraboic equation can be soved by marching techniques and has the forma soution as [4] ikδx( Q ) u( x+δ xz, ) = e u( xz, ) () We assume the wave function of pane wave is of the form as beow (, ) exp( ( cosα ) sin α) u x z = ikx + ikz () Equation () can be approximated by using first-order Tayor expansions of the square root Q, which yieds the Standard Paraboic Equation (SPE) when set n=, one wi get u u + ik + k ( n ) u = 0 z x (3) u u + ik = 0 z x we referred equation (4) as SPE in vacuum. Using the properties of the Fourier transform, we can write the Fourier transform of Equation (4) as U( x, f) 4 π f U ( x, f ) + ik = 0 x This is now an ordinary differentia equation which can be soved in cosed form, giving f x U( x, f) e k U(0, f) We can write the inverse Fourier transform of Equation (4) as = (6) (4) (5)
3 Kong Meng et a. / Procedia Engineering 9 (0) where u(0,z) is initia fied. The forma soution of SPE is given as f x uxz (, ) e k = I I( u(0, z)) ik x Q Δ u x+δ xz = e u xz (8) (, ) (, ) Then one wi get the compete spit-step Fourier transform soution of equation (8) Where Δx is step. 3.FD agorithm of Paraboic Equation 3. Spit-step Padé PE f Δx ux ( xz, ) e k +Δ =I I( uxz (, )) In this paper, we approximate the wide-ange/exponentia pseudo-differentia operator directy with a sum of Padé(,) functions as ( Q ) N ( ) ~ b ( Q ) ikδx Q ikδx + + a Q e = e + (0) = + in which ccompex coefficients (a, b) are determined by the constraints of Tayor series expansion and stabiity conditions, whie its accuracy is aso dependent on step Δ x size. The approximation accuracy equivaent Padé (N, N), which has N-order accuracy [5]. Substituting (0) into equation (), we wi get the form of SSPPE: where N u( x+δ xz, ) = u( xz, ) + v( x xz, ) +Δ () = ( +Δ, ) = + ( ) ( ) (, ) v x xz b Q a Q u xz Discretizing in z, we obtain finite-difference equations for the where j j v = v( x+δx, jδ z), s = kδ z, u = u( x, jδ z) v s in the form: { } j+ j s j ( ( ) ) a j+ j j v + v + + s n v = u + u + ( s ( n ) ) u b b Combined Equation () and (3), fied can be obtained recursivey. To simuate the infinite space of the eectromagnetic scattering probem, we added Non-oca boundary conditions [6] (NLBC) as the truncation boundary conditions. (7) (9) () (3)
4 48 Kong Meng et a. / Procedia Engineering 9 (0) Near-Fied/Far-Fied Transformation Recaed that the -D bistatic RCS of an object in direction θ is given as ( cos, sin ) r s ψ r θ θ σ( θ) = im πr i r ψ ( rcos θ, rsin θ) (4) If the incident fied is a pane wave with unit ampitude, after compicated mathematics derives, equation (4) can be expressed as s ikzsinθ σ ( θ ) = k cos θ u ( x, z) e dz (5) In which x 0 can usuay be chosen as 0λ (λ is the waveength of incident wave). 4.Numerica Resuts 4. Large-scae radio wave propagation + We assume the incident fied is of the form in equation (), and the frequency of wave is 0.3GHz in a vacuum. The rea and imaginary parts of the fied u(x, z) were cacuated using the SSFT method,and set α = 0 o,x=000m, Δ x =00 m. The computation is accompished on a Pentium IV 3.0GHz computer using Matab7.0, the CPU time for this simuation is 0. seconds, and the resuts is shown in Figure. It can be seen that the resut of SSFT method agrees we with that of anaytica resuts. When the FD method is used to hande the probem, it is difficut to compete the operation because of the imited step and arge-scae matrix operations. Due to arge step desirabe, SSFT method is more appropriate for deaing with the probem. 0 Figure. (a) Figure. (b) Figure. Radio waves propagation in cacuating the vacuum vaue with SSFT :(a)the rea part of fied,(b) the imaginary part of fied 4. Eectromagnetic scattering of eectricay arge objects In this numerica exampe, the incident fied is a pane wave with its ampitude set to be.0, the incident source with its waveength of 0. m is horizonta poarized and propagating in a vacuum, which iuminated a PEC circuar cyinder of radius 5λ. The computation domain is set to be X Z=[0, 0λ] [0, 40λ], we use PML absorbing boundary conditions at the top and bottom of the domain. As shown in Fig. the ampitude of the scattered near fied is computed with the SPE and SSPPE (N=4) methods. In Fig.3 bistatic RCS resuts obtained by SSPPE method is compared with the resuts of method of moments (MOM), they agree we from 0 to 50. And the CPU consumed by this exampe is 5.8 seconds. When we use the SPE method [4] to hande the same probem, the CPU time is 3.0s, but
5 Kong Meng et a. / Procedia Engineering 9 (0) the anges accuracy up to 5.To obtain the fu bistatic RCS, we need ony two rotated PE runs for SSPPE whie in SPE method the number of rotation is six [4]. The resuts which show in Figure 4 proved the efficiency of the present method. Figure.(a) Figure.(b) Figure.3 Figure.4 Figure. The ampitude of the scattered near fied computed with the SPE (a) and SSPPE(b). Figure.3 Bistatic RCS of PEC radius 5λ Cyinder. Figure.4 Bistatic RCS of PEC radius 5λ cyinder obtained with rotation paraboic axis. 5.Concusions Two agorithms for soving Paraboic Equation, Sip-step Fourier Transform and Finite Difference, are presented and compared. Meanwhie the appication of the two agorithms to radio propagation and eectromagnetic scattering probems is discussed. According to numerica resuts, we suggest to hande arge-scae radio wave propagation probems with SSFT method to obtain good resuts, whie FD method is suitabe for deaing with eectromagnetic scattering probems. Acknowedgement This work is supported by the Key Program of Nationa Natura Science Foundation of China (609300) and Science and Technoogica Fund of Anhui Province for Outstanding Youth ( Y08) and the Key Scientific Research Base program of Hefei Norma University(0jd06) and partiay by the Natura Science Foundation of the Anhui Higher Education Institution of China (KJ0A40). References [] Levy M F,Borsboom P P.Radar cross-section computations using the paraboic equation method.eectron Lett;996,3(3),p [] A. Tafove and S.C. Hagness. Computationa eectrodynamics: the finite-difference time-domain method. Norwood MA, Artech House, 000. [3] Zaporozhets A A, Levy M F.Bistatic RCS cacuations with the vector paraboic equation method.ieee Trans on AP,999,47(),p [4] M. F. Levy. Paraboic equation methods for eectromagnetic wave propagation. London: The Institution of Eectrica Engineers, 000. [5] Coins M D, Generaization of the spit-step padé soution.acoust.soc.amer,994,96,p [6] Dan Givoi. Non-refecting Boundary Conditions. Journa of Computationa Physics, 99; 94, p.-9.
A Simple and Efficient Algorithm of 3-D Single-Source Localization with Uniform Cross Array Bing Xue 1 2 a) * Guangyou Fang 1 2 b and Yicai Ji 1 2 c)
A Simpe Efficient Agorithm of 3-D Singe-Source Locaization with Uniform Cross Array Bing Xue a * Guangyou Fang b Yicai Ji c Key Laboratory of Eectromagnetic Radiation Sensing Technoogy, Institute of Eectronics,
More informationHYDROGEN ATOM SELECTION RULES TRANSITION RATES
DOING PHYSICS WITH MATLAB QUANTUM PHYSICS Ian Cooper Schoo of Physics, University of Sydney ian.cooper@sydney.edu.au HYDROGEN ATOM SELECTION RULES TRANSITION RATES DOWNLOAD DIRECTORY FOR MATLAB SCRIPTS
More informationHigh-order approximations to the Mie series for electromagnetic scattering in three dimensions
Proceedings of the 9th WSEAS Internationa Conference on Appied Mathematics Istanbu Turkey May 27-29 2006 (pp199-204) High-order approximations to the Mie series for eectromagnetic scattering in three dimensions
More informationGauss Law. 2. Gauss s Law: connects charge and field 3. Applications of Gauss s Law
Gauss Law 1. Review on 1) Couomb s Law (charge and force) 2) Eectric Fied (fied and force) 2. Gauss s Law: connects charge and fied 3. Appications of Gauss s Law Couomb s Law and Eectric Fied Couomb s
More informationResearch on liquid sloshing performance in vane type tank under microgravity
IOP Conference Series: Materias Science and Engineering PAPER OPEN ACCESS Research on iquid soshing performance in vane type tan under microgravity Reated content - Numerica simuation of fuid fow in the
More informationInternational Journal of Mass Spectrometry
Internationa Journa of Mass Spectrometry 280 (2009) 179 183 Contents ists avaiabe at ScienceDirect Internationa Journa of Mass Spectrometry journa homepage: www.esevier.com/ocate/ijms Stark mixing by ion-rydberg
More informationDISTRIBUTION OF TEMPERATURE IN A SPATIALLY ONE- DIMENSIONAL OBJECT AS A RESULT OF THE ACTIVE POINT SOURCE
DISTRIBUTION OF TEMPERATURE IN A SPATIALLY ONE- DIMENSIONAL OBJECT AS A RESULT OF THE ACTIVE POINT SOURCE Yury Iyushin and Anton Mokeev Saint-Petersburg Mining University, Vasiievsky Isand, 1 st ine, Saint-Petersburg,
More informationMultiple Beam Interference
MutipeBeamInterference.nb James C. Wyant 1 Mutipe Beam Interference 1. Airy's Formua We wi first derive Airy's formua for the case of no absorption. ü 1.1 Basic refectance and transmittance Refected ight
More informationLecture 9. Stability of Elastic Structures. Lecture 10. Advanced Topic in Column Buckling
Lecture 9 Stabiity of Eastic Structures Lecture 1 Advanced Topic in Coumn Bucking robem 9-1: A camped-free coumn is oaded at its tip by a oad. The issue here is to find the itica bucking oad. a) Suggest
More informationA Fictitious Time Integration Method for a One-Dimensional Hyperbolic Boundary Value Problem
Journa o mathematics and computer science 14 (15) 87-96 A Fictitious ime Integration Method or a One-Dimensiona Hyperboic Boundary Vaue Probem Mir Saad Hashemi 1,*, Maryam Sariri 1 1 Department o Mathematics,
More informationModule 22: Simple Harmonic Oscillation and Torque
Modue : Simpe Harmonic Osciation and Torque.1 Introduction We have aready used Newton s Second Law or Conservation of Energy to anayze systems ike the boc-spring system that osciate. We sha now use torque
More informationApproximation and Fast Calculation of Non-local Boundary Conditions for the Time-dependent Schrödinger Equation
Approximation and Fast Cacuation of Non-oca Boundary Conditions for the Time-dependent Schrödinger Equation Anton Arnod, Matthias Ehrhardt 2, and Ivan Sofronov 3 Universität Münster, Institut für Numerische
More information4 Separation of Variables
4 Separation of Variabes In this chapter we describe a cassica technique for constructing forma soutions to inear boundary vaue probems. The soution of three cassica (paraboic, hyperboic and eiptic) PDE
More informationA Novel Learning Method for Elman Neural Network Using Local Search
Neura Information Processing Letters and Reviews Vo. 11, No. 8, August 2007 LETTER A Nove Learning Method for Eman Neura Networ Using Loca Search Facuty of Engineering, Toyama University, Gofuu 3190 Toyama
More informationPublished in: Proceedings of the Twenty Second Nordic Seminar on Computational Mechanics
Aaborg Universitet An Efficient Formuation of the Easto-pastic Constitutive Matrix on Yied Surface Corners Causen, Johan Christian; Andersen, Lars Vabbersgaard; Damkide, Lars Pubished in: Proceedings of
More informationPhysics 506 Winter 2006 Homework Assignment #6 Solutions
Physics 506 Winter 006 Homework Assignment #6 Soutions Textbook probems: Ch. 10: 10., 10.3, 10.7, 10.10 10. Eectromagnetic radiation with eiptic poarization, described (in the notation of Section 7. by
More informationCONCHOID OF NICOMEDES AND LIMACON OF PASCAL AS ELECTRODE OF STATIC FIELD AND AS WAVEGUIDE OF HIGH FREQUENCY WAVE
Progress In Eectromagnetics Research, PIER 30, 73 84, 001 CONCHOID OF NICOMEDES AND LIMACON OF PASCAL AS ELECTRODE OF STATIC FIELD AND AS WAVEGUIDE OF HIGH FREQUENCY WAVE W. Lin and Z. Yu University of
More informationSupporting Information for Suppressing Klein tunneling in graphene using a one-dimensional array of localized scatterers
Supporting Information for Suppressing Kein tunneing in graphene using a one-dimensiona array of ocaized scatterers Jamie D Was, and Danie Hadad Department of Chemistry, University of Miami, Cora Gabes,
More informationIdentification of macro and micro parameters in solidification model
BULLETIN OF THE POLISH ACADEMY OF SCIENCES TECHNICAL SCIENCES Vo. 55, No. 1, 27 Identification of macro and micro parameters in soidification mode B. MOCHNACKI 1 and E. MAJCHRZAK 2,1 1 Czestochowa University
More informationFirst-Order Corrections to Gutzwiller s Trace Formula for Systems with Discrete Symmetries
c 26 Noninear Phenomena in Compex Systems First-Order Corrections to Gutzwier s Trace Formua for Systems with Discrete Symmetries Hoger Cartarius, Jörg Main, and Günter Wunner Institut für Theoretische
More informationTunnel Geological Prediction Radar Alternating Electromagnetic Field Propagation Attenuation in Lossy Inhomogeneous Medium
Tunne Geoogica Prediction Radar Aternating Eectromagnetic Fied Propagation Attenuation in Lossy Inhomogeneous Medium Si Yang Chen,a, Yan Peng Zhu,b, Zhong Li,c, Tian Yu Zhang Coage of civi engineering,lanhou
More informationPaper presented at the Workshop on Space Charge Physics in High Intensity Hadron Rings, sponsored by Brookhaven National Laboratory, May 4-7,1998
Paper presented at the Workshop on Space Charge Physics in High ntensity Hadron Rings, sponsored by Brookhaven Nationa Laboratory, May 4-7,998 Noninear Sef Consistent High Resoution Beam Hao Agorithm in
More informationSolution of Wave Equation by the Method of Separation of Variables Using the Foss Tools Maxima
Internationa Journa of Pure and Appied Mathematics Voume 117 No. 14 2017, 167-174 ISSN: 1311-8080 (printed version); ISSN: 1314-3395 (on-ine version) ur: http://www.ijpam.eu Specia Issue ijpam.eu Soution
More informationA Solution to the 4-bit Parity Problem with a Single Quaternary Neuron
Neura Information Processing - Letters and Reviews Vo. 5, No. 2, November 2004 LETTER A Soution to the 4-bit Parity Probem with a Singe Quaternary Neuron Tohru Nitta Nationa Institute of Advanced Industria
More informationAvailable online at ScienceDirect. Procedia Computer Science 96 (2016 )
Avaiabe onine at www.sciencedirect.com ScienceDirect Procedia Computer Science 96 (206 92 99 20th Internationa Conference on Knowedge Based and Inteigent Information and Engineering Systems Connected categorica
More information$, (2.1) n="# #. (2.2)
Chapter. Eectrostatic II Notes: Most of the materia presented in this chapter is taken from Jackson, Chap.,, and 4, and Di Bartoo, Chap... Mathematica Considerations.. The Fourier series and the Fourier
More informationThe influence of temperature of photovoltaic modules on performance of solar power plant
IOSR Journa of Engineering (IOSRJEN) ISSN (e): 2250-3021, ISSN (p): 2278-8719 Vo. 05, Issue 04 (Apri. 2015), V1 PP 09-15 www.iosrjen.org The infuence of temperature of photovotaic modues on performance
More informationMethods for Ordinary Differential Equations. Jacob White
Introduction to Simuation - Lecture 12 for Ordinary Differentia Equations Jacob White Thanks to Deepak Ramaswamy, Jaime Peraire, Micha Rewienski, and Karen Veroy Outine Initia Vaue probem exampes Signa
More informationIntroduction to Simulation - Lecture 14. Multistep Methods II. Jacob White. Thanks to Deepak Ramaswamy, Michal Rewienski, and Karen Veroy
Introduction to Simuation - Lecture 14 Mutistep Methods II Jacob White Thans to Deepa Ramaswamy, Micha Rewiensi, and Karen Veroy Outine Sma Timestep issues for Mutistep Methods Reminder about LTE minimization
More informationA Brief Introduction to Markov Chains and Hidden Markov Models
A Brief Introduction to Markov Chains and Hidden Markov Modes Aen B MacKenzie Notes for December 1, 3, &8, 2015 Discrete-Time Markov Chains You may reca that when we first introduced random processes,
More informationLecture 8 February 18, 2010
Sources of Eectromagnetic Fieds Lecture 8 February 18, 2010 We now start to discuss radiation in free space. We wi reorder the materia of Chapter 9, bringing sections 6 7 up front. We wi aso cover some
More informationDYNAMIC RESPONSE OF CIRCULAR FOOTINGS ON SATURATED POROELASTIC HALFSPACE
3 th Word Conference on Earthquake Engineering Vancouver, B.C., Canada August -6, 4 Paper No. 38 DYNAMIC RESPONSE OF CIRCULAR FOOTINGS ON SATURATED POROELASTIC HALFSPACE Bo JIN SUMMARY The dynamic responses
More informationLECTURE 10. The world of pendula
LECTURE 0 The word of pendua For the next few ectures we are going to ook at the word of the pane penduum (Figure 0.). In a previous probem set we showed that we coud use the Euer- Lagrange method to derive
More informationPhysics 235 Chapter 8. Chapter 8 Central-Force Motion
Physics 35 Chapter 8 Chapter 8 Centra-Force Motion In this Chapter we wi use the theory we have discussed in Chapter 6 and 7 and appy it to very important probems in physics, in which we study the motion
More informationLecture 17 - The Secrets we have Swept Under the Rug
Lecture 17 - The Secrets we have Swept Under the Rug Today s ectures examines some of the uirky features of eectrostatics that we have negected up unti this point A Puzze... Let s go back to the basics
More information(This is a sample cover image for this issue. The actual cover is not yet available at this time.)
(This is a sampe cover image for this issue The actua cover is not yet avaiabe at this time) This artice appeared in a journa pubished by Esevier The attached copy is furnished to the author for interna
More informationD. Prémel, J.M. Decitre and G. Pichenot. CEA, LIST, F Gif-sur-Yvette, France
SIMULATION OF EDDY CURRENT INSPECTION INCLUDING MAGNETIC FIELD SENSOR SUCH AS A GIANT MAGNETO-RESISTANCE OVER PLANAR STRATIFIED MEDIA COMPONENTS WITH EMBEDDED FLAWS D. Préme, J.M. Decitre and G. Pichenot
More informationSTA 216 Project: Spline Approach to Discrete Survival Analysis
: Spine Approach to Discrete Surviva Anaysis November 4, 005 1 Introduction Athough continuous surviva anaysis differs much from the discrete surviva anaysis, there is certain ink between the two modeing
More informationOn the Number of Limit Cycles for Discontinuous Generalized Liénard Polynomial Differential Systems
Internationa Journa of Bifurcation and Chaos Vo. 25 No. 10 2015 1550131 10 pages c Word Scientific Pubishing Company DOI: 10.112/S02181271550131X On the Number of Limit Cyces for Discontinuous Generaized
More informationCombining reaction kinetics to the multi-phase Gibbs energy calculation
7 th European Symposium on Computer Aided Process Engineering ESCAPE7 V. Pesu and P.S. Agachi (Editors) 2007 Esevier B.V. A rights reserved. Combining reaction inetics to the muti-phase Gibbs energy cacuation
More informationTheory of Generalized k-difference Operator and Its Application in Number Theory
Internationa Journa of Mathematica Anaysis Vo. 9, 2015, no. 19, 955-964 HIKARI Ltd, www.m-hiari.com http://dx.doi.org/10.12988/ijma.2015.5389 Theory of Generaized -Difference Operator and Its Appication
More informationSTABILITY OF A PARAMETRICALLY EXCITED DAMPED INVERTED PENDULUM 1. INTRODUCTION
Journa of Sound and Vibration (996) 98(5), 643 65 STABILITY OF A PARAMETRICALLY EXCITED DAMPED INVERTED PENDULUM G. ERDOS AND T. SINGH Department of Mechanica and Aerospace Engineering, SUNY at Buffao,
More informationSEMINAR 2. PENDULUMS. V = mgl cos θ. (2) L = T V = 1 2 ml2 θ2 + mgl cos θ, (3) d dt ml2 θ2 + mgl sin θ = 0, (4) θ + g l
Probem 7. Simpe Penduum SEMINAR. PENDULUMS A simpe penduum means a mass m suspended by a string weightess rigid rod of ength so that it can swing in a pane. The y-axis is directed down, x-axis is directed
More informationBP neural network-based sports performance prediction model applied research
Avaiabe onine www.jocpr.com Journa of Chemica and Pharmaceutica Research, 204, 6(7:93-936 Research Artice ISSN : 0975-7384 CODEN(USA : JCPRC5 BP neura networ-based sports performance prediction mode appied
More informationResearch Article Solution of Point Reactor Neutron Kinetics Equations with Temperature Feedback by Singularly Perturbed Method
Science and Technoogy of Nucear Instaations Voume 213, Artice ID 261327, 6 pages http://dx.doi.org/1.1155/213/261327 Research Artice Soution of Point Reactor Neutron Kinetics Equations with Temperature
More informationSelf Inductance of a Solenoid with a Permanent-Magnet Core
1 Probem Sef Inductance of a Soenoid with a Permanent-Magnet Core Kirk T. McDonad Joseph Henry Laboratories, Princeton University, Princeton, NJ 08544 (March 3, 2013; updated October 19, 2018) Deduce the
More informationTracking Control of Multiple Mobile Robots
Proceedings of the 2001 IEEE Internationa Conference on Robotics & Automation Seou, Korea May 21-26, 2001 Tracking Contro of Mutipe Mobie Robots A Case Study of Inter-Robot Coision-Free Probem Jurachart
More informationCS229 Lecture notes. Andrew Ng
CS229 Lecture notes Andrew Ng Part IX The EM agorithm In the previous set of notes, we taked about the EM agorithm as appied to fitting a mixture of Gaussians. In this set of notes, we give a broader view
More informationAn H 2 type Riemannian metric on the space of planar curves
An H 2 type Riemannian metric on the space of panar curves Jayant hah Mathematics Department, Northeastern University, Boston MA emai: shah@neu.edu Abstract An H 2 type metric on the space of panar curves
More informationSimplified analysis of EXAFS data and determination of bond lengths
Indian Journa of Pure & Appied Physics Vo. 49, January 0, pp. 5-9 Simpified anaysis of EXAFS data and determination of bond engths A Mishra, N Parsai & B D Shrivastava * Schoo of Physics, Devi Ahiya University,
More informationMATH 172: MOTIVATION FOR FOURIER SERIES: SEPARATION OF VARIABLES
MATH 172: MOTIVATION FOR FOURIER SERIES: SEPARATION OF VARIABLES Separation of variabes is a method to sove certain PDEs which have a warped product structure. First, on R n, a inear PDE of order m is
More informationCopyright information to be inserted by the Publishers. Unsplitting BGK-type Schemes for the Shallow. Water Equations KUN XU
Copyright information to be inserted by the Pubishers Unspitting BGK-type Schemes for the Shaow Water Equations KUN XU Mathematics Department, Hong Kong University of Science and Technoogy, Cear Water
More information(Refer Slide Time: 2:34) L C V
Microwave Integrated Circuits Professor Jayanta Mukherjee Department of Eectrica Engineering Indian Intitute of Technoogy Bombay Modue 1 Lecture No 2 Refection Coefficient, SWR, Smith Chart. Heo wecome
More informationMath 124B January 31, 2012
Math 124B January 31, 212 Viktor Grigoryan 7 Inhomogeneous boundary vaue probems Having studied the theory of Fourier series, with which we successfuy soved boundary vaue probems for the homogeneous heat
More informationLegendre Polynomials - Lecture 8
Legendre Poynomias - Lecture 8 Introduction In spherica coordinates the separation of variabes for the function of the poar ange resuts in Legendre s equation when the soution is independent of the azimutha
More informationData Search Algorithms based on Quantum Walk
Data Search Agorithms based on Quantum Wak Masataka Fujisaki, Hiromi Miyajima, oritaka Shigei Abstract For searching any item in an unsorted database with items, a cassica computer takes O() steps but
More informationA NOVEL WAVELET-GALERKIN METHOD FOR MOD- ELING RADIO WAVE PROPAGATION IN TROPO- SPHERIC DUCTS
Progress In Eectromagnetics Research B, Vo 36, 35 52, 2012 A NOVEL WAVELET-GALERKIN METHOD FOR MOD- ELING RADIO WAVE PROPAGATION IN TROPO- SPHERIC DUCTS A Iqba * and V Jeoti Department of Eectrica and
More informationMath 124B January 17, 2012
Math 124B January 17, 212 Viktor Grigoryan 3 Fu Fourier series We saw in previous ectures how the Dirichet and Neumann boundary conditions ead to respectivey sine and cosine Fourier series of the initia
More informationParallel-Axis Theorem
Parae-Axis Theorem In the previous exampes, the axis of rotation coincided with the axis of symmetry of the object For an arbitrary axis, the paraeaxis theorem often simpifies cacuations The theorem states
More informationJackson 4.10 Homework Problem Solution Dr. Christopher S. Baird University of Massachusetts Lowell
Jackson 4.10 Homework Probem Soution Dr. Christopher S. Baird University of Massachusetts Lowe PROBLEM: Two concentric conducting spheres of inner and outer radii a and b, respectivey, carry charges ±.
More informationMathematical Scheme Comparing of. the Three-Level Economical Systems
Appied Mathematica Sciences, Vo. 11, 2017, no. 15, 703-709 IKAI td, www.m-hikari.com https://doi.org/10.12988/ams.2017.7252 Mathematica Scheme Comparing of the Three-eve Economica Systems S.M. Brykaov
More informationLecture Notes for Math 251: ODE and PDE. Lecture 34: 10.7 Wave Equation and Vibrations of an Elastic String
ecture Notes for Math 251: ODE and PDE. ecture 3: 1.7 Wave Equation and Vibrations of an Eastic String Shawn D. Ryan Spring 212 ast Time: We studied other Heat Equation probems with various other boundary
More informationIntroduction. Figure 1 W8LC Line Array, box and horn element. Highlighted section modelled.
imuation of the acoustic fied produced by cavities using the Boundary Eement Rayeigh Integra Method () and its appication to a horn oudspeaer. tephen Kirup East Lancashire Institute, Due treet, Bacburn,
More informationSection 6: Magnetostatics
agnetic fieds in matter Section 6: agnetostatics In the previous sections we assumed that the current density J is a known function of coordinates. In the presence of matter this is not aways true. The
More informationNumerical Simulation of Launch Tube Based on Container type Missile Launch Technology
Avaiabe onine at www.scienceirect.com Proceia Engineering 31 (2012) 302 307 Internationa Conference on Avances in Computationa oeing an Simuation Numerica Simuation of Launch Tube Base on Container type
More informationScalable Spectrum Allocation for Large Networks Based on Sparse Optimization
Scaabe Spectrum ocation for Large Networks ased on Sparse Optimization innan Zhuang Modem R&D Lab Samsung Semiconductor, Inc. San Diego, C Dongning Guo, Ermin Wei, and Michae L. Honig Department of Eectrica
More informationPhysics 505 Fall Homework Assignment #4 Solutions
Physics 505 Fa 2005 Homework Assignment #4 Soutions Textbook probems: Ch. 3: 3.4, 3.6, 3.9, 3.0 3.4 The surface of a hoow conducting sphere of inner radius a is divided into an even number of equa segments
More information6.434J/16.391J Statistics for Engineers and Scientists May 4 MIT, Spring 2006 Handout #17. Solution 7
6.434J/16.391J Statistics for Engineers and Scientists May 4 MIT, Spring 2006 Handout #17 Soution 7 Probem 1: Generating Random Variabes Each part of this probem requires impementation in MATLAB. For the
More informationSymbolic models for nonlinear control systems using approximate bisimulation
Symboic modes for noninear contro systems using approximate bisimuation Giordano Poa, Antoine Girard and Pauo Tabuada Abstract Contro systems are usuay modeed by differentia equations describing how physica
More informationPHYS 110B - HW #1 Fall 2005, Solutions by David Pace Equations referenced as Eq. # are from Griffiths Problem statements are paraphrased
PHYS 110B - HW #1 Fa 2005, Soutions by David Pace Equations referenced as Eq. # are from Griffiths Probem statements are paraphrased [1.] Probem 6.8 from Griffiths A ong cyinder has radius R and a magnetization
More informationPHYSICS LOCUS / / d dt. ( vi) mass, m moment of inertia, I. ( ix) linear momentum, p Angular momentum, l p mv l I
6 n terms of moment of inertia, equation (7.8) can be written as The vector form of the above equation is...(7.9 a)...(7.9 b) The anguar acceeration produced is aong the direction of appied externa torque.
More informationScienceDirect. Numerical modelling of debris bed water quenching
Avaiabe onine at www.sciencedirect.com ScienceDirect Procedia IUTAM 15 2015 64 71 IUTAM Symposium on Mutiphase fows with phase change: chaenges and opportunities, Hyderabad, India December 08 December
More informationMeasurement of acceleration due to gravity (g) by a compound pendulum
Measurement of acceeration due to gravity (g) by a compound penduum Aim: (i) To determine the acceeration due to gravity (g) by means of a compound penduum. (ii) To determine radius of gyration about an
More informationIntroduction to Simulation - Lecture 13. Convergence of Multistep Methods. Jacob White. Thanks to Deepak Ramaswamy, Michal Rewienski, and Karen Veroy
Introduction to Simuation - Lecture 13 Convergence of Mutistep Methods Jacob White Thans to Deepa Ramaswamy, Micha Rewiensi, and Karen Veroy Outine Sma Timestep issues for Mutistep Methods Loca truncation
More informationarxiv: v1 [math.ca] 6 Mar 2017
Indefinite Integras of Spherica Besse Functions MIT-CTP/487 arxiv:703.0648v [math.ca] 6 Mar 07 Joyon K. Boomfied,, Stephen H. P. Face,, and Zander Moss, Center for Theoretica Physics, Laboratory for Nucear
More informationGaussian Curvature in a p-orbital, Hydrogen-like Atoms
Advanced Studies in Theoretica Physics Vo. 9, 015, no. 6, 81-85 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/10.1988/astp.015.5115 Gaussian Curvature in a p-orbita, Hydrogen-ike Atoms Sandro-Jose Berrio-Guzman
More informationFEMAG: A High Performance Parallel Finite Element Toolbox for Electromagnetic Computations
Internationa Journa of Energy and Power Engineering 2016; 5(1-1): 57-64 Pubished onine November 28, 2015 (http://www.sciencepubishinggroup.com/j/ijepe) doi: 10.11648/j.ijepe.s.2016050101.19 ISSN: 2326-957X
More informationSymmetric Beauty: Three-steps Force Amplified Clamping Device Based on Lever-Toggle Mechanism Driven by Pneumatic Muscle
06 Internationa Conference on Computer, Mechatronics and Eectronic Engineering (CMEE 06) ISBN: 978--60595-406-6 Symmetric Beauty: Three-steps orce Ampified Camping Device Based on ever-togge Mechanism
More informationNumerical solution of one dimensional contaminant transport equation with variable coefficient (temporal) by using Haar wavelet
Goba Journa of Pure and Appied Mathematics. ISSN 973-1768 Voume 1, Number (16), pp. 183-19 Research India Pubications http://www.ripubication.com Numerica soution of one dimensiona contaminant transport
More informationGeneral Certificate of Education Advanced Level Examination June 2010
Genera Certificate of Education Advanced Leve Examination June 2010 Human Bioogy HBI6T/P10/task Unit 6T A2 Investigative Skis Assignment Task Sheet The effect of temperature on the rate of photosynthesis
More informationTM Electromagnetic Scattering from 2D Multilayered Dielectric Bodies Numerical Solution
TM Eectromagnetic Scattering from D Mutiayered Dieectric Bodies Numerica Soution F. Seydou,, R. Duraiswami, N.A. Gumerov & T. Seppänen. Department of Eectrica and Information Engineering University of
More informationTIME DEPENDENT TEMPERATURE DISTRIBUTION MODEL IN LAYERED HUMAN DERMAL PART
VOL. 8, No. II, DECEMBER, 0, pp 66-76 TIME DEPENDENT TEMPERATURE DISTRIBUTION MODEL IN LAYERED HUMAN DERMAL PART Saraswati Acharya*, D. B. Gurung, V. P. Saxena Department of Natura Sciences (Mathematics),
More informationMinimizing Total Weighted Completion Time on Uniform Machines with Unbounded Batch
The Eighth Internationa Symposium on Operations Research and Its Appications (ISORA 09) Zhangiaie, China, September 20 22, 2009 Copyright 2009 ORSC & APORC, pp. 402 408 Minimizing Tota Weighted Competion
More informationUnconditional security of differential phase shift quantum key distribution
Unconditiona security of differentia phase shift quantum key distribution Kai Wen, Yoshihisa Yamamoto Ginzton Lab and Dept of Eectrica Engineering Stanford University Basic idea of DPS-QKD Protoco. Aice
More informationHomework #04 Answers and Hints (MATH4052 Partial Differential Equations)
Homework #4 Answers and Hints (MATH452 Partia Differentia Equations) Probem 1 (Page 89, Q2) Consider a meta rod ( < x < ), insuated aong its sides but not at its ends, which is initiay at temperature =
More informationRELUCTANCE The resistance of a material to the flow of charge (current) is determined for electric circuits by the equation
INTRODUCTION Magnetism pays an integra part in amost every eectrica device used today in industry, research, or the home. Generators, motors, transformers, circuit breakers, teevisions, computers, tape
More informationExperimental Investigation and Numerical Analysis of New Multi-Ribbed Slab Structure
Experimenta Investigation and Numerica Anaysis of New Muti-Ribbed Sab Structure Jie TIAN Xi an University of Technoogy, China Wei HUANG Xi an University of Architecture & Technoogy, China Junong LU Xi
More informationHILBERT? What is HILBERT? Matlab Implementation of Adaptive 2D BEM. Dirk Praetorius. Features of HILBERT
Söerhaus-Workshop 2009 October 16, 2009 What is HILBERT? HILBERT Matab Impementation of Adaptive 2D BEM joint work with M. Aurada, M. Ebner, S. Ferraz-Leite, P. Godenits, M. Karkuik, M. Mayr Hibert Is
More informationMathematical Model for Potassium Release from Polymer-coated Fertiliser
ARTICLE IN PRESS Biosystems Engineering (2004) 88 (3), 395 400 Avaiabe onine at www.sciencedirect.com doi:10.1016/j.biosystemseng.2004.03.004 SW}Soi and Water Mathematica Mode for Potassium Reease from
More information1.2 Partial Wave Analysis
February, 205 Lecture X.2 Partia Wave Anaysis We have described scattering in terms of an incoming pane wave, a momentum eigenet, and and outgoing spherica wave, aso with definite momentum. We now consider
More informationSeparation of Variables and a Spherical Shell with Surface Charge
Separation of Variabes and a Spherica She with Surface Charge In cass we worked out the eectrostatic potentia due to a spherica she of radius R with a surface charge density σθ = σ cos θ. This cacuation
More informationFormulas for Angular-Momentum Barrier Factors Version II
BNL PREPRINT BNL-QGS-06-101 brfactor1.tex Formuas for Anguar-Momentum Barrier Factors Version II S. U. Chung Physics Department, Brookhaven Nationa Laboratory, Upton, NY 11973 March 19, 2015 abstract A
More informationA GENERAL METHOD FOR EVALUATING OUTAGE PROBABILITIES USING PADÉ APPROXIMATIONS
A GENERAL METHOD FOR EVALUATING OUTAGE PROBABILITIES USING PADÉ APPROXIMATIONS Jack W. Stokes, Microsoft Corporation One Microsoft Way, Redmond, WA 9852, jstokes@microsoft.com James A. Ritcey, University
More informationUI FORMULATION FOR CABLE STATE OF EXISTING CABLE-STAYED BRIDGE
UI FORMULATION FOR CABLE STATE OF EXISTING CABLE-STAYED BRIDGE Juan Huang, Ronghui Wang and Tao Tang Coege of Traffic and Communications, South China University of Technoogy, Guangzhou, Guangdong 51641,
More informationProbabilistic Assessment of Total Transfer Capability Using SQP and Weather Effects
J Eectr Eng Techno Vo. 9, No. 5: 52-526, 24 http://dx.doi.org/.537/jeet.24.9.5.52 ISSN(Print) 975-2 ISSN(Onine) 293-7423 Probabiistic Assessment of Tota Transfer Capabiity Using SQP and Weather Effects
More informationAPARTIALLY covered trench located at the corner of
IEEE TRANSACTIONS ON ANTENNAS AND PROPAGATION, Penetration, Radiation and Scattering for a Cavity-backed Gap in a Corner Danio Erricoo, Senior Member, IEEE, Piergiorgio L. E. Usenghi, Feow, IEEE Abstract
More informationEXISTENCE OF SOLUTIONS FOR ONE-DIMENSIONAL WAVE EQUATIONS WITH NONLOCAL CONDITIONS
Eectronic Journa of Differentia Equations, Vo. 21(21), No. 76, pp. 1 8. ISSN: 172-6691. URL: http://ejde.math.swt.edu or http://ejde.math.unt.edu ftp ejde.math.swt.edu (ogin: ftp) EXISTENCE OF SOLUTIONS
More informationSolving Maxwell s Equations Using the Ultra Weak Variational Formulation
Soving Maxwe s Equations Using the Utra Weak Variationa Formuation T. Huttunen, M. Mainen and P. Monk Department of Appied Physics, University of Kuopio, P.O.Box 1627, 7211 Kuopio, Finand Department of
More informationDIGITAL FILTER DESIGN OF IIR FILTERS USING REAL VALUED GENETIC ALGORITHM
DIGITAL FILTER DESIGN OF IIR FILTERS USING REAL VALUED GENETIC ALGORITHM MIKAEL NILSSON, MATTIAS DAHL AND INGVAR CLAESSON Bekinge Institute of Technoogy Department of Teecommunications and Signa Processing
More informationMigration of Ground Penetrating Radar data in heterogeneous and dispersive media
New Strategies for European Remote Sensing, Oui (ed.) 25 Mipress, Rotterdam, ISBN 9 5966 3 X Migration of Ground Penetrating Radar data in heterogeneous and dispersive media Armando R. Sena, Pau L. Stoffa
More information