Email faizjsofeeceutacir M OJAGHI Zanjan Regional Electric Company Zanjan Iran There are different numerical techniques for computing electric fields These numerical techniques enable the designer to study the problems that are difficult to be solve ID: 25491
Download Pdf The PPT/PDF document "Instructive Review of Computation of Ele..." is the property of its rightful owner. Permission is granted to download and print the materials on this web site for personal, non-commercial use only, and to display it on your personal computer provided you do not modify the materials and that you retain all copyright notices contained in the materials. By downloading content from our website, you accept the terms of this agreement.
InstructiveReviewofComputationofElectricFieldsusingDifferentNumericalJAWADFAIZDepartmentofElectricalandComputerEngineering,FacultyofEngineering,UniversityofTehran,Tehran,Iran.E-mail:faiz-j@sofe.ece.ut.ac.irM.OJAGHIZanjanRegionalElectricCompany,Zanjan,IranTherearedifferentnumericaltechniquesforcomputingelectricfields.Thesenumericaltechniquesenablethedesignertostudytheproblemsthataredifficulttobesolvedbyanalyticalmethods.Thispaperattemptstogiveaninstructivereviewofdifferentnumericaltechniquesinelectricfieldanalysis.Thesetechniquesinvolvefinitedifference,finiteelement,boundaryelement,chargesimulation,finiteelementwithvariablefieldintensityandMonteCarlomethods.Themeritsandlimitsofthevariousmethodsareoutlined.Someexamplesaregiveninwhichthefieldcomputationsusingdifferentnumericaltechniquesarecompared.Asanexamplethemostconvenienttechniqueapplicabletoelectricfieldcomputationwithinthetankofpowertransformersisintroduced.INTRODUCTIONUNTILNOWelectricandelectromagneticfieldshavebeenregardedasprobablythemostabstractanddifficultpartoftheundergraduateelectricalengineeringcurriculum.ThisislargelyduetothefactthatsuchfieldscannotbevisualizeddirectlydirectlyElectricfieldscanbecomputedusingvariousmethodswithdifferentprecision.However,forinsulationofelectricalequipmentamoreaccurateelectricfieldpredictionisrequired.Withtheadventofcomputingpowerandnumericaltech-niquesinrecentyears,ithasbecomepracticaltousedifferenttechniquestocomputetheelectricfields.Suchnumericaltechniquesenablethedesignertosolveproblemsthataredifficult,anduseofanalyticalapproachwithmanyempiricalfactorsisimpossible.Theaimofthispaperistogiveareviewoftheapplicationofdifferentnumericaltechniquesintheelectricfieldscomputationthatisusefulforanundergraduatecourse.Thismaterialshouldenablethegraduatestudenttousetheseanalysesingraduateresearchaswellaslaterintheworkplace.Theauthorswillalsogivetheirpersonalviewsonthefieldevaluationoftransformers.Intwo-dimensional(2D)analysis,ifthefieldisnottimevarying,theelectricpotentialVintheactualspaceissatisfiedbyPoisson'sequationquationr2Vÿ= r0 1whereistheelectricchargedensity,relativeandabsolutepermittivityofthefreespace.InordertosolvepartialdifferentialEqn.1,anumberofboundaryconditionsmustbeimposed.Analyticalsolutionsforpracticalboundarycon-ditionsaredifficult,ifnotimpossibleandtherefore,numericaltechniquesarenecessary.Eachtechniquehasitsownmeritsanddraw-backsandonetechniquecannotbegenerallypreferredtoothers.Basedontheproposedproblem,themostconvenienttechniquemustbeselected.Differentnumericaltechniques,sofarusedfortheelectricfieldanalysis,arebrieflyreviewedandthemostconvenienttechniqueisthensuggestedforelectricfieldevaluationwithintheinteriorspaceofpowertransformers.BOUNDARYCONDITIONSTherearetwofollowingtypesofboundaryconditionsthatareconsideredintheelectricfield1.Boundarybetweenconductinganddielectricmaterials.2.Boundarybetweendifferentdielectricmaterials.Fromtheelectricalpointofview,oneofthefollowingconditionsmaybesatisfiedontheboundariesofthefirsttype:a.Theelectricpotentialsofallpointsontheboundaryareknown.Thisispossibleiftheavailableconductorisconnectedtoafixedpotentialsource(Dirchletcondition). *AcceptedInt.J.EngngEd.Vol.18,No.3,pp.344±356,20020949-149X/91$3.00+0.00PrintedinGreatBritain.2002TEMPUSPublications. Thewholechargeontheboundarysurfaceisknownwhiletheelectricpotentialoverdifferentpointsisunknown.Thisisthecasewhentheavailableconductorisnotconnectedtoafixedpotential(floating).Normallyinsuchacasethewholechargeonthesurfaceboundaryisequaltozero(Neumanncondition).FINITEDIFFERENCEMETHODFinitedifferencemethod(FDM)istheoldesttechniqueinthefieldcomputationsthatwasintro-ducedbyGauss.ThenBoltzmanpublisheditinhisnotesin1892.However,becauseofthelargeamountofcomputationsrequired,theextensiveuseoftheFDMdatesbackonlytotheeventofthecomputer.ElectricpotentialVoveranyregiondependsonthe(x,yandz)coordinatesanditsderivatives.PotentialofanypointmaybegivenversustheelectricpotentialoftheadjacentpointsusingTylor'sexpansion.Forinstance,considerFig.1havingpoint0andsixadjacentpoints.Thecoor-dinatesofthesepointsare:P0 x0;y0;z0;P1 x0;y0;z0h;P2 x0;y0ÿh;z0;P3 x0;y0;z0ÿh;P4 x0;y0h;z0;P5 x0h;y0;z0;P6 x0ÿh;y0;z0Forsimplicity,thedifferencebetweentheoriginandtheadjacentpointsistakentobeh.PotentialofanarbitrarypointP x;y;zmaybecalculatedversuspotentialofpointP0asfollows:V x;y;zV x0;y0;z0 xÿx0Vx x0;y0;z0 yÿy0Vy x0;y0;z0 zÿz0Vz x0;y0;z0=1! xÿx02Vxx x0;y0;z02 xÿx0 yÿy0Vxy x0;y0;z02 xÿx0 zÿz0Vxz x0;y0;z0 yÿy02Vyy x0;y0;z02 yÿy0 zÿz0Vyz x0;y0;z0 zÿz02Vzz x0;y0;z0=2!... 2whereVi@V=@i;Vii@2V=@i2;Vij@2V=@i@j;i;jx;y;z;h xÿx0 yÿy0 zÿz0Ifhtendstoaverysmallvalue,thetermscontain-ingthethird-andhigher-orderderivativesmaybeignoredandEqn.2canberewrittenasfollows:V1V0hVz P00:5h2Vzz P0V2V0ÿhVy P00:5h2Vyy P0V3V0ÿhVz P00:5h2Vzz P0V4V0hVy P00:5h2Vyy P0V5V0hVx P00:5h2Vxx P0V6V0ÿhVx P00:5h2Vzz P0 3X6i1Vi6V0h2Vxx P0Vyy P0Vzz P0 4SolutionofEqn.4andsubstitutingfromEqns.1±2yields:V0V1V2V3V4V5V6h20= 5where0istheelectricchargedensityatpoint0.AsshowninEqn.5,thereisalinearrelationshipbetweenthepotentialofpoint0andthepotentialsoftheadjacentpoints.Forthecaseswherediffer-encebetweenthepoint0andtheadjacentpointsisnotthesameorthesepointsarewithindifferentinsulatingmaterials,andalsoforrotatingfieldsinelectricalmachines,similarequationscanbederived[3,4].IntheFDM,theproposedregionisdiscretizedusingtheequationssimilartoEqn.5.Dimensionsofthemeshesmustbesuchthattheapproximationisacceptable.Thevertexesofthemeshesarenodesontheboundaryoftheregionandtheirpotentialsareknownortheycorrespondtopoint0inFig.1,enclosedbytheothernodes.Forthelatternote(notontheboundary),equationssimilartoEqn.5canbewrittenversusthepotentialsoftheadjacentnodes.Ifthenumberofsuchnodesisequalton,nlinearalgebraicequationswithnunknownvaluesofthenodepotentialsareobtained.Potentialsofthepointsinsidethemeshesmaybedeterminedusingdifferentinterpolationtechniques.TheFDMisnotcapableofcalculatingelectricfielddirectlyatdifferentpointsontheproposedregion.Whenthepotentialsofthenodesareobtained,anumericalderivativeevaluationtech-niqueisusedtocalculatetheelectricfieldintensity:EÿrV. Fig.1.Point0andsixadjacentpoints.InstructiveReviewofComputationofElectricFieldusingDifferentNumericalTechniques345 ELEMENTMETHODSAmongthevariousnumericaltechniques,thefiniteelementmethod(FEM)hasadominantpositionbecauseitisversatile,havingastronginterchangeabilityandcanbeincorporatedintostandardprograms[5,6].FEMisbasedonthisfactthatthephysicalsystemsstabilizesattheminimumlevelofenergy.Thegeneralequationofenergyinanelectricfieldis:W0:5 vrv rvdvÿ vdvÿ S0svds 6wherevisthevolumeoftheproposedregionistheelectricchargedensitywithinvolumevS0istheboundarysurfacewithNummanconditionssisthesurfaceelectricchargedensitywithinS0.SinceEqn.1describestheelectricpotentialdis-tributionintherealsystems,basedonthemini-mumenergyleveltheorem,itisconcludedthatEqn.1minimizestheenergypresentedbyEqn.6.IntheFEM,thevolumeoftheproposedregionisdividedintoMsmallpolyhedronelementswheretheirsidesformagridwithNnodes.Thepotentialfunctionisthenapproximatedby:V rXNi1fi rVi 7whererisanypointontheproposedregion.fi riscalledtheshapefunctionhavingthefollowingfeatures:a)fi risequaltozeroanywhere,exceptonthesubregionWi.Thesub-regionWiconsistsoftheelementswherenodeiisoneoftheirvertices.b)fi riscontinuousontheboundariesduetoWi,andapolyhedroninsideeachelement.c)anyfi risequaltounitatthelocationofnodeiandzeroattheothernodes:fi ri1forijfi rj0fori6jViinEqn.7isequaltothepotentialofnodei.SubstitutingEqn.7intoEqn.6,theapproximateenergyispresentedbyWwhichisminimizedunderthefollowingconditions:@W=@Vi0;i1;2;...;N 8SinceWisaquadraticfunctionofVi,applyingconditions8leadstothefollowinglinearalgebraicequations:GVA 9whereVistheknownvectorwithelementsVi;Aistheknownvectorobtainedfromthevolumechargedensityintheproposedregionandbound-aryconditions;Gisthenon-singularsquaresymmetricalmatrix.SolutionofthissystemofequationsgivesthevaluesofVisandhenceanapproximatedistribu-tionofthepotentialcanbedeterminedbasedonEqn.7.Electricfieldintensitywithineachelementisobtainedusingthegradientexpressionasfollows:Em ÿrVmÿXNi1Vi:fmi 10Oftenthefirstderivativeoffiisnon-continuous.Therefore,reductionofthemaximumsizeoftheelementsandtendingtozero,leadsEqn.7totherealdistributionofthepotential.Inspiteofthis,nocontinuityofthefieldintensityonthebound-ariesoftheelementsremainsinforce.Iffunctionfiisconsideredasacompleten-orderpolynomial,betterresultscanbeobtained.Ifhpresentsthemaximumsizeoftheelements,reductionofhcanreducethepotentialerrorwithratio hn1;andelectricfieldiscontinuousanditserrorisreducedbyratio hn[7].Fig.2showsthemeshingandequipotentiallinesdeterminedusingtheFEM.TheFEMcouldbealsousedwherethepermittivityoftheproposedregionisnotconstant.Insuchacase,itisnecessarytoreplace rwith,whichshowsthepositiondependency.Methodshavebeenintroducedtoconsiderthefloatingelectrodeswithunknownpotentialordifferentinsulatingmaterials[7].Thereareseveralreportsforautomaticmeshingoftheproposedregion[9±11].BOUNDARYELEMENTMETHODIfdistributionofelectricchargeforeveryregion(includingboundarysurfaces)isknown,electricpotentialandfieldintensityforeachpointcanbecomputedusingCoulomb'slawofGauss'slaw[2].Inpractice,Laplacianequationisnormallyused.Thismeansthatelectricchargeisenclosedonlyinsidetheboundariesoftheproposedregionandvolumechargedensityinsidetheregionisequaltozeroornegligible.Theelectricpotentialandfieldintensityare:V I ss= 4Rds 11E I ss 4R2ds:^aR 12whereIistheproposedpoint,Sisthesummaryofallboundarysurfaces,sisthesurfacechargedensityoversurfaces,RisthedistanceofthedifferentialelementsfrompointIand^aRisunitvectoralongRdirectedfromdifferentialelementdstopointI.InpracticesisunknownanditseemsthatsolvingEqns.11±12isimpossible.However,toJawadFaizandM.Ojaghi346 thisdifficulty,theBEMmaybeemployed.Inthismethod,theboundarysurfaceisdividedintoNelements.Figure3showsatypicalboundaryelementsduetoaflatboundarysurfaceonxy-plane.Thentakingintoaccountthesmalldimensionsoftheboundaryelements,ageneralformmaybeconsid-eredforthesurfacechargedensityonelementj sj.Thisgeneralformisoftenapolynomialwithunknowncoefficients.Forinstance,ifitistakentobeaquadraticpolynomialasfollows:sja1ja2jxa3jya4jxy 13substitutingsjinEqns.11and12gives:V IXNj1 sjs= 4RdsXNj1 sja1ja2jxa3jya4jxy= 4Rds 14E IXNj1 sjs= 4R2ds^aR Fig.2.ComputationofelectricfieldusingFEM:a)meshing,b)equipotentiallines[4].InstructiveReviewofComputationofElectricFieldusingDifferentNumericalTechniques347 XNj1 sja1ja2jxa3jya4jxy= 4R2ds^aR 15ItisclearthatV(I)andE(I)arelinearfunctionsofthecoefficientsofthepolynomials.Byusingthefollowingdefinitions:P1j sjds= 4RP2j sj xds= 4RP3j sj yds= 4RP4j sj xyds= 4Rf1j sjds= 4R2^aRf2j sj xds= 4R2^aRf3j sj yds= 4R2^aRf4j sj xyds= 4R2^aREqns.14and15become:V IXNj1X4k1Pkjakj 16E IXNj1X4k1fkjakj 17Inthenextstage,thenumberofunknownsfunc-tionPsjisselectedoneachboundaryelement.Then,basedontheboundarytypewheresjispartofit,theboundaryconditionsequationforanyselectedpointisformedusingEqns.16and17.Therefore,asystemofalgebraiclinearequationsisobtainedwhichfinallyproducesthecoefficientsofthepolynomial.Pkjandfkjarenumericallyoranalyticallyobtainedbyintegration.Hencethesurfacechargedensitydistributiononallboundarysurfacesisknown.FinallytheelectricpotentialandfieldintensitycanbedeterminedusingEqns.16and17. Fig.3.Aflatquadrangularboundaryelement. Fig.4.ElectricfieldcalculationusingBEM:a)actualelectricalsystem,b)electricfielddensityvectorsonthesphere,c)equipotentialcountors[12].JawadFaizandM.Ojaghi348 4presentsatypicalproblemsolvedbytheBEM,wheretheelectricfieldbetweenthetwosphereshavingthesameradiusr1isobtained.Electricfieldintensityhasbeenshownwithasuit-ablevector.Inaddition,theequipotentialsbetweentwosphereshavebeencalculatedbycomputationoftheelectricpotentialindifferentpoints.Amethodhasbeengivenin[13]forthecurved-shapeboundaryelementsinordertomodelprac-ticalsurfaceswiththedesirableaccuracy.Bound-aryelementswithfullaxialsymmetryhavebeenpresentedintheliterature[14±16].InputandoutputdataprocessingmethodsintheBEMhavebeenintroducedin[12,17,18].CHARGESIMULATIONMETHODThismethodissimilartotheBEM.Thediffer-encebetweenchargesimulationmethod(CSM)andtheBEMisthesimulationofthesurfacechargeexistingontheboundarysurfaces.IntheBEMthesurfacechargedensityfunctiononthedifferentsurfaceboundariesareestimated,whileintheCSMthesurfacechargedensityissubstitutedbyasetofdiscretizedlinearchargedistribution.Thesubstitutedlinearchargedistributionissuchthattheelectricpotentialandfieldintensityversustheirchargesareanalyticknownfunctions.Chargedistributiononaninfinitelengthlinewithaconstantdensity,onafiniteline,onacircleetc.,areexamplesofthechargedistribution.Differentdistributiontypesandtheirequationshavebeengivenin[4]and[19].Electricpotentialandfieldintensityequationsduetotheabovementionedchargedistribution,forthepointsonthecharges,havesingularity.Toovercomethissingularity,thepositionofthereplacedchargesimulationisconsideredoutsidetheproposedspaceandnormallyinsidetheelec-trodes.Theexactpositionandthereplacedchargedistributiontypearearbitrarilyselectedbasedontheexperience.Chargevalueortheirlinearchargedensityiscomputedsuchthattheboundarycon-ditionsaresatisfiedonsomesurfaceboundaries,asdescribedinthefollowingpart.Electricpotentialandfieldintensityduetothereplacedchargesimulationindifferentpointsarelinearfunctionsofthechargevalueorchargedensity:Vj rPj rqj 18Ej rfj rqj 19wherePj risthepotentialfactor,fj risthefieldintensityfactor,ristherelativepositionoftheproposedpointandqiisthejthreplacedchargedistribution.Pj randfj raredifferentfordiffer-entchargedistributions.Thereisthefollowingrelationshipbetweenthem:fj rÿrPj r 20Sincethereisalinearrelationshipbetweentheelectricpotential(andfieldintensity)duetothereplacedchargedistributionindifferentpointsandchargevalueorchargedensity,forasetofsuchchargedistribution,thesuperpositiontheoremcanbeappliedtocalculatetheelectricpotentialandfield:Vj rXNj1Pj rqj 21Er rXNj1fj rqj 22IntheCSM,thenumberofpointsselectedonthesurfaceboundaryisequaltothereplacedsimula-tioncharges.Dependingontheselectedpoints,Eqns.21and22areusedandtheboundaryconditionsforindividualpointsareconsidered.Theseequationsarelinearfunctionsofqj.Sincetheexactpositionofthechargesandtheselectedpointsareknown,PjsandfjsinEqns.21and22areexactlycalculated.Thenonlyqisaretheunknownvaluesoftheabovelinearequations.Therefore,computationsoftheelectricpotentialandfieldintensityarepossibleusingEqns.21and22.ItisclearthatintheCSM,boundarycondi-tionsaresatisfiedonlyinthepointsselectedtowritetheequations.Beforeusingqjsduetothesolutionoftheequations,itisnecessarytostudytheboundaryconditionsonotherpointsoftheboundarysurfaces.WhenqjsintheCSMiscalculated,theboundaryconditionsondifferentsurfacepointsmustbedetermined.Iftheaccu-racyisnotenough,number,positionandtypeofthereplacedsimulationchargesandalsopositionoftheselectedpointsforwritingtheequationsmustbevariedinordertoobtainsufficientaccuracy.Figure5showsatypicalproblemsolvedusingtheCSM.TheproblemwascomputationoftheelectricfieldbetweentwospheresshowninFig.5a.Figure5bpresentstheequipotentialslinesbysubstitutingthesurfacechargeofeachspherewithtwopoint-charges.Figure5cindicatesthecorrespondingresultwhenthreepoint-chargesaresubstitutedforeachsphere,inwhichtheaccuracyishigher.Inaddition,twopoint-chargeshavebeenusedinFig.5d,butthepositionofchargesisdifferentwiththatofFig.5b.Accuracyofthelattercaseisbetterthantheothertwo.FINITEELEMENTMETHODWITHVARIABLEFIELDDENSITYIntheFEM,themainvariableistheelectricpotentialV,whiletheelectricfieldintensityisnormallytherequiredquantity.Forexample,inthedesignofhighvoltagedeviceinsulation,itisInstructiveReviewofComputationofElectricFieldusingDifferentNumericalTechniques349 rytohavetheamplitudeanddirectionoftheelectricfieldonthesurfacesoftheelectrodesandinsulation.Instudyofthedischargephenom-enon,thepathoftheforcelinesisrequiredwhichcanbedeterminediftheelectricfielddistributionintheproposedspaceisknown.IntheFEManumericalintegrationtechniqueisusedinordertocalculateEthatnormallyhaserror.ThereisamoreaccurateFEMtechniqueinwhichthemainvariablehasbeentakenEinsteadofV[20].Ifthevolumechargedensityintheproposedspaceiszeroand Fig.5.CalculationofelectricfieldusingCSM:a)electricsystem,b)equipotentialwhentwochargepointsusedforeachsphere,c)equipotentialwhenthreechargepointsusedforeachsphere,d)equipotentialwhentwochargepointsindifferentpositionsusedforeachsphere.JawadFaizandM.Ojaghi350 ttivityEisconstantandisotropic,thentheMaxwell'sequationsforelectrostaticfieldsareasfollows:divD0 23CurlE0 24whereDistheelectricchargedensityandEistheelectricfieldintensityand:DE 25ApplyingEqn.24andusingEqn.23and25leadsto:r2E0 26whichistheLaplacianequation.For2Dfields: @2Ex=@x2@2Ex=@y2i @2Ey=@x2@2Ey=@y2j0 27Butavectorwillbezeroifallitscomponentsarezero:@2Ex=@x2@2Ex=@y20@2Ey=@x2@2Ey=@y20 28Foreachcomponentofthefieldintensity,theFEMisusedonceandfinallytheelectricfieldintensityoverthewholeproposedregionisobtained.However,withdefinedboundaryconditions,theequationwillhaveauniquesolution.Inpractice,theboundaryconditionsonthefirst-typeboundaryversuselectricpoten-tialareknownbutthereisnoideaconcerningtheirelectricfieldintensity.In[20],theBEMhasbeenusedinordertoovercomethisdifficult;itmeansthattheBEMisemployedonthementionedboundariesthenfiniteelementmethodwithvariablefielddensity(FEMVFD)isapplied.AnalysishasbeencarriedoutbasedontheFEMandFEMVFDandtheresultshavebeenpresentedinFig.6.Theproblemwascomputationoftheelectricfieldintheregionbetweenthetwocylindershavingpotentialsof100Vand200V.Becauseoftheaxialsymmetryanditsboundariesandinfinitelength,thefieldhasonlyaradialcomponentvaryinginradialdirection.Therefore,analysisofthefieldispossibleintwodimensions.Duetothesymmetryonlyone-quarterofthecylinderisusedforanaly-sis.Figure6bshowsthemeshingforbothtech-niquesandFig.6cpresentstheequi-fieldintensitylines.Asseen,inspiteofthelowernumberofmeshesintheFEMVFD,theequi-fieldintensitylineshaveabettercontinuityandtheiraccuracyisalsohigher.Inthisexample,themaximumerrorinsolutionbytheFEMVFDisabout2%whilethiserroris9.3%whentheFEMisused.MONTECARLOMETHODIntheFDM,thepotentialofeverynodeisequaltothemeanvalueofthepotentialsoftheadjacentnodes.Whenthedistancebetweentheproposednodeandadjacentnodesarethesame,sameweightpotentialsareinvolvedinthecomputationofthemeanvalue.Otherwise,theweightswillnotbeequal[4].Buttheclosernodeswillbeheavier.However,thesumofallweightsisunity.GenerallyintheFDM,potentialofeachnode(V0versuspotentialofnadjacentnodesisasfollows:V0XNI1WiVi 29whereWi1,V1isthepotentialofthenodeadjacenttoi-thnodeandWiisitsweight.ValueofWidependsontherelativedistanceofthei-thnodefromtheproposednode.ThebasicequationofMonteCarloMethod(MCM)issimilartoEqn.29.Therefore,thesetwomethodsaresimilar,exceptthatdeterminationoftheadjacentnodesandmethodofcalculationofWidiffers.IntheFDM,allnodesaredefinedaftermeshingprocessandatthesametimeadjacentnodesofeverynodearedetermined.InWicom-putation,analyticalrelationshipsareused.ButintheMCM,adjacentnodesarealwaysontheboundariesandtheirexactpositionsaredeter-minedusingarandomprocessandWisisobtainedusingprobabilitytechniques.InFig.7,Sisshowboundarysurfacesoftheproblem.TheseboundarysurfacesareDrichlettypewithpotentialVi.Calculationofpotentialatanypointsuchasr0isrequiredusingequationssimilartoEqn.29.IntheMCM,simulationofarandommovementisusedinordertodetermineeachinode.Anymovementbeginsfromr0andaftersuccessivejumpswithvariablelengthsandinrandomdirections,leadstoaSj.Conditionsgoverningoneveryrandommove-mentareasfollows:a)Allmovementsbeginfromr0.b)Lengthofanyjumpisequaltotheminimumdistanceofthebeginningpointwithboundarysurfaces(Sis).c)Directionofveryjumpisrandom.d)Theendofanyrandommovementwillreachwhentheminimumdistancementionedinitembissmallerthanthatofthepredefinedvaluesuchas0.Attheendofarandommovement,apointononeoftheboundarieswiththeclosestdistancetotheendpointisselectedasanadjacentpoint.Finally,potentialofpointr0iscalculatedasfollows:V r0XNi1V ri=N 30whereNisthetotalnumberoftherandomInstructiveReviewofComputationofElectricFieldusingDifferentNumericalTechniques351 entsandriistheadjacentpointduetothej-thrandommovement.Inordertoeliminatethestatisticalerrors,itisnecessarythatNbecomeslargeenoughandoftenseveralthousands.InEqn.30,allweightsareapparentlythesameandequalto1/N.Infact,itisnotso,becausetakingintoaccountthelargenumberofrandommovements(N),probabilityofseveralreputationsofonepointexistsandforrelativelyshorter Fig.6.ComparisonofcomputedelectricfieldusingFEMandFEMVFD:a)electricsystem,b)meshingbasedonthetwotechniques,c)equi-fieldintensitylinesusingbothmethods.JawadFaizandM.Ojaghi352 cefromr0,thisprobabilityislarger.SupposeconstantpotentialViforanyboundarysurface,Eqn.30canbewrittenasfollows:V r0Xni1NiVi 31wherenisthenumberofboundarysurfacesandNiisthenumberofrandommovementsendedtoapointatSi.AlthoughtheMCMwasdescribedusingtheFDM,itisnecessarytonotethattheMCMisitselfanindependentmethodhavingspecial Fig.7.Representationofthepossiblerandommovementpathsforreachingfromtheproposedpoint(r0)tooneoftheboundarysurfaces(Si). Fig.8.Geometryoflaboratoryhighvoltageelectrodes[7].InstructiveReviewofComputationofElectricFieldusingDifferentNumericalTechniques353 ntals[21].Computationofelectricpoten-tialusingMCMissimilartothemicroscopicanalysisofgaspressurewheretherandommove-mentofthemoleculesissimulated.Atechniqueforcomputationoftheelectricfieldoverthespacesconsistingofdifferentisolatingmaterialshasbeenpresentedin[22].In[7,22],anumberoftechniqueshavebeenintroducedtospeedupthecalculationusingtheMCM.Figure8showsthegeometryofthelaboratoryhighvoltageelectrodesinwhichthe Fig.9.ComparisonofthecomputedelectricfieldsofthesystemshowninFig.8usingthreetechniques:a)alongpathAB,b)alongpathCDandc)alongpathEF[7].JawadFaizandM.Ojaghi354 potentialandfieldalonglineAB,CDandEFhavebeenestimatedusingthethreemethods:MCM,CSMandFEM.TheresultsarepresentedinFig.9.ItisclearthattheaccuracyoftheMCMissimilartotheothertechniques.Thereareabout24000randommovementsinthissimulation.COMPARISONOFMETHODSANDCONCLUSIONSAnynumericaltechniqueforelectricfieldevaluationhasitsownmeritsanddrawbacksanditisnotgenerallypossibletopreferonetechniqueovertheothers.Structureofthepowertransformerindicatesthatnormally3Dcomputationoftheelectricfieldisrequiredbecausethereisnosymmetryinordertoignoreonedimension.Inaddition,numer-icaltechniquemustbecapableofdealingwiththenarrowinsulatingorconductinglayers.FDMandFEMintheelectricfieldcomputationhavetwodrawbacks:a)Indeterminationofpotentialdistributionusingtwomethods,numericalderivativetechniquesmustbeusedinordertoobtaintheelectricfieldintensity.Thishasconsiderableerrorthatleadstoalargeerrorintheelectricfieldcomputation.Inmanycases,thefieldisrequiredforthedesignofinsulationofelectricalequipment.Meanwhile,studyofsomephenomenasuchaselectricaldischargeispossiblebyelectricfieldcomputation.b)Inordertopreventalargeelectricfieldanditsdrawbacks,thesharpedgesonthedifferentsurfacesareavoided.Therefore,thecurvedsurfacesareoftenpreferred.FDM,FEMandFEMVFDhavedifficultyinmodelingsuchcurvedsurfaces;butCSMandBEMcanbeeasilyadaptedforsuchcases.Ontheotherhand,althoughFDM,FEMandFEMVFDcantheoreticallycompute3Dfields,therearemanyproblemsindealingwiththismatter.Oneseriousproblemis3Dmeshgenerationanditsmodificationtoapproachtherequiredaccuracy.Generallymanualcalculationiscumbersomeandtimeconsumingandalsocomputerprogrammingisreallycomplicated.Anotherdifficultyisthelargesizeofthecoeffi-cientmatrixofthesystemofequations.Inthesemethodsnumberofequationsisproportionaltothememoryrequiredforcomputationwhileinothermethods(CSMandBEM)thisnumberisproportionalwithareaoftheboundarysurfaces.ItmeansthatinFDM,FEMandFEMVFD,thecoefficientmatrixhasonemoredimensionthanthecoefficientmatrixduetotheotherabove-mentionedmethods.Therefore,morecomputermemoryandlongercomputationtimearerequired.HenceFDM,FEMandFEMVFDmaynotbecon-sideredconvenienttechniquesforelectricfieldcomputation.MCM,atleastinelectricfieldcomputation,isnotsocommonandhasnoconsiderableprogressinrecentyears.Atthepresent,applicationofthismethodinnarrowlayerproblemsisdifficult,ifnotimpossible.Inspiteofthesimplicityofcomputerprogram-mingandhighaccuracyofthemethod,incompu-tationof3DelectricfieldshavingnarrowlayersCSMisconfrontedwithamajordifficulty.Inthismethoditisnecessarytoconsiderchargeswithinthementionedlayerssuchthattheyhaveenoughdistancefromtwosidesofthelayer.Butthethicknessistoonarrowandsuchanassumptionmaynotbecorrect.FinallyBEMiscapableofanalyzing3Dfieldsandtherearesomereportsshowingitsapplicationstonarrowlayers.Inconclusion,BEMmaybeconsideredthemostconvenienttechniqueforelec-tricfieldcomputationwithintheinteriorofthepowertransformertank[23,24].REFERENCES1.N.N.Rao,PC-assistedinstructionintroductoryelectromagnetics,IEEETrans.onEducation,33(1),pp.51±55,Feb.1990.2.W.H.HaytJr,EngineeringElectromagnetics,4thEdition,McGraw-HillBookCo.,1981.3.BharatHeavyElectricals,Transformers,TataMcGraw-HillPublishingCompany,NewDelhi,1987.4.E.KuffelandW.S.Zeungl,HighVoltageEngineeringFundamentals,PergamonPress,1984.5.M.V.K.Chari,G.BedrosianandI.D.Angelo,Finiteelementapplicationsinelectricalengineering,IEEETrans.Magnetics,29,1993,pp.1306±1314.6.C.R.I.Emson,J.SimkinandC.W.Trowbridge,Astatusreportelectromagneticfieldcomputation,IEEETrans.Magnetics,30,1994,pp.1533±1540.7.M.D.R.Beasley,Comparativestudyofthreemethodsforcomputingelectricfields,IEEProc.,126(1),January1979,pp.126±134.8.O.W.Anderson,Finiteelementsolutionofcomplexpotentialelectricfields,IEEETrans.PowerApparatusandSystems,PAS-96(4),July/August1977,pp.1156±1161.9.O.W.Anderson,Laplacianelectrostaticfieldcalculationsbyfiniteelementswithautomaticgridgeneration,IEEETrans.PowerApparatusandSystems,PAS-92,Sept./Oct.1973,pp.1485±1492.10.A.KamitaniandM.Miyauchi,Automaticnumericalgenerationoffiniteelementsfortwo-dimensionalregionofarbitraryshape,Proc.Int.Symp.Electromagnetics(ISEM),Sapporo,Japan,1993,pp.89±92.InstructiveReviewofComputationofElectricFieldusingDifferentNumericalTechniques355 J.FaizandE.Shafagh,AutomaticmeshgenerationusingAuto-Cad,IEEETrans.Education,41(4),November1998,pp.325±330.12.H.Tsuboi,etal.Postprocessingforboundaryelementanalysisinelectromagneticfieldproblems,Proc.Int.Symp.Electromagnetics(ISEM),Sapporo,Japan,1993,pp.117±120.13.T.Misakietal.Computationofthree-dimensionalelectricfieldproblemsbyasurfacechargemethodanditsapplicationtooptimuminsulatordesign,IEEETrans.PowerApparatusandSystems,PAS-101(3),March1982,pp.627±634.14.S.SutoandB.Bachmann,Athreedimensionalhighspeedsurfacechargesimulationmethod(3D-HSSSM),Int.Symp.HighVoltageEngineering,Athens,Greece,Sept.1983,p.11.08.15.S.SatoandW.W.S.Zaengle,Effective3-dimensionalelectricfieldcalculationbysurfacechargesimulationmethod,IEEProc.,133,Pt.a,No.2,March1986,pp.77±83.16.M.ReisterandP.Weib,Computationofelectricfieldsbyuseofsurfacechargesimulationmethod,Int.Symp.HighVoltageEngineering,Athens,Greece,Sept.1983,p.11.06.17.M.NishinoandT.Takedo,Preprocessorforthe3Dboundaryelementmethodbasedontheconstructivesolid-geometry,Proc.Int.Symp.Electromagnetics(ISEM),Sapporo,Japan,1993,pp.93±96.18M.NakahiraandT.Takeda,Ameshrefinementforboundaryelementmethodusinganerrorestimatingindex,Proc.Int.Symp.Electromagnetics(ISEM),Sapporo,Japan,1993,pp.151±154.19.H.Singer,H.SteinbiglerandP.Weise,Achargesimulationmethodforthecalculationofhighvoltagefields,IEEETrans.PowerApparatusandSystems,PAS-93,Sept./Oct.1974,pp.1660±1668.20.H.Yamashita,K.ShinozakiandE.Nakamae,Aboundaryfiniteelementmethodtocomputedirectlyelectricfieldintensitywithhighaccuracy,IEEETrans.PowerDelivery,3(4),October1988,pp.1754±1760.21.R.M.Brvensee,Probabilisticpotentialtheoryappliedtoelectricalengineeringproblems,Proc.IEEE,61(4),April1973,pp.423±437.22.M.KrauseandK.Muller,AMonteCarlomethodfortwo-andthree-dimensionalelectrostaticfieldcalculationinmaterialsofdifferentpermittivity,Int.Symp.HighVoltageEngineering,Athens,Greece,Sept.1983,p.11.04.23.J.FaizandM.Ojaghi,Afastboundaryelementmethodtoelectricfieldcomputationwithinthetankofpowertransformers,Int.J.ComputationandMathematicsinElectricalandElectronicEngineering(COMPEL),17(1/2/3),1998,pp.69±77.24.J.FaizandM.Ojaghi,Noveltechniquesfortreatingsingularityproblemsintheboundaryelementmethodofevaluationwithinthetankofapowertransformers,IEEETrans.PowerDelivery,15(2),April2000,pp.592±598.JawadFaizreceivedtheBachelorandMastersdegreesinElectricalEngineeringfromTabrizUniversityinIranin1974and1975respectivelygraduatingwithFirstClassHonours.HereceivedthePh.D.inElectricalEngineeringfromtheUniversityofNewcastle-upon-Tyne,England,in1988.Earlyinhiscareer,heservedasafacultymemberinTabrizUniversityfor10years.AfterobtaininghisPh.D.degreeherejoinedTabrizUniversitywhereheheldthepositionofAssistantProfessorfrom1988to1992,AssociateProfessorfrom1992to1997,andhasbeenaProfessorsince1998.SinceFebruary1999hehasbeenworkingasaProfessoratDepartmentofElectricalandComputerEngineering,FacultyofEngineering,UniversityofTehran.Heistheauthorof50publicationsininternationaljournalsand60publicationsinconferenceproceedings.Dr.FaizisaSeniorMemberofIEEEandmemberofIranAcademyofSciences.HisteachingandresearchinterestsareswitchedreluctanceandVRmotordesign;designandmodellingofelectricalmachinesanddrives.MansourOjaghireceivedtheB.Sc.degreeinElectricalEngineeringfromtheShahidChamranUniversity,Ahwaz,Iranin1993andtheM.Sc.degreeinElectricalPowerEngineeringfromtheUniversityofTabriz,Iran,in1997.HeisnowworkingasaseniorengineerofZanjanRegionalElectricCompany,Zanjan,Iran.Hisareaofinterestisonthemodelingandfieldanalysisoftransformers,protectiverelayingandelectricalmachinesmodeling.JawadFaizandM.Ojaghi356