A constant recursive convolution technique for frequency dependent scalar wave equation based FDTD algorithm
|
|
- Simon Rogers
- 6 years ago
- Views:
Transcription
1 J Comput Electron (213) 12: DOI 1.17/s A constant recursve convoluton technque for frequency dependent scalar wave equaton bed FDTD algorthm M. Burak Özakın Serkan Aksoy Publshed onlne: 8 June 213 Sprnger ScenceBusness Meda New York 213 Abstract A scalar wave equaton bed recursve convoluton fnte-dfference tme-doman algorthm s developed for a frequency-dependent Debye medum n ths paper. Ths algorthm s bed on a recursve evaluaton of a convoluton ntegral n the tme doman. A numercal example s presented for a problem of wde-band reflecton from an ar-water nterface. The obtaned results are compared wth an analytcal soluton. The excellent agreement s observed between the numercal results. Memory and computatonal tme advantages of the proposed method over Maxwell s equatons bed soluton are also shown. Keywords Dspersve medum Fnte dfference tme doman method Recursve convoluton 1 Introducton The Fnte Dfference Tme Doman (FDTD) method s an effectve numercal tool for soluton of complex electromagnetc problems 1]. In partcular, the FDTD method s n rapd progress for modelng the dspersve materals such Debye, Lorentz types, etc. Fundamentally, three man Maxwell s equatons bed approaches n the dspersve applcaton of the FDTD method are Recursve Convoluton (RC), Auxlary Dfferental Equaton and Z-Transform technques. The RC technque h dfferent forms of Constant Recursve Convoluton (CRC) 2], Pecewse Lnear Recursve Convoluton 3], Pecewse Constant Recursve Convoluton 4], Trapezodal Recursve Convoluton 5], JE M.B. Özakın ( ) S. Aksoy Department of Electroncs Engneerng, Gebze Insttute of Technology, Kocael, Turkey e-mal: bozakn@gyte.edu.tr Convoluton 6] and Pecewse Lnear JE Recursve Convoluton 7]. Soluton of wave propagaton problems n a dspersve medum s a crucal topc because the wave velocty s a frequency dependent parameter n the dspersve medum. The materal dsperson strongly affects propagated wdeband sgnals. Ths leads to dstorton of ths sgnal n dspersve meda. In partcular, electromagnetc wave propagaton n the dspersve medum can be modeled ether wth the Maxwell s equatons or the wave equaton. Although the dfferent Maxwell s equatons bed approaches for the FDTD dspersve modelng are present mentoned before, the wave equaton bed approaches for the dspersve modelng are not wdespread n the lterature. However, the wave equaton s used n many dfferent scentfc dscplnes such acoustcs 8], flud dynamcs 9] and so on. It also smplfes an electromagnetc problem because of usng only one electromagnetc feld component n the soluton. Ths s especally mportant for applcaton of numercal methods because both calculaton tme and memory requrements become less. In ths paper, the electromagnetc wave propagaton n a lnear dspersve Debye medum s nvestgated by the scalar wave equaton bed recursve convoluton fnte-dfference tme-doman method. For ths am, frst, a recursve accumulator s extracted. Then, an FDTD update equaton bed on the CRC algorthm (generally known the RC technque) for the scalar wave equaton of the Debye medum s formulated and valdated wth an analytcal soluton of an example problem. Ths paper s organzed follows: In Sect. 2, a scalar wave equaton bed formulaton of the CRC-FDTD method s gven. In Sect. 3, a numercal example of an ar-water nterface for the Debye medum s shown by a comparson of an analytcal soluton. Memory and computatonal tme
2 J Comput Electron (213) 12: comparsons between the proposed and Maxwell s equatons bed solutons are also gven. Secton 4 gves summary and fnal remarks. 2 Formulaton of the problem A Fourer transformed Helmholtz equaton n a onedmensonal Cartesan coordnate s u(x, ω) k 2 (ω)u(x, ω) (1) where s the Laplace operator, u(x, ω) s a feld component and ω s an angular frequency. The key pont n ths equaton about treatment of the wave number k k(ω) for dspersve meda k 2 (ω) ω 2 ε(ω)μ. (2) Here, only, the frequency-dependent delectrc permttvty ε(ω) s used for dspersve medum modelng. The magnetc permeablty of the medum μ s sumed to be equal to a free space permeablty μ. In ths ce, ε(ω) can be wrtten n a more general form ε(ω) ε ε χ e (ω) ] (3) where ε and ε are the delectrc permttvty of the free space and the delectrc permttvty at a hgh-frequency lmt, respectvely. χ e (ω) shows the delectrc susceptblty of the dspersve medum. Substtutng Eq. (3) over Eq. (2) nto Eq. (1), the scalar wave equaton becomes u(x, ω) ε ε μ ω 2 u(x, ω) ε μ ω 2 χ e (ω)u(x, ω). (4) Here, one can note n Eq. (4) that the multplcaton of two functons n the frequency doman s present. It means that a convoluton ntegral wll appear n the tme doman. Applyng the nverse Fourer transform through Eq. (4) leads to the followng equaton u(x, t) ε 2 c 2 u(x, t) t2 1 2 { u(x, t) χe c 2 t 2 (t) } (5) where c 1/ ε μ s the speed of lght n the free space and the convoluton term can be wrtten u(x, t) ε 2 c 2 u(x, t) t2 1 2 { t } c 2 t 2 u(x, t τ)χ e (τ)dτ. (6) Now, let us defne the convoluton ntegral for smplcty ξ(x,t) t u(x, t τ)χ e (τ)dτ. (7) Therefore Eq. (5) can be wrtten u(x, t) ε c 2 2 t 2 u(x, t) 1 c 2 2 ξ(x,t). (8) t2 In order to extract the FDTD update equaton for the dspersve medum, the feld must be dscretzed over the tme and space varables u(x, t) u( x, n t). Here, and n are the ndces of dscrete space and tme varables, respectvely. x and t show the unt spatal and tme steps, respectvely. In ths ce, the scalar wave equaton becomes n ε 2 c 2 t 2 n 1 2 c 2 t 2 ξ n (9) where the man concern s to treat the convoluton ntegral n a dscrete form. In order to overcome ths problem, the feld component n the convoluton ntegral s sumed to be constant (the CRC technque) over the each tme nterval t. Therefore, the convoluton ntegral n the dscrete form wll be rearranged ξ n n t n1 n tτ χ e τ dτ (m1) t nm n1 χ e τ dτ nm χ e m. (1) Now, central fnte dfferences for the spatal and tme dervatves can be appled nto Eq. (9). Then, the FDTD update equaton for the scalar wave equaton becomes n 1 2n n 1 x 2 1 ξ n1 c 2 where ε n1 2 n n1 c 2 t 2 2ξ n ξ n1 t 2 (11) ξ n1 n1m χ e m, (12) n1 ξ n nm χ e m n χ e ξ n1 n n1m χ e m ( n1 χ e n χ e 1 χ e m1, (13) ) χ e m2. (14)
3 754 J Comput Electron (213) 12: Expresson of ξ n1 2ξ n ξ n1 n Eq. (11) can be wrtten by usng Eq. (12), Eq. (13) and Eq. (14) ξ n1 2ξ n ξ n1 ( n1 χ e n χ e 1 2n χ e ) n1m ( χe m2 2χ e m1 χ e m ). (15) Let us defne a parameter of χ e for smplcty χ e m χ e m2 2χ e m1 χ e m. (16) Then Eq. (15) wll be ξ n1 2ξ n ξ n1 n1 χ e n χ e 1 2n χ e χ e m. (17) Substtutng Eq. (17) nto Eq. (11), the dscrete scalar wave equaton becomes n 1 2n n 1 x 2 1 ( n1 ε n1 c 2 2 n n1 t 2 χ e n χ e 1 2n χ e ) 1 n1m χ e m. (18) After more smplfcatons n Eq. (18), the FDTD update equaton can be rearranged n1 x 2 (ε χ e )n 1 x 2 (ε χ e )n 1 x2 (2ε χ e 1 2χ e ) 2(c t)2 x 2 (ε χ e ) n ε ε χ e n1 1 ε χ e n1m χ e m (19) where, at a frst glance, t seems that the teraton progress n the summaton term wll requre storng a large number of pt feld tme values. However, n order to crcumvent ths problem, t s shown that a recursve relaton for χ e m can be obtanable for the Debye medum. Let us consder the Debye medum wth ts susceptblty functon χ e (ω) ε s ε (2) 1 jω where ε s s the statc delectrc permttvty at the low frequences and s the characterstc relaxaton tme of the medum (j 1). After applcaton of the nverse Fourer transform to Eq. (2), the tme doman representaton of χ e (t) wll be ( ) εs ε χ e (t) e t H(t) (21) where H(t) s the unt step functon, whch shows the satsfacton of the causalty prncple. χ e (t) can be used for the calculaton of χ e m n Eq. (16) χ e m (m1) t (m1) t χ e τ dτ ( ε s ε (ε s ε )e )e τ dτ Thus, by usng Eq. (22), Eq. (16) becomes χ e m (ε s ε )e t 1 3e (ε s ε )e 1 e t ]. (22) 3e 2 t e 3 t ] 1 e t ] 3. (23) Now, let us try to fnd the recursve relaton for the convoluton term n Eq. (1). For ths am, defne a new parameter ψ n ψ n χ e m. (24) One can examne a few tme steps for ψ n startng wth n 2 ψ 2 1m χ e m 1 χ e (25) and for n 3 1 ψ 3 2m χ e m 2 χ e 1 χ e 1 (26) where χ e m1 can be lnked to χ e m by usng Eq. (23) χ e m1 e t χ e m. (27) Then, the recursve relaton for ψ 3 wll be fnally obtaned ψ 3 2 χ e e t 1 χ e 2 χ e e t ψ 2 (28) where t s clear that a next tme value of ψ s found from ts prevous tme value. Therefore, a recursve update equaton can be reformulated n a more general form for n 3 ψ n n1 χ e e t ψ n1 (29)
4 J Comput Electron (213) 12: Fg. 2 Demonstraton of the ncdent Gaussan pulse whch propagates toward to the water regon at the 8th tme step Fg. 1 The flowchart of the scalar wave equaton bed CRC-FDTD method where ψ ψ 1 ψ 2. Here, t s worthy to note that the recursve update equaton accumulator for ψ n s extracted specally for an exponental form of χ e (t). As a result, substtutng ths recursve equaton nto Eq. (19), the fnal form of the FDTD update equaton for the scalar wave equaton becomes n1 x 2 (ε χ e )n 1 x 2 (ε χ e )n 1 x2 (2ε χ e 1 2χ e ) 2(c t)2 x 2 (ε χ e ) n ε ε χ e n1 1 ε χ e ψ n (3) where t s of nterest to note that the two-tme steps of χ e are present n the formulaton. However, χ e 1 can be calculated from χ e by usng Eq. (22). A flowchart of the scalar wave equaton bed CRC- FDTD algorthm can be drawn n Fg Numercal example wth valdaton In ths secton, the proposed method s valdated by an example problem of the wde-band reflecton from an ar-water nterface. The ar regon s modeled the free-space medum wth the parameters of ε 1, ε s 1, ε 1/36π 1 9 F/m, μ 4π 1 7 H/m and s. The water s modeled the Debye type dspersve medum wth the parameters of ε 1.8, ε s 81 and s. The well known CFL lmt s used n order to keep the stablty of the FDTD soluton. In addton, the problem space s termnated wth a frst order Mur type absorbng boundary condton. A one-dmensonal numercal example s consdered n 2]. It conssts of 1 cells separated by 499 cells and 51 cells for the free-space and water, respectvely. Each cell h equal length of x 37.5 µm and t.625 ps. A normalzed Gaussan pulse shown n Fg. 2 s used an exctaton source wth a spatal wdth of 4 cells and located at 1th cell ( ) πn t H(n t) n1 1 sn 2 ] H(n t T) (31) T where T s a spatal wdth of the pulse. After mpngng the electromagnetc wave on the arwater nterface, the pulse s separated nto the reflected and propagated parts. Ths phenomenon can be clearly seen n Fg. 3. In Fg. 4, a comparson between the proposed FDTD soluton and the analytcal soluton 2] of ths problem s shown for the reflecton coeffcent at 2th tme step. The reflecton coeffcent of the numercal example w obtaned by applyng a Ft Fourer Transform to the tme doman values of an observaton pont. An excellent agreement s observed between the results. In Table 1, the memory and the computatonal tme advantages of the proposed method over the Maxwell s equatons bed soluton are shown for the one-dmensonal example gven above. The comparsons are performed by an Intel Core2 Quad computer of 2.4 GHz CPU speed. The proposed formulaton needs % 33 less memory wth less computaton tme. Ths s due to the usng one feld component
5 756 J Comput Electron (213) 12: Summary and fnal remarks Fg. 3 Demonstraton of the pulses whch reflected from the nterface and propagates nto the water at the 2th tme step In ths study, the scalar wave equaton bed recursve convoluton formulaton of the frequency dependent fnte dfference tme doman method s gven for the lnear dspersve Debye medum. In the frst step, the scalar wave equaton n the tme doman for the dspersve medum s extracted from the Helmholtz equaton n the frequency doman. Relaton of the frequency-dependent parameters n that equaton s clearly revealed. Then, the nverse Fourer transform s employed n order to obtan the tme doman form. Ths leads to the convoluton ntegral n the formulaton. A tme doman recursve accumulator for ths convoluton ntegral s extracted and lnked wth the FDTD algorthm. The proposed formulaton s valdated wth an analytcal result over the numercal example of the wde-band reflecton coeffcent from the ar regon to water. The excellent agreement s observed between the results. The memory and the computatonal tme advantages of the proposed method over the Maxwell s equatons bed soluton are also shown. For the future work, we are plannng to extend the proposed scalar wave equaton bed FDTD formulaton for other types of the dspersve medums such Drude, Lorentz and so on wth ncludng the lossy term. References Fg. 4 Comparson of the frequency-dependent pulse reflecton coeffcent between the analytcal and the FDTD solutons Table 1 Comparson of the memory and the computatonal tme between the proposed and Maxwell s equatons bed solutons n onedmensonal ce Method (1 D) Memory (kb) Computatonal tme (s) Proposed Maxwell s eq. bed n the soluton, although the Maxwell s equatons bed soluton needs two feld components. These advantages wll be domnant n the ce of two and three dmensonal problems. 1. Taflove, A., Hagness, S.C.: Computatonal Electrodynamcs: The Fnte-Dfference Tme-Doman Method, 3rd edn. Artech House, London (25) 2. Luebbers, R., Hunsberger, F.P., Kunz, K.S., Standler, R.B., Schneder, M.: A frequency-dependent fnte-dfference tme-doman formulaton for dspersve materals. IEEE Trans. Electromagn. Compat. 32(3), (199) 3. Kelley, D.F., Luebbers, R.J.: Pecewse lnear recursve convoluton for dspersve meda usng FDTD. IEEE Trans. Antenn Propag. 44(6), (1996) 4. Schuster, J.W., Luebbers, R.J.: An accurate FDTD algorthm for dspersve meda usng a pecewse constant recursve convoluton technque. In: IEEE Antenn and Propagaton Socety Internatonal Symposum, Atlanta, GA, USA, June (1998) 5. Sushansan, R., LoVetr, J.: A comparson of numercal technques for modelng electromagnetc dspersve meda. IEEE Mcrow. Guded Wave Lett. 5(12), (1995) 6. Sushansan, R., LoVetr, J.: An FDTD formulaton for dspersve meda usng a current densty. IEEE Trans. Antenn Propag. 46(1), (1998) 7. Sushansan, R., LoVetr, J.: A Novel FDTD Formulaton for Dspersve Meda. IEEE Mcrow. Wrel. Compon. Lett. 13(5), (23) 8. Knsler, L.E., Frey, A.R., Coppens, A.B., Sanders, J.V.: Fundamentals of Acoustcs, 4th edn. Wley, New York (2) 9. Cebec, T., Shao, J.P., Kafyeke, F., Laurendeau, E.: Computatonal Flud Dynamcs for Engneers. Sprnger, Berln (25)
B. H. Jung Department of Information and Communication Engineering Hoseo University Asan, Chungnam , Korea
Progress In Electromagnetcs Research, PIER 79, 339 352, 2008 SOLVING TIME DOMAIN HELMHOLTZ WAVE EQUATION WITH MOD-FDM B. H. Jung Department of Informaton and Communcaton Engneerng Hoseo Unversty Asan,
More informationCHAPTER 5 NUMERICAL EVALUATION OF DYNAMIC RESPONSE
CHAPTER 5 NUMERICAL EVALUATION OF DYNAMIC RESPONSE Analytcal soluton s usually not possble when exctaton vares arbtrarly wth tme or f the system s nonlnear. Such problems can be solved by numercal tmesteppng
More informationNote: Please use the actual date you accessed this material in your citation.
MIT OpenCourseWare http://ocw.mt.edu 6.13/ESD.13J Electromagnetcs and Applcatons, Fall 5 Please use the followng ctaton format: Markus Zahn, Erch Ippen, and Davd Staeln, 6.13/ESD.13J Electromagnetcs and
More information(Online First)A Lattice Boltzmann Scheme for Diffusion Equation in Spherical Coordinate
Internatonal Journal of Mathematcs and Systems Scence (018) Volume 1 do:10.494/jmss.v1.815 (Onlne Frst)A Lattce Boltzmann Scheme for Dffuson Equaton n Sphercal Coordnate Debabrata Datta 1 *, T K Pal 1
More informationHomework 4. 1 Electromagnetic surface waves (55 pts.) Nano Optics, Fall Semester 2015 Photonics Laboratory, ETH Zürich
Homework 4 Contact: frmmerm@ethz.ch Due date: December 04, 015 Nano Optcs, Fall Semester 015 Photoncs Laboratory, ETH Zürch www.photoncs.ethz.ch The goal of ths problem set s to understand how surface
More informationAppendix B. The Finite Difference Scheme
140 APPENDIXES Appendx B. The Fnte Dfference Scheme In ths appendx we present numercal technques whch are used to approxmate solutons of system 3.1 3.3. A comprehensve treatment of theoretcal and mplementaton
More informationA new Approach for Solving Linear Ordinary Differential Equations
, ISSN 974-57X (Onlne), ISSN 974-5718 (Prnt), Vol. ; Issue No. 1; Year 14, Copyrght 13-14 by CESER PUBLICATIONS A new Approach for Solvng Lnear Ordnary Dfferental Equatons Fawz Abdelwahd Department of
More informationModule 3 LOSSY IMAGE COMPRESSION SYSTEMS. Version 2 ECE IIT, Kharagpur
Module 3 LOSSY IMAGE COMPRESSION SYSTEMS Verson ECE IIT, Kharagpur Lesson 6 Theory of Quantzaton Verson ECE IIT, Kharagpur Instructonal Objectves At the end of ths lesson, the students should be able to:
More informationECE 107: Electromagnetism
ECE 107: Electromagnetsm Set 8: Plane waves Instructor: Prof. Vtaly Lomakn Department of Electrcal and Computer Engneerng Unversty of Calforna, San Dego, CA 92093 1 Wave equaton Source-free lossless Maxwell
More informationNumerical Heat and Mass Transfer
Master degree n Mechancal Engneerng Numercal Heat and Mass Transfer 06-Fnte-Dfference Method (One-dmensonal, steady state heat conducton) Fausto Arpno f.arpno@uncas.t Introducton Why we use models and
More informationChapter 4 The Wave Equation
Chapter 4 The Wave Equaton Another classcal example of a hyperbolc PDE s a wave equaton. The wave equaton s a second-order lnear hyperbolc PDE that descrbes the propagaton of a varety of waves, such as
More information1-Dimensional Advection-Diffusion Finite Difference Model Due to a Flow under Propagating Solitary Wave
014 4th Internatonal Conference on Future nvronment and nergy IPCB vol.61 (014) (014) IACSIT Press, Sngapore I: 10.776/IPCB. 014. V61. 6 1-mensonal Advecton-ffuson Fnte fference Model ue to a Flow under
More informationComputational Electromagnetics in Antenna Analysis and Design
Computatonal Electromagnetcs n Antenna Analyss and Desgn Introducton It s rare for real-lfe EM problems to fall neatly nto a class that can be solved by the analytcal methods presented n the precedng lectures.
More informationOne-sided finite-difference approximations suitable for use with Richardson extrapolation
Journal of Computatonal Physcs 219 (2006) 13 20 Short note One-sded fnte-dfference approxmatons sutable for use wth Rchardson extrapolaton Kumar Rahul, S.N. Bhattacharyya * Department of Mechancal Engneerng,
More informationEffect of Losses in a Layered Structure Containing DPS and DNG Media
PIERS ONLINE, VOL. 4, NO. 5, 8 546 Effect of Losses n a Layered Structure Contanng DPS and DNG Meda J. R. Canto, S. A. Matos, C. R. Pava, and A. M. Barbosa Insttuto de Telecomuncações and Department of
More informationImplementation of the Matrix Method
Computatonal Photoncs, Prof. Thomas Pertsch, Abbe School of Photoncs, FSU Jena Computatonal Photoncs Semnar 0 Implementaton of the Matr Method calculaton of the transfer matr calculaton of reflecton and
More informationThree-dimensional eddy current analysis by the boundary element method using vector potential
Physcs Electrcty & Magnetsm felds Okayama Unversty Year 1990 Three-dmensonal eddy current analyss by the boundary element method usng vector potental H. Tsubo M. Tanaka Okayama Unversty Okayama Unversty
More informationNUMERICAL MODELING OF ACTIVE DEVICES CHARACTERIZED BY MEASURED S-PARAMETERS IN FDTD
Progress In Electromagnetcs Research, PIER 80, 381 392, 2008 NUMERICAL MODELING OF ACTIVE DEVICES CHARACTERIZED BY MEASURED S-PARAMETERS IN FDTD D. Y. Su, D. M. Fu, and Z. H. Chen Natonal Key Lab of Antennas
More informationSUPPLEMENTARY INFORMATION
do: 0.08/nature09 I. Resonant absorpton of XUV pulses n Kr + usng the reduced densty matrx approach The quantum beats nvestgated n ths paper are the result of nterference between two exctaton paths of
More informationNanoplasmonics FDTD Simulations Using a Generalized Dispersive Material Model
Nanoplasmoncs FDTD Smulatons Usng a Generalzed Dspersve Materal Model Ludmla J. Prokopeva, Alexander V. Kldshev, * Jeran Fang, Joshua Borneman, Mark D. Thoreson,,3 Vladmr M. Shalaev, and Vladmr P. Drachev
More informationResearch Article Green s Theorem for Sign Data
Internatonal Scholarly Research Network ISRN Appled Mathematcs Volume 2012, Artcle ID 539359, 10 pages do:10.5402/2012/539359 Research Artcle Green s Theorem for Sgn Data Lous M. Houston The Unversty of
More informationApplication of B-Spline to Numerical Solution of a System of Singularly Perturbed Problems
Mathematca Aeterna, Vol. 1, 011, no. 06, 405 415 Applcaton of B-Splne to Numercal Soluton of a System of Sngularly Perturbed Problems Yogesh Gupta Department of Mathematcs Unted College of Engneerng &
More informationVisco-Rubber Elastic Model for Pressure Sensitive Adhesive
Vsco-Rubber Elastc Model for Pressure Senstve Adhesve Kazuhsa Maeda, Shgenobu Okazawa, Koj Nshgch and Takash Iwamoto Abstract A materal model to descrbe large deformaton of pressure senstve adhesve (PSA
More informationA Hybrid Variational Iteration Method for Blasius Equation
Avalable at http://pvamu.edu/aam Appl. Appl. Math. ISSN: 1932-9466 Vol. 10, Issue 1 (June 2015), pp. 223-229 Applcatons and Appled Mathematcs: An Internatonal Journal (AAM) A Hybrd Varatonal Iteraton Method
More informationTransfer Functions. Convenient representation of a linear, dynamic model. A transfer function (TF) relates one input and one output: ( ) system
Transfer Functons Convenent representaton of a lnear, dynamc model. A transfer functon (TF) relates one nput and one output: x t X s y t system Y s The followng termnology s used: x y nput output forcng
More informationImplementation of the Matrix Method
Computatonal Photoncs, Prof. Thomas Pertsch, Abbe School of Photoncs, FSU Jena Computatonal Photoncs Semnar 0 Implementaton of the Matr Method calculaton of the transfer matr calculaton of reflecton and
More informationELASTIC WAVE PROPAGATION IN A CONTINUOUS MEDIUM
ELASTIC WAVE PROPAGATION IN A CONTINUOUS MEDIUM An elastc wave s a deformaton of the body that travels throughout the body n all drectons. We can examne the deformaton over a perod of tme by fxng our look
More informationNumerical Solutions of a Generalized Nth Order Boundary Value Problems Using Power Series Approximation Method
Appled Mathematcs, 6, 7, 5-4 Publshed Onlne Jul 6 n ScRes. http://www.scrp.org/journal/am http://.do.org/.436/am.6.77 umercal Solutons of a Generalzed th Order Boundar Value Problems Usng Power Seres Approxmaton
More informationLecture 13 APPROXIMATION OF SECOMD ORDER DERIVATIVES
COMPUTATIONAL FLUID DYNAMICS: FDM: Appromaton of Second Order Dervatves Lecture APPROXIMATION OF SECOMD ORDER DERIVATIVES. APPROXIMATION OF SECOND ORDER DERIVATIVES Second order dervatves appear n dffusve
More informationThe Exact Formulation of the Inverse of the Tridiagonal Matrix for Solving the 1D Poisson Equation with the Finite Difference Method
Journal of Electromagnetc Analyss and Applcatons, 04, 6, 0-08 Publshed Onlne September 04 n ScRes. http://www.scrp.org/journal/jemaa http://dx.do.org/0.46/jemaa.04.6000 The Exact Formulaton of the Inverse
More informationCOMPOSITE BEAM WITH WEAK SHEAR CONNECTION SUBJECTED TO THERMAL LOAD
COMPOSITE BEAM WITH WEAK SHEAR CONNECTION SUBJECTED TO THERMAL LOAD Ákos Jósef Lengyel, István Ecsed Assstant Lecturer, Professor of Mechancs, Insttute of Appled Mechancs, Unversty of Mskolc, Mskolc-Egyetemváros,
More informationEEE 241: Linear Systems
EEE : Lnear Systems Summary #: Backpropagaton BACKPROPAGATION The perceptron rule as well as the Wdrow Hoff learnng were desgned to tran sngle layer networks. They suffer from the same dsadvantage: they
More informationDigital Signal Processing
Dgtal Sgnal Processng Dscrete-tme System Analyss Manar Mohasen Offce: F8 Emal: manar.subh@ut.ac.r School of IT Engneerng Revew of Precedent Class Contnuous Sgnal The value of the sgnal s avalable over
More informationAnalytical Gradient Evaluation of Cost Functions in. General Field Solvers: A Novel Approach for. Optimization of Microwave Structures
IMS 2 Workshop Analytcal Gradent Evaluaton of Cost Functons n General Feld Solvers: A Novel Approach for Optmzaton of Mcrowave Structures P. Harscher, S. Amar* and R. Vahldeck and J. Bornemann* Swss Federal
More informationFormal solvers of the RT equation
Formal solvers of the RT equaton Formal RT solvers Runge- Kutta (reference solver) Pskunov N.: 979, Master Thess Long characterstcs (Feautrer scheme) Cannon C.J.: 970, ApJ 6, 55 Short characterstcs (Hermtan
More informationEVALUATION OF THE VISCO-ELASTIC PROPERTIES IN ASPHALT RUBBER AND CONVENTIONAL MIXES
EVALUATION OF THE VISCO-ELASTIC PROPERTIES IN ASPHALT RUBBER AND CONVENTIONAL MIXES Manuel J. C. Mnhoto Polytechnc Insttute of Bragança, Bragança, Portugal E-mal: mnhoto@pb.pt Paulo A. A. Perera and Jorge
More informationModule 1 : The equation of continuity. Lecture 1: Equation of Continuity
1 Module 1 : The equaton of contnuty Lecture 1: Equaton of Contnuty 2 Advanced Heat and Mass Transfer: Modules 1. THE EQUATION OF CONTINUITY : Lectures 1-6 () () () (v) (v) Overall Mass Balance Momentum
More informationNON-CENTRAL 7-POINT FORMULA IN THE METHOD OF LINES FOR PARABOLIC AND BURGERS' EQUATIONS
IJRRAS 8 (3 September 011 www.arpapress.com/volumes/vol8issue3/ijrras_8_3_08.pdf NON-CENTRAL 7-POINT FORMULA IN THE METHOD OF LINES FOR PARABOLIC AND BURGERS' EQUATIONS H.O. Bakodah Dept. of Mathematc
More informationInductance Calculation for Conductors of Arbitrary Shape
CRYO/02/028 Aprl 5, 2002 Inductance Calculaton for Conductors of Arbtrary Shape L. Bottura Dstrbuton: Internal Summary In ths note we descrbe a method for the numercal calculaton of nductances among conductors
More informationarxiv: v1 [math.co] 12 Sep 2014
arxv:1409.3707v1 [math.co] 12 Sep 2014 On the bnomal sums of Horadam sequence Nazmye Ylmaz and Necat Taskara Department of Mathematcs, Scence Faculty, Selcuk Unversty, 42075, Campus, Konya, Turkey March
More informationON A DETERMINATION OF THE INITIAL FUNCTIONS FROM THE OBSERVED VALUES OF THE BOUNDARY FUNCTIONS FOR THE SECOND-ORDER HYPERBOLIC EQUATION
Advanced Mathematcal Models & Applcatons Vol.3, No.3, 2018, pp.215-222 ON A DETERMINATION OF THE INITIAL FUNCTIONS FROM THE OBSERVED VALUES OF THE BOUNDARY FUNCTIONS FOR THE SECOND-ORDER HYPERBOLIC EUATION
More informationNUMERICAL DIFFERENTIATION
NUMERICAL DIFFERENTIATION 1 Introducton Dfferentaton s a method to compute the rate at whch a dependent output y changes wth respect to the change n the ndependent nput x. Ths rate of change s called the
More informationImplicit Integration Henyey Method
Implct Integraton Henyey Method In realstc stellar evoluton codes nstead of a drect ntegraton usng for example the Runge-Kutta method one employs an teratve mplct technque. Ths s because the structure
More informationHybrid FDTD/FETD Technique Using Parametric Quadratic Programming for Nonlinear Maxwell s Equations
Progress In Electromagnetcs Research M, Vol. 54, 113 123, 217 Hybrd FDTD/FETD Technque Usng Parametrc Quadratc Programmng for Nonlnear Maxwell s Equatons Hongxa L 1,BaoZhu 2 and Jefu Chen 3, * Abstract
More informationProfessor Terje Haukaas University of British Columbia, Vancouver The Q4 Element
Professor Terje Haukaas Unversty of Brtsh Columba, ancouver www.nrsk.ubc.ca The Q Element Ths document consders fnte elements that carry load only n ther plane. These elements are sometmes referred to
More informationA MODIFIED METHOD FOR SOLVING SYSTEM OF NONLINEAR EQUATIONS
Journal of Mathematcs and Statstcs 9 (1): 4-8, 1 ISSN 1549-644 1 Scence Publcatons do:1.844/jmssp.1.4.8 Publshed Onlne 9 (1) 1 (http://www.thescpub.com/jmss.toc) A MODIFIED METHOD FOR SOLVING SYSTEM OF
More informationn α j x j = 0 j=1 has a nontrivial solution. Here A is the n k matrix whose jth column is the vector for all t j=0
MODULE 2 Topcs: Lnear ndependence, bass and dmenson We have seen that f n a set of vectors one vector s a lnear combnaton of the remanng vectors n the set then the span of the set s unchanged f that vector
More informationHongyi Miao, College of Science, Nanjing Forestry University, Nanjing ,China. (Received 20 June 2013, accepted 11 March 2014) I)ϕ (k)
ISSN 1749-3889 (prnt), 1749-3897 (onlne) Internatonal Journal of Nonlnear Scence Vol.17(2014) No.2,pp.188-192 Modfed Block Jacob-Davdson Method for Solvng Large Sparse Egenproblems Hongy Mao, College of
More informationThe Finite Element Method
The Fnte Element Method GENERAL INTRODUCTION Read: Chapters 1 and 2 CONTENTS Engneerng and analyss Smulaton of a physcal process Examples mathematcal model development Approxmate solutons and methods of
More information) is the unite step-function, which signifies that the second term of the right-hand side of the
Casmr nteracton of excted meda n electromagnetc felds Yury Sherkunov Introducton The long-range electrc dpole nteracton between an excted atom and a ground-state atom s consdered n ref. [1,] wth the help
More informationSTATIC ANALYSIS OF TWO-LAYERED PIEZOELECTRIC BEAMS WITH IMPERFECT SHEAR CONNECTION
STATIC ANALYSIS OF TWO-LERED PIEZOELECTRIC BEAMS WITH IMPERFECT SHEAR CONNECTION Ákos József Lengyel István Ecsed Assstant Lecturer Emertus Professor Insttute of Appled Mechancs Unversty of Mskolc Mskolc-Egyetemváros
More informationChapter 5. Solution of System of Linear Equations. Module No. 6. Solution of Inconsistent and Ill Conditioned Systems
Numercal Analyss by Dr. Anta Pal Assstant Professor Department of Mathematcs Natonal Insttute of Technology Durgapur Durgapur-713209 emal: anta.bue@gmal.com 1 . Chapter 5 Soluton of System of Lnear Equatons
More informationPhysics 5153 Classical Mechanics. D Alembert s Principle and The Lagrangian-1
P. Guterrez Physcs 5153 Classcal Mechancs D Alembert s Prncple and The Lagrangan 1 Introducton The prncple of vrtual work provdes a method of solvng problems of statc equlbrum wthout havng to consder the
More information16 Reflection and transmission, TE mode
16 Reflecton transmsson TE mode Last lecture we learned how to represent plane-tem waves propagatng n a drecton ˆ n terms of feld phasors such that η = Ẽ = E o e j r H = ˆ Ẽ η µ ɛ = ˆ = ω µɛ E o =0. Such
More informationChapter - 2. Distribution System Power Flow Analysis
Chapter - 2 Dstrbuton System Power Flow Analyss CHAPTER - 2 Radal Dstrbuton System Load Flow 2.1 Introducton Load flow s an mportant tool [66] for analyzng electrcal power system network performance. Load
More informationThe Geometry of Logit and Probit
The Geometry of Logt and Probt Ths short note s meant as a supplement to Chapters and 3 of Spatal Models of Parlamentary Votng and the notaton and reference to fgures n the text below s to those two chapters.
More informationErratum: A Generalized Path Integral Control Approach to Reinforcement Learning
Journal of Machne Learnng Research 00-9 Submtted /0; Publshed 7/ Erratum: A Generalzed Path Integral Control Approach to Renforcement Learnng Evangelos ATheodorou Jonas Buchl Stefan Schaal Department of
More informationCanonical transformations
Canoncal transformatons November 23, 2014 Recall that we have defned a symplectc transformaton to be any lnear transformaton M A B leavng the symplectc form nvarant, Ω AB M A CM B DΩ CD Coordnate transformatons,
More informationSupplemental document
Electronc Supplementary Materal (ESI) for Physcal Chemstry Chemcal Physcs. Ths journal s the Owner Socetes 01 Supplemental document Behnam Nkoobakht School of Chemstry, The Unversty of Sydney, Sydney,
More informationCHAPTER II THEORETICAL BACKGROUND
3 CHAPTER II THEORETICAL BACKGROUND.1. Lght Propagaton nsde the Photonc Crystal The frst person that studes the one dmenson photonc crystal s Lord Raylegh n 1887. He showed that the lght propagaton depend
More informationON MECHANICS WITH VARIABLE NONCOMMUTATIVITY
ON MECHANICS WITH VARIABLE NONCOMMUTATIVITY CIPRIAN ACATRINEI Natonal Insttute of Nuclear Physcs and Engneerng P.O. Box MG-6, 07725-Bucharest, Romana E-mal: acatrne@theory.npne.ro. Receved March 6, 2008
More informationFinite Element Modelling of truss/cable structures
Pet Schreurs Endhoven Unversty of echnology Department of Mechancal Engneerng Materals echnology November 3, 214 Fnte Element Modellng of truss/cable structures 1 Fnte Element Analyss of prestressed structures
More informationAn efficient algorithm for multivariate Maclaurin Newton transformation
Annales UMCS Informatca AI VIII, 2 2008) 5 14 DOI: 10.2478/v10065-008-0020-6 An effcent algorthm for multvarate Maclaurn Newton transformaton Joanna Kapusta Insttute of Mathematcs and Computer Scence,
More informationModule 3: Element Properties Lecture 1: Natural Coordinates
Module 3: Element Propertes Lecture : Natural Coordnates Natural coordnate system s bascally a local coordnate system whch allows the specfcaton of a pont wthn the element by a set of dmensonless numbers
More informationChapter 12. Ordinary Differential Equation Boundary Value (BV) Problems
Chapter. Ordnar Dfferental Equaton Boundar Value (BV) Problems In ths chapter we wll learn how to solve ODE boundar value problem. BV ODE s usuall gven wth x beng the ndependent space varable. p( x) q(
More informationA PROCEDURE FOR SIMULATING THE NONLINEAR CONDUCTION HEAT TRANSFER IN A BODY WITH TEMPERATURE DEPENDENT THERMAL CONDUCTIVITY.
Proceedngs of the th Brazlan Congress of Thermal Scences and Engneerng -- ENCIT 006 Braz. Soc. of Mechancal Scences and Engneerng -- ABCM, Curtba, Brazl,- Dec. 5-8, 006 A PROCEDURE FOR SIMULATING THE NONLINEAR
More informationLecture 12: Discrete Laplacian
Lecture 12: Dscrete Laplacan Scrbe: Tanye Lu Our goal s to come up wth a dscrete verson of Laplacan operator for trangulated surfaces, so that we can use t n practce to solve related problems We are mostly
More informationStatistics Chapter 4
Statstcs Chapter 4 "There are three knds of les: les, damned les, and statstcs." Benjamn Dsrael, 1895 (Brtsh statesman) Gaussan Dstrbuton, 4-1 If a measurement s repeated many tmes a statstcal treatment
More informationInterconnect Modeling
Interconnect Modelng Modelng of Interconnects Interconnect R, C and computaton Interconnect models umped RC model Dstrbuted crcut models Hgher-order waveform n dstrbuted RC trees Accuracy and fdelty Prepared
More informationUsing T.O.M to Estimate Parameter of distributions that have not Single Exponential Family
IOSR Journal of Mathematcs IOSR-JM) ISSN: 2278-5728. Volume 3, Issue 3 Sep-Oct. 202), PP 44-48 www.osrjournals.org Usng T.O.M to Estmate Parameter of dstrbutons that have not Sngle Exponental Famly Jubran
More information8.592J: Solutions for Assignment 7 Spring 2005
8.59J: Solutons for Assgnment 7 Sprng 5 Problem 1 (a) A flament of length l can be created by addton of a monomer to one of length l 1 (at rate a) or removal of a monomer from a flament of length l + 1
More informationImplementation of the Matrix Method
Computatonal Photoncs, Summer Term 01, Abbe School of Photoncs, FSU Jena, Prof. Thomas Pertsch Computatonal Photoncs Semnar 03, 7 May 01 Implementaton of the Matr Method calculaton of the transfer matr
More informationMathematical Preparations
1 Introducton Mathematcal Preparatons The theory of relatvty was developed to explan experments whch studed the propagaton of electromagnetc radaton n movng coordnate systems. Wthn expermental error the
More informationGeorgia Tech PHYS 6124 Mathematical Methods of Physics I
Georga Tech PHYS 624 Mathematcal Methods of Physcs I Instructor: Predrag Cvtanovć Fall semester 202 Homework Set #7 due October 30 202 == show all your work for maxmum credt == put labels ttle legends
More informationCIVL 8/7117 Chapter 10 - Isoparametric Formulation 42/56
CIVL 8/77 Chapter 0 - Isoparametrc Formulaton 4/56 Newton-Cotes Example Usng the Newton-Cotes method wth = ntervals (n = 3 samplng ponts), evaluate the ntegrals: x x cos dx 3 x x dx 3 x x dx 4.3333333
More informationAPPENDIX A Some Linear Algebra
APPENDIX A Some Lnear Algebra The collecton of m, n matrces A.1 Matrces a 1,1,..., a 1,n A = a m,1,..., a m,n wth real elements a,j s denoted by R m,n. If n = 1 then A s called a column vector. Smlarly,
More informationProcedia Computer Science
Avalable onlne at www.scencedrect.com Proceda Proceda Computer Computer Scence Scence 1 (01) 00 (009) 589 597 000 000 Proceda Computer Scence www.elsever.com/locate/proceda Internatonal Conference on Computatonal
More informationThe equation of motion of a dynamical system is given by a set of differential equations. That is (1)
Dynamcal Systems Many engneerng and natural systems are dynamcal systems. For example a pendulum s a dynamcal system. State l The state of the dynamcal system specfes t condtons. For a pendulum n the absence
More informationTHE SMOOTH INDENTATION OF A CYLINDRICAL INDENTOR AND ANGLE-PLY LAMINATES
THE SMOOTH INDENTATION OF A CYLINDRICAL INDENTOR AND ANGLE-PLY LAMINATES W. C. Lao Department of Cvl Engneerng, Feng Cha Unverst 00 Wen Hwa Rd, Tachung, Tawan SUMMARY: The ndentaton etween clndrcal ndentor
More informationProblem 1: To prove that under the assumptions at hand, the group velocity of an EM wave is less than c, I am going to show that
PHY 387 K. Solutons for problem set #7. Problem 1: To prove that under the assumptons at hand, the group velocty of an EM wave s less than c, I am gong to show that (a) v group < v phase, and (b) v group
More informationIntroduction to Antennas & Arrays
Introducton to Antennas & Arrays Antenna transton regon (structure) between guded eaves (.e. coaxal cable) and free space waves. On transmsson, antenna accepts energy from TL and radates t nto space. J.D.
More informationBeyond Zudilin s Conjectured q-analog of Schmidt s problem
Beyond Zudln s Conectured q-analog of Schmdt s problem Thotsaporn Ae Thanatpanonda thotsaporn@gmalcom Mathematcs Subect Classfcaton: 11B65 33B99 Abstract Usng the methodology of (rgorous expermental mathematcs
More informationUncertainty in measurements of power and energy on power networks
Uncertanty n measurements of power and energy on power networks E. Manov, N. Kolev Department of Measurement and Instrumentaton, Techncal Unversty Sofa, bul. Klment Ohrdsk No8, bl., 000 Sofa, Bulgara Tel./fax:
More informationLecture 10 Support Vector Machines II
Lecture 10 Support Vector Machnes II 22 February 2016 Taylor B. Arnold Yale Statstcs STAT 365/665 1/28 Notes: Problem 3 s posted and due ths upcomng Frday There was an early bug n the fake-test data; fxed
More informationComputational Astrophysics
Computatonal Astrophyscs Solvng for Gravty Alexander Knebe, Unversdad Autonoma de Madrd Computatonal Astrophyscs Solvng for Gravty the equatons full set of equatons collsonless matter (e.g. dark matter
More informationNumerical Simulation of One-Dimensional Wave Equation by Non-Polynomial Quintic Spline
IOSR Journal of Matematcs (IOSR-JM) e-issn: 78-578, p-issn: 319-765X. Volume 14, Issue 6 Ver. I (Nov - Dec 018), PP 6-30 www.osrournals.org Numercal Smulaton of One-Dmensonal Wave Equaton by Non-Polynomal
More informationANALYSIS OF ELECTROMAGNETIC FIELD USING THE CONSTRAINED INTERPOLATION PROFILE METHOD PHÂN TÍCH TRƯỜNG ĐIỆN TỪ SỬ DỤNG PHƯƠNG PHÁP CIP
ANALYSIS OF ELECTROMAGNETIC FIELD USING THE CONSTRAINED INTERPOLATION PROFILE METHOD PHÂN TÍCH TRƯỜNG ĐIỆN TỪ SỬ DỤNG PHƯƠNG PHÁP CIP LÊ VŨ HƯNG Cao đẳng kỹ thuật quốc ga Kushro, Nhật Bản Kushro Natonal
More informationSIO 224. m(r) =(ρ(r),k s (r),µ(r))
SIO 224 1. A bref look at resoluton analyss Here s some background for the Masters and Gubbns resoluton paper. Global Earth models are usually found teratvely by assumng a startng model and fndng small
More informationTensor Smooth Length for SPH Modelling of High Speed Impact
Tensor Smooth Length for SPH Modellng of Hgh Speed Impact Roman Cherepanov and Alexander Gerasmov Insttute of Appled mathematcs and mechancs, Tomsk State Unversty 634050, Lenna av. 36, Tomsk, Russa RCherepanov82@gmal.com,Ger@npmm.tsu.ru
More informationSTUDY ON TWO PHASE FLOW IN MICRO CHANNEL BASED ON EXPERI- MENTS AND NUMERICAL EXAMINATIONS
Blucher Mechancal Engneerng Proceedngs May 0, vol., num. www.proceedngs.blucher.com.br/evento/0wccm STUDY ON TWO PHASE FLOW IN MICRO CHANNEL BASED ON EXPERI- MENTS AND NUMERICAL EXAMINATIONS Takahko Kurahash,
More informationAdvanced Circuits Topics - Part 1 by Dr. Colton (Fall 2017)
Advanced rcuts Topcs - Part by Dr. olton (Fall 07) Part : Some thngs you should already know from Physcs 0 and 45 These are all thngs that you should have learned n Physcs 0 and/or 45. Ths secton s organzed
More informationColor Rendering Uncertainty
Australan Journal of Basc and Appled Scences 4(10): 4601-4608 010 ISSN 1991-8178 Color Renderng Uncertanty 1 A.el Bally M.M. El-Ganany 3 A. Al-amel 1 Physcs Department Photometry department- NIS Abstract:
More informationWeek3, Chapter 4. Position and Displacement. Motion in Two Dimensions. Instantaneous Velocity. Average Velocity
Week3, Chapter 4 Moton n Two Dmensons Lecture Quz A partcle confned to moton along the x axs moves wth constant acceleraton from x =.0 m to x = 8.0 m durng a 1-s tme nterval. The velocty of the partcle
More informationTHE EFFECT OF TORSIONAL RIGIDITY BETWEEN ELEMENTS ON FREE VIBRATIONS OF A TELESCOPIC HYDRAULIC CYLINDER SUBJECTED TO EULER S LOAD
Journal of Appled Mathematcs and Computatonal Mechancs 7, 6(3), 7- www.amcm.pcz.pl p-issn 99-9965 DOI:.75/jamcm.7.3. e-issn 353-588 THE EFFECT OF TORSIONAL RIGIDITY BETWEEN ELEMENTS ON FREE VIBRATIONS
More informationSupporting Information
Supportng Informaton Water structure at the ar-aqueous nterface of dvalent caton and ntrate solutons Man Xu, Rck Spnney, Heather C. Allen* allen@chemstry.oho-state.edu Fresnel factors and spectra normalzaton
More informationConsistency & Convergence
/9/007 CHE 374 Computatonal Methods n Engneerng Ordnary Dfferental Equatons Consstency, Convergence, Stablty, Stffness and Adaptve and Implct Methods ODE s n MATLAB, etc Consstency & Convergence Consstency
More informationDifference Equations
Dfference Equatons c Jan Vrbk 1 Bascs Suppose a sequence of numbers, say a 0,a 1,a,a 3,... s defned by a certan general relatonshp between, say, three consecutve values of the sequence, e.g. a + +3a +1
More informationCOMPARISON OF SOME RELIABILITY CHARACTERISTICS BETWEEN REDUNDANT SYSTEMS REQUIRING SUPPORTING UNITS FOR THEIR OPERATIONS
Avalable onlne at http://sck.org J. Math. Comput. Sc. 3 (3), No., 6-3 ISSN: 97-537 COMPARISON OF SOME RELIABILITY CHARACTERISTICS BETWEEN REDUNDANT SYSTEMS REQUIRING SUPPORTING UNITS FOR THEIR OPERATIONS
More informationNew Method for Solving Poisson Equation. on Irregular Domains
Appled Mathematcal Scences Vol. 6 01 no. 8 369 380 New Method for Solvng Posson Equaton on Irregular Domans J. Izadan and N. Karamooz Department of Mathematcs Facult of Scences Mashhad BranchIslamc Azad
More informationBézier curves. Michael S. Floater. September 10, These notes provide an introduction to Bézier curves. i=0
Bézer curves Mchael S. Floater September 1, 215 These notes provde an ntroducton to Bézer curves. 1 Bernsten polynomals Recall that a real polynomal of a real varable x R, wth degree n, s a functon of
More information