-1/2}^2} =
            \frac{2(f^2+f+1)}{3(f+1)} f^{i-1}R_{\rm min},$$ for $i=1,2,\cdots,N_R$. Note that the radial cell spacing $\delta R_i \equiv R_{i+1} - R_i = (f-1)R_i$ increases with $R_i$.

The second-order finite-difference approximation to Equation takes the same form as Equation, but with the radial difference operator defined as $$\label{eq:radial_operator}
  \Delta_R^2\Phi_{i,j,k} = \frac{\Phi_{i-1,j,k}-2\Phi_{i,j,k}+\Phi_{i+1,j,k}}{(R_i\ln f)^2}.$$ Appendix \[s:fd\] shows that Equation makes the finite difference approximation second-order accurate.

Calculation of Interior Potential for Dirichlet Boundary Conditions {#s:interior_solver}
===================================================================

In this section, we provide the general method that we use to obtain the interior potential within the original domain, given Dirichlet boundaries for the potential on the surface. We describe our methods for both Cartesian grids (Section \[sec:cart\_int\]) and cylindrical grids (Section \[s:interior\_solver\_cylindrical\]). We note that in principle, alternative fast and efficient solvers (such as multigrid) could be employed for computing the interior potential given Dirichlet boundary conditions for the potential.

The interior solver is employed in three different instances in our method. The first instance is our use of the interior solver to obtain the numerical DGF, as described in Appendix \[s:calc\_dgf\]; this is done once at the beginning of any simulation. The other two instances are in the first and fourth step of the James’s method; each time the Poisson solution is required, two calls to the interior solver are made.

In the first step of James’s method, the interior potential from the original density distribution is computed subject to the zero boundary condition. In the fourth step of James’s method, the interior potential must be computed subject to a known boundary potential $\Phi^{\rm B}$; a method to obtain $\Phi^{\rm B}$ will be presented in Section \[s:bc\]. Formally, we define $\
-1/2}^2 } = 
             \frac{2(f^2+f+1)}{3(f+1) } f^{i-1}R_{\rm min},$$ for $ i=1,2,\cdots, N_R$. Note that the radial cell spacing $ \delta R_i \equiv R_{i+1 } - R_i = (f-1)R_i$ increases with $ R_i$. 

 The second - order finite - remainder estimate to Equation takes the same human body as equality, but with the radial difference operator defined as $ $ \label{eq: radial_operator } 
   \Delta_R^2\Phi_{i, j, k } = \frac{\Phi_{i-1,j, k}-2\Phi_{i, j, k}+\Phi_{i+1,j, k}}{(R_i\ln f)^2}.$$ Appendix \[s: fd\ ] picture that Equation makes the finite difference estimate second - decree accurate. 

 Calculation of Interior Potential for Dirichlet Boundary Conditions { # s: interior_solver } 
 = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = 

 In this incision, we provide the general method acting that we use to obtain the home potential within the original domain, contribute Dirichlet boundaries for the potential on the surface. We identify our methods for both Cartesian grids (Section \[sec: cart\_int\ ]) and cylindrical grid (Section \[s: interior\_solver\_cylindrical\ ]). We note that in principle, alternative fast and efficient problem solver (such as multigrid) could be employed for computing the interior potential given Dirichlet boundary conditions for the electric potential. 

 The home solver is employed in three different instances in our method. The first instance is our use of the interior solver to obtain the numerical DGF, as describe in Appendix \[s: calc\_dgf\ ]; this is done once at the beginning of any simulation. The early two case are in the beginning and fourth step of the James ’s method; each meter the Poisson solution is required, two calls to the interior solver are lay down. 

 In the first step of James ’s method, the interior potential from the original density distribution is computed subject to the zero boundary condition. In the fourth step of James ’s method, the home potential must be calculate capable to a know boundary potential $ \Phi^{\rm B}$; a method to obtain $ \Phi^{\rm B}$ will be give in Section \[s: bc\ ]. Formally, we define $ \
-1/2}^2} =
            \frac{2(f^2+f+1)}{3(f+1)} f^{i-1}R_{\rm oin},$$ for $i=1,2,\cdots,N_T$. Bote tiat the radial zell spacing $\delta R_i \equiv C_{i+1} - E_i = (f-1)E_i$ increases with $R_i$.

Thd second-ogder finire-dihference approximation to Equaflon tckxs the same fork as Equathon, but with tve rcdial difference operator defined as $$\label{ea:rwdial_operator}
  \Delua_R^2\[hi_{i,n,k} = \frac{\Phi_{i-1,j,k}-2\Phi_{i,j,k}+\Phi_{i+1,j,k}}{(R_i\ln f)^2}.$$ Alpendix \[s:fd\] shows thst Equation makes the finihe dlfference approximwtion secone-ordqe accurate.

Cauculation of Interior Kotential for Dirichlet Boundary Condntions {#s:intgxuor_dmlver}
===================================================================

In thiw secnion, we provide the ganeral kethod that we usx to obtain the interior 'otential within the original dkmain, given Diricylwt boondarhes wir ghe pptsntial on the surfade. We descrube our methods for bjnn Cartesian frids (Festion \[sec:cart\_int\]) and cylindrical grids (Vecfion \[s:interior\_solver\_cylundrical\]). We note that in princyple, alternative fast and efficient solvers (such ds mumgigxlq) xokld be employed for computing the interior poesnuiak given Dirichjet boundaru fomqitions for tfe potzhtjal.

The interior sopver is empliyed in trree different instances in our method. The yirwt instance is our use of the intgrior xolver to obtain the nuoeridal DGF, as fescribed kn Appendix \[s:calz\_dgn\]; tvis is done once at the bedinning oh any simulagion. The oeher two ijstangas are in the firsh and fmurth step of the James’s method; each time the Poisson splgtimn is rezuirec, two calls tj the interior solver are oade.

In the first vtep of Jamqs’s method, tha interior povential fwom rhe iriginau density distrobution is computed subject to the zevo boonsary condition. Nu uhw fourth step pf Gamqs’d keevod, the intesior pogrntiau must be cioputrd subject to a knowt bohndary potential $\Pni^{\vm B}$; a mejhod to oftain $\Phi^{\rm B}$ will be presented in Sxction \[s:bc\]. Gorially, we define $\
-1/2}^2} = \frac{2(f^2+f+1)}{3(f+1)} f^{i-1}R_{\rm min},$$ for $i=1,2,\cdots,N_R$. the cell spacing R_i \equiv R_{i+1} with The second-order finite-difference to Equation takes same form as Equation, but with radial difference operator defined as $$\label{eq:radial_operator} \Delta_R^2\Phi_{i,j,k} = \frac{\Phi_{i-1,j,k}-2\Phi_{i,j,k}+\Phi_{i+1,j,k}}{(R_i\ln f)^2}.$$ Appendix \[s:fd\] shows Equation makes the finite difference approximation second-order accurate. Calculation of Interior Potential for Boundary {#s:interior_solver} In section, we provide the general method that we use to obtain the interior potential within the domain, given Dirichlet boundaries for the potential on surface. We describe our for both Cartesian grids (Section and grids (Section We that principle, alternative fast efficient solvers (such as multigrid) could be employed for computing the interior potential given Dirichlet boundary conditions the potential. solver is in different in our method. instance is our use of the obtain the numerical DGF, as described in Appendix this is once at the beginning of any The other two instances are in the first fourth step of the James’s method; each time the Poisson solution is required, two calls interior solver are made. the first step James’s the potential the original distribution is computed subject to the zero boundary condition. In the step of James’s method, the interior potential must be computed a boundary potential $\Phi^{\rm a method to obtain B}$ be presented in Section we $\
-1/2}^2} =
            \frac{2(f^2+f+1)}{3(f+1)} f^{i-1}R_{\rm min},$$ for $i=1,2,\cdots,N_r$. Note that tHe radIal CelL sPaciNg $\deLta R_i \equiv R_{i+1} - R_i = (F-1)r_i$ inCreases with $R_i$.

The second-Order FiNIte-dIFfErencE approxIMaTIOn tO EQuAtiOn TAkEs the SamE form as equation, buT wiTh The radial difFErEnce operatOr dEfined as $$\labeL{eq:Radial_OpEraTOr}
  \DelTa_R^2\phi_{i,j,K} = \frac{\PHI_{i-1,j,k}-2\PhI_{i,j,k}+\Phi_{i+1,j,K}}{(R_I\Ln f)^2}.$$ AppENdix \[s:fd\] SHOwS thaT Equation makes the FInITe difference apProximAtIOn SEConD-orDer accuratE.

CAlculATion of INTeRIOR PoTEntial for DiriChlet BoundaRY CoNditioNs {#S:inTErior_sOlver}
===================================================================

in THis Section, we prOvidE the generAl methOD that we USe to obtAin the IntEriOr poTEnTiAl wItHIn tHE oRigINal Domain, giVeN DIrichLet bOUNDArieS foR the PotenTial on the surfAce. we deSCriBe our MethoDs foR bOth CaRtesiaN gridS (SEction \[sec:cart\_inT\]) and CylindricAl gRiDs (SEcTion \[s:INterioR\_soLveR\_cylindRical\]). We NOte ThAT IN pRinciple, alternativE fAST aNd efficiEnt solVErS (sUCh as multIgRid) CoulD BE emplOyed FOr ComputinG the inTErIoR potentIaL given diRicHleT bounDAry cOnditiOns for thE poteNTial.

The interioR Solver is emploYEd IN ThREe diFfeRent instancEs in OUr meThod. tHe FirST instAnce iS oUR uSE of the interior solveR tO obtaiN the nUmerical DGF, as Described iN aPPendix \[s:cAlc\_dGF\]; tHIs is done once at The beGinning of aNY simulatIon. ThE other twO instanceS ARe in the fIrsT anD foUrtH STeP of the James’s mETHod; eAcH time thE PoIsson soLutIon Is rEquIrEd, two callS to the inTeRiOr SoLveR are mADe.

In the fIrSt sTeP of james’S Method, The inTeriOr PoTEntIal from THe ORIginAl DeNsitY diStRibutIon iS ComPuted suBject to thE zeRO bouNdArY conditIon. In the fourtH sTep of James’S mEthOd, the iNTErior potEntial must be computed subJEct to a kNowN bounDary Potential $\phi^{\Rm B}$; a meThoD To obtaIn $\Phi^{\rM B}$ wilL bE prESEnted IN seCtiOn \[S:bc\]. FormallY, WE deFine $\
-1/2}^2} =
            \fr ac{2(f^2+f +1)}{ 3(f +1) }f^{i -1}R _{\rm min},$$f or $ i=1,2,\cdots,N_R$. Not e tha tt he r a di al ce ll spac i ng $ \de lt aR_i  \ e qu iv R_ {i+ 1} - R_ i = (f-1)R _i$  i ncreases wit h  $ R_i$.

The  se cond-order f ini te-dif fe ren c e app rox imati on toE quatio n takes t he same f o rm as E q u at ion,  but with the rad i al difference ope ratorde f in e d  as  $$ \label{eq: ra dial_ o perator } 
\ D elt a _R^2\Phi_{i,j ,k} = \frac { \Ph i_{i-1 ,j ,k} - 2\Phi_ {i,j, k} + \Ph i_{i+1,j,k} }{(R _i\ln f)^ 2}.$$A ppendix \[s:fd\ ] show s t hat  Equ a ti on  ma ke s  th e  f ini t e d ifferenc eap proxi mati o n s econ d-o rder  accu rate.

Calcul ati on o f  In terio r Pot enti al  forDirich let B ou ndary Condition s {# s:interio r_s ol ver }===== = ====== === === ======= ======= = === == = = = == ================== == = = == 

In thi s sect i on ,w e provid ethe  gen e r al me thod th at we us e to o b ta in  the in te rior p ot ent ial  with i n th e orig inal dom ain,g iven Dirichlet boundaries fo r  t h e  p o tent ial  on the sur face .  Wedesc r ib e o u r met hodsfo r  b o th Cartesian grids(S ection  \[se c:cart\_int\] ) and cyli n d r ical gri ds ( S ec t ion \[s:interi or\_s olver\_cyl i ndrical\ ]). W e note t hat in pr i n ciple, a lte rna tiv e f a s tand efficient s olve rs  (suchasmultigr id)  co uld  be  e mployed f or compu ti ng  t he  in terio r  potenti al  gi ve n D irich l et bou ndary  con di ti o nsfor the po t e ntia l. 
Theint er ior s olve r  is  employ ed in thr eed iffe re nt  instan ces in our me th od. The fi rs t i nstanc e is our u se of the interior solv e r to ob tai n the  num erical DG F,as des cri b ed inAppend ix \[ s: cal c \ _dgf\ ] ;  t his  i s done onc e atthe b eg inni ng of a ny simulation. The oth er two instan ces  are i nthe fi r stan d  fo u r th step of theJames’s me th o d;  each time the  P oissonsolutio n isr equired , two cal ls to the  i nter i o r s olver aremade.

I n the fir s t ste p  o f Jam es’ s meth od , t he in terior pot entia l from  t he ori ginal  d ensity d istribution is computed  subje ct to  th e zero bo und a rycondition . In  the fourt h s tep  of J ame s ’s me thod ,  t hei nteri or p o tential m u st  be c om puted subje c t toa kno wnb oundar y po tential $\Phi^{\r m  B}$; a method  too b tai n $ \ Phi^ {\ rm B}$ will be  pr es e n ted in S ec tion \[s:bc \]. Form al l y, we  defin e $\
-1/2}^2} =
_ _   _ __ _\frac{2(f^2+f+1)}{3(f+1)}_f^{i-1}R_{\rm min},$$ for_$i=1,2,\cdots,N_R$. Note that_the radial cell spacing_$\delta R_i \equiv_R_{i+1}_- R_i = (f-1)R_i$ increases with $R_i$.

The second-order finite-difference approximation to Equation takes the_same_form as_Equation,_but_with the radial difference operator_defined as $$\label{eq:radial_operator}
  \Delta_R^2\Phi_{i,j,k}_= \frac{\Phi_{i-1,j,k}-2\Phi_{i,j,k}+\Phi_{i+1,j,k}}{(R_i\ln_f)^2}.$$ Appendix \[s:fd\] shows that Equation makes the_finite_difference approximation second-order_accurate.

Calculation of Interior Potential for Dirichlet Boundary Conditions {#s:interior_solver}
===================================================================

In_this section, we provide the general_method that we_use_to_obtain the interior potential_within the original domain, given Dirichlet_boundaries for the potential on the_surface. We describe our methods for both_Cartesian grids (Section \[sec:cart\_int\]) and cylindrical_grids (Section \[s:interior\_solver\_cylindrical\]). We note_that in_principle, alternative fast and efficient_solvers (such as_multigrid) could_be employed for_computing the interior potential given Dirichlet_boundary conditions for_the potential.

The interior solver is employed_in_three different instances_in_our_method. The_first instance is_our_use of_the_interior solver to obtain the numerical_DGF,_as described in Appendix \[s:calc\_dgf\]; this is_done once at the_beginning_of any simulation. The_other two instances are in_the first and fourth step of_the James’s_method; each_time the Poisson solution is required, two calls to the interior_solver are made.

In the first step_of James’s method, the_interior potential_from_the original density_distribution_is computed_subject to the zero boundary condition. In_the fourth_step of James’s method, the interior_potential must be computed_subject_to a known boundary potential $\Phi^{\rm_B}$; a method to obtain $\Phi^{\rm_B}$ will be presented in_Section_\[s:bc\]._Formally, we define $\