 ${\text{\textit{St}}}_\theta$, and for simulations S1 and S2. The results show that as ${\text{\textit{St}}}_\theta$ is increased, $\chi_p$ increases. This is mainly because as ${\text{\textit{St}}}_\theta$ is increased, the thermal memory of the particle increases, and the particle temperature depends strongly on its encounter with the fluid temperature field along its trajectory history for times up to ${\textit{O}\left( \tau_\theta \right)}$ in the past. As a result, the particle temperature can differ strongly from the local fluid temperature. The results also show that $\chi_p$ is dramatically suppressed when two-way thermal coupling is accounted for. One reason for this is that as shown earlier, two-way thermal coupling leads to a suppression in the fluid temperature gradients. As these gradients are suppressed, the fluid temperature along the particle trajectory history differs less from the local fluid temperature than it would have in the absence of two-way thermal coupling, and as a result $\chi_p$ is decreased.

The results for kurtosis of $\dot{\theta}_p$, as a function of ${\text{\textit{St}}}_\theta$ and for various ${\text{\textit{St}}}$ are shown in figure \[fig:PDFdeT\](f). The results show that the kurtosis decreases with increasing ${\text{\textit{St}}}_\theta$. This is mainly due to the filtering effect mentioned earlier, wherein as ${\text{\textit{St}}}_\theta$ is increased, the particles are less able to respond to rapid fluctuations in the fluid temperature along their trajectory. Further, the kurtosis is typically larger when the two-way thermal coupling is taken into account (simulations S1), and is maximum for ${\text{\textit{St}}}=1$. This is due to the particle clustering on the fronts of the fluid temperature field, as will be discussed in section \[sec:FATT\].

Our results for the PDF of $\dot{\theta}_p$ and its moments differ somewhat from those in [@Bec2014]. This is in part due to the difference in the forcing methods employed by [@Bec2014] and that in our study. The solution of may be written as [@Bec2014] $${\left\langle \dot{\theta_p}^2 \right\rangle} = \frac{1}{2\tau_\theta^3}\int_0^\infty {\left\langle \
$ { \text{\textit{St}}}_\theta$, and for simulations S1 and S2. The results show that as $ { \text{\textit{St}}}_\theta$ is increased, $ \chi_p$ increase. This is chiefly because as $ { \text{\textit{St}}}_\theta$ is increased, the thermal memory of the particle increase, and the atom temperature depends strongly on its encounter with the fluent temperature playing field along its trajectory history for times up to $ { \textit{O}\left (\tau_\theta \right)}$ in the past. As a consequence, the particle temperature can differ powerfully from the local fluid temperature. The results also show that $ \chi_p$ is dramatically suppress when two - way thermal yoke is accounted for. One rationality for this is that as read earlier, two - way thermal coupling lead to a suppression in the fluid temperature gradients. As these gradients are suppress, the fluid temperature along the particle trajectory history differs less from the local fluid temperature than it would have in the absence of two - direction thermal coupling, and as a result $ \chi_p$ is decreased. 

 The results for kurtosis of $ \dot{\theta}_p$, as a function of $ { \text{\textit{St}}}_\theta$ and for various $ { \text{\textit{St}}}$ are shown in figure \[fig: PDFdeT\](f). The results show that the kurtosis decreases with increasing $ { \text{\textit{St}}}_\theta$. This is mainly due to the filtering effect mentioned sooner, wherein as $ { \text{\textit{St}}}_\theta$ is increased, the particles are less able to answer to rapid fluctuation in the fluid temperature along their trajectory. Further, the kurtosis is typically larger when the two - direction thermal coupling is taken into account (simulations S1), and is maximal for $ { \text{\textit{St}}}=1$. This is due to the particle clustering on the fronts of the fluid temperature field, as will be discussed in section \[sec: FATT\ ]. 

 Our results for the PDF of $ \dot{\theta}_p$ and its moments differ somewhat from those in [ @Bec2014 ]. This is in function due to the deviation in the forcing methods employed by [ @Bec2014 ] and that in our study. The solution of may be written as [ @Bec2014 ] $ $ { \left\langle \dot{\theta_p}^2 \right\rangle } = \frac{1}{2\tau_\theta^3}\int_0^\infty { \left\langle \
 ${\tedt{\textit{St}}}_\theta$, and for rimulations S1 aue S2. Thx resulfs show ghat as ${\text{\textit{St}}}_\theta$ is ibcreawed, $\chi_p$ increases. Thir is mainpy becauwe aw ${\text{\textiv{At}}}_\theta$ is indveaseb, vhe thermal mempry of the particle incraared, and the particle temperature depegds strpnhly on its encjunttr rith nht fluid temperature field along ifs trajtctory history for times up to ${\textit{O}\left( \twu_\thfta \right)}$ in the pwst. As a rewult, rhe particle temperatugz can diffet strongly from the local fluid gempexature. The tzwulhv also show that $\chi_p$ is dramaticallf supprrssed when two-eay thwrmal coupling is accmunted for. One reafon for tvia is that as showb warligr, twm-way rhefmam roulling peavs to a suplression in the fluid temperatirq gradients. As these gwadients are suppressed, the fluid tempegatude along the particle teajectory history difvers less from the local fluid temperature than it would hdve ii ghe ansenzw lf two-way thermal coupling, and as a result $\chy_l$ os decreased.

The vesults for kurtosos og $\dot{\theta}_p$, as a funefikn of ${\text{\textit{St}}}_\hheta$ agd foe various ${\texy{\textit{St}}}$ are shown in figuee \[fig:PDFdeT\](y). Tye results show thct the kurtoris cecresses with increasing ${\tert{\textjt{St}}}_\theta$. Tjis is majvly due to the fkltvrinc effect mentioned earlier, wherein es ${\tert{\textit{Rt}}}_\thgta$ is yncreased, hhe particles are less ahle tl sespond to rapid fluctuations in the fluiv temperature anonc their crajecbory. Further, thq kurtosis is jypically largdr when thv two-way vhermal coupjing is taken lnto account (simulatyons S1), abd is mxbimum for ${\text{\yextit{St}}}=1$. Nhns is due to the particle ckusjeding on the frouus if the fluid tekpefattrv fmeld, wv will be divcusred on seztion \[ste:FABT\].

Ojr rrsults for the PDF ox $\dof{\theta}_p$ and its mokekts diffet somewhae from those on [@Bec2014]. This is in kart dne to vhe digfetence in the forcing methods emlloyed by [@Beg2014] and that in our study. The splution of may be written as [@Bec2014] $${\left\laigle \dot{\theta_p}^2 \right\rantle} = \frac{1}{2\tau_\theta^3}\inj_0^\innty {\left\langke \
${\text{\textit{St}}}_\theta$, and for simulations S1 and S2. show as ${\text{\textit{St}}}_\theta$ increased, $\chi_p$ increases. ${\text{\textit{St}}}_\theta$ increased, the thermal of the particle and the particle temperature depends strongly its encounter with the fluid temperature field along its trajectory history for times to ${\textit{O}\left( \tau_\theta \right)}$ in the past. As a result, the particle temperature differ from local temperature. The results also show that $\chi_p$ is dramatically suppressed when two-way thermal coupling is accounted One reason for this is that as shown two-way thermal coupling leads a suppression in the fluid gradients. these gradients suppressed, fluid along the particle history differs less from the local fluid temperature than it would have in the absence of two-way coupling, and result $\chi_p$ decreased. results kurtosis of $\dot{\theta}_p$, function of ${\text{\textit{St}}}_\theta$ and for various in figure \[fig:PDFdeT\](f). The results show that the decreases with ${\text{\textit{St}}}_\theta$. This is mainly due to filtering effect mentioned earlier, wherein as ${\text{\textit{St}}}_\theta$ is the particles are less able to respond to rapid fluctuations in the fluid temperature along Further, the kurtosis is larger when the thermal is into (simulations S1), is maximum for ${\text{\textit{St}}}=1$. This is due to the particle clustering the fronts of the fluid temperature field, as will be section Our results for PDF of $\dot{\theta}_p$ and moments somewhat from those in is part difference the methods employed by [@Bec2014] that in our study. The of may be written = \frac{1}{2\tau_\theta^3}\int_0^\infty {\left\langle \
 ${\text{\textit{St}}}_\theta$, and for siMulations S1 And S2. THe rEsuLtS shoW thaT as ${\text{\textit{ST}}}_\ThetA$ is increased, $\chi_p$ increaSes. ThIs IS maiNLy BecauSe as ${\texT{\TeXTIt{ST}}}_\tHeTa$ iS iNCrEased, The Thermal Memory of thE paRtIcle increaseS, AnD the particLe tEmperature dePenDs stroNgLy oN Its enCouNter wIth the FLuid teMperature FiELd alonG Its trajECToRy hiStory for times up to ${\TExTIt{O}\left( \tau_\thetA \right)}$ In THe PASt. AS a rEsult, the paRtIcle tEMperatuRE cAN DIffER strongly from The local fluID teMperatUrE. ThE ResultS also ShOW thAt $\chi_p$ is draMatiCally suppRessed WHen two-wAY thermaL couplIng Is aCcouNTeD fOr. ONe REasON fOr tHIs iS that as sHoWn EarliEr, twO-WAY TherMal CoupLing lEads to a suppreSsiOn in THe fLuid tEmperAturE gRadieNts. As tHese gRaDients are suppreSsed, The fluid tEmpErAtuRe Along THe partIclE trAjectorY historY DifFeRS LEsS from the local fluid TeMPErAture thaN it wouLD hAvE In the absEnCe oF two-WAY therMal cOUpLing, and aS a resuLT $\cHi_P$ is decrEaSed.

The ReSulTs fOr kurTOsis Of $\dot{\tHeta}_p$, as a FunctIOn of ${\text{\textit{sT}}}_\theta$ and for vARiOUS ${\tEXt{\teXtiT{St}}}$ are shown In fiGUre \[fIg:PDfDet\](f). THE resuLts shOw THaT The kurtosis decreaseS wIth incReasiNg ${\text{\textit{ST}}}_\theta$. This IS MAinly due To thE FiLTering effect meNtionEd earlier, wHErein as ${\tExt{\teXtit{St}}}_\thEta$ is incrEASed, the paRtiCleS arE leSS AbLe to respond to RAPid fLuCtuatioNs iN the fluId tEmpEraTurE aLong their TrajectoRy. fuRtHeR, thE kurtOSis is typIcAllY lArgEr wheN The two-Way thErmaL cOuPLinG is takeN InTO AccoUnT (sImulAtiOnS S1), and Is maXImuM for ${\texT{\textit{St}}}=1$. thiS Is duE tO tHe partiCle clustering On The fronts oF tHe fLuid teMPErature fIeld, as will be discussed in SEction \[sEc:FaTT\].

OuR resUlts for thE PDf of $\dot{\TheTA}_p$ and iTs momeNts diFfEr sOMEwhat FROm ThoSe In [@Bec2014]. This iS IN paRt due To The dIfferenCe in the forcing methODs eMployed by [@Bec2014] aNd tHat iN OUr StuDY. THE soLuTIon OF May be written as [@BEc2014] $${\left\langLe \DOt{\Theta_p}^2 \righT\RanGlE} = \frac{1}{2\taU_\theta^3}\iNt_0^\infTY {\left\laNgle \
 ${\text{\textit{St}}}_\th eta$, andfor s imu lat io ns S 1 an d S2. The resu l ts s how that as ${\text{\t extit {S t }}}_ \ th eta$is incr e as e d , $ \c hi _p$  i n cr eases . T his ismainly bec aus eas ${\text{\ t ex tit{St}}}_ \th eta$ is incr eas ed, th ethe r mal m emo ry of  the p a rticle  increase s, and th e  partic l e  t empe rature depends st r on g ly on its enco unterwi t ht h e f lui d temperat ur e fie l d along it s t raj e ctory history  for timesu p t o ${\t ex tit { O}\lef t( \t au _ \th eta \right) }$ i n the pas t. Asa  result ,  the pa rticle  te mpe ratu r eca n d if f ers tr ong l y f rom thelo ca l flu id t e m p e ratu re.  The  resu lts also show  th at $ \ chi _p$ i s dra mati ca lly s uppres sed w he n two-way therm al c oupling i s a cc oun te d for .  One r eas onfor thi s is th a t a ss h o wn  earlier, two-wayth e r ma l coupli ng lea d sto a suppre ss ion  int h e flu id t e mp eraturegradie n ts .As thes egradie nt s a resuppr e ssed , thefluid te mpera t ure along thep article traje c to r y  h i stor y d iffers less  fro m  the  loc a lflu i d tem perat ur e  t h an it would have in  t he abs enceof two-way th ermal coup l i n g, and a s ar es u lt $\chi_p$ is  decr eased.

Th e  results  forkurtosis  of $\dot { \ theta}_p $,asa f unc t i on  of ${\text{\ t e xtit {S t}}}_\t het a$ andfor  va rio us${ \text{\te xtit{St} }} $ar esho wn in figure \ [f ig: PD Fde T\](f ) . Theresul ts s ho wt hat  the ku r to s i s de cr ea seswit hincre asin g  ${ \text{\ textit{St }}} _ \the ta $.  This i s mainly dueto  the filte ri ngeffect m entioned  earlier, wherein as ${ \ text{\t ext it{St }}}_ \theta$ i s i ncreas ed, the pa rticle s are  l ess a ble t o re spo nd  to rapidf l uct uatio ns  inthe flu id temperature alo n g t heir trajecto ry.  Fur t h er , t h ek urt os i s i s typically large r when the  t w o- way therma l  co up ling is  takenintoa ccount(simulati ons S1),an d is m axi mum for ${ \text{\t extit{St} } }=1$. Th is is  du e to t he  pa rticl e clus t eri ng on  the f ro nts of  thefl uid temp erature field, as willbe dis cusse d i n section  \[ s ec: FATT\].
Ourresults fo r t hePDF o f $ \ dot{\ thet a }_ p$a nd it s mo m ents diff e rsom e w ha t from thos e i n [ @Bec2 014 ] . This  isin part due to th e  difference in  the f orc ing meth od s employed by[@B ec 2 0 14] andth at in our s tudy. Th es oluti on ofmay be  writte n as [@Bec2 014]  $$ {\left\la ngl e\ dot{\th et a_ p }^2 \r ight \r angle}  = \fr a c{1} { 2 \tau_\theta^3}\i nt_0^ \ i nfty{ \le ft\la ng le \
 ${\text{\textit{St}}}_\theta$,_and for_simulations S1 and S2._The results_show_that as_${\text{\textit{St}}}_\theta$_is increased, $\chi_p$_increases. This is_mainly because as ${\text{\textit{St}}}_\theta$_is increased, the_thermal_memory of the particle increases, and the particle temperature depends strongly on its encounter_with_the fluid_temperature_field_along its trajectory history for_times up to ${\textit{O}\left( \tau_\theta_\right)}$ in_the past. As a result, the particle temperature_can_differ strongly from_the local fluid temperature. The results also show that_$\chi_p$ is dramatically suppressed when two-way_thermal coupling is_accounted_for._One reason for this_is that as shown earlier, two-way_thermal coupling leads to a suppression_in the fluid temperature gradients. As these_gradients are suppressed, the fluid temperature_along the particle trajectory history_differs less_from the local fluid temperature_than it would_have in_the absence of_two-way thermal coupling, and as a_result $\chi_p$ is_decreased.

The results for kurtosis of $\dot{\theta}_p$,_as_a function of_${\text{\textit{St}}}_\theta$_and_for various_${\text{\textit{St}}}$ are shown_in_figure \[fig:PDFdeT\](f)._The_results show that the kurtosis decreases_with_increasing ${\text{\textit{St}}}_\theta$. This is mainly due to_the filtering effect mentioned_earlier,_wherein as ${\text{\textit{St}}}_\theta$ is_increased, the particles are less_able to respond to rapid fluctuations_in the_fluid temperature_along their trajectory. Further, the kurtosis is typically larger when the_two-way thermal coupling is taken into_account (simulations S1), and_is maximum_for_${\text{\textit{St}}}=1$. This is_due_to the_particle clustering on the fronts of the_fluid temperature_field, as will be discussed in_section \[sec:FATT\].

Our results for_the_PDF of $\dot{\theta}_p$ and its moments_differ somewhat from those in [@Bec2014]._This is in part due_to_the_difference in the forcing methods_employed by [@Bec2014] and that in_our study. The_solution of may be written as [@Bec2014]_$${\left\langle_\dot{\theta_p}^2 \right\rangle} = \frac{1}{2\tau_\theta^3}\int_0^\infty {\left\langle \