planar\], 2.1 $\%$ error in the value of $u(0, 2.5)$ is observed on considering only the first two roots i.e., the steady state result and residue for the two non zero poles. The errors arising due to non inclusion of higher order terms is more initially as the higher order poles contribute only at very small times because of the exponential term. The error falls sharply to a negligible value (0.005$\%$) on considering the contribution from the first 6 roots i.e., first 11 poles. More accurate results can be obtained by adding residues from higher order poles.

![\[convergence-planar\]Percentage error in the radiation energy density $u(x,\tau)$ in the slab as a function of number of roots considered (N).](figure10.eps){width="8cm"}

Fig. \[uxt-planar-eps0\] shows the plot of radiation energy density $u(x,\tau)$ as a function of space and time for $\epsilon =0$. Contrary to the results for finite $\epsilon$, the radiation energy density attains a finite value even at very early times due to the absence of retardation effects. However, the material energy density shows the same trend as for finite $\epsilon$.

![\[uxt-planar-eps0\]Scaled radiation energy density $u(x,\tau)$ vs position (x) in the slab of scaled thickness $b=1$ at different times for $ \epsilon = 0$.](figure11.eps){width="8cm"}

Spherical shell {#spherical-shell-1}
---------------

For the spherical shell, initially ($\tau$=0.01) the radiation energy density falls rapidly from the inner surface (scaled radius $X_1=1$) where radiation is incident towards the outer surface (scaled radius $X_2=2$) as shown in Fig. \[uxt-sp\]. Though the trend is similar to the planar slab, the values of the scaled energy densities are less. Also, contrary to the planar case, the variation in energy densities remain sharper in the inner meshes compared to the outer ones and the variations in energy densities are not linear with position even after attaining steady state. This is evident because the mass of the material to be heated in the radially outward direction increases. Similar to the planar slab, the material energy density lags behind the radiation energy densities at early stages and finally reaches equilibrium (
planar\ ], 2.1 $ \%$ error in the value of $ u(0, 2.5)$ is observed on considering only the inaugural two root i.e., the steady state result and remainder for the two non zero poles. The errors arising due to non inclusion body of high order terminus is more initially as the higher order pole contribute only at very small times because of the exponential term. The error hang sharply to a negligible value (0.005$\%$) on considering the contribution from the first 6 root i.e., first 11 poles. More accurate results can be obtained by add residues from higher order poles. 

! [ \[convergence - planar\]Percentage error in the radiation sickness energy density $ u(x,\tau)$ in the slab as a function of number of roots considered (N).](figure10.eps){width="8 cm " } 

 Fig. \[uxt - planar - eps0\ ] shows the plot of radiation energy density $ u(x,\tau)$ as a function of space and time for $ \epsilon = 0$. Contrary to the results for finite $ \epsilon$, the radiation energy concentration achieve a finite value even at very early times due to the absence of slowdown effects. However, the material energy density indicate the same trend as for finite $ \epsilon$. 

! [ \[uxt - planar - eps0\]Scaled radiation energy concentration $ u(x,\tau)$ vs position (x) in the slab of scaled thickness $ b=1 $ at different times for $ \epsilon = 0$.](figure11.eps){width="8 cm " } 

 Spherical shell { # spherical - shell-1 } 
 --------------- 

 For the spherical shell, initially ($ \tau$=0.01) the radiation energy density fall rapidly from the inside open (scaled radius $ X_1=1 $) where radiation is incident towards the outer open (scaled radius $ X_2=2 $) as shown in Fig. \[uxt - sp\ ]. Though the drift is similar to the planar slab, the values of the scaled energy concentration are less. Also, adverse to the planar case, the variation in energy densities remain sharp in the inner mesh compared to the out ones and the variations in energy densities are not linear with military position even after attaining firm state. This is evident because the mass of the material to be heated in the radially outbound direction increases. alike to the planar slab, the material energy density lags behind the radiation energy densities at early stage and finally reaches equilibrium (
plajar\], 2.1 $\%$ error in the value of $u(0, 2.5)$ is obserrwd on ronsideding onlh the first two roots i.e., the sready state result and resiaue for tje two nin ztro poles. The errors arislug dus to uoi inclusion of migher ordes terms is mora knntially as the higher order poles cogtributr lnly at very siall eimea because of the exponential term. Fhe errmr falls sharlly to a negligible value (0.005$\%$) on fonsidering the cojtribution drom rhe first 6 ruots i.e., figvt 11 poles. Jore accurate results can be obgaineb by adding rwsifoes from higier orqer poles.

![\[conyvrgence-[lanar\]Prrcentage errov in vhe eadiation energy densmty $u(x,\tau)$ in the slaf as a futccion of number of roors conshderad (N).](wugufe10.eks){wmdtg="8cm"}

Fig. \[uxv-planar-eps0\] ahows the poot of radiation entrgr density $u(x,\tah)$ as a ftnction of space and time for $\epsilon =0$. Bontdary to the results for finite $\epsilon$, the rwdiation qnergy density attains a finite value even at verf earmh tnnes djw ho the absence of retardation effects. However, fht msterial energy density shpwd yre same trend as fox fjnite $\epsilon$.

![\[uxt-plwnar-eps0\]Fcalee radiatijn emergy density $u(x,\tau)$ vs posirion (x) in thv slqb of scaled thickuess $b=1$ at diyferenj timex for $ \epsilon = 0$.](figure11.e's){widtg="8cm"}

Sphericap shell {#slferical-shell-1}
---------------

For ghe s[herical shell, initially ($\twu$=0.01) the raviatipn enerey dgnsity salls rapifly fvmm the inner surfafe (sccled sadius $X_1=1$) wjere radiation is incident towacvs the outer xusfabe (scaled radims $X_2=2$) as shown ig Fig. \[uxt-sp\]. Thpugh thz trena is similzr to tie planar slwb, the values lf the scalev energy qensutiew are ldrs. Also, contraty to the planar cawe, the variation ik enetgg densities remcnn sharper in the inver mvshxs coi[ared to the outdr umes avd the varictkons in energy densities are not linear with ppsltion eveb after wttaining stesdy state. This is tvidenv becanse thr mwss of the material to be heatsd in the raqially outwawd dlrecjion increaxes. Similar to the planar slab, the matecial energy density lagw behind the radiatnok energy denxitief at earlf stages and finally reaches equilibrlum (
planar\], 2.1 $\%$ error in the value 2.5)$ observed on only the first state and residue for two non zero The errors arising due to non of higher order terms is more initially as the higher order poles contribute at very small times because of the exponential term. The error falls sharply a value on the contribution from the first 6 roots i.e., first 11 poles. More accurate results can be by adding residues from higher order poles. ![\[convergence-planar\]Percentage in the radiation energy $u(x,\tau)$ in the slab as function number of considered Fig. shows the plot radiation energy density $u(x,\tau)$ as a function of space and time for $\epsilon =0$. Contrary to the for finite radiation energy attains finite even at very due to the absence of retardation material energy density shows the same trend as finite $\epsilon$. radiation energy density $u(x,\tau)$ vs position in the slab of scaled thickness $b=1$ at times for $ \epsilon = 0$.](figure11.eps){width="8cm"} Spherical shell {#spherical-shell-1} --------------- For the spherical shell, initially radiation energy density falls from the inner (scaled $X_1=1$) radiation incident towards outer surface (scaled radius $X_2=2$) as shown in Fig. \[uxt-sp\]. Though trend is similar to the planar slab, the values of energy are less. Also, to the planar case, variation energy densities remain sharper inner compared ones the in energy densities are linear with position even after steady state. This is the material to be heated in the radially direction increases. Similar to the planar slab, material energy density lags behind the radiation energy densities at early stages finally reaches
planar\], 2.1 $\%$ error in the value of $u(0, 2.5)$ iS observed oN consIdeRinG oNly tHe fiRst two roots i.e., tHE steAdy state result and residUe for ThE Two nON zEro poLes. The eRRoRS AriSiNg Due To NOn IncluSioN of highEr order terMs iS mOre initially AS tHe higher orDer Poles contribUte Only at VeRy sMAll tiMes BecauSe of thE ExponeNtial term. thE Error fALls sharPLY tO a neGligible value (0.005$\%$) on coNSiDEring the contriBution FrOM tHE FirSt 6 rOots i.e., firsT 11 pOles. MORe accurATe RESUltS Can be obtained By adding resIDueS from hIgHer ORder poLes.

![\[coNvERgeNce-planar\]PeRcenTage error In the rADiation ENergy deNsity $u(X,\taU)$ in The sLAb As A fuNcTIon OF nUmbER of Roots conSiDeRed (N).](fIgurE10.EPS){WidtH="8cm"}

fig. \[uXt-plaNar-eps0\] shows thE plOt of RAdiAtion EnergY denSiTy $u(x,\tAu)$ as a fUnctiOn Of space and time fOr $\epSilon =0$. ContRarY tO thE rEsultS For finIte $\EpsIlon$, the RadiatiON enErGY DEnSity attains a finite VaLUE eVen at verY early TImEs DUe to the aBsEncE of rETArdatIon eFFeCts. HowevEr, the mATeRiAl energY dEnsity ShOws The Same tREnd aS for fiNite $\epsiLon$.

![\[uxT-Planar-eps0\]ScaleD Radiation enerGY dENSiTY $u(x,\tAu)$ vS position (x) iN the SLab oF scaLEd ThiCKness $B=1$ at diFfEReNT times for $ \epsilon = 0$.](figUrE11.eps){wiDth="8cm"}

spherical shelL {#spherical-SHELl-1}
---------------

For the SpheRIcAL shell, initiallY ($\tau$=0.01) tHe radiatioN Energy deNsity Falls rapIdly from tHE Inner surFacE (scAleD raDIUs $x_1=1$) where radiatiON Is inCiDent towArdS the outEr sUrfAce (ScaLeD radius $X_2=2$) aS shown in fiG. \[uXt-Sp\]. thoUgh thE Trend is sImIlaR tO thE planAR slab, tHe valUes oF tHe SCalEd energY DeNSItieS aRe Less. alsO, cOntraRy to THe pLanar caSe, the variAtiON in eNeRgY densitIes remain sharPeR in the inneR mEshEs compAREd to the oUter ones and the variationS In energY deNsitiEs arE not lineaR wiTh posiTioN Even afTer attAininG sTeaDY State. tHIs Is eViDent becausE THe mAss of ThE matErial to Be heated in the radiaLLy oUtward directiOn iNcreASEs. simILaR To tHe PLanAR Slab, the material Energy densItY LaGs behind thE RadIaTion eneRgy densIties AT early sTages and fInally reaChEs eqUILibRium (
planar\], 2.1 $\%$ error i n the valu e of$u( 0,2. 5)$is o bserved on con s ider ing only the first two  root si .e., th e ste ady sta t er e sul tan d r es i du e for  th e two n on zero po les .The errors a r is ing due to  no n inclusionofhigher  o rde r  term s i s mor e init i ally a s the hig he r  order poles c o n tr ibut e only at very sm a ll times becauseof the  e x po n e nti alterm. Theer ror f a lls sha r pl y t o a negligible va lue (0.005$ \ %$)  on co ns ide r ing th e con tr i but ion from th e fi rst 6 roo ts i.e . , first 11 pole s. Mor e a ccu rate re su lts  c a n b e  o bta i ned  by addi ng  r esidu es f r o m high erorde r pol es.

![\[conv erg ence - pla nar\] Perce ntag eerror  in th e rad ia tion energy den sity  $u(x,\ta u)$  i n t he  slab as a f unc tio n of nu mber of roo ts c o ns idered (N).](figur e1 0 . ep s){width ="8cm" } 
Fi g . \[uxt- pl ana r-ep s 0 \] sh owst he  plot of  radia t io nenergyde nsity$u (x, \ta u)$ a s  a f unctio n of spa ce an d  time for $\ep s ilon =0$. Con t ra r y  t o  the  re sults for f init e  $\e psil o n$ , t h e rad iatio ne ne r gy density attainsafinite  valu e even at ver y early ti m e s  due tothea bs e nce of retarda tioneffects. H o wever, t he ma terial e nergy den s i ty shows  th e s ame  tr e n das for finite $ \eps il on$.

! [\[ uxt-pla nar -ep s0\ ]Sc al ed radiat ion ener gy  d en si ty$u(x, \ tau)$ vs  p osi ti on(x) i n  the s lab o f sc al ed thi ckness$ b= 1 $  atdi ff eren t t im es fo r $\ eps ilon =0$.](figu re1 1 .eps ){ wi dth="8c m"}

Spherica lshell {#sp he ric al-she l l -1}
---- -----------

For the sp h ericalshe ll, i niti ally ($\t au$ =0.01)  th e  radia tion e nergy  d ens i t y fal l s  r api dl y from the i nne r sur fa ce ( scaledradius $X_1=1$) wh e reradiation isinc iden t to war d st heou t ers u rface (scaled r adius $X_2 =2 $ )as shown i n  Fi g.  \[uxt- sp\]. T hough the tre nd is sim ilar to t he  pla n a r s lab, the v alues of  the scal e d ene r gy  dens iti es are  l ess . Als o, con t rar y tothe pl an ar cas e, th evariatio n in energy densities r emainsharp erin the in ner mes hes compa redto the out erone s and  th e  vari atio n sine nergy  den s ities are no t l i n ea r with posi t i o n e ven a fte r  attai ning  steady state. Th i s is evident b ecau s e  th e m a ss o fthe material t o b eh e ated inth e radiallyoutwarddi r ectio n incr eases.  Simila r to the pl anar  sl ab, the m ate ri a l energ yde n sity l agsbe hind t he rad i atio n energy densities  at e a r ly st a ges  andfi nally r e ache s equilibr ium (
planar\], 2.1_$\%$ error_in the value of_$u(0, 2.5)$_is_observed on_considering_only the first_two roots i.e.,_the steady state result_and residue for_the_two non zero poles. The errors arising due to non inclusion of higher order_terms_is more_initially_as_the higher order poles contribute_only at very small times_because of_the exponential term. The error falls sharply to_a_negligible value (0.005$\%$)_on considering the contribution from the first 6 roots_i.e., first 11 poles. More accurate_results can be_obtained_by_adding residues from higher_order poles.

![\[convergence-planar\]Percentage error in the radiation_energy density $u(x,\tau)$ in the slab_as a function of number of roots_considered (N).](figure10.eps){width="8cm"}

Fig. \[uxt-planar-eps0\] shows the plot_of radiation energy density $u(x,\tau)$_as a_function of space and time_for $\epsilon =0$._Contrary to_the results for_finite $\epsilon$, the radiation energy density_attains a finite_value even at very early times_due_to the absence_of_retardation_effects. However,_the material energy_density_shows the_same_trend as for finite $\epsilon$.

![\[uxt-planar-eps0\]Scaled radiation_energy_density $u(x,\tau)$ vs position (x) in the_slab of scaled thickness_$b=1$_at different times for_$ \epsilon = 0$.](figure11.eps){width="8cm"}

Spherical shell_{#spherical-shell-1}
---------------

For the spherical shell, initially ($\tau$=0.01)_the radiation_energy density_falls rapidly from the inner surface (scaled radius $X_1=1$) where radiation_is incident towards the outer surface_(scaled radius $X_2=2$) as_shown in_Fig._\[uxt-sp\]. Though the_trend_is similar_to the planar slab, the values of_the scaled_energy densities are less. Also, contrary_to the planar case,_the_variation in energy densities remain sharper_in the inner meshes compared to_the outer ones and the_variations_in_energy densities are not linear_with position even after attaining steady_state. This is_evident because the mass of the material_to_be heated in the radially outward_direction_increases. Similar to the planar slab,_the_material_energy density lags behind the_radiation energy densities at early stages_and finally reaches equilibrium (