 D/L$, where $L$ is the height of the sample and $D$ its width, $g$ the gravitational acceleration, $\beta$ is the thermal expansion coefficient and $\nu$ and $\kappa$ the kinematic viscosity and the thermal diffusivity, respectively.

RB convection is used as a simple model for all types of natural convection, regardless of geometry. For 2D RB, at a moderate ${\text{Ra}}= 10^8$, the heat flux is expected to be independent of the aspect-ratio, and thus of the sidewall boundary conditions, beyond $\Gamma \gtrsim 26$ [@poe12]. This value of $\Gamma$ is expected to be lower for 3D RB [@bai10]. Low aspect-ratio applications with a no-slip BC on the lateral walls, such as cooling fins and domestic heating, have moderate Rayleigh numbers (${\text{Ra}}< 10^{10}$) that are accessible with contemporary direct numerical simulations (DNS). However, Nature is dominated by flows with large aspect-ratios and inaccessible high Ra [@chi12].

In both experiments and numerics, the aspect-ratio $\Gamma$ is commonly reduced to decrease the volume of the setup, as Ra depends only on the height of the system. In simulations, reducing the aspect-ratio for fixed Ra decreases the computational complexity by $\mathcal{O}(\Gamma^2)$ in three-dimensions, while in experiments it means a smaller setup volume. In both cases it makes higher Rayleigh numbers more accesible. Despite the additional effort required for large aspect-ratio, research in this direction has been done in simulations [@bai10] and experiments [@fun05; @pui13; @hog13]. However, these endeavours were all limited in Ra due the large aspect-ratios used. In simulations one can approximate infinitely large $\Gamma$ by using periodic boundary conditions at the sidewalls. This leads to the question how large the simulated domain has to be to sufficiently approximate infinite $\Gamma$, as a truly infinite domain can in principle support any wavenumber while a periodic domain allows only for multiples of the fundamental one imposed by its size. This is addressed in this study in 2D RB, which is expected to require a larger domain than 3D and therefore is more demanding in terms of computational resources.

Different temperature boundary conditions at the plates were studied by [@joh09], where no significant difference between constant flux and isothermal plates was found in the higher Ra
D / L$, where $ L$ is the height of the sample and $ D$ its width, $ g$ the gravitational acceleration, $ \beta$ is the thermal expansion coefficient and $ \nu$ and $ \kappa$ the kinematic viscosity and the thermal diffusivity, respectively. 

 RB convection is used as a simple exemplar for all type of natural convection, regardless of geometry. For 2D RB, at a moderate $ { \text{Ra}}= 10 ^ 8 $, the heat flux is ask to be independent of the aspect - proportion, and thus of the sidewall boundary conditions, beyond $ \Gamma \gtrsim 26 $ [ @poe12 ]. This value of $ \Gamma$ is expected to be broken for 3D RB [ @bai10 ]. Low aspect - proportion applications with a no - slip BC on the lateral walls, such as cooling flipper and domestic heating, get moderate Rayleigh numbers ($ { \text{Ra } } < 10^{10}$) that are accessible with contemporary lineal numerical simulations (DNS). However, Nature is dominated by flows with big aspect - ratios and inaccessible eminent Ra [ @chi12 ]. 

 In both experiments and numerics, the view - ratio $ \Gamma$ is commonly reduced to decrease the volume of the setup, as Ra depend only on the height of the system. In simulations, reducing the aspect - ratio for fixed Ra decreases the computational complexity by $ \mathcal{O}(\Gamma^2)$ in three - dimensions, while in experiments it means a smaller setup volume. In both cases it makes high Rayleigh numbers more accesible. Despite the additional campaign required for large aspect - ratio, inquiry in this direction has been suffice in simulations [ @bai10 ] and experiments [ @fun05; @pui13; @hog13 ]. However, these enterprise were all limited in Ra due the large aspect - ratios use. In simulations one can approximate infinitely large $ \Gamma$ by use periodic boundary condition at the sidewall. This leads to the question how large the simulate domain has to be to sufficiently approximate infinite $ \Gamma$, as a in truth infinite domain can in principle support any wavenumber while a periodic domain allow only for multiples of the fundamental one imposed by its size. This is addressed in this survey in 2D RB, which is expected to necessitate a larger knowledge domain than 3D and therefore is more demanding in term of computational resources. 

 Different temperature boundary conditions at the plates were studied by [ @joh09 ], where no significant difference between constant flux and isothermal plates was found in the high Ra
 D/L$, where $L$ is the height on the sample and $D$ its xidth, $g$ the graxitational acceleration, $\beta$ iw the thermal expansion coewficient wnd $\nu$ abd $\keppa$ the kinematmd viscosity ahf thz vhermal diffusiyity, respecdively.

RB convewtkou is used as a simple model for all eypes og jatural convecjion, gedardmvsw of geometry. For 2D RB, at a mkderate ${\text{Ra}}= 10^8$, the neat flux is expected to bf infependent of the adpect-ratio, qnd eyus of the skdewall boludary condijions, beyond $\Gamma \gtrsim 26$ [@poe12]. Tfis vclue of $\Gamna$ is gxpected to ue lowvr for 3D RB [@bai10]. Low avpect-rayio applicatioks wivh a no-slip BC on the latxral walls, such as cjoling fits and domestic hearibg, haee mmderxre Faymemgh numbegs (${\vext{Ra}}< 10^{10}$) thaf are acceswible with contempotawj direct numedical fiiulations (DNS). However, Nature is dominattd by flows with large aspecr-ratios and inaccessihle high Wa [@chi12].

In both experiments and numerics, the aspect-satio $\Eamna$ kw fommonly reduced to decrease the volume of thq sttui, as Ra depends okly on the height pf tng system. In sioulations, reducing the aspeft-ratio for dixed Ra qecrrases the computational complexity by $\mcthxal{O}(\Gamma^2)$ in three-bimensions, wkile im exprriments it means a smauler setup volule. In botg cases it makes figmer Rayleigh numbers more accqsible. Dewpitz the adaitipnal esfort requlred nmr large aspect-ratlo, rededrch in thls direction has been done in smkulations [@bao10] dnd experimznts [@fmn05; @pui13; @hog13]. Howeder, these endecvours wzre alu limited pn Ra due the large wspect-ratios gded. In simuletions onq cab approximagd infinitely lsrge $\Gammc$ by ysing periodic boukdary donditions at tkt sudewalls. This lradr tj nhe qtastion how ldrge thd simuuated domaik hxs tp be to sufficiently appdoximate infinite $\Baima$, as a truly igfinite domaim can in principle suppmrt any eavgnumber while a periodic domain allows ojly for multiplqs on thg fundamentcl one imposed by its size. This is addrewsed in this studi in 2D RB, which is exkecued to requirx a lawger domahn than 3D and therefire is more demancing in terms of comphtatiotal rfsources.

Different temperature boundary conditions at the plates were studued by [@joh09], where nk sibnifibanc dnfferense bxtxeen constant fluq and isothermal plates was founv in the hhgker Ra
D/L$, where $L$ is the height of and its width, the gravitational acceleration, coefficient $\nu$ and $\kappa$ kinematic viscosity and thermal diffusivity, respectively. RB convection is as a simple model for all types of natural convection, regardless of geometry. 2D RB, at a moderate ${\text{Ra}}= 10^8$, the heat flux is expected to independent the and of the sidewall boundary conditions, beyond $\Gamma \gtrsim 26$ [@poe12]. This value of $\Gamma$ is expected be lower for 3D RB [@bai10]. Low aspect-ratio with a no-slip BC the lateral walls, such as fins domestic heating, moderate numbers 10^{10}$) that are with contemporary direct numerical simulations (DNS). However, Nature is dominated by flows with large aspect-ratios and inaccessible Ra [@chi12]. experiments and the $\Gamma$ commonly reduced to volume of the setup, as Ra the height of the system. In simulations, reducing aspect-ratio for Ra decreases the computational complexity by in three-dimensions, while in experiments it means a setup volume. In both cases it makes higher Rayleigh numbers more accesible. Despite the additional for large aspect-ratio, research this direction has done simulations and [@fun05; @pui13; However, these endeavours were all limited in Ra due the large used. In simulations one can approximate infinitely large $\Gamma$ by boundary at the sidewalls. leads to the question large simulated domain has to sufficiently infinite truly domain in principle support any while a periodic domain allows for multiples of the size. This is addressed in this study in RB, which is expected to require a domain than 3D and therefore is more demanding in terms of computational Different temperature at the plates were studied by [@joh09], where significant difference between constant and isothermal plates was found in the higher Ra
 D/L$, where $L$ is the height of the sAmple and $D$ iTs widTh, $g$ The GrAvitAtioNal acceleratioN, $\Beta$ Is the thermal expansion cOeffiCiENt anD $\Nu$ And $\kaPpa$ the kINeMATic ViScOsiTy ANd The thErmAl diffuSivity, respEctIvEly.

RB convectIOn Is used as a sImpLe model for alL tyPes of nAtUraL ConveCtiOn, regArdlesS Of geomEtry. For 2D Rb, aT A moderATe ${\text{RA}}= 10^8$, THe Heat Flux is expected to bE InDEpendent of the aSpect-rAtIO, aND ThuS of The sidewalL bOundaRY conditIOnS, BEYonD $\gamma \gtrsim 26$ [@poE12]. This value oF $\gamMa$ is exPeCteD To be loWer foR 3D rb [@baI10]. Low aspect-rAtio ApplicatiOns witH A no-slip bc on the lAteral WalLs, sUch aS CoOlIng FiNS anD DoMesTIc hEating, haVe MoDeratE RayLEIGH numBerS (${\texT{Ra}}< 10^{10}$) thAt are accessibLe wIth cONteMporaRy dirEct nUmEricaL simulAtionS (DnS). However, Nature Is doMinated by FloWs WitH lArge aSPect-raTioS anD inacceSsible hIGh RA [@cHI12].

iN bOth experiments and nUmERIcS, the aspeCt-ratiO $\gaMmA$ Is commonLy RedUced TO DecreAse tHE vOlume of tHe setuP, As ra Depends OnLy on thE hEigHt oF the sYStem. in simuLations, rEduciNG the aspect-ratiO For fixed Ra decREaSES tHE comPutAtional compLexiTY by $\mAthcAL{O}(\gamMA^2)$ in thRee-diMeNSiONs, while in experimentS iT means A smalLer setup volumE. In both casES IT makes hiGher rAyLEigh numbers morE acceSible. DespiTE the addiTionaL effort rEquired foR LArge aspeCt-rAtiO, reSeaRCH iN this directioN HAs beEn Done in sImuLations [@Bai10] And ExpEriMeNts [@fun05; @pui13; @Hog13]. HowevEr, ThEsE eNdeAvourS Were all lImIteD iN Ra Due thE Large aSpect-RatiOs UsED. In SimulatIOnS ONe caN aPpRoxiMatE iNfiniTely LArgE $\Gamma$ bY using perIodIC bouNdArY conditIons at the sideWaLls. This leaDs To tHe quesTIOn how larGe the simulated domain has TO be to suFfiCientLy apProximate InfInite $\GAmmA$, As a truLy infiNite dOmAin CAN in prINCiPle SuPport any waVENumBer whIlE a peRiodic dOmain allows only for MUltIples of the funDamEntaL ONe ImpOSeD By iTs SIze. tHIs is addressed in This study iN 2D rb, wHich is expeCTed To Require A larger DomaiN Than 3D anD thereforE is more deMaNdinG IN teRms of compuTational Resources.

dIfferENt TempeRatUre bouNdAry CondiTions aT The PlateS were sTuDied by [@Joh09], whErE no signiFicant difference between ConstaNt fluX anD isothermAl pLAteS was found In thE higher Ra
 D/L$, where $L$ is the he ight of th e sam ple  an d$D$itswidth, $g$ the grav itational acceleration , $\b et a $ is th e the rmal ex p an s i onco ef fic ie n tand $ \nu $ and $ \kappa$ th e k in ematic visco s it y and thethe rmal diffusi vit y, res pe cti v ely.
RB  conv ection is use d as a si mp l e mode l  for al l ty pesof natural convec t io n , regardless o f geom et r y. F or2DRB, at a m od erate ${\text { Ra } } =  10 ^ 8$, the heatflux is exp e cte d to b eind e penden t ofth e  as pect-ratio,  and  thus ofthe si d ewall b o undarycondit ion s,beyo n d$\ Gam ma \gt r si m 2 6 $ [ @poe12].  T hi s val ue o f $ \ Gamm a$is e xpect ed to be lowe r f or 3 D  RB  [@ba i10].  Low  a spect -ratio  appl ic ations with a n o-sl ip BC onthe  l ate ra l wal l s, suc h a s c oolingfins an d  do me s t i cheating, have mode ra t e  R ayleighnumber s  ( ${ \ text{Ra} }<  10 ^{10 } $ ) tha t ar e  a ccessibl e with co nt emporar ydirect  n ume ric al si m ulat ions ( DNS). Ho wever ,  Nature is dom i nated by flow s  w i t hl arge  as pect-ratios  and inac cess i bl e h i gh Ra  [@ch i1 2 ]. 
In both experiment sand nu meric s, the aspect -ratio $\G a m m a$ is co mmon l yr educed to decr easethe volume of the s etup,  as Ra d epends on l y  on thehei ght  of  th e sy stem. In simu l a tion s,  reduci ngthe asp ect -ra tio  fo rfixed Radecrease sth eco mpu tatio n al compl ex ity  b y $ \math c al{O}( \Gamm a^2) $in thr ee-dime n si o n s, w hi le  inexp er iment s it mea ns a sm aller set upv olum e.  I n bothcases it make shigher Ray le igh  numbe r s  more ac cesible. Despite the ad d itional  ef fortrequ ired forlar ge asp ect - ratio,  resea rch i nthi s direc t i on  ha sbeen donei n  si mulat io ns [ @bai10]  and experiments [ @ fun 05; @pui13; @ hog 13]. H ow eve r ,t hes ee nde a v ours were all l imited inRa du e the larg e  as pe ct-rati os used . Ins imulati ons one c an approx im atei n fin itely larg e $\Gamm a$ by usi n g per i od ic bo und ary co nd iti ons a t thes ide walls . This  l eads t o the  q uestionhow large the simulated  domai n has  to  be to su ffi c ien tly appro xima te infinit e $ \Ga mma$,  as a tru ly i n fi nit e  doma in c a n in prin c ip les u pp ort any wav e n u mbe r whi lea  perio dicdomain allows onl y  for multiples  oft h e f und a ment al  one imposed b y i ts s ize. Thi sis addresse d in thi ss tudyin 2DRB, wh ich ise x pe c ted to  req uir e a large r d om a in than  3 Da nd the refo re  is mo re dem a ndin g in terms of comp utati o n al re s our ces.
D ifferen t  tem perature b oundary con dition s at  theplateswe re stu die dby [@joh09 ] , where n o sig nifican tdiff ere nce be twee n const antfl uxand isoth e r ma l  p la t eswasfound  i n th e higherR a
 D/L$,_where $L$_is the height of_the sample_and_$D$ its_width,_$g$ the gravitational_acceleration, $\beta$ is_the thermal expansion coefficient_and $\nu$ and_$\kappa$_the kinematic viscosity and the thermal diffusivity, respectively.

RB convection is used as a simple_model_for all_types_of_natural convection, regardless of geometry._For 2D RB, at a_moderate ${\text{Ra}}=_10^8$, the heat flux is expected to be_independent_of the aspect-ratio,_and thus of the sidewall boundary conditions, beyond $\Gamma_\gtrsim 26$ [@poe12]. This value of_$\Gamma$ is expected_to_be_lower for 3D RB_[@bai10]. Low aspect-ratio applications with a_no-slip BC on the lateral walls,_such as cooling fins and domestic heating,_have moderate Rayleigh numbers (${\text{Ra}}< 10^{10}$)_that are accessible with contemporary_direct numerical_simulations (DNS). However, Nature is_dominated by flows_with large_aspect-ratios and inaccessible_high Ra [@chi12].

In both experiments and_numerics, the aspect-ratio_$\Gamma$ is commonly reduced to decrease_the_volume of the_setup,_as_Ra depends_only on the_height_of the_system._In simulations, reducing the aspect-ratio for_fixed_Ra decreases the computational complexity by $\mathcal{O}(\Gamma^2)$_in three-dimensions, while in_experiments_it means a smaller_setup volume. In both cases_it makes higher Rayleigh numbers more_accesible. Despite_the additional_effort required for large aspect-ratio, research in this direction has been_done in simulations [@bai10] and experiments_[@fun05; @pui13; @hog13]. However,_these endeavours_were_all limited in_Ra_due the_large aspect-ratios used. In simulations one can_approximate infinitely_large $\Gamma$ by using periodic boundary_conditions at the sidewalls._This_leads to the question how large_the simulated domain has to be_to sufficiently approximate infinite $\Gamma$,_as_a_truly infinite domain can in_principle support any wavenumber while a_periodic domain allows_only for multiples of the fundamental one_imposed_by its size. This is addressed_in_this study in 2D RB, which_is_expected_to require a larger domain_than 3D and therefore is more_demanding in terms of computational resources.

Different temperature boundary conditions_at the plates_were studied by [@joh09], where_no_significant_difference between constant flux and isothermal plates was found in_the higher_Ra