Combined RKDG and LDG Method for the Simulation of the Bipolar Charge Transport in Solid Dielectrics
|
|
- Horace Quinn
- 6 years ago
- Views:
Transcription
1 PIERS ONLINE, VOL. 5, NO. 1, Combined RKDG nd LDG Method for the Simution of the Bipor Chrge Trnsport in Soid Dieectrics Jihun Tin, Jun Zou, nd Jinsheng Yun Deprtment of Eectric Engineering, Tsinghu University, Beijing , Chin Abstrct The spce chrge trnsport described by the convection-rection equtions is simuted by using the Runge-Kutt discontinuous Gerkin RKDG method. Compred to the trdition finite voume method, the RKDG method hs higher ccurcy nd better performnce of cpturing the steep front in the chrge profie. Combined with the RKDG method, the oc discontinuous Gerkin LDG method ws dopted to sove the Poisson s eqution insted of the previous boundry eement method BEM, in order tht the numeric qudrture for the soution of the convection-rection equtions is repced by nytic formus. The bipor chrge trnsport under dc votge is simuted by the RKDG+LDG method, which is more efficient nd produces identic resuts with those of the RKDG+BEM method. 1. INTRODUCTION The spce chrge dynmics pys n importnt roe on the degrdtion nd the brekdown of soid dieectrics under the high votge due to the intensifiction of the oc eectric fied. The trnsport of the different species of spce chrges under the eectric fied cn be described by set of convection-rection equtions couped with the Poisson s eqution. In our previous work, the RKDG method ws utiized to sove the convection-rection equtions [1]. Unike the trdition QUICKEST method with the ULTIMATE fux imiter [, 3], the RKDG method cn hnde the convection nd the rection terms simutneousy without incorporting the spit procedure [4], nd the high-order ccurcy cn be obtined without incorporting wide stenci, which is more suitbe for the shock cpturing thn the QUICKEST method [5]. For the soution of the Poisson s eqution, the oc discontinuous Gerkin LDG method proposed by B. Cockburn et. [5, 6] is dopted in this pper insted of the genery used boundry eement method BEM [7], in order to repce the numeric qudrture with nytic formus in the RKDG wek formution of the convection-rection equtions. This combined RKDG+LDG method is verified by the simution of the spce chrge trnsport in ow-density poyethyene LDPE smpe under high dc votge. The resuts show tht the proposed method is more efficient nd produces identic spce chrge profies with those of the RKDG+BEM method.. MATHEMATICAL MODEL DESCRIPTIONS The fuid mode for the spce chrge trnsport in this pper ws initiy proposed by Aison nd Hi [8]. It hs redy been widey dopted by Le Roy nd Begroui et. in [, 3, 7]. The mode is given s foows: n t + f x = S convection-rection eqution, 1 E x = ρ ε 0 ε r Poisson s eqution, f = ±µ n E trnsport eqution, 3 where n is the number of crriers per unit voume nd the subscript represents different crrier types. Usuy, four species of crriers re tken into ccount to describe the chrge dynmics in soid dieectrics, i.e., free eectrons nd hoes, which represent the crriers in the show trps; trpped eectrons nd hoes, which represent the immobie crriers residing in the deep trps. The source term S in 1 comprises the effects of trpping nd recombintion. Both the source term S nd the grdient of the crrier convection fux f contribute to the vrition of the crrier concentrtion ccording to 1. At the eectrodes, the inet boundry condition, i.e., crrier injection, for the convection-rection eqution is usuy prescribed s the Schottky injection w, whie the Ohm s w is dopted s the outet boundry condition, i.e., crrier extrction. For the definitions nd the physic menings of the mode prmeters, reders my refer to [1, 3, 7, 8] for detis.
2 PIERS ONLINE, VOL. 5, NO. 1, RKDG METHOD FOR CONVECTION-REACTION EQUATION The uthors of [3, 7] used the QUICKEST method with the ULTIMATE fux imiter to sove the convection-rection Eqution 1. However, it cnnot ccurtey simute the spce chrge profie with very steep wvefront, i.e., shock wve, due to tht wide stenci is dopted to interpote the numeric fux. Moreover, the spit tretment of the originy couped convection nd rection terms brings ddition errors [1]. In our previous work [1], the RKDG method ws ppied to the convection-rection eqution, in which the high-order piecewise continuous Legendre poynomis were chosen s the bsis functions to expnd the spce chrge profie. A speci Runge-Kutt method combined with sope imiter ws used for the time discretiztion, which ensures the numeric stbiity. The source term S nd the convection term f / x were hnded simutneousy by the RKDG method without the spit procedure. High-order ccurcy cn be obtined without incorporting wide stenci, which is more suitbe for the shock cpturing thn the QUICKEST method. The bsic steps of the RKDG method re s foows: 1. Discontinuous Gerkin DG spce discretiztion: The crrier concentrtion n of the jth eement = x j 1, x j+ 1 cn be pproximted by the Legendre poynomis {P x} =0,...,k s: n = k =0 P x xj, 4 where is the order of the poynomi nd is the ength of. By tking {P x} =0,...,k s the weighting functions, the convection-rection Eqution 1 cn be weky enforced. After the simpifiction by utiizing the orthogonity of the Legendre poynomis, the semi-discrete form of the convection-rection eqution cn be obtined: d dt α t = + 1 where ˆf is the numeric fux. P x xj f P x xj ˆf x j+1/ x j 1/ + S P x xj, 5. Numeric fux: The Lx-Friedrichs numeric fux ws dopted in the RKDG method ccording to [9]. For the right boundry of the jth eement, it is given s: where ˆf j+1/ n j+1/, n + j+1/ = 1 [ ] f j n j+1/ + f j n + j+1/ Cn + j+1/ n j+1/, 6 C = mx v t = 0, x, 7 x n j+ nd n + 1 j+ re the eft nd right imit vues of n 1 t x j+ 1, v t = 0 is the initi crrier speed. 3. RK time discretiztion: According to [9], the RK time discretiztion shoud be performed in k + 1 intermedite steps to mtch the ccurcy degree of the spce domin. With the chrge distribution n m t the time point m t, the RK time stepping is s foows: Set n 0 = n m ; b For the s intermedite steps s = 1,..., k + 1: n s = ΛΠ s 1 α s ω s, 8 =0
3 PIERS ONLINE, VOL. 5, NO. 1, where ω s = n + β s α s tl j {n }, 9 L j {n } represents the right-hnd side of 5. ΛΠ is sope imiter proposed by Bisws et., which ensures the numeric stbiity of the RKDG method [10]. The determintion of the coefficients α s nd β s is introduced by Gottieb et. in [11]. c Set n m+1 = n k LDG METHOD COMBINED WITH RKDG METHOD The purpose of dopting the LDG method for the soution of the Poisson s eqution is to promote the efficiency of the RKDG method by nyticy evuting the integrs insted of the numeric qudrture. In the RKDG semi-discrete formution of the convection-rection Eqution 5, there exists the integrnd which is the product of the crrier concentrtion n nd the eectric fied E: f P x xj = x xj. 10 If the eectric fied E is soved by the BEM, this integr shoud be evuted by numeric methods, such s the Guss-Legendre qudrture. However, for the combined RKDG+LDG method, the Legendre poynomis re seected s the bsis functions for both the crrier concentrtion n nd the eectric fied E. This fct enbes us to nyticy ccute the bove integr LDG Method for Poisson s Eqution For the one-dimension Poisson s eqution defined in domin Ω with the first kind condition on the domin boundry Γ D : dq = f in Ω, 11 q = E = du, u = g D on Γ D, its wek soution stisfies the foowing two equtions: x j+ 1 qr = rû u dr, 1 x j 1 q dv = x j+ 1 fv + vˆq, 13 x j 1 where r nd v beongs to the spce of Legendre poynomis defined on with the degree t most k, i.e., r, v P k = {P x : = 0,..., k}. The numeric fux û in 1 nd ˆq in 13 re defined s: ˆq {q} C11 C = + 1 [u] û {u} C 1 C [q] ˆq q = û + C 11 u + g D n on Γ g D, n = D in Ω, 14 { 1 eft boundry 1 right boundry, 15 where C 11 = C = 1, C 1 = C 1 = 0, {q} = 1 q+ + q, [q] = q + q. {u} nd [u] dopt j+ 1 j+ 1 j+ 1 j+ 1 the simir definitions s {q} nd [q]. By expnding u nd its grdient q with Legendre poynomis, the coefficients of the bsis functions cn be obtined by soving 1 nd 13.
4 PIERS ONLINE, VOL. 5, NO. 1, Anytic Evution of the Integr In the RKDG+LDG method, both the chrge density n nd the eectric fied E re expnded by the Legendre poynomis with the form s in 4. Therefore, the integr in 10 cn be directy given by nytic formus. Tke the Legendre poynomi spce P k where k = 1 nd k = s exmpe, the nytic soution of the integr cn be given s: x xj = { 0 = 0 ±µ 0 b b 1 = 1 x xj 0 =0 = ±µ 0 b b b =1 ±µ 0 b b b b 1 = for k = 1, 16 for k =, 17 where {b } 0 represent the coefficients of the Legendre poynomis for E. The bove nytic formus re fexibe nd cn be ppied to vrious situtions, due to tht they re derived from the integrtion with respect to the product of the Legendre poynomi series which re used to pproximte the rbitrry distribution of the eectric fied nd the crrier density. 5. RESULTS AND COMPARISONS A 118 µm LDPE smpe under extern dc votge 1 kv is used s test mode for the bipor chrge trnsport. The dopted prmeter vues re the sme s those in [3]. To verify the combined RKDG+LDG method with different bsis function orders k = 1 nd k =, its simution resuts t t = 0, 30, 50, 100 s re compred with those obtined from the RKDG+BEM method. In the impementtion of the RKDG+BEM method, five-point Guss-Legendre numeric qudrture ws dopted in ech spci eement in. The retive errors of the eectric fied nd the net chrge density re given in Tbe 1. The retive errors re defined s foows: ε E = E BEM E LDG E BEM, ε ρ = ρ /BEM ρ /LDG ρ /BEM 18 It cn be seen from these dt tht for the cse k = 1, the RKDG+LDG method produces most identic resuts s those of the RKDG+BEM method. For k =, the retive error of the eectric fied is sti very sm. Athough the retive error of the net chrge density increses to bout 1%, the simuted chrge profies of the two methods re sti very cose to ech other see Figure 1. Moreover, by using the RKDG+LDG method, of the numeric qudrtures impemented in the RKDG+BEM method re repced by nytic soutions, which sves 40% of the computtion time. From our numeric tests, the difference in the computtion time between the BEM nd the LDG method for soving the Poisson s eqution is negigibe nd the incresed efficiency of the RKDG+LDG method is miny due to repcing the numeric qudrtures by nytic formus in the soution of the convection-rection eqution. Therefore, the combined RKDG+LDG method is verified nd shows to be more efficient thn the RKDG+BEM method. 8 Net Chrge Density C/m s 50s 30s 0s RKDG+LDG k= RKDG+BEM k= Xm x 10-4 Figure 1: Comprison of the simuted spce chrge profies by using the rd order RKDG+LDG nd RKDG+BEM methods.
5 PIERS ONLINE, VOL. 5, NO. 1, Tbe 1: Retive errors of the eectric fied nd the net chrge density by using the RKDG+BEM nd the RKDG+LDG methods. Timess ε E k = 1 ε E k = ε ρ k = 1 ε ρ k = e e e e e e e e e e e e CONCLUSIONS The combined RKDG nd LDG method proposed in this pper hs high-order ccurcy nd is suitbe for simuting the spce chrge trnsport with steep fronts. Both the RKDG nd the LDG methods re bsed upon the discontinuous Gerkin spce discretiztion with the Legendre poynomis s the bsis functions, which enbes us to repce the numeric qudrtures in the RKDG+BEM method with nytic formus. The RKDG+LDG method is verified by the simution of the bipor chrge trnsport in n LDPE smpe, which produces identic resuts with those of the RKDG+BEM method. Due to the empoyment of the nytic integr, the RKDG+LDG method cn promote the efficiency of the numeric scheme by sving 40% of the computtion time tht needed by the RKDG+BEM method. ACKNOWLEDGMENT This work is supported by the Ntion Ntur Science Foundtion of Chin REFERENCES 1. Tin, J., J. Zou, Y. Wng, J. Liu, J. Yun, nd Y. Zhou, Simution of bipor chrge trnsport with trpping nd recombintion in poymeric insutors using Runge-Kutt discontinuous gerkin method, Journ of Physics D: Appied Physics, Vo. 41, No. 19, 10, October Le Roy, S., Numeric methods in the simution of chrge trnsport in soid dieectrics, IEEE Trnsctions on Dieectrics nd Eectric Insution, Vo. 13, No., 39 46, Begroui, E., I. Boukhris, A. Ke, G. Teyssedre, nd C. Lurent, A new numeric mode ppied to bipor chrge trnsport, trpping nd recombintion under ow nd high dc votges, Journ of Physics D: Appied Physics, Vo. 40, No. 1, , Toro, E. F., Riemnn Sovers nd Numeric Methods for Fuid Dynmics, Chpter 15, Spitting Schemes for PDEs with Source Terms, Springer, Cockburn, B. nd C. W. Shu, Runge-kutt discontinuous gerkin methods for convectiondominted probems, J. Sci. Comput, Vo. 16, No. 3, , Cockburn, B. nd C. W. Shu, The oc discontinuous gerkin method for time-dependent convection-diffusion systems, SIAM Journ on Numeric Anysis, Vo. 35, No. 6, , Le Roy, S., P. Segur, G. Teyssedre, nd C. Lurent, Description of bipor chrge trnsport in poyethyene using fuid mode with constnt mobiity: Mode prediction, Journ of Physics D: Appied Physics, Vo. 37, No., , Aison, J. M. nd R. M. Hi, A mode for bipor chrge trnsport, trpping nd recombintion in degssed crossinked poyethene, Journ of Physics D: Appied Physics, Vo. 7, , Cockburn, B. nd C.-W. Shu, TVB Runge-Kutt oc projection discontinuous Gerkin finite eement method for conservtion ws II: Gener frmework, Mthemtics of Computtion, Vo. 5, No. 186, , Bisws, R., K. D. Devine, nd J. E. Fherty, Pre, dptive finite eement methods for conservtion ws, App. Numer. Mth., Vo. 14, No. 1 3, 55 83, Gottieb, S. nd C.-W. Shu, Tot Vrition diminishing Runge-Kutt schemes, Mthemtics of Computtion, Vo. 67, No. 1, 73 85, 1998.
Development of the Sinc Method for Nonlinear Integro-Differential Eequations
Austrin Journ of Bsic nd Appied Sciences, 4(): 558-555, ISS 99-878 Deveopment of the Sinc Method for oniner Integro-Differenti Eequtions K. Jei, M. Zrebni, 3 M. Mirzee Chi,3 Ismic Azd University Brnch
More information1. The vibrating string problem revisited.
Weeks 7 8: S eprtion of Vribes In the pst few weeks we hve expored the possibiity of soving first nd second order PDEs by trnsforming them into simper forms ( method of chrcteristics. Unfortuntey, this
More informationMath 124B January 24, 2012
Mth 24B Jnury 24, 22 Viktor Grigoryn 5 Convergence of Fourier series Strting from the method of seprtion of vribes for the homogeneous Dirichet nd Neumnn boundry vue probems, we studied the eigenvue probem
More informationIn this appendix, we evaluate the derivative of Eq. 9 in the main text, i.e., we need to calculate
Supporting Tet Evoution of the Averge Synptic Updte Rue In this ppendi e evute the derivtive of Eq. 9 in the min tet i.e. e need to ccute Py ( ) Py ( Y ) og γ og. [] P( y Y ) P% ( y Y ) Before e strt et
More informationMAGIC058 & MATH64062: Partial Differential Equations 1
MAGIC58 & MATH646: Prti Differenti Equtions 1 Section 4 Fourier series 4.1 Preiminry definitions Definition: Periodic function A function f( is sid to be periodic, with period p if, for, f( + p = f( where
More informationFluid Flow through a Tube
. Theory through Tube In this experiment we wi determine how we physic retionship (so ced w ), nmey Poiseue s eqution, ppies. In the suppementry reding mteri this eqution ws derived s p Q 8 where Q is
More informationIntroduction to statically indeterminate structures
Sttics of Buiding Structures I., EASUS Introduction to stticy indeterminte structures Deprtment of Structur echnics Fcuty of Civi Engineering, VŠB-Technic University of Ostrv Outine of Lecture Stticy indeterminte
More informationStatistical Physics. Solutions Sheet 5.
Sttistic Physics. Soutions Sheet 5. Exercise. HS 04 Prof. Mnfred Sigrist Ide fermionic quntum gs in hrmonic trp In this exercise we study the fermionic spiness ide gs confined in three-dimension hrmonic
More informationStability condition of finite difference solution for viscoelastic wave equations
Erthq Sci (009): 79 85 79 Doi: 0007/s589-009-079- Stbiity condition of finite difference soution for viscoestic wve equtions Chengyu Sun, Yunfei Xio Xingyo Yin nd Hongcho Peng Chin University of Petroeum,
More informationDynamic Analysis of the Turnout Diverging Track for HSR with Variable Curvature Sections
Word Journ of Engineering nd Technoogy, 07, 5, 4-57 http://www.scirp.org/journ/wjet ISSN Onine: -449 ISSN Print: -4 Dynmic Anysis of the Turnout Diverging Trck for HSR with Vribe Curvture Sections Wdysw
More informationCIT and QIT with a Periodic Impulsional Potential
Adv. Studies Theor. Phys., Vol. 6, 2012, no. 6, 285-291 nd with Periodic Impulsionl Potentil S. Seddighi Chhrborj srkhosh@mth.upm.edu.my Islmic Azd University, Science nd Reserch Brnch Deprtment of Mthemtics,
More information3.4 Numerical integration
3.4. Numericl integrtion 63 3.4 Numericl integrtion In mny economic pplictions it is necessry to compute the definite integrl of relvlued function f with respect to "weight" function w over n intervl [,
More informationWhere does Oxygen Extinction Occur in a Soil Profile?
2st Interntion Congress on Modeing nd Simution, God Cost, Austri, 29 Nov to 4 ec 25 www.mssnz.org.u/modsim25 Where does Oxygen Extinction Occur in Soi Profie? Freemn J Cook nd John H Knightb Freemn Cook
More informationarxiv: v1 [math.co] 5 Jun 2015
First non-trivi upper bound on the circur chromtic number of the pne. Konstnty Junosz-Szniwski, Fcuty of Mthemtics nd Informtion Science, Wrsw University of Technoogy, Pond Abstrct rxiv:1506.01886v1 [mth.co]
More informationComplete Description of the Thelen2003Muscle Model
Compete Description o the he23usce ode Chnd John One o the stndrd musce modes used in OpenSim is the he23usce ctutor Unortuntey, to my knowedge, no other pper or document, incuding the he, 23 pper describing
More informationTheoretical foundations of Gaussian quadrature
Theoreticl foundtions of Gussin qudrture 1 Inner product vector spce Definition 1. A vector spce (or liner spce) is set V = {u, v, w,...} in which the following two opertions re defined: (A) Addition of
More informationANALYZING EFFECT OF TEMPERATURE VARIATION ON SKIN SUB-LAYERS USING FINITE ELEMENT APPROACH
Origin Reserch Artice Aied sciences Interntion Journ of Phrm nd Bio Sciences ISSN 0975-699 ANALYZING EFFECT OF TEMPERATURE VARIATION ON SKIN SUB-LAYERS USING FINITE ELEMENT APPROACH MONA DWIVEDI 1 * AND
More informationTravelling Profile Solutions For Nonlinear Degenerate Parabolic Equation And Contour Enhancement In Image Processing
Applied Mthemtics E-Notes 8(8) - c IN 67-5 Avilble free t mirror sites of http://www.mth.nthu.edu.tw/ men/ Trvelling Profile olutions For Nonliner Degenerte Prbolic Eqution And Contour Enhncement In Imge
More informationEnergy Balance of Solar Collector
Renewbe Energy Grou Gret Ides Grow Better Beow Zero! Wecome! Energy Bnce of Sor Coector Mohmd Khrseh E-mi:m.Khrseh@gmi.com Renewbe Energy Grou Gret Ides Grow Better Beow Zero! Liuid Ft Pte Coectors. Het
More informationMutipy by r sin RT P to get sin R r r R + T sin (sin T )+ P P = (7) ffi So we hve P P ffi = m (8) choose m re so tht P is sinusoi. If we put this in b
Topic 4: Lpce Eqution in Spheric Co-orintes n Mutipoe Expnsion Reing Assignment: Jckson Chpter 3.-3.5. Lpce Eqution in Spheric Coorintes Review of spheric por coorintes: x = r sin cos ffi y = r sin sin
More informationFRACTURE OF PIEZOELECTRIC MATERIALS
CF000OR FRACUR OF PZOCRC MARAS ong-yi Zhng* nd Mingho Zho eprtment of Mechnic ngineering Hong ong University of Science nd echnoogy Cer Wter By owoon Hong ong Chin * -mi: mezhngt@ust.hk ABSRAC he present
More informationIdentification of Axial Forces on Statically Indeterminate Pin-Jointed Trusses by a Nondestructive Mechanical Test
Send Orders of Reprints t reprints@benthmscience.net 50 The Open Civi Engineering Journ, 0, 7, 50-57 Open Access Identifiction of Axi Forces on Stticy Indeterminte Pin-Jointed Trusses by Nondestructive
More information4. Calculus of Variations
4. Clculus of Vritions Introduction - Typicl Problems The clculus of vritions generlises the theory of mxim nd minim. Exmple (): Shortest distnce between two points. On given surfce (e.g. plne), nd the
More informationAdvanced Calculus: MATH 410 Notes on Integrals and Integrability Professor David Levermore 17 October 2004
Advnced Clculus: MATH 410 Notes on Integrls nd Integrbility Professor Dvid Levermore 17 October 2004 1. Definite Integrls In this section we revisit the definite integrl tht you were introduced to when
More informationSuggested Solution to Assignment 5
MATH 4 (5-6) prti diferenti equtions Suggested Soution to Assignment 5 Exercise 5.. () (b) A m = A m = = ( )m+ mπ x sin mπx dx = x mπ cos mπx + + 4( )m 4 m π. 4x cos mπx dx mπ x cos mπxdx = x mπ sin mπx
More informationThermal Relaxation Times in Biological Tissues Subjected to Pulsed Laser Irradiation
9th AIAA/ASME Joint Thermophysics nd Het Trnsfer Conference 5-8 June 006, Sn Frncisco, Ciforni AIAA 006-938 Therm Rextion Times in Bioogic Tissues Subjected to Pused Lser Irrdition Kyunghn Kim * nd Zhixiong
More informationA Study on Prediction of Cavitation for Centrifugal Pump Myung Jin Kim, Hyun Bae Jin, and Wui Jun Chung
Word Acdemy of Science, Engineering nd Technoogy Interntion Journ of Mechnic nd Mechtronics Engineering Vo:6, No:12, 2012 A Study on Prediction of Cvittion for Centrifug Pump Myung Jin Kim, Hyun Be Jin,
More informationA Slipping and Buried Strike-Slip Fault in a Multi-Layered Elastic Model
Geosciences 7, 7(): 68-76 DOI:.59/j.geo.77. A Sipping nd Buried Strike-Sip Fut in Muti-Lyered Estic Mode Asish Krmkr,*, Snjy Sen Udirmpur Pisree Sikshytn (H.S.), Udirmpur, P.O. Knyngr, Pin, Indi Deprtment
More informationMath 5440 Problem Set 3 Solutions
Mth 544 Mth 544 Problem Set 3 Solutions Aron Fogelson Fll, 25 1: Logn, 1.5 # 2) Repet the derivtion for the eqution of motion of vibrting string when, in ddition, the verticl motion is retrded by dmping
More informationHomework Assignment #5 Solutions
Physics 506 Winter 008 Textbook probems: Ch. 9: 9., 9., 9.4 Ch. 10: 10.1 Homework Assignment #5 Soutions 9. A spheric hoe of rdius in conducting medium cn serve s n eectromgnetic resonnt cvity. Assuming
More informationP R A S I C ' 02. Simpozionul na.ional cu participare interna.ional/ PRoiectarea ASIstat/ de Calculator
UNIVERSITTE TRNSILVNI DIN R0OV tedr Design de Produs i Robotic Simpozionu n.ion cu prticipre intern.ion/ PRoiectre SIstt/ de cutor P R S I ' 0 Vo. I Mecnisme i Triboogie 7-8 Noiembrie rov Români ISN 97-65-064-9
More informationThe Regulated and Riemann Integrals
Chpter 1 The Regulted nd Riemnn Integrls 1.1 Introduction We will consider severl different pproches to defining the definite integrl f(x) dx of function f(x). These definitions will ll ssign the sme vlue
More information1 The Lagrange interpolation formula
Notes on Qudrture 1 The Lgrnge interpoltion formul We briefly recll the Lgrnge interpoltion formul. The strting point is collection of N + 1 rel points (x 0, y 0 ), (x 1, y 1 ),..., (x N, y N ), with x
More informationGEOMETRIC PROBABILITY MODELS TO ANALYZE STRATEGIES FOR FINDING A BURIED CABLE
Journ of the Opertions Reserch Society of Jpn Vo. 6, No. 3, Juy 17, pp. 4 417 c The Opertions Reserch Society of Jpn GEOMETRIC PROBABILITY MODELS TO ANALYZE STRATEGIES FOR FINDING A BURIED CABLE Ken-ichi
More informationJournal of Computational and Applied Mathematics. A perturbation solution for bacterial growth and bioremediation in a porous medium with bio-clogging
Journ of Computtion nd Appied Mthemtics 234 (2) 279 2723 Contents ists vibe t ScienceDirect Journ of Computtion nd Appied Mthemtics journ homepge: www.esevier.com/octe/cm A perturbtion soution for bcteri
More informationMath 5440 Problem Set 3 Solutions
Mth 544 Mth 544 Problem Set 3 Solutions Aron Fogelson Fll, 213 1: (Logn, 1.5 # 2) Repet the derivtion for the eqution of motion of vibrting string when, in ddition, the verticl motion is retrded by dmping
More informationLecture 6 Notes, Electromagnetic Theory I Dr. Christopher S. Baird University of Massachusetts Lowell
Lecture 6 Notes, Eectrognetic Theory I Dr. Christopher S. Bird University of Msschusetts Lowe. Associted Legendre Poynois - We now return to soving the Lpce eqution in spheric coordintes when there is
More informationSPECTRAL APPROXIMATION OF TIME-HARMONIC MAXWELL EQUATIONS IN THREE-DIMENSIONAL EXTERIOR DOMAINS. 1. Introduction
PECTRAL APPROXIMATION OF TIME-HARMONIC MAXWELL EQUATION IN THREE-DIMENIONAL EXTERIOR DOMAIN LINA MA 1, JIE HEN 1 AND LI-LIAN WANG 2 Abstrct. We deveop in this pper n efficient nd robust spectr-gerkin method
More informationChapter 0. What is the Lebesgue integral about?
Chpter 0. Wht is the Lebesgue integrl bout? The pln is to hve tutoril sheet ech week, most often on Fridy, (to be done during the clss) where you will try to get used to the ides introduced in the previous
More information8. Modeling of Rotary Inverted Pendulum
NCU Deprtment of Eectric n Computer Engineering 5 Spring Course by Prof. Yon-Ping Chen 8. oeing of otry Inverte Penuum z m (, y, z,) g y L Figure 8- he ynmic moe
More informationMatrices, Moments and Quadrature, cont d
Jim Lmbers MAT 285 Summer Session 2015-16 Lecture 2 Notes Mtrices, Moments nd Qudrture, cont d We hve described how Jcobi mtrices cn be used to compute nodes nd weights for Gussin qudrture rules for generl
More informationConservation Law. Chapter Goal. 5.2 Theory
Chpter 5 Conservtion Lw 5.1 Gol Our long term gol is to understnd how mny mthemticl models re derived. We study how certin quntity chnges with time in given region (sptil domin). We first derive the very
More informationA REVIEW OF CALCULUS CONCEPTS FOR JDEP 384H. Thomas Shores Department of Mathematics University of Nebraska Spring 2007
A REVIEW OF CALCULUS CONCEPTS FOR JDEP 384H Thoms Shores Deprtment of Mthemtics University of Nebrsk Spring 2007 Contents Rtes of Chnge nd Derivtives 1 Dierentils 4 Are nd Integrls 5 Multivrite Clculus
More informationBending Analysis of Castellated Beams. By Sahar Elaiwi Boksun Kim Long-Yuan Li
Athens Journ of Technoogy nd Engineering X Y Bending Anysis of Csteted Bems By Shr Eiwi Boksun Kim Long-Yun Li Existing studies hve shown tht the od-crrying cpcity of csteted bems cn be infuenced by the
More informationLecture 6: Singular Integrals, Open Quadrature rules, and Gauss Quadrature
Lecture notes on Vritionl nd Approximte Methods in Applied Mthemtics - A Peirce UBC Lecture 6: Singulr Integrls, Open Qudrture rules, nd Guss Qudrture (Compiled 6 August 7) In this lecture we discuss the
More informationKRASNOSEL SKII TYPE FIXED POINT THEOREM FOR NONLINEAR EXPANSION
Fixed Point Theory, 13(2012), No. 1, 285-291 http://www.mth.ubbcluj.ro/ nodecj/sfptcj.html KRASNOSEL SKII TYPE FIXED POINT THEOREM FOR NONLINEAR EXPANSION FULI WANG AND FENG WANG School of Mthemtics nd
More informationNEW INEQUALITIES OF SIMPSON S TYPE FOR s CONVEX FUNCTIONS WITH APPLICATIONS. := f (4) (x) <. The following inequality. 2 b a
NEW INEQUALITIES OF SIMPSON S TYPE FOR s CONVEX FUNCTIONS WITH APPLICATIONS MOHAMMAD ALOMARI A MASLINA DARUS A AND SEVER S DRAGOMIR B Abstrct In terms of the first derivtive some ineulities of Simpson
More informationMAC-solutions of the nonexistent solutions of mathematical physics
Proceedings of the 4th WSEAS Interntionl Conference on Finite Differences - Finite Elements - Finite Volumes - Boundry Elements MAC-solutions of the nonexistent solutions of mthemticl physics IGO NEYGEBAUE
More informationNumerical quadrature based on interpolating functions: A MATLAB implementation
SEMINAR REPORT Numericl qudrture bsed on interpolting functions: A MATLAB implementtion by Venkt Ayylsomyjul A seminr report submitted in prtil fulfillment for the degree of Mster of Science (M.Sc) in
More informationLecture 14: Quadrature
Lecture 14: Qudrture This lecture is concerned with the evlution of integrls fx)dx 1) over finite intervl [, b] The integrnd fx) is ssumed to be rel-vlues nd smooth The pproximtion of n integrl by numericl
More informationTransport Calculations. Tseelmaa Byambaakhuu, Dean Wang*, and Sicong Xiao
A Locl hp Adptive Diffusion Synthetic Accelertion Method for Neutron Trnsport Clcultions Tseelm Bymbkhuu, Den Wng*, nd Sicong Xio University of Msschusetts Lowell, 1 University Ave, Lowell, MA 01854 USA
More informationNUMERICAL INTEGRATION
NUMERICAL INTEGRATION How do we evlute I = f (x) dx By the fundmentl theorem of clculus, if F (x) is n ntiderivtive of f (x), then I = f (x) dx = F (x) b = F (b) F () However, in prctice most integrls
More informationA Fully-coupled Thermo-Hydro-Mechanical Model for the Description of the Behavior of Swelling Porous Media
Act Poytechnic Hungric Vo. 8, No. 4, 211 A Fuy-couped Thermo-Hydro-Mechnic Mode for the Description of the Behvior of Sweing Porous Medi Hnifi Missoum, Ndi Lredj, rim Bendni, Mustph Miki Construction,
More informationQUADRATURE is an old-fashioned word that refers to
World Acdemy of Science Engineering nd Technology Interntionl Journl of Mthemticl nd Computtionl Sciences Vol:5 No:7 011 A New Qudrture Rule Derived from Spline Interpoltion with Error Anlysis Hdi Tghvfrd
More informationMinimum Energy State of Plasmas with an Internal Transport Barrier
Minimum Energy Stte of Plsms with n Internl Trnsport Brrier T. Tmno ), I. Ktnum ), Y. Skmoto ) ) Formerly, Plsm Reserch Center, University of Tsukub, Tsukub, Ibrki, Jpn ) Plsm Reserch Center, University
More informationAPPROXIMATE INTEGRATION
APPROXIMATE INTEGRATION. Introduction We hve seen tht there re functions whose nti-derivtives cnnot be expressed in closed form. For these resons ny definite integrl involving these integrnds cnnot be
More information5.7 Improper Integrals
458 pplictions of definite integrls 5.7 Improper Integrls In Section 5.4, we computed the work required to lift pylod of mss m from the surfce of moon of mss nd rdius R to height H bove the surfce of the
More informationCMDA 4604: Intermediate Topics in Mathematical Modeling Lecture 19: Interpolation and Quadrature
CMDA 4604: Intermedite Topics in Mthemticl Modeling Lecture 19: Interpoltion nd Qudrture In this lecture we mke brief diversion into the res of interpoltion nd qudrture. Given function f C[, b], we sy
More informationCalculus 2: Integration. Differentiation. Integration
Clculus 2: Integrtion The reverse process to differentition is known s integrtion. Differentition f() f () Integrtion As it is the opposite of finding the derivtive, the function obtined b integrtion is
More informationWeek 10: Line Integrals
Week 10: Line Integrls Introduction In this finl week we return to prmetrised curves nd consider integrtion long such curves. We lredy sw this in Week 2 when we integrted long curve to find its length.
More informationThe Basic Functional 2 1
2 The Bsic Functionl 2 1 Chpter 2: THE BASIC FUNCTIONAL TABLE OF CONTENTS Pge 2.1 Introduction..................... 2 3 2.2 The First Vrition.................. 2 3 2.3 The Euler Eqution..................
More informationu t = k 2 u x 2 (1) a n sin nπx sin 2 L e k(nπ/l) t f(x) = sin nπx f(x) sin nπx dx (6) 2 L f(x 0 ) sin nπx 0 2 L sin nπx 0 nπx
Chpter 9: Green s functions for time-independent problems Introductory emples One-dimensionl het eqution Consider the one-dimensionl het eqution with boundry conditions nd initil condition We lredy know
More informationLecture 23: Interpolatory Quadrature
Lecture 3: Interpoltory Qudrture. Qudrture. The computtion of continuous lest squres pproximtions to f C[, b] required evlutions of the inner product f, φ j = fxφ jx dx, where φ j is polynomil bsis function
More informationNumerical integration
2 Numericl integrtion This is pge i Printer: Opque this 2. Introduction Numericl integrtion is problem tht is prt of mny problems in the economics nd econometrics literture. The orgniztion of this chpter
More informationLinear Systems COS 323
Liner Systems COS 33 Soving Liner Systems of Equtions Define iner system Singurities in iner systems Gussin Eimintion: A gener purpose method Nïve Guss Guss with pivoting Asymptotic nysis Tringur systems
More informationS. S. Dragomir. 1. Introduction. In [1], Guessab and Schmeisser have proved among others, the following companion of Ostrowski s inequality:
FACTA UNIVERSITATIS NIŠ) Ser Mth Inform 9 00) 6 SOME COMPANIONS OF OSTROWSKI S INEQUALITY FOR ABSOLUTELY CONTINUOUS FUNCTIONS AND APPLICATIONS S S Drgomir Dedicted to Prof G Mstroinni for his 65th birthdy
More informationMath 115 ( ) Yum-Tong Siu 1. Lagrange Multipliers and Variational Problems with Constraints. F (x,y,y )dx
Mth 5 2006-2007) Yum-Tong Siu Lgrnge Multipliers nd Vritionl Problems with Constrints Integrl Constrints. Consider the vritionl problem of finding the extremls for the functionl J[y] = F x,y,y )dx with
More informationAPPROXIMATE LIMIT CYCLES FOR THE RAYLEIGH MODEL
ROMAI J, 4, 228, 73 8 APPROXIMATE LIMIT CYCLES FOR THE RAYLEIGH MODEL Adelin Georgescu, Petre Băzăvn, Mihel Sterpu Acdemy of Romnin Scientists, Buchrest Deprtment of Mthemtics nd Computer Science, University
More informationf(x)dx . Show that there 1, 0 < x 1 does not exist a differentiable function g : [ 1, 1] R such that g (x) = f(x) for all
3 Definite Integrl 3.1 Introduction In school one comes cross the definition of the integrl of rel vlued function defined on closed nd bounded intervl [, b] between the limits nd b, i.e., f(x)dx s the
More information221B Lecture Notes WKB Method
Clssicl Limit B Lecture Notes WKB Method Hmilton Jcobi Eqution We strt from the Schrödinger eqution for single prticle in potentil i h t ψ x, t = [ ] h m + V x ψ x, t. We cn rewrite this eqution by using
More informationCandidates must show on each answer book the type of calculator used.
UNIVERSITY OF EAST ANGLIA School of Mthemtics My/June UG Exmintion 2007 2008 ELECTRICITY AND MAGNETISM Time llowed: 3 hours Attempt FIVE questions. Cndidtes must show on ech nswer book the type of clcultor
More informationHomework 4 , (1) 1+( NA +N D , (2)
Homework 4. Problem. Find the resistivity ρ (in ohm-cm) for piece of Si doped with both cceptors (N A = 9 cm 3 ) nd donors (N D = 6 cm 3 ). Since the electron nd hole mobilities depend on the concentrtion
More informationNumerical Integration
Chpter 5 Numericl Integrtion Numericl integrtion is the study of how the numericl vlue of n integrl cn be found. Methods of function pproximtion discussed in Chpter??, i.e., function pproximtion vi the
More information1 Linear Least Squares
Lest Squres Pge 1 1 Liner Lest Squres I will try to be consistent in nottion, with n being the number of dt points, nd m < n being the number of prmeters in model function. We re interested in solving
More informationLecture 14 Numerical integration: advanced topics
Lecture 14 Numericl integrtion: dvnced topics Weinn E 1,2 nd Tiejun Li 2 1 Deprtment of Mthemtics, Princeton University, weinn@princeton.edu 2 School of Mthemticl Sciences, Peking University, tieli@pku.edu.cn
More informationReview of Calculus, cont d
Jim Lmbers MAT 460 Fll Semester 2009-10 Lecture 3 Notes These notes correspond to Section 1.1 in the text. Review of Clculus, cont d Riemnn Sums nd the Definite Integrl There re mny cses in which some
More informationUndergraduate Research
Undergrdute Reserch A Trigonometric Simpson s Rule By Ctherine Cusimno Kirby nd Sony Stnley Biogrphicl Sketch Ctherine Cusimno Kirby is the dughter of Donn nd Sm Cusimno. Originlly from Vestvi Hills, Albm,
More informationLYAPUNOV-TYPE INEQUALITIES FOR NONLINEAR SYSTEMS INVOLVING THE (p 1, p 2,..., p n )-LAPLACIAN
Electronic Journl of Differentil Equtions, Vol. 203 (203), No. 28, pp. 0. ISSN: 072-669. URL: http://ejde.mth.txstte.edu or http://ejde.mth.unt.edu ftp ejde.mth.txstte.edu LYAPUNOV-TYPE INEQUALITIES FOR
More informationApplication of Exp-Function Method to. a Huxley Equation with Variable Coefficient *
Interntionl Mthemticl Forum, 4, 9, no., 7-3 Appliction of Exp-Function Method to Huxley Eqution with Vrible Coefficient * Li Yo, Lin Wng nd Xin-Wei Zhou. Deprtment of Mthemtics, Kunming College Kunming,Yunnn,
More informationMA Handout 2: Notation and Background Concepts from Analysis
MA350059 Hndout 2: Nottion nd Bckground Concepts from Anlysis This hndout summrises some nottion we will use nd lso gives recp of some concepts from other units (MA20023: PDEs nd CM, MA20218: Anlysis 2A,
More informationReversing the Chain Rule. As we have seen from the Second Fundamental Theorem ( 4.3), the easiest way to evaluate an integral b
Mth 32 Substitution Method Stewrt 4.5 Reversing the Chin Rule. As we hve seen from the Second Fundmentl Theorem ( 4.3), the esiest wy to evlute n integrl b f(x) dx is to find n ntiderivtive, the indefinite
More informationADVANCEMENT OF THE CLOSELY COUPLED PROBES POTENTIAL DROP TECHNIQUE FOR NDE OF SURFACE CRACKS
ADVANCEMENT OF THE CLOSELY COUPLED PROBES POTENTIAL DROP TECHNIQUE FOR NDE OF SURFACE CRACKS F. Tkeo 1 nd M. Sk 1 Hchinohe Ntionl College of Technology, Hchinohe, Jpn; Tohoku University, Sendi, Jpn Abstrct:
More informationMath 270A: Numerical Linear Algebra
Mth 70A: Numericl Liner Algebr Instructor: Michel Holst Fll Qurter 014 Homework Assignment #3 Due Give to TA t lest few dys before finl if you wnt feedbck. Exercise 3.1. (The Bsic Liner Method for Liner
More informationV. More General Eigenfunction Series. A. The Sturm-Liouville Family of Problems
V. More Gener Eigenfunction Series c 214, Phiip D. Loewen A. The Sturm-Liouvie Fmiy of Probems Given n re interv, b], functions p(x), q(x), r(x), nd constnts c, c 1, d, d 1, consider (ODE) (p(x)y (x))
More informationPhysics 116C Solution of inhomogeneous ordinary differential equations using Green s functions
Physics 6C Solution of inhomogeneous ordinry differentil equtions using Green s functions Peter Young November 5, 29 Homogeneous Equtions We hve studied, especilly in long HW problem, second order liner
More information1 1D heat and wave equations on a finite interval
1 1D het nd wve equtions on finite intervl In this section we consider generl method of seprtion of vribles nd its pplictions to solving het eqution nd wve eqution on finite intervl ( 1, 2. Since by trnsltion
More informationBEAM DIAGRAMS AND FORMULAS. Nomenclature
BEA DIAGAS AND FOULAS Nomencture E = moduus of esticity of stee t 9,000 ksi I = moment of inerti of em (in. 4 ) L = tot ength of em etween rection points (ft) m = mimum moment (kip-in.) = mimum moment
More informationS. S. Dragomir. 2, we have the inequality. b a
Bull Koren Mth Soc 005 No pp 3 30 SOME COMPANIONS OF OSTROWSKI S INEQUALITY FOR ABSOLUTELY CONTINUOUS FUNCTIONS AND APPLICATIONS S S Drgomir Abstrct Compnions of Ostrowski s integrl ineulity for bsolutely
More informationu = 0 in Ω, u = f on B,
Preprint April 1991 Colordo Stte University Deprtment of Mthemtics COMPUTATION OF WEAKLY AND NEARLY SINGULAR INTEGRALS OVER TRIANGLES IN R 3 EUGENE L. ALLGOWER 1,2,4, KURT GEORG 1,3,4 nd KARL KALIK 3,5
More informationAn iterative method for solving nonlinear functional equations
J. Mth. Anl. Appl. 316 (26) 753 763 www.elsevier.com/locte/jm An itertive method for solving nonliner functionl equtions Vrsh Dftrdr-Gejji, Hossein Jfri Deprtment of Mthemtics, University of Pune, Gneshkhind,
More informationSIMULATION OF WOOD DRYING STRESSES USING CVFEM
Ltin Americn Appied Reserch 4:2- (2) SIMULATION OF WOOD DRYING STRESSES USING CVFEM C. H. SALINAS, C. A. CHAVEZ, Y. GATICA nd R.A. ANANIAS Deprtmento de Ingenierí Mecánic, Universidd de Bío-Bío, Av. Coo
More informationA Modified ADM for Solving Systems of Linear Fredholm Integral Equations of the Second Kind
Applied Mthemticl Sciences, Vol. 6, 2012, no. 26, 1267-1273 A Modified ADM for Solving Systems of Liner Fredholm Integrl Equtions of the Second Kind A. R. Vhidi nd T. Dmercheli Deprtment of Mthemtics,
More informationDynamic analysis of an underwater suspension system in ocean currents with method of averaging
Indin Journ of Geo Mrine Sciences Vo. (),November 7, pp. -7 Dynmic nysis of n underwter suspension system in ocen currents with method of verging Xing Li,, Jifeng ui,, Min Zho,,* & Tong Ge, Stte Key Lbortory
More informationThe Riemann Integral
Deprtment of Mthemtics King Sud University 2017-2018 Tble of contents 1 Anti-derivtive Function nd Indefinite Integrls 2 3 4 5 Indefinite Integrls & Anti-derivtive Function Definition Let f : I R be function
More informationAN INTEGRAL INEQUALITY FOR CONVEX FUNCTIONS AND APPLICATIONS IN NUMERICAL INTEGRATION
Applied Mthemtics E-Notes, 5(005), 53-60 c ISSN 1607-510 Avilble free t mirror sites of http://www.mth.nthu.edu.tw/ men/ AN INTEGRAL INEQUALITY FOR CONVEX FUNCTIONS AND APPLICATIONS IN NUMERICAL INTEGRATION
More informationGeneralizations of the Basic Functional
3 Generliztions of the Bsic Functionl 3 1 Chpter 3: GENERALIZATIONS OF THE BASIC FUNCTIONAL TABLE OF CONTENTS Pge 3.1 Functionls with Higher Order Derivtives.......... 3 3 3.2 Severl Dependent Vribles...............
More informationMonte Carlo method in solving numerical integration and differential equation
Monte Crlo method in solving numericl integrtion nd differentil eqution Ye Jin Chemistry Deprtment Duke University yj66@duke.edu Abstrct: Monte Crlo method is commonly used in rel physics problem. The
More informationLecture 19: Continuous Least Squares Approximation
Lecture 19: Continuous Lest Squres Approximtion 33 Continuous lest squres pproximtion We begn 31 with the problem of pproximting some f C[, b] with polynomil p P n t the discrete points x, x 1,, x m for
More informationNUMERICAL INTEGRATION. The inverse process to differentiation in calculus is integration. Mathematically, integration is represented by.
NUMERICAL INTEGRATION 1 Introduction The inverse process to differentition in clculus is integrtion. Mthemticlly, integrtion is represented by f(x) dx which stnds for the integrl of the function f(x) with
More informationWe partition C into n small arcs by forming a partition of [a, b] by picking s i as follows: a = s 0 < s 1 < < s n = b.
Mth 255 - Vector lculus II Notes 4.2 Pth nd Line Integrls We begin with discussion of pth integrls (the book clls them sclr line integrls). We will do this for function of two vribles, but these ides cn
More information