 the construction of stable PI methods for the full BGK equation with ${\nu}= 1$ and stable TPI methods for with ${\nu}= \rho$ is practically identical to that in [@LafitteMelisSamaey2017] and [@MelisSamaey2017], respectively. The choice of method parameters for the full BGK equation will be discussed more closely in section \[sec:method\_param\].

Linearized Boltzmann equation and its spectrum {#subsec:lin_boltzmann_equation}
----------------------------------------------

### Linearization of the Boltzmann equation {#subsubsec:lin_boltzmann_equation}

To simplify the analysis of the Boltzmann equation, it is customary to linearize the collision operator ${{\mathcal{Q}}{\mathopen{}\mathclose\bgroup\originalleft}({f^{\varepsilon}}{\aftergroup\egroup\originalright})({\mathbf{v}})}$ around the global Maxwellian distribution ${{\mathcal{M}}^{{\rho^{\infty}},{{\mathbf{{\bar{v}}}}^{\infty}},{T^{\infty}}}_{{\mathbf{v}}}}= {\mathcal{M}}_{\mathbf{v}}^{1,0,1}$, which, for $D_v = 2$, is given by: $$\label{eq:global_maxwellian}
    {\mathcal{M}}_{\mathbf{v}}^{1,0,1} = \frac{1}{2\pi} \exp{\mathopen{}\mathclose\bgroup\originalleft}(-\frac{{{\mathopen{}\mathclose\bgroup\originalleft}\vert{\mathbf{v}}{\aftergroup\egroup\originalright}\vert}^2}{2}{\aftergroup\egroup\originalright}),$$ see, for instance, [@Cercignani1988; @Golse2005; @Saint-Raymond2009]. Subsequently, we consider small fluctuations of ${f^{\varepsilon}}$ around the global equilibrium, that is: $$\label{eq:linearization_f}
    {f^{\varepsilon}}({\mathbf{v}}) = {\mathcal{M}}_{\mathbf{v}}^{1,0,1}\big(1 + {g^{{\varepsilon}}}({\mathbf{v}})\big).$$ Using a similar shorthand notation as before, that is, $g' = {g^{{\varepsilon}}}({\mathbf{v}}')$, ${\mathcal{M}}' = {\mathcal{M}}_{{\mathbf{v}}'}^{1,0,1}$, and so on, we linearize the quadratic terms in as: $$\label{eq:linearizing_
the construction of stable PI methods for the full BGK equality with $ { \nu}= 1 $ and static TPI method acting for with $ { \nu}= \rho$ is practically identical to that in [ @LafitteMelisSamaey2017 ] and [ @MelisSamaey2017 ], respectively. The option of method parameter for the full BGK equation will be discourse more closely in section   \[sec: method\_param\ ]. 

 linearize Boltzmann equality and its spectrum { # subsec: lin_boltzmann_equation } 
 ---------------------------------------------- 

 # # # Linearization of the Boltzmann equation { # subsubsec: lin_boltzmann_equation } 

 To simplify the analysis of the Boltzmann equality, it is customary to linearize the collision operator $ { { \mathcal{Q}}{\mathopen{}\mathclose\bgroup\originalleft}({f^{\varepsilon}}{\aftergroup\egroup\originalright})({\mathbf{v}})}$ around the ball-shaped Maxwellian distribution $ { { \mathcal{M}}^{{\rho^{\infty}},{{\mathbf{{\bar{v}}}}^{\infty}},{T^{\infty}}}_{{\mathbf{v}}}}= { \mathcal{M}}_{\mathbf{v}}^{1,0,1}$, which, for $ D_v = 2 $, is given by: $ $ \label{eq: global_maxwellian } 
     { \mathcal{M}}_{\mathbf{v}}^{1,0,1 } = \frac{1}{2\pi } \exp{\mathopen{}\mathclose\bgroup\originalleft}(-\frac{{{\mathopen{}\mathclose\bgroup\originalleft}\vert{\mathbf{v}}{\aftergroup\egroup\originalright}\vert}^2}{2}{\aftergroup\egroup\originalright}),$$ see, for instance, [ @Cercignani1988; @Golse2005; @Saint - Raymond2009 ]. Subsequently, we consider small fluctuations of $ { f^{\varepsilon}}$ around the global equilibrium, that is: $ $ \label{eq: linearization_f } 
     { f^{\varepsilon}}({\mathbf{v } }) = { \mathcal{M}}_{\mathbf{v}}^{1,0,1}\big(1 + { g^{{\varepsilon}}}({\mathbf{v}})\big).$$ use a similar shorthand notation as before, that is, $ thousand' = { g^{{\varepsilon}}}({\mathbf{v}}')$, $ { \mathcal{M } }' = { \mathcal{M}}_{{\mathbf{v}}'}^{1,0,1}$, and so on, we linearize the quadratic terms in as: $ $ \label{eq: linearizing _
 thf construction of stable PI methods for the funl BGK equatiov with ${\nu}= 1$ and stable TPI mevhodw for with ${\nu}= \rho$ is practizally idejtical ti thet in [@LafitteMelmaSamaey2017] and [@MspisScmeey2017], respectiveli. The choice of method pardmdtzrs for the full BGK equation will bq discuxsfd more closeli in xqctikn \[sec:method\_param\].

Linearized Boltzmahn equauion and its specttum {#subsec:lin_boltzmann_equahion}
----------------------------------------------

### Linearization of hhe Boltzmabn ezyation {#subsucsec:lin_boluzkann_equatikn}

To simplify the analysis of tfe Boktzmann eqocrioj, it is custimary to linearizc the cmllisiom operator ${{\matmcal{Q}}{\katyopen{}\mathclose\bgroup\ociginalleft}({f^{\varepsiljn}}{\aftergrmu'\egroup\originalright})({\mqtybf{v}})}$ droutd tfw guobzl Mzxwelllan distributjon ${{\mathcal{N}}^{{\rho^{\infty}},{{\mathbf{{\bar{v}}}}^{\onsnu}},{T^{\infty}}}_{{\mathbr{v}}}}= {\matrcwl{M}}_{\mathbf{v}}^{1,0,1}$, which, for $D_v = 2$, is given by: $$\ldbem{eq:global_maxwellian}
    {\nathcal{M}}_{\mathbf{v}}^{1,0,1} = \frac{1}{2\ki} \exp{\mathjpen{}\mathclose\bgroup\originalleft}(-\frac{{{\mathopen{}\mathclmse\bgcojp\oxlninauoevt}\vert{\mathbf{v}}{\aftergroup\egroup\originalright}\vere}^2}{2}{\zfuernroup\egroup\origikalright}),$$ see, for imshamse, [@Cercignani1988; @Golse2005; @Sajnt-Raymond2009]. Subsequfntly, wg consuder smalj flictuations of ${f^{\varepsilon}}$ aeound the gljval equilibrium, thct is: $$\label{ed:lingarizayion_f}
    {f^{\varepsilon}}({\matkbf{v}}) = {\jathcal{M}}_{\matjbf{v}}^{1,0,1}\big(1 + {f^{{\xarepsilon}}}({\mathbf{x}})\bin).$$ Uving a similar shorthand njtation aw beyore, thag is, $g' = {g^{{\vwrepsilon}}}({\mwthbf{y}}')$, ${\mathcal{M}}' = {\mathcal{L}}_{{\mathyf{v}}'}^{1,0,1}$, atd so on, wf linearize the quadratic terms in as: $$\label{eq:kiteagizing_
the construction of stable PI methods for BGK with ${\nu}= and stable TPI is identical to that [@LafitteMelisSamaey2017] and [@MelisSamaey2017], The choice of method parameters for full BGK equation will be discussed more closely in section \[sec:method\_param\]. Linearized Boltzmann and its spectrum {#subsec:lin_boltzmann_equation} ---------------------------------------------- ### Linearization of the Boltzmann equation {#subsubsec:lin_boltzmann_equation} To the of Boltzmann it is customary to linearize the collision operator ${{\mathcal{Q}}{\mathopen{}\mathclose\bgroup\originalleft}({f^{\varepsilon}}{\aftergroup\egroup\originalright})({\mathbf{v}})}$ around the global Maxwellian distribution ${{\mathcal{M}}^{{\rho^{\infty}},{{\mathbf{{\bar{v}}}}^{\infty}},{T^{\infty}}}_{{\mathbf{v}}}}= {\mathcal{M}}_{\mathbf{v}}^{1,0,1}$, for $D_v = 2$, is given by: $$\label{eq:global_maxwellian} = \frac{1}{2\pi} \exp{\mathopen{}\mathclose\bgroup\originalleft}(-\frac{{{\mathopen{}\mathclose\bgroup\originalleft}\vert{\mathbf{v}}{\aftergroup\egroup\originalright}\vert}^2}{2}{\aftergroup\egroup\originalright}),$$ see, instance, [@Cercignani1988; @Golse2005; @Saint-Raymond2009]. Subsequently, consider fluctuations of around global that is: $$\label{eq:linearization_f} = {\mathcal{M}}_{\mathbf{v}}^{1,0,1}\big(1 + {g^{{\varepsilon}}}({\mathbf{v}})\big).$$ Using a similar shorthand notation as before, that is, $g' = {g^{{\varepsilon}}}({\mathbf{v}}')$, ${\mathcal{M}}' {\mathcal{M}}_{{\mathbf{v}}'}^{1,0,1}$, and we linearize quadratic in $$\label{eq:linearizing_
 the construction of stable PI Methods for The fuLl BgK eQuAtioN witH ${\nu}= 1$ and stable TPi MethOds for with ${\nu}= \rho$ is practIcallY iDEntiCAl To thaT in [@LafiTTemELissaMaEy2017] aNd [@mElIsSamAey2017], RespectIvely. The chOicE oF method paramETeRs for the fuLl BgK equation wiLl bE discuSsEd mORe cloSelY in seCtion \[sEC:methoD\_param\].

LinEaRIzed BoLTzmann eQUAtIon aNd its spectrum {#subsEC:lIN_boltzmann_equaTion}
----------------------------------------------

### LiNeARiZATioN of The BoltzmaNn EquatIOn {#subsuBSeC:LIN_boLTzmann_equatioN}

To simplify THe aNalysiS oF thE boltzmAnn eqUaTIon, It is customaRy to Linearize The colLIsion opERator ${{\maThcal{Q}}{\MatHopEn{}\maTHcLoSe\bGrOUp\oRIgInaLLefT}({f^{\varepsIlOn}}{\AfterGrouP\EGROup\oRigInalRight})({\Mathbf{v}})}$ around The GlobAL MaXwellIan diStriBuTion ${{\mAthcal{m}}^{{\rho^{\iNfTy}},{{\mathbf{{\bar{v}}}}^{\infTy}},{T^{\iNfty}}}_{{\mathbF{v}}}}= {\mAtHcaL{M}}_{\MathbF{V}}^{1,0,1}$, which, For $d_v = 2$, iS given bY: $$\label{eQ:GloBaL_MAXwEllian}
    {\mathcal{M}}_{\mathBf{V}}^{1,0,1} = \FRaC{1}{2\pi} \exp{\maThopen{}\MAtHcLOse\bgrouP\oRigInalLEFt}(-\fraC{{{\matHOpEn{}\mathclOse\bgrOUp\OrIginallEfT}\vert{\mAtHbf{V}}{\afTergrOUp\egRoup\orIginalriGht}\veRT}^2}{2}{\aftergroup\egrOUp\originalrigHT}),$$ sEE, FoR InstAncE, [@Cercignani1988; @golsE2005; @sainT-RayMOnD2009]. SuBSequeNtly, wE cONsIDer small fluctuationS oF ${f^{\varePsiloN}}$ around the gloBal equilibRIUM, that is: $$\lAbel{EQ:lINearization_f}
    {f^{\vArepsIlon}}({\mathbf{V}}) = {\Mathcal{M}}_{\MathbF{v}}^{1,0,1}\big(1 + {g^{{\vaRepsilon}}}({\mATHbf{v}})\big).$$ USinG a sImiLar SHOrThand notation AS BefoRe, That is, $g' = {G^{{\vaRepsiloN}}}({\maThbF{v}}')$, ${\mAthCaL{M}}' = {\mathcal{m}}_{{\mathbf{v}}'}^{1,0,1}$, AnD sO oN, wE liNeariZE the quadRaTic TeRms In as: $$\lABel{eq:lInearIzinG_
 the construction of stabl e PI metho ds fo r t hefu ll B GK e quation with $ { \nu} = 1$ and stable TPI me thods  f o r wi t h${\nu }= \rho $  i s pra ct ic all yi de ntica l t o thatin [@Lafit teM el isSamaey2017 ]  a nd [@Melis Sam aey2017], re spe ctivel y.  Th e  choi ceof me thod p a ramete rs for th ef ull BG K  equati o n  w illbe discussed more cl o sely in sectio n \[se c: m et h o d\_ par am\].

Lin ea rized Boltzma n ne q u ati o n and its spe ctrum {#sub s ec: lin_bo lt zma n n_equa tion} 
- - --- ----------- ---- --------- ------ - ------- - --

###  Linea riz ati on o f  t he  Bo lt z man n  e qua t ion  {#subsu bs ec :lin_ bolt z m a n n_eq uat ion} 

Tosimplify theana lysi s  of  theBoltz mann  e quati on, it  is c us tomary to linea rize  the coll isi on  op er ator$ {{\mat hca l{Q }}{\mat hopen{} \ mat hc l o s e\ bgroup\originallef t} ( { f^ {\vareps ilon}} { \a ft e rgroup\e gr oup \ori g i nalri ght} ) ({ \mathbf{ v}})}$ ar ou nd thegl obal M ax wel lia n dis t ribu tion $ {{\mathc al{M} } ^{{\rho^{\inft y }},{{\mathbf{ { \b a r {v } }}}^ {\i nfty}},{T^{ \inf t y}}} _{{\ m at hbf { v}}}} = {\m at h ca l {M}}_{\mathbf{v}}^{ 1, 0,1}$,  whic h, for $D_v =  2$, is gi v e n  by: $$\ labe l {e q :global_maxwel lian} 
    {\mat h cal{M}}_ {\mat hbf{v}}^ {1,0,1} = \ frac{1}{ 2\p i}\ex p{\ m a th open{}\mathcl o s e\bg ro up\orig ina lleft}( -\f rac {{{ \ma th open{}\ma thclose\ bg ro up \o rig inall e ft}\vert {\ mat hb f{v }}{\a f tergro up\eg roup \o ri g ina lright} \ ve r t }^2} {2 }{ \aft erg ro up\eg roup \ ori ginalri ght}),$$see ,  for  i ns tance,[@Cercignani1 98 8; @Golse2 00 5;@Saint - R aymond20 09]. Subsequently, we c o nsidersma ll fl uctu ations of  ${ f^{\va rep s ilon}} $ arou nd th eglo b a l equ i l ib riu m,  that is:$ $ \la bel{e q: line arizati on_f}
    {f^{\var e psi lon}}({\mathb f{v }})= {\ mat h ca l {M} }_ { \ma t h bf{v}}^{1,0,1}\ big(1 + {g ^{ { \v arepsilon} } }({ \m athbf{v }})\big ).$$U sing asimilar s horthandno tati o n  as  before, t hat is,$g' = {g^ { {\var e ps ilon} }}( {\math bf {v} }')$,  ${\ma t hca l{M}} ' = {\ ma thcal{ M}}_{ {\ mathbf{v }}'}^{1,0,1}$, and so o n, welinea riz e the qua dra t icterms inas:$$\label{e q:l ine arizi ng_ 
 the_construction of_stable PI methods for_the full_BGK_equation with_${\nu}=_1$ and stable_TPI methods for_with ${\nu}= \rho$ is_practically identical to_that_in [@LafitteMelisSamaey2017] and [@MelisSamaey2017], respectively. The choice of method parameters for the full BGK_equation_will be_discussed_more_closely in section \[sec:method\_param\].

Linearized Boltzmann equation_and its spectrum {#subsec:lin_boltzmann_equation}
----------------------------------------------

### Linearization_of the_Boltzmann equation {#subsubsec:lin_boltzmann_equation}

To simplify the analysis of the_Boltzmann_equation, it is_customary to linearize the collision operator ${{\mathcal{Q}}{\mathopen{}\mathclose\bgroup\originalleft}({f^{\varepsilon}}{\aftergroup\egroup\originalright})({\mathbf{v}})}$ around the_global Maxwellian distribution ${{\mathcal{M}}^{{\rho^{\infty}},{{\mathbf{{\bar{v}}}}^{\infty}},{T^{\infty}}}_{{\mathbf{v}}}}= {\mathcal{M}}_{\mathbf{v}}^{1,0,1}$, which,_for $D_v =_2$,_is_given by: $$\label{eq:global_maxwellian}
 _  {\mathcal{M}}_{\mathbf{v}}^{1,0,1} = \frac{1}{2\pi} \exp{\mathopen{}\mathclose\bgroup\originalleft}(-\frac{{{\mathopen{}\mathclose\bgroup\originalleft}\vert{\mathbf{v}}{\aftergroup\egroup\originalright}\vert}^2}{2}{\aftergroup\egroup\originalright}),$$_see, for instance, [@Cercignani1988; @Golse2005; @Saint-Raymond2009]._Subsequently, we consider small fluctuations of ${f^{\varepsilon}}$_around the global equilibrium, that is:_$$\label{eq:linearization_f}
    {f^{\varepsilon}}({\mathbf{v}})_= {\mathcal{M}}_{\mathbf{v}}^{1,0,1}\big(1_+ {g^{{\varepsilon}}}({\mathbf{v}})\big).$$ Using a similar_shorthand notation as_before, that_is, $g' =_{g^{{\varepsilon}}}({\mathbf{v}}')$, ${\mathcal{M}}' = {\mathcal{M}}_{{\mathbf{v}}'}^{1,0,1}$, and so_on, we linearize_the quadratic terms in as: $$\label{eq:linearizing_