 the hydrogen atom, respectively. The gas is assumed to be a perfect gas with $\gamma = 5/3$ and $\mu = 1.4 m_H$.

The cooling function of @AudHen05 includes the cooling by fine-structure lines of CII and OI, further the cooling by H (Ly$\alpha$-line) and the electron recombination onto positively charged grains. The heating is due to the photo-electric effect on small grains and polycyclic aromatic hydrocarbons (PAH) caused by the far-ultraviolet galactic background radiation. For more information about this cooling function see @Wol95 [@WolMcK03; @Spi78; @BakTie94] and @Hab68. The pressure-equillibrium curve resulting from the cooling function is plotted as black curve in Figure \[fig:phasediag\_static5\]. For the numerical solution of the fluid dynamical equations, we used the radiative cooling routine implemented into Enzo. For each hydrodynamical time step, the state variables are iterated over several subcycles, and the resulting total energy increment for the whole time step is added.

For our numerical study, the two-dimensional setup of @AudHen05 and [@HenAud07] was adopted with small modifications. The initial gas content corresponds to the warm neutral material (WNM) in the ISM, i. e., the temperature is $T=7100$ K, the pressure is $P=7\times10^{-13}$ erg cm$^{-3}$ and the number density of neutral hydrogen is $n=0.71$ cm$^{-3}$. From the left and the right boundaries, warm gas with identical thermodynamic properties flows into the computational domain, where the cosine-shaped inflow velocity profile is modulated with small perturbations, realised by randomly shifted phases. These phase shifts are kept constant for the different simulations, so that initial conditions are exactly the same for all runs to ensure comparability. Following @HenBan08, the top and bottom boundary conditions are periodic. The physical dimensions of the computational domain are $20\times20$ pc. The two inflows of gas collide in the middle of the domain. The supersonic collision causes a steep rise in the gas density that triggers the thermal instability, and gas undergoes a transition into the phase of the cold neutral material (CNM) in the ISM. In this phase, the gas has temperatures in the range $30$–$100$ K and number densities within
the hydrogen atom, respectively. The gas is assumed to be a perfect gas with $ \gamma = 5/3 $ and $ \mu = 1.4 m_H$. 

 The cool affair of @AudHen05 includes the cooling by fine - social organization lines of CII and OI, further the cooling by H (Ly$\alpha$-line) and the electron recombination onto positively blame grains. The heating is ascribable to the photograph - electric effect on small grain and polycyclic aromatic hydrocarbons (PAH) caused by the far - ultraviolet galactic background radiation. For more information about this cool function see @Wol95 [ @WolMcK03; @Spi78; @BakTie94 ] and @Hab68. The pressure - equillibrium bend resulting from the cooling function is plotted as bootleg curve in Figure \[fig: phasediag\_static5\ ]. For the numeral solution of the fluid dynamical equations, we use the radiative cooling routine implemented into Enzo. For each hydrodynamical time step, the department of state variables are iterated over several subcycles, and the resulting total energy increment for the whole time step is add. 

 For our numerical study, the two - dimensional frame-up of @AudHen05 and [ @HenAud07 ] was assume with small modifications. The initial gas contentedness corresponds to the warm neutral material (WNM) in the ISM, i.   e., the temperature is $ T=7100 $ kilobyte, the pressure is $ P=7\times10^{-13}$ erg cm$^{-3}$ and the number density of neutral hydrogen is $ n=0.71 $ cm$^{-3}$. From the leftover and the right boundaries, warm gas with identical thermodynamic properties flows into the computational domain, where the cosine - determine inflow velocity profile is modulated with small perturbations, recognize by randomly shifted phases. These phase shifts are kept constant for the different simulations, so that initial conditions are precisely the same for all runs to guarantee comparability. be @HenBan08, the top and bottom boundary conditions are periodic. The forcible dimensions of the computational domain are $ 20\times20 $ pc. The two inflows of gas collide in the center of the domain. The supersonic collision causes a steep raise in the gas density that triggers the thermal instability, and gas undergo a transition into the phase of the cold neutral fabric (CNM) in the ISM. In this phase, the gas has temperatures in the stove $ 30$–$100 $ K and number densities within
 thf hydrogen atom, respectiyely. The gas is cwsumed to be a perfezt gas with $\gamma = 5/3$ and $\mu = 1.4 m_Y$.

The xooling function of @AuaHen05 inclldes the xoolmng by fine-strucvhre lincf of GII aud OI, further thg cooling by H (Ly$\alpha$-line) avd the electron recombination onto pofitivelu fharged grains. The reatjng is due to the photo-electric efrect on small grains and polycyclic aromatic hjdrofarbons (PAH) caused by the far-omtrwciolet galacgic background radiatikn. For more information about tfis cpoling funxtuon vee @Wol95 [@WolNcK03; @Sii78; @BakTie94] and @Hab68. The pressute-equillibrium cucve eesulting from the comling function is klotted as bmack curve in Figyrw \[fig:khasegiag\_rratkc5\]. Roc tge numfriral solutioh of the flyid dynamical equatoogw, we used the radiaeide cooling routine implemented into Enzm. Fkr each hydrodynamical rime step, the state vwriables wre iterated over several subcycles, and the resulding vogal ekergh ijcrement for the whole time step is added.

For jhr nlmerical study, thc two-dimensional srtkp jf @AudHen05 and [@HenAub07] wzs adopted with smwll modyficarions. The iniyial gas content correspondw to the wari neutral material (CNM) in the IRM, i. g., the yemperature is $T=7100$ K, the 'ressude is $P=7\timed10^{-13}$ erg cm$^{-3}$ zvd the number devsiny ox neutral hydrogen is $n=0.71$ cm$^{-3}$. From the lefc and thd ribht botndaries, wwrm gas with identical tjermobynamhc propertles flows into the computational domain, where tve bosine-sha'ed innlow velocity pwofile is modukated wnth smxll perturgations, realised br randomly shhvted phases. Vhese phafe syiftw are kdot constant fot the difytrent simuoations, so that inltial donditions are zraxtly the same fpr xll rlns tj ensure compdrabkligu. Foluowing @HenBcn08, the top and bottom boungary conditions are petijdic. The physicaj dimensions pf the computationwl dokaii are $20\yimgs20$ pc. The two inflows of gas comlide in hhe middle of tre dimain. The su'ersonic collision causes a steep rise ii the gas density that rriggers the thermak instabilitb, and das undercoes a transition inro the phase of tme cold neutral materizl (CNM) in tje ISM. In this phase, the gas has temperatures in the range $30$–$100$ K and number eensitmef within
the hydrogen atom, respectively. The gas is be perfect gas $\gamma = 5/3$ The function of @AudHen05 the cooling by lines of CII and OI, further cooling by H (Ly$\alpha$-line) and the electron recombination onto positively charged grains. The is due to the photo-electric effect on small grains and polycyclic aromatic hydrocarbons caused the galactic radiation. For more information about this cooling function see @Wol95 [@WolMcK03; @Spi78; @BakTie94] and @Hab68. The curve resulting from the cooling function is plotted black curve in Figure For the numerical solution of fluid equations, we the cooling implemented into Enzo. each hydrodynamical time step, the state variables are iterated over several subcycles, and the resulting total energy for the step is For numerical the two-dimensional setup and [@HenAud07] was adopted with small gas content corresponds to the warm neutral material in the i. e., the temperature is $T=7100$ the pressure is $P=7\times10^{-13}$ erg cm$^{-3}$ and the density of neutral hydrogen is $n=0.71$ cm$^{-3}$. From the left and the right boundaries, warm identical thermodynamic properties flows the computational domain, the inflow profile modulated with perturbations, realised by randomly shifted phases. These phase shifts are kept for the different simulations, so that initial conditions are exactly for runs to ensure Following @HenBan08, the top bottom conditions are periodic. The of computational pc. two of gas collide in middle of the domain. The collision causes a steep that triggers the thermal instability, and gas undergoes transition into the phase of the cold material (CNM) in the ISM. In this phase, the gas has temperatures the range and number densities within
 the hydrogen atom, respectiveLy. The gas is AssumEd tO be A pErfeCt gaS with $\gamma = 5/3$ and $\mU = 1.4 M_H$.

ThE cooling function of @AudHEn05 incLuDEs thE CoOling By fine-sTRuCTUre LiNeS of cIi AnD OI, fuRthEr the coOling by H (Ly$\AlpHa$-Line) and the elECtRon recombiNatIon onto positIveLy charGeD grAIns. ThE heAting Is due tO The phoTo-electriC eFFect on SMall graINS aNd poLycyclic aromatic hYDrOCarbons (PAH) causEd by thE fAR-uLTRavIolEt galactic BaCkgroUNd radiaTIoN. fOR moRE information aBout this cooLIng FunctiOn See @wOl95 [@WolMCK03; @Spi78; @baKtie94] And @Hab68. The prEssuRe-equilliBrium cURve resuLTing froM the coOliNg fUnctIOn Is PloTtED as BLaCk cURve In Figure \[FiG:pHasedIag\_sTATIC5\]. For The NumeRical Solution of the FluId dyNAmiCal eqUatioNs, we UsEd the RadiatIve coOlIng routine impleMentEd into EnzO. FoR eAch HyDrodyNAmical TimE stEp, the stAte variABleS aRE ITeRated over several suBcYCLeS, and the rEsultiNG tOtAL energy iNcRemEnt fOR The whOle tIMe Step is adDed.

For OUr NuMerical StUdy, the TwO-diMenSionaL SetuP of @Audhen05 and [@HeNAud07] wAS adopted with smALl modificatioNS. THE InITial Gas Content corrEspoNDs to The wARm NeuTRal maTeriaL (Wnm) iN The ISM, i. e., the temperatUrE is $T=7100$ K, tHe preSsure is $P=7\times10^{-13}$ Erg cm$^{-3}$ and thE NUMber densIty oF NeUTral hydrogen is $N=0.71$ cm$^{-3}$. FrOm the left aND the righT bounDaries, waRm gas with IDEntical tHerModYnaMic PROpErties flows inTO The cOmPutatioNal Domain, wHerE thE coSinE-sHaped inflOw velociTy PrOfIlE is ModulATed with sMaLl pErTurBatioNS, realiSed by RandOmLy SHifTed phasES. THESe phAsE sHiftS arE kEpt coNstaNT foR the difFerent simUlaTIons, So ThAt initiAl conditions aRe Exactly the SaMe fOr all rUNS to ensurE comparability. Following @hEnBan08, thE toP and bOttoM boundary ConDitionS arE PeriodIc. The pHysicAl DimENSions OF ThE coMpUtational dOMAin Are $20\tiMeS20$ pc. THe two inFlows of gas collide iN The Middle of the doMaiN. The SUPeRsoNIc COllIsIOn cAUSes a steep rise in The gas densItY ThAt triggers THe tHeRmal insTabilitY, and gAS undergOes a transItion into ThE phaSE Of tHe cold neutRal materIal (CNM) in tHE ISM. IN ThIs phaSe, tHe gas hAs TemPeratUres in THe rAnge $30$–$100$ K And numBeR densiTies wItHin
 the hydrogen atom, respec tively. Th e gas  is  as su medto b e a perfect ga s  wit h $\gamma = 5/3$ and $ \mu =  1 . 4 m_ H $. 

The  coolin g  f u n cti on  o f @ Au d He n05 i ncl udes th e coolingbyfi ne-structure li nes of CII  an d OI, furthe r t he coo li ngb y H ( Ly$ \alph a$-lin e ) andthe elect ro n  recom b ination o nt o po sitively chargedg ra i ns. The heatin g is d ue to t hepho to-electri ceffec t  on sma l lg r a ins and polycycli c aromatich ydr ocarbo ns  (P A H) cau sed b yt hefar-ultravi olet  galactic  backg r ound ra d iation.  For m ore  in form a ti on  ab ou t  th i scoo l ing  functio nse e @Wo l95[ @ W o lMcK 03;  @Sp i78;@BakTie94] an d @ Hab6 8 . T he pr essur e-eq ui llibr ium cu rve r es ulting from the  coo ling func tio nispl otted as bla ckcur ve in F igure \ [ fig :p h a s ed iag\_static5\]. Fo rt h enumerica l solu t io no f the fl ui d d ynam i c al eq uati o ns , we use d ther ad ia tive co ol ing ro ut ine  im pleme n tedinto E nzo. For  each hydrodynamical time step, th e  s t a te vari abl es are iter ated over  sev e ra l s u bcycl es, a nd th e  resulting total en er gy inc remen t for the who le time st e p is added .

F o ro ur numerical s tudy,  the two-d i mensiona l set up of @A udHen05 a n d  [@HenAu d07 ] w asado p t ed  with small m o d ific at ions. T heinitial  ga s c ont ent  c orrespond s to the  w ar mne utr al ma t erial (W NM ) i nthe  ISM, i. e.,  thetemp er at u reis $T=7 1 00 $ K, t he  p ress ure  i s $P= 7\ti m es1 0^{-13} $ erg cm$ ^{- 3 }$ a nd  t he numb er density of  n eutral hyd ro gen  is $n = 0 .71$ cm$ ^{-3}$. From the left a n d the r igh t bou ndar ies, warm  ga s with  id e ntical  therm odyna mi c p r o perti e s  f low sinto the c o m put ation al  dom ain, wh ere the cosine-sha p edinflow veloci typrof i l eism od u lat ed wit h small perturbat ions, real is e dby randoml y  sh if ted pha ses. Th ese p h ase shi fts are k ept const an t fo r the  different  simulat ions, sot hat i n it ial c ond itions  a reexact ly the sam e for  all r un s to e nsure  c omparabi lity. Following @HenBan 08, th e top  an d bottombou n dar y conditi onsare period ic.  Th e phy sic a l dim ensi o ns  of the c ompu t ational d o ma ina r e$20\times20 $ p c.The t woi nflows  ofgas collide in th e  middle of the  dom a i n.The supe rs onic collision  ca us e s  a steep  r ise in thegas dens it y  that  trigg ers th e therm a l  i n stabil ity,  an d gas und erg oe s  a tran si ti o n into  the  p hase o f thec oldn e utral material ( CNM)i n  theI SM.  In t hi s phase ,  the  gas has t emperatures  in th e ra nge $ 30$–$10 0$  K and  nu mb er densiti e s within
 the_hydrogen atom,_respectively. The gas is_assumed to_be_a perfect_gas_with $\gamma =_5/3$ and $\mu_= 1.4 m_H$.

The cooling_function of @AudHen05_includes_the cooling by fine-structure lines of CII and OI, further the cooling by H_(Ly$\alpha$-line)_and the_electron_recombination_onto positively charged grains. The_heating is due to the_photo-electric effect_on small grains and polycyclic aromatic hydrocarbons (PAH)_caused_by the far-ultraviolet_galactic background radiation. For more information about this cooling_function see @Wol95 [@WolMcK03; @Spi78; @BakTie94]_and @Hab68. The_pressure-equillibrium_curve_resulting from the cooling_function is plotted as black curve_in Figure \[fig:phasediag\_static5\]. For the numerical_solution of the fluid dynamical equations, we_used the radiative cooling routine implemented_into Enzo. For each hydrodynamical_time step,_the state variables are iterated_over several subcycles,_and the_resulting total energy_increment for the whole time step_is added.

For our_numerical study, the two-dimensional setup of_@AudHen05_and [@HenAud07] was_adopted_with_small modifications._The initial gas_content_corresponds to_the_warm neutral material (WNM) in the_ISM,_i. e., the temperature is $T=7100$ K, the_pressure is $P=7\times10^{-13}$ erg_cm$^{-3}$_and the number density_of neutral hydrogen is $n=0.71$_cm$^{-3}$. From the left and the_right boundaries,_warm gas_with identical thermodynamic properties flows into the computational domain, where the_cosine-shaped inflow velocity profile is modulated_with small perturbations, realised_by randomly_shifted_phases. These phase_shifts_are kept_constant for the different simulations, so that_initial conditions_are exactly the same for all_runs to ensure comparability._Following_@HenBan08, the top and bottom boundary_conditions are periodic. The physical dimensions_of the computational domain are_$20\times20$_pc._The two inflows of gas_collide in the middle of the_domain. The supersonic_collision causes a steep rise in the_gas_density that triggers the thermal instability,_and_gas undergoes a transition into the_phase_of_the cold neutral material (CNM)_in the ISM. In this phase,_the gas has temperatures in the range $30$–$100$ K_and number densities_within