there are basically two approaches to modeling the edge and scrape - off layer ( sol ) regions in magnetically confined plasmas : mean field transport models and turbulence models . in the first approach@xcite , mean field
fluid equations are solved in a realistic geometry typically including the divertor and x - point regions .
these models have very detailed descriptions of collisions , neutral particles , impurities , and material surfaces , whereas turbulent transport is not calculated self - consistently but modeled by effective turbulent diffusion terms .
turbulence models@xcite take a different approach . here ,
a detailed description of the turbulence is the main goal , whereas the level of detail of collisions , neutral particles , impurities and material surfaces , if included in the models at all , is crude compared with the highly detailed descriptions in mean field transport models .
furthermore , in order to reduce the computational requirements the magnetic field geometry is often simplified .
both approaches have provided useful results in their respective regimes of validity .
nevertheless , several open questions in contemporary fusion plasma research require that these complementary approaches are merged .
examples are : the formation of a elevated density shoulder@xcite as the particle density approaches the greenwald density limit and the associated changes of filament properties@xcite , the high particle density front@xcite on the high field side of single null diverted plasmas , and eventually the possible influence of turbulent transport on the transition to the detached@xcite divertor regime .    here
, we present a step towards more detailed edge and sol turbulence models by showing how classical collisional transport across the magnetic field is included in drift - fluid models@xcite based on the braginskii fluid closure@xcite . compared with neoclassical transport ,
classical collisional transport across the magnetic field is small in most regions of a magnetically confined fusion plasma . in turbulence dominated systems
collisional transport is important because of its diffusive nature which inevitably determines the smallest possible scale length of the turbulence below which the turbulent energy is dissipated .
it is well - known that turbulence behaves remarkably different when collisional dissipation is completely neglected compared with the situation where a very small but finite dissipation is present .
the results of neoclassical transport theory are concerned with the steady - state and does not take turbulence into account nor resolve the much faster turbulent time - scale . in the pfirsch - schl " uter regime ,
neoclassical transport is a result of collisional diffusive transport along the magnetic field in combination with guiding - center drifts due to magnetic field inhomogeneity . on the turbulent time - scale
, it is therefore not evident that neoclassical effects through parallel collisional transport is the dominant regulator of turbulence in particular when the spatial scale of the turbulence is small as in e.g. , drift wave turbulence .
a consistent model for perpendicular collisional transport may therefore be of importance for the interplay between collisions and turbulence and certainly provide a valuable validation tool for models which do not take classical perpendicular collisional transport into account .
moreover , classical perpendicular transport is believed to play a role for the perpendicular convection of filaments in the sol region .
a filament at the outboard midplane typically has a size which is on the order of @xmath0 ion gyroradii , but when then filament approaches the x - point region the perpendicular scale length can easily be smaller than the ion gyroradius due the strong magnetic shear near the x - point@xcite . in that case
the perpendicular conductivity and viscosity are important mechanisms for closing the filament current loop@xcite and thereby modify the filament convection and possibly decouple perturbations above and below the x - point .
finally , classical perpendicular transport are important in other fully ionized plasmas without neoclassical transport , e.g. cylindrical plasmas .    in drift - fluid models ,
the components of the momentum density equations perpendicular to the magnetic field are solved iteratively essentially under the assumption that the dynamics is evolving much slower than the ion gyration time scale . in this limit , every force density , including resistivity and collisional viscosity , gives rise to a perpendicular drift@xcite .
most of the collisional drifts are included in transport models@xcite but not in a way which provides energy conservation . in turbulence models , on the other hand , the collisional drifts are usually neglected or do only include selected collisional effects@xcite . in a model
including ion temperature dynamics , we show how these collisional terms in combination with heat fluxes , heat exchange terms , and viscous heating terms in the electron and ion pressure equations are included in drift - fluid turbulence models in a consistent way .
consistency here mainly refers to energy conservation . in this paper
we show that collisional transport across the magnetic field can be included in drift - fluid models without altering the energy theorem@xcite .
for example , we show that the collisional viscous damping of turbulence and mean flows gives rise to a conservative energy transfer between the kinetic energy and the ion thermal energy . in other words
, collisions can be added to drift - fluid models without introducing energy sinks or sources .
the inclusion of perpendicular collisional transport in the drift - fluid equation gives rise to diffusion of particle density , momentum , and ion and electron pressure , and hence potentially renders the common use of artificial diffusion terms in turbulence models redundant .
the collisonal terms are complex functions of the fluid and the electromagnetic field variables
. thus when included in turbulence models , they place significant demands on computing powers .
therefore , we also present a partly linearized model for the collisional transport across the magnetic field which allows significantly faster computing times .
the model is embedded in the hot edge - sol electrostatic ( hesel ) model@xcite , the successor of the esel model@xcite .
hesel is a two - dimensional ( 2d ) model describing interchange turbulence in a plane perpendicular to the magnetic field at the outboard midplane of a tokamak .
it evolves the particle density , the vorticity , and the electron and ion pressures .
we show that our partly linearised collision terms , including inter species energy exchange , do not alter the energy theorem , and we show by means of 1d numerical simulations without turbulence that the simplified and full models give very similar results .
the article is organized as follows . in section [ sec : perpcolltrans ] we describe how collisional effects influence the dynamics perpendicular to the magnetic field in a low - frequency drift - fluid turbulence model .
we derive the global energy theorem and discuss energy exchange mechanisms due to collisions . in section [ sec : heselmodel ] a simplified model for collisional effects is derived which is included in the hesel model also presented here .
numerical simulations of the full and simplified models for the perpendicular collisional effects are presented , and finally in sec .
[ sec : conclusions ] we summarise our findings and draw conclusions .
in this section we show how collisional transport across the magnetic field is consistently incorporated into drift - fluid models .
our derivation is based on the braginskii@xcite fluid equations describing the time evolution of particle density @xmath1 , velocity @xmath2 and temperature @xmath3 for a collisional , quasi - neutral @xmath4 , electrostatic , simple electron - ion plasma .
the index @xmath5 is a species label , which refers to ions and electrons , respectively .
we start our derivation from the braginskii momentum equation @xmath6\bm{u}_a = -\nabla p_a - \nabla \cdot \bm{\pi}_a + q_an_a(\bm{e } + \bm{u}_a\times \bm{b } ) + \bm{r}_a .
\label{eq : momentum}\end{aligned}\ ] ] here , @xmath7 denotes mass , @xmath8 is charge , @xmath9 is the scalar pressure , @xmath10 is the stress tensor , and @xmath11 denotes the resistive force .
the resistivity in the electron momentum equation consists of a frictional force @xmath12          \label{eq : ru}\end{aligned}\ ] ] and a thermal force @xmath13 where the electron - ion collision frequency is defined as @xmath14 the magnitude of the electron gyrofrequency is @xmath15 , and @xmath16 denotes the coulomb logarithm .
we also introduced a unit vector parallel to the magnetic field @xmath17 and the magnetic field aligned component of the gradient operator @xmath18 . if @xmath19 is an arbitrary vector , then we introduced the notation @xmath20 for the projection onto the magnetic field unit vector @xmath21 , and @xmath22 for the perpendicular part of the vector .
momentum conservation implies @xmath23 .
the origin of the thermal force @xmath24 is the velocity dependence of the particle collision frequency ( @xmath25 ) .
the stress tensor @xmath10 consists of three parts : i ) a part @xmath26 describing viscosity along the magnetic field due to like - particle collisons , ii ) a gyro - frequency dependent part @xmath27 describing collisional momentum transport across the magnetic field due to like particle collisions , and iii ) a gyro - viscous part @xmath28 . here
we consider the @xmath29 and @xmath30 parts of the ion stress tensor .
the electron stress tensor is neglected because it is smaller than the ion contributions by the electron - ion mass ratio . in the ion stress tensor
we assume that the perpendicular flow is incompressible and assume a constant magnetic field .
the pure perpendicular ion viscous tensor@xcite written in a local coordinate system @xmath31 , where @xmath32 is aligned with @xmath21 , then reads :    @xmath33    [ eq : visc_perp_coll ]    here the ion viscosity coefficient is given as @xmath34 and the ion - ion collision frequency is @xmath35 the parallel - perpendicular parts , e.g. , @xmath36 , do not contribute to the perpendicular dynamics to lowest order and are therefore not treated here since we do not consider the dynamics parallel to the magnetic field . in the parallel direction the parallel - perpendicular part of the stress tensor give rise to perpendicular viscosity . in the drift ordering @xcite the fundamental assumptions are that the characteristic frequencies are much smaller than the ion gyrofrequency @xmath37 and that perpendicular gradient length scales @xmath38 are larger than the ion gyro - radius @xmath39 we further assume @xmath40 , here @xmath41 denotes a characteristic parallel gradient length scale .
this ordering permits iterative solutions of the momentum equations for the perpendicular fluid drifts .
the first order drifts are given as @xmath42 which are the @xmath43  and diamagnetic drifts , respectively .
the second order drifts are @xmath44 in all second order drift we take @xmath45 .
the first term @xmath46 is the polarization drift . due to the mass dependence
only the ion polarization drift is retained .
the second term @xmath47 is the drift@xcite associated with resistivity , which is identical for electrons and ions due to momentum conservation .
the resistive drift consists of a friction force drift and a thermal force drift , respectively : @xmath48 here @xmath49 denotes the total scalar pressure .
note that only the diamagnetic current appears in the friction force drift @xmath50 because the electron and ion @xmath43-drift contributions cancel .
the diamagnetic drift is the fluid representation of the particle grad - b and curvature drifts , and the magnetization current associated with the gyrating charged particles@xcite .
therefore , the diamagnetic drift does not represent transport of guiding - centers unless the magnetic field is inhomogeneous . in that sense the friction force drift @xmath50 is a fluid representation of friction between electrons and ions in opposite directed larmor orbits .
the direction of thermal force drift @xmath51 is parallel to the electron temperature gradient .
as already mentioned , this drift arises because fast particles experience less collisions than slow particles , which result in unbalanced fluxes , and hence implies up - gradient transport .
the last first order drift is the viscous drift @xmath52 .
as earlier mentioned this term is mass dependent and hence only the ion drift is retained . in the local coordinate system @xmath31 the purely perpendicular parts can be written as @xmath53 being a second order drift only the first order drifts @xmath54 , see eq .
( [ eq : vperp0 ] ) , are retained in the stress tensor @xmath55 .      with the algebraic expressions for the perpendicular drifts at hand
, we can write down the resulting drift - fluid equations . we will omit the magnetic field aligned parts of the momentum equations since we are here mainly concerned with perpendicular dynamics . by inserting
the perpendicular drifts given in equations ( [ eq : vperp0 ] ) and ( [ eq : vperp1 ] ) into the braginskii@xcite particle density and pressure equations , we get :    @xmath56\big )        + p_e \nabla \cdot [ \bm{u}_e + \bm{u}_{de } + { \bm{\hat{b}}}u_{\|e}+\bm{u}_r ] \notag    \\
+ \nabla \cdot \bm{q}_{e }         + \nabla_{\perp } \cdot \bm{q}_{e}^ * = -q_{\delta } + \bm{r } \cdot ( \bm{u}_i-\bm{u}_e ) \label{eq : df_pe}\\      \frac{3}{2}{\frac{\partial}{\partialt}}p_i        + \frac{3}{2}\nabla \cdot \big(p_i[\bm{u}_e + \bm{u}_{di } + { \bm{\hat{b}}}u_{\|i } + \bm{u}_{pi } + \bm{u}_{\pi i } + \bm{u}_r]\big)\notag \\
+ p_i \nabla \cdot [ \bm{u}_e + \bm{u}_{di } + { \bm{\hat{b}}}u_{\|i } + \bm{u}_{pi } + \bm{u}_{\pi i } + \bm{u}_r ]               + \nabla \cdot \bm{q } _ { i }         + \nabla_{\perp } \cdot \bm{q}_{\perp i}^ *        + \bm{\pi}_i^{\perp } : \nabla \bm{u}_{i\perp,1 }   = q_{\delta}.   \label{eq : df_pi}\end{aligned}\ ] ]    equation ( [ eq : df_w ] ) is the vorticity equation which results from subtracting the ion and electron particle density equations assuming quasi - neutrality . in all equations higher order terms
are retained even though some of these terms are formally small .
the motivation for not neglecting these higher order terms is that they , e.g. , guarantee energy conservation , provide diffusion or are responsible for energy exchange between different plasma species .
drift - fluid equations have been derived by several authors@xcite but without self - consistently including perpendicular collisional transport . in the remainder of this section we describe and discuss the collisional terms entering the drift - fluid moment equations .
thorough descriptions and discussions of classical transport outside the drift - fluid context is found in e.g. refs .  .
we start with the electron particle density equation ( [ eq : df_n ] ) where the divergence of the resistive flux is the only collisional term . using equations ( [ eq : uru ] ) and ( [ eq : urt ] ) , the divergence of the resistive flux can be written as @xmath57 \bigg ) , \label{eq : resflux}\end{aligned}\ ] ] demonstrating the origin of the resistive drifts , namely a perpendicular random - walk process with step length @xmath58 and frequency @xmath59 .
we see that the perpendicular friction force @xmath60 gives rise to particle density diffusion with a diffusion coefficient @xmath61 in drift - fluid models , similar to the explicit appearance of spatial diffusion in low - frequency ordered collision operators@xcite .
the thermal force drift @xmath51 is partly canceled by the electron temperature gradient dependent part of the friction force drift @xmath50 , leaving only the last term on the right hand side . in order to highlight the @xmath62 and @xmath63 dependency of the electron - ion collision frequency @xmath64 we define @xmath65 where @xmath66 and @xmath67 denote a constant reference particle density and a constant reference electron temperature , respectively . equation ( [ eq : resflux ] ) then reads : @xmath68\bigg ) , \label{eq : resfluxii}\end{aligned}\ ] ] showing that the particle density diffusion coefficient has a @xmath69 dependence .
next , consider the electron pressure equation ( [ eq : df_pe ] ) . here , collisional effects enter through classical heat fluxes , heat transfer terms , and resistive drift terms .
the collisional electron heat flux@xcite consists of two parts : @xmath70 where the heat conductivities@xcite are @xmath71 and where only the first order perpendicular drifts given in eq.([eq : vperp0 ] ) were used in @xmath72 .
the origin of the heat flux @xmath72 is the same as the thermal force @xmath73 given in eq .
( [ eq : rt ] ) , which are both associated with the velocity dependence of the electron - ion collision frequency .
note the exact onsager@xcite symmetric coefficients in @xmath73 and @xmath72 .
it is convenient to express the perpendicular parts of the electron heat flux in terms of the perpendicular resistive drifts given in equations ( [ eq : uru ] ) and ( [ eq : urt ] ) @xmath74 when written in this form it is evident that the divergence of the frictional electron pressure flux cancels the perpendicular part of @xmath72 : @xmath75 on the right hand side of the electron pressure equation ( [ eq : df_pe ] ) , the @xmath76 term transfers energy between the electron and ion fluid kinetic energies an the electron thermal energy . in the drift - fluid electron pressure equation
the energy transfer term can be combined with the remaining perpendicular heat flux terms and terms which depend on the resistive drift @xmath77 the two first terms on the right hand side , @xmath78 \bigg ) ,        \label{eq : pe_diffii }   \end{aligned}\ ] ] among others , yield diffusion of electron pressure and in general relax gradients , and hence the up - gradient flux contained in @xmath79 is canceled by perpendicular heat conduction in @xmath80 . the last term in eq .
( [ eq : pe_diff ] ) transfers energy between the electron and ion thermal energy due to resistivity .
hence , no explicit energy transfer channel between kinetic and thermal energy due to resistivity remains in the drift - fluid equations .
energy exchange between comoving electron and ion fluids is described by the heat exchange term @xmath81 entering the electron and ion pressure equations ( [ eq : df_pe ] ) and ( [ eq : df_pi ] ) , respectively .
the equilibration occurs on the slower @xmath82 collision time - scale because the energy transfer in each scattering event is proportional to the mass ratio of the scattering particles @xmath83 .
we elaborate further on inter - species energy exchange in section [ sec : energyconservation ] .
we now turn to the ion pressure equation ( [ eq : df_pi ] ) where the heat conduction is given as @xmath84 and the thermal conductivities@xcite are @xmath85 the ion heat flux is solely driven by ion temperature gradients and is independent of ion - electron collisions .
the ratio between the perpendicular electron and ion heat conductivities is of the order @xmath86 for @xmath87 , whereas @xmath88 .
perpendicular collisional heat transport is therefore dominated by the ions , and parallel heat transport is electron dominated .
the terms involving the resistive drift in the ion pressure equation can be written as @xmath89\bigg ) - \bm{u}_r \cdot { \nabla_{\perp}}p_i ,      \label{eq : pi_diff}\end{aligned}\ ] ] which shows that electron - ion collisions through the resistive drift give rise to ion pressure diffusion , but also that the thermal force drift @xmath51 potentially drives the ion pressure up the electron temperature gradient .
note that the corresponding effect in the electron pressure equation was canceled by the electron heat conduction .
the diffusion coefficient is proportional to the electron diffusion coefficient times the ratio of the ion to the electron temperature .
the ratio of the resistive diffusion and the perpendicular ion heat conduction can be estimated as @xmath90 where @xmath91 and @xmath92 are characteristic electron and ion temperature gradient length scales .
perpendicular heat conduction therefore dominates resistive diffusion except for plasmas where the ion temperature is significantly higher than the electron temperature or where the electron temperature has much steeper gradients than the ion temperature .    in both pressure equations the diamagnetic heat flux
is defined as @xmath93
lastly , we consider viscous effects in the vorticity and ion pressure equations ( [ eq : df_w ] ) and ( [ eq : df_pi ] ) , respectively . in both equations
viscous effects enter through the viscous drift @xmath94 defined in eq .
( [ eq : vperp1 ] ) , but the viscous tensor only explicitly appears in the ion pressure equation .
notice that only the divergence of the viscous particle density flux remains in the vorticity equation .
the resistive electron and ion fluxes cancel as a consequence of momentum conservation .
the divergence of the perpendicular viscous flux results in diffusion of the magnetic field aligned @xmath43  and ion diamagnetic vorticities due to ion - ion collisions , as we will show in section [ sec : heselmodel ] . in the ion pressure equation
the corresponding ion pressure flux gives rise to hyperviscosity `` @xmath95 '' of the ion pressure .
however , compared with the perpendicular ion heat conduction the viscous ion pressure flux term is in most cases small @xmath96 the viscous drift and the viscous tensor provide an energy transfer channel between ion thermal energy and kinetic energy through the @xmath97 and the @xmath98 terms .
the underlying energy transfer mechanism is the randomization of ordered perpendicular fluid motion due to ion - ion collision .
as we will show in the next section , both terms are important for the conservation of energy .
the gyro - viscous part vanishes@xcite exactly @xmath99 , when the gyro - viscous tensor @xmath100 is evaluated with the first order drift velocities @xmath101 .      in this section explicit expressions for the polarization and gyro - viscous drifts
are given .
these drifts enter the _ gyro - viscous _ cancellation which eliminates the advection of vorticity by the diamagnetic drift in the polarization equation ( [ eq : df_w ] ) , but also the advection of ion parallel momentum and ion heat fluxes by the diamagnetic drift in their corresponding moment equations@xcite .
these equations are significantly altered by the gyro - viscous tensor .
the gyro - viscous cancellation is not complete in the sense that numerous small correction terms remains in the vorticity equation , in the parallel momentum equation , and in the heat flux equation .
the actual derivation is cumbersome and the level of complication rises with the level of detail included in the polarization drift terms , e.g. , the inclusion of the polarization heat flux and anisotropic pressure@xcite .
here , we aim at formulating a workable model for the use in numerical turbulence simulations . therefore , we do not include the polarization heat flux and we do not account for anisotropic pressures .
furthermore , we bring the magnetic unit vector under the material derivative in the polarization drift @xmath102 and neglect the corresponding correction terms in the gyro - viscous drift . this approximation yields simpler expressions compared to the full expressions in e.g. refs .  .
as already discussed , we leave out parallel momentum for the purpose of exposition . after carrying out the gyro - viscous cancellation@xcite ,
the combined divergences of the ion polarization and gyro - viscous fluxes are : @xmath103 where we introduced the short - hand notation @xmath104 and where the material derivative is defined as @xmath105 .
the advecting velocities are @xmath106 here the drift @xmath107 represents the remainder of the gyro - viscous cancellation .
@xmath108 is derived from the scalar function @xmath109 given in eq .
37 in ref .
the scalar function @xmath109 enters the momentum equation as a correction to the scalar pressure and hence only contributes when the magnetic field is inhomogeneous .
the contributions originating from @xmath109 are therefore small . here
, we choose to keep only the terms in @xmath109 necessary for energy conservation .
furthermore , terms in @xmath109 originating from the purely perpendicular heat flux are neglected because energy conservation requires that these contributions are retained together with the polarization heat flux which is not kept here .    as a rule of thumb@xcite
all ion drifts retained in the ion continuity equation , except the ion diamagnetic drift , must be kept in the advection part of the polarization drift in order to conserve energy .
unfortunately , this leads to a recursive definition of the polarization drift since the polarization drift includes a polarization drift advection term .
if exact energy conservation is a crucial demand , this feature makes the model unsuitable for numerical calculations , but the model is convenient when discussing energy conservation in drift - fluid models in general .      in this section
we present the energy theorem for the drift - fluid model eqs .
( [ eq : df_n])-([eq : df_pi ] ) . without taking collisional effects into account
it has previously been shown , e.g. , refs .  , that drift - fluid models conserve energy . here , the drift - fluid energy theorem is presented including the perpendicular collisional transport terms and the collisional energy exchange terms described in section [ sec : moments ] .
the global energy is obtained by multiplying the vorticity equation ( [ eq : df_w ] ) by @xmath110 and integrate over all space neglecting surface terms @xmath111\notag \\
=     \int dv \ , & nm_i \frac{{\nabla_{\perp}}\phi}{b } \cdot   d_t^\ #   { \nabla_{\perp}}\phi^ * +    p_i { \nabla \times \frac{{\bm{\hat{b}}}}{q_ib } } \cdot \nabla { \nabla_{\perp}}\phi^ * \cdot \frac{{\nabla_{\perp}}\phi}{b }      -\bm{u}_e \cdot \nabla p   + ( \nabla \cdot \bm{\pi}^{\perp}_i ) \cdot \bm{u}_e .
\label{eq : energy_kin}\end{aligned}\ ] ] terms involving parallel dynamics are everywhere left out .
detailed descriptions of energy conservation and energy transfer channels in the parallel direction in the absence of perpendicular collisional effects is thoroughly described in refs .   and .
next , the ion and electron pressure equations are integrated over space using eqs .
( [ eq : electron_heatflux_cancel ] ) , ( [ eq : pe_diff ] ) , and ( [ eq : pi_diff ] ) , again neglecting surface terms @xmath112 where @xmath113 and @xmath108 were inserted into the @xmath114 term , resulting in the second and third terms on the right hand side of eq .
( [ eq : energy_pi ] ) , respectively .
the energy theorem @xmath115 is obtained by adding eqs .
( [ eq : energy_kin])-([eq : energy_pi ] ) and by bringing the ion particle density under the time derivative in terms involving the ion polarization drift using the particle density equation ( [ eq : df_n ] ) and the vorticity equation ( [ eq : df_w ] ) .
as expected , the energy theorem is not altered by the inclusion of perpendicular collisional effects . in other words
, the collisional effects do not give rise to sinks no sources in the energy theorem , but solely provide energy transfer channels .
in that respect the energy theorem is a sanity check of the model and can be used as a verification tool in numerical simulations .
the perpendicular kinetic energy @xmath116 is a peculiar quantity .
the @xmath43-drift @xmath117 part is easy to interpret , because the fluid moments in eqs .
( [ eq : df_n])-([eq : df_pi ] ) are all advected by the @xmath43-drift , and hence this part describes the kinetic energy associated with the @xmath43-drift . on the contrary , the terms in the kinetic energy @xmath116 involving the diamagnetic drift are more difficult to interpret because the diamagnetic drift to lowest order does not advect fluid elements .
the diamagnetic part of the kinetic energy is more appropriately described as a finite larmor radius ( flr ) correction to the ion pressure and to the ion @xmath43-kinetic energy@xcite . to emphasize this property of the kinetic energy , it is denoted the _ modified kinetic energy_.    to elucidate how energy is transferred between the constituents of the energy theorem eq .
( [ eq : cons_energy ] ) it is instructive to consider the time evolution of each term separately @xmath118 terms which enter with opposite signs represent conservative energy transfer channels .
the interchange transfer terms , first terms on the right hand sides of eqs .
( [ eq : ddt_ekin])-([eq : ddt_epe ] ) , transfer energy between ion and electron thermal energies and the modified kinetic energy . the compression of the polarization drift , second term on the right hand sides of eqs .
( [ eq : ddt_ekin ] ) and ( [ eq : ddt_epi ] ) , provides an energy transfer between ion thermal energy and the modified kinetic energy through the @xmath119 terms .
this energy transfer channel may play an important role in generating and sustaining @xmath43  mean flows@xcite .    the heat exchange term @xmath120 , entering eqs .
( [ eq : ddt_epe ] ) and ( [ eq : ddt_epi ] ) with opposite signs , describes energy exchange between electrons and ions due to elastic electron - ion collisions .
resistivity provides an additional energy transfer channel through the resistive drift transfer terms @xmath121 entering eqs .
( [ eq : ddt_epe ] ) and ( [ eq : ddt_epi ] ) with opposite signs .
the parallel resistivity , on the contrary , gives rise to an energy exchange between the ion and electron parallel kinetic energies @xmath122 and the electron thermal energy . here
, resistivity opposes the motion of the charge carrying particles along the magnetic field . in the perpendicular direction
, resistivity opposes the diamagnetic current , which to lowest order describes particle gyration around magnetic field lines .
since the kinetic energy of the gyrating motion is included in the thermal energies , perpendicular resistivity in the drift - fluid model accordingly give rise to a conservative energy between the ion and electron thermal energies .
note that a finite resistive energy transfer requires a finite ion pressure gradient .
when @xmath123 , ions are homogeneously distributed and therefore the probability of an electron colliding with a gyrating ion is the same in all directions . only when @xmath124 , electrons do on average experience an average effect of the gyrating ions .
the direction of the energy transfer can be deduced by writing out the transfer term in the electron pressure equation explicitly @xmath125 \cdot { \nabla_{\perp}}p_i ,            \label{eq : r_e_trans}\end{aligned}\ ] ] where the thermal force and the electron temperature gradient dependent part of the friction force were combined .
the second and third terms on the right hand side heat the electrons whenever @xmath126 .
the heating of electrons due to friction is known as joule heating .
however , it can be shown that the second and third terms are also capable of cooling the electrons .
necessary conditions are @xmath127 and that the ion temperature gradient length scale exceeds the particle density gradient length scale .
the electron temperature gradient dependent term transfers thermal energy from the electron to the ions when @xmath128 , an effect known as the ettingshausen effect@xcite .
it is important to keep in mind that both the heat exchange term @xmath120 and the resistive transfer terms occur on the slow time - scale @xmath129 .
comparing the magnitudes of these energy exchange terms @xmath130 shows that unless the temperatures are very different , the resistive energy transfer is smaller than the heat exchange when considering profiles , but also that the resistive can be important when pressure gradients are steep , e.g. , in the edge region and in general when turbulent transport is strong .
viscosity on the perpendicular flows due to ion - ion collisions gives rise to the energy transfer term @xmath131 between the modified kinetic energy and the ion thermal energy .
we should expect that the energy transfer goes from kinetic energy to thermal energy because collisions increase entropy .
if the perpendicular electric field vanishes the energy transfer term vanishes accordingly .
it can be shown that the @xmath132 dependent part of @xmath133 entering the transfer term gives rise to an unidirectional energy transfer from modified kinetic energy to the ion thermal energy .
the direction of the energy transfer by the remaining part , which depends on the ion diamagnetic drift , is not unidirectional .
however , recall that the modified kinetic energy does not solely describe kinetic energy . therefore , in order to investigate the direction of the energy transfer one must instead consider the energy transfer from the ( flr ) corrected ion @xmath43  kinetic energy .
disentangling the true kinetic energy from the modified kinetic energy is beyond the scope of this work .
turbulent transport in the sol region is intermittent and is carried by coherent structures@xcite .
these coherent structures of elevated pressure , often referred to as filaments or blobs , are elongated along the magnetic field .
they are mainly born in the vicinity of the last closed flux surface ( lcfs ) at the outboard midplane due to ballooning@xcite . once these plasma outbreaks are expelled into the the open field line region
, they expand@xcite into the rare sol plasma in the parallel direction with velocities comparable to the ion sound speed @xmath134 , while traversing the sol region radially at velocities reaching a fraction of @xmath135 .
typically , the blobs imply fluctuation levels which exceed unity .    in this section we derive a reduced model for investigations of blob formation and transport across the sol region of magnetically confined toroidal plasmas .
the model is named hesel@xcite and is the successor of the esel model@xcite . like the esel model
, the hesel model considers the dynamics in a plane perpendicular to the magnetic field located at the outboard midplane .
perpendicular collisional transport is described by partly linearized expressions derived from the more general results presented in the preceding section .
parallel transport in the sol region is represented by damping terms .
parallel drift - wave dynamics and magnetic perturbations are ignored .
turbulent fluctuations are hence driven by the interchange instability due to an inhomogenious toroidal magnetic field .
these reductions facilitate long but also computationally inexpensive computer simulations which can provide sufficient data for statistical analysis and comparison with experiment . compared to the esel model
the hesel model has the following enhancements : 1 ) a dynamical description of ion pressure dynamics , 2 ) the inertial response to changes in the ion pressure is included , 3 ) sheath dissipation , and 4 ) an improved description of collisional processes , including inter - species energy exchange , based on the results presented in the previous section .
the model domain is a plane perpendicular to the magnetic field at the outboard midplane .
the magnetic field geometry is approximated by local slab coordinates @xmath31 .
@xmath136 and @xmath137 coordinates denote the radial and the poloidal azimuthal positions in the domain and @xmath138 denotes the position along the magnetic field . the toroidal magnetic field magnitude is approximated as @xmath139 , where @xmath140 denotes the inverse aspect ratio , and @xmath141 and @xmath142 are the minor and major radii , respectively .
the hesel model consists of two regions joined at the lcfs ( @xmath143 ) modeling the edge and sol regions . besides radial boundary conditions ,
these regions are only set apart by damping terms due to losses along open field lines which are only active in the sol region . in this section
we present and discuss the simplifications of the drift - fluid model equations ( [ eq : df_n])-([eq : df_pi ] ) that lead to the hesel model .
first we turn our attention to the vorticity equation ( [ eq : df_w ] ) .
here we invoke an approximation commonly known as the thin - layer approximation .
this approximation shares many features with the boussinesq approximation often used in the description of incompressible , buoyancy - driven flows in neutral fluid dynamics , but is fundamentally different since the mechanism driving the dynamics in plasma drift - fluid equations is finite compressibility of the plasma fluid drifts .
this rather crude approximation neglects particle density variations in the polarization flux entering the vorticity equation ( [ eq : df_w ] ) , and hence assumes a constant inertia of all fluid parcels irrespective of the local particle density . the motivation for this approximation , which in particular is debatable when applied to edge plasma models , is the desire for computational efficiency and that energy conservation requires polarization drift advection in the polarization drift itself as discused in section [ sec : pola_gyr_canc ] .
this recursive definition is very inconvenient in numerical computations .
we note that gyrofluid models@xcite do not suffer from these difficulties .    specifically , we approximate the ion polarization flux defined in eq .
( [ eq : gv_cancel ] ) as @xmath144)=           -\nabla \cdot \bigg(\frac{n}{\omega_i } d_t^ \ # { \nabla_{\perp}}\phi^ * \bigg )          + \nabla \cdot ( n\bm{u}_{\chi } )      \simeq    -\frac{n_0}{\omega_{i0 } } \nabla \cdot ( d_t^0 { \nabla_{\perp}}\phi^*\big ) ,         \label{eq : hesel_tl}\end{aligned}\ ] ] where we have introduced @xmath145 and where subscripts `` @xmath146 '' on quantities on the right hand side refer to constant characteristic reference values .
the material derivative @xmath147 only includes the @xmath43-drift at constant magnetic field @xmath148 . dictated by energy conservation , all other advection terms in the original polarization flux
are neglected .
energy conservation requires that the divergence of the polarization drift in the ion pressure equation ( [ eq : df_pi ] ) is approximated in the same way @xmath149 in the final equations the divergence of the polarization drift does not enter the ion pressure equation , because the vorticity equation is substituted into the ion pressure equation .
the divergence of the polarization flux @xmath150)$ ] is formally small and does not affect energy conservation .
it is therefore neglected .
we start by considering the particle density equation ( [ eq : df_n ] ) where the only perpendicular collisional term is the divergence of the resistive flux written out in eq .
( [ eq : resfluxii ] ) .
as previously noted , the divergence of the resistive flux gives rise to particle density diffusion with a diffusion coefficient @xmath151 .
naturally , the particle density and electron temperature fields in the edge and sol regions of a tokamak plasma will not be completely correlated .
however , we expect that the profiles of both fields share features such as being flat in the sol region and having steep gradients near the separatrix .
furthermore , measurements in the sol region at the outboard midplane show that @xmath62 and @xmath63 are approximately proportional in blobs@xcite in which case @xmath152 .
therefore , we choose to fix the diffusion coefficient in eq .
( [ eq : resfluxii ] ) evaluating it using characteristic constant reference values @xmath66 and @xmath67 for the particle density and the electron temperature , respectively .
recall that the thermal force was partly canceled leaving only the last term in eq .
( [ eq : resfluxii ] ) .
this residual thermal force is expected to be partly canceled by the ion temperature gradient part whenever the ion and electron temperatures are correlated .
therefore , we neglect the temperature gradient terms in the resistive flux eq .
( [ eq : resfluxii ] ) .
the divergence of the perpendicular resistive flux in the particle density equation ( [ eq : df_n ] ) then reads @xmath153 where @xmath154 is the constant diffusion coefficient calculated using @xmath66 and @xmath67 , and @xmath155 is the reference ion to electron temperature ratio .
usage of the same approximate resistive drift @xmath156 in both the electron and ion pressure equations ( [ eq : df_pe])-([eq : df_pi ] ) avoids undesired temperature dynamics such as temperature anti - diffusion .
next we consider the electron pressure equation ( [ eq : df_pe ] ) . as was shown in eq .
( [ eq : pe_diff ] ) , the resistive fluxes and the heat fluxes partially cancel .
apart from the energy exchange terms , the remaining terms were given in eq .
( [ eq : pe_diffii ] ) .
these terms are approximated as @xmath157 \bigg )     \notag   \\      & \simeq \nabla \cdot ( p_e \bm{u}_{r0 } )   - \nabla \cdot \bigg(d_e \frac{11}{12 } n{\nabla_{\perp}}t_e \bigg ] \bigg )   .
\end{aligned}\ ] ] the resistive energy transfer terms entering the electron and ion pressure equations are approximated as @xmath158    in the ion pressure equation ( [ eq : df_pi ] ) the collisional transport has a different character .
here , the perpendicular ion heat conduction dominate the divergence of the resistive flux .
however , the resitive flux is retained in order to handle situations with flat ion temperature gradients properly , see eq .
( [ eq : ion_cond_vs_resi ] ) . as in the electron pressure equation the divergence of the ion pressure resistive flux
is approximated as @xmath159 the perpendicular ion heat conduction eq .
( [ eq : ion_heatflux ] ) @xmath160 is approximated by taking the ion diffusion coefficient @xmath161 to be constant .    in the electron and ion pressure equations ( [ eq : df_pe])-([eq : df_pi ] ) the heat exchange term defined in eq .
( [ eq : energy_equilibration ] ) is evaluated with a constant electron - ion collision frequency @xmath162    finally , we consider approximations of terms related to the perpendicular collisional ion viscosity tensor @xmath29 , given in equation ( [ eq : visc_perp_coll ] ) , which enters the vorticity and ion pressure equations . in the vorticity equation ( [ eq : df_w ] ) ,
@xmath29 gives rise to diffusion of vorticity via the divergence of the flux associated with the viscous drift @xmath163 .
recall that in order to reduce the computational costs , the thin - layer approximation to the polarization flux was invoked as described in section [ sec : thinlayer ] .
the thin - layer approximation implied that the first order perpendicular ion drifts entering the polarization drift itself were approximated as @xmath164 accordingly , only @xmath165 enters @xmath29 .
expressed in local coordinates @xmath31 we get @xmath166 furthermore , the thin - layer approximation takes the particle density in the ion polarization flux to be constant . making the same approximation to the viscous drift flux
, we get @xmath167 which demonstrates that the perpendicular viscous drift gives rise to diffusion of the modified vorticity @xmath168 .
vorticity diffuses due to ion - ion collisions which randomize the ions and hence transfer energy from the modified kinetic energy , see eq .
( [ eq : ddt_ekin ] ) , to the ion thermal energy through the terms @xmath169 and @xmath170 entering the ion pressure equation ( [ eq : df_pi ] ) .
energy conservation requires that the perpendicular drifts entering @xmath29 and @xmath163 are approximated according to equation ( [ eq : uperpzero ] ) : @xmath171,\\        p_i\nabla \cdot \bm{u}_{\pi^{\perp } }         \simeq p_i\nabla \cdot \bm{u}^0_{\pi^{\perp}_i }         & = p_i\frac{3}{10}\frac{d_i}{\omega_{i0}}\nabla_{\perp}^2 \nabla_{\perp}^2 \phi^*.        \label{eq : visc_flux_comp_pi}\end{aligned}\ ] ] the @xmath172 term in the ion pressure equation is formally small , see eq .
( [ eq : viscpiflux_to_qi ] ) .
this term does not change the energy theorem and is therefore neglected .      in a toroidal plasma perpendicular collisional diffusion
is enhanced by neoclassical transport .
specifically , in a collisional plasma , @xmath173 , the transport along the magnetic field is diffusive ; here @xmath174 denotes the mean free path and @xmath41 is the parallel length scale . despite being interrupted by collisions ,
guiding - centers drifts , such as the grad - b drift , give rise to transport across magnetic flux surfaces . in a tokamak ,
the direction , inwards or outwards , depends on whether particles are above or below the outboard midplane .
when the diffusive transport along the magnetic field and the perpendicular guiding center drifts are superimposed they give rise to an enhanced pfirsch - schlter@xcite perpendicular diffusion @xmath175 , where @xmath176 is the safety factor .
neoclassical transport theory assumes that the plasma evolves on a time - scale much slower than the typical time - scale of turbulent transport .
therefore , neoclassical estimates of perpendicular transport do not immediately apply to intermittent transport in the sol region .
furthermore , existing literature is only concerned with neoclassical transport on closed field lines . here , we are interested in local transport on closed _ and _ open field - lines .
however , it is evident that the physical mechanism driving pfirsch - schlter transport is also present on open field lines and that it will influence profiles as well as transients , but since a rigorous derivation of these effects is beyond the scope of this work , we resort to the existing estimates keeping in mind that these will only hold approximately .
specifically , diffusion coefficients are multiplied by a common factor @xmath177 which then include variations of the inverse aspect ratio @xmath178 at the outboard midplane . since the safety factor becomes infinite at the lcfs and is not defined on open field lines we estimate the effect of a finite magnetic field line pitch by the safety factor at @xmath179 of the poloidal magnetic flux @xmath180 .
a more detailed discussion is given in fundamenskii et al ref .  .
it is important to note that since the perpendicular diffusive particle flux is enhanced by neoclassical effects , so is the resistive energy transfer , second last term in eq .
( [ eq : ddt_epe ] ) , between the electron and ion thermal energies@xcite .
the hesel model is 2d and therefore parallel losses on open field lines are parametrized .
the parametrization is in line with the parametrization in esel model given in ref .  , but it is augmented by parametrization of parallel ion pressure dynamics and the addition of a damping term representing sheath currents at material surfaces such as divertors and limiters . as already mentioned , the primary aim of the hesel model is to describe interchange driven turbulence transport at the outboard midplane . as a consequence of ballooning the turbulent plasma source mainly resides on the high field side at the outboard midplane .
for this reason , we assume that all parametrized parallel terms in the hesel model act as sinks .
blobs expand@xcite in the parallel direction with a velocity comparable to the ion sound speed @xmath181 .
the expansion depletes the particle density in the model - domain located at the outboard midplane . under the assumption that both ends of the blob expand with the same velocity , we get the characteristic parallel particle density damping rate @xmath182 where @xmath183 denotes the mach number and @xmath184 is the parallel blob size .
we assume that blobs are predominantly born on the high field side in a region centred around the outboard midplane with a @xmath185 poloidal extend@xcite , and hence approximate the filament blob size as @xmath186    in the sol region vorticity losses are due to two mechanisms : first , parallel advection of vorticity enters through the divergence of the ion polarization drift particle density flux @xmath187 .
since vorticity is mainly transported and generated by blobs in the sol region , we assume that the parallel gradient length scales of vorticity and particle density are approximately the same . in that case , the damping rate parametrising parallel advection of vorticity is then taken to be identical to the particle density damping rate : @xmath188 secondly , the parallel current in the vorticity equation provides an additional path for current density generated by the diamagnetic current in blobs . the role played by the parallel current
is roughly determined by the parallel conductivity , the magnetic topology particularly in the vicinity of the x - point , the sheath in front of material surfaces e.g. the divertor in most tokamaks , and the ratio of the characteristic time @xmath189 it takes a blob to move a distance in the radial direction comparable to its own size to the characteristic time @xmath190 of the parallel communication of current perturbations between the outboard midplane and the sheath . here , we consider blobs in the inertial regime where radial blob convection is barely inhibited by parallel currents .
this regime is typically relevant for plasmas with a high particle density , a long connection length @xmath191 between the outboard midplane and the divertor , or a low conductivity e.g. , due to high neutral gas densities in the sol region@xcite .
nevertheless , we assume that the divertor sheath influences the background profile of the electrostatic potential .
an average in the parallel direction@xcite of the parallel current leaves us with the value of the current at the sheath entrance given by the bohm criterion@xcite . for this reason we only apply the sheath damping to the averaged fields @xmath192
,      \label{eq : sheathdamping}\end{aligned}\ ] ] where @xmath193 is the bohm potential and the azimuthal average is defined as @xmath194 here , @xmath195 is the domain size in the @xmath137-direction and @xmath196 is an arbitrary function .
it is straightforward to use other closures of the parallel current e.g. , the sheath - connected regime@xcite can be investigated by applying the sheath damping not only to the averaged fields but to the fluctuating parts as well .
we now turn to the electron pressure equation ( [ eq : df_pe ] ) . here , the parallel heat fluxes given in eqs .
( [ eq : electron_heatflux_qu])-([eq : electron_heatflux_qt ] ) and the parallel advection term are parametrized .
quasi - neutrality only allows small deviations between the ion and electron parallel velocities .
consequently , the thermal heat flux @xmath197 given in eq .
( [ eq : electron_heatflux_qu ] ) is therefore small in comparison with the divergence of the parallel pressure flux and is therefore neglected in the model .
the spitzer - hrm@xcite parallel electron heat conduction @xmath198 given eq .
( [ eq : electron_heatflux_qt ] ) , on the other hand , is not negligible .
it has a strong electron temperature dependence @xmath199 , which attenuates parallel electron temperature gradients .
keep in mind that the braginskii closure is only valid in the collisional regime .
in particular @xmath198 is not well behaved as we approach the `` collisionless '' regime .
this is readily seen when writing the heat conductivity in terms of the mean free path @xmath200 @xmath201 as a rule of thumb the braginskii heat conduction is therefore only valid as long as the electron collisionality is sufficiently high ( @xmath202 ) . at lower collisionalities the spitzer - hrm heat conduction exceeds the free - streaming heat flux@xcite .
a remedy for this shortcoming is to introduce _ heat - flux - limiters _ , which in an ad - hoc manner limit the parallel heat flux to a fraction of the free - streaming value .
here we restrict ourselves to collisional sol conditions where the spitzer - hrm heat conductivity is valid .
conduction is then parametrized as @xmath203 where we introduced the spitzer - hrm damping rate @xmath204 in the parametrization of the electron heat conduction we take the connection length @xmath205 as the parallel gradient scale length and not the parallel filament size @xmath184 to avoid an unphysical damping of the electron temperature background .
parallel advection of electron and ion pressure are approximated as a one - dimensional adiabatic expansion ( @xmath206 ) @xmath207    in the ion pressure equation parallel heat conduction is neglected . for the parameters used here
( ion collisionality @xmath208 and @xmath209 , see fig.14 in ref . ) , parallel ion heat transport is dominated by parallel advection .
we are now ready to write down the full set of equations which constitute the hesel model . in order to highlight characteristic quantities and for convenience , the model equations are gyro - bohm normalized @xmath210 where @xmath211 is the characteristic ion gyro frequency , @xmath212 is the cold - ion hybrid thermal gyro - radius , @xmath66 , @xmath67 , and @xmath213 are characteristic reference values for the particle density and the temperatures .
the full system of equations is :    @xmath214    where the advective derivatives are defined as @xmath215 @xmath43-advection is written in terms of the anti - symmetric bracket @xmath216 the curvature operator is defined as @xmath217 and the normalized modified potential reads @xmath218 in the ion pressure equation , the divergences of the polarization and viscous drifts were eliminated by substitution of the vorticity equation .
diamagnetic advection in the electron and ion pressure equations is no longer present as a consequence of the diamagnetic cancellation @xmath219 all terms related to diffusion and parallel damping are grouped on the right hand sides of eqs .
( [ eq : n_norm])-([eq : ionpressure_thinlayer ] ) and are given as @xmath220\\            \lambda_{p_i } =         & \frac{5}{2}d_e ( 1+\tau ) \nabla \cdot ( t_i{\nabla_{\perp}}n )      -d_e ( 1+\tau){\nabla_{\perp}}\ln n \cdot { \nabla_{\perp}}p_i       + 2d_i\nabla \cdot ( n{\nabla_{\perp}}t_i ) \notag \\      & \frac{3}{10}d_i \big [           ( \partial^2_{xx}\phi^ * - \partial^2_{yy}\phi^*)^2 + 4 ( \partial^2_{xy } \phi^*)^2      \big ]      + \frac{3 m_e}{m_i}\nu_{ei0}(p_e- p_i )         + \sigma(x)\bigg (       p_i \mathcal{s }          - \frac{9}{2}\frac{p_i}{\tau_n}\bigg )    , \end{aligned}\ ] ] where @xmath155 and @xmath221\end{aligned}\ ] ] is a `` smooth step function '' which defines the transition from closed to open field line regions in the model .
the normalized diffusion coefficients are @xmath222 and the damping rates due to parallel advection and parallel spitzer - h " arm heat conduction are @xmath223 where @xmath224 is the mach number .
the sheath damping term entering the vorticity and ion pressure equations is given as @xmath225.\end{aligned}\ ] ]      the thin - layer approximation alters the conserved energy compared to conserved energy of the full drift - fluid model described in sec .
[ sec : energyconservation ] .
the derivation paths are very similar .
first , the pressure equations ( [ eq : electronpressure ] ) and ( [ eq : ionpressure_thinlayer ] ) are integrated neglecting surface terms .
next , the vorticity equation ( [ eq : w_norm ] ) is multiplied by `` @xmath226 '' ( [ eq : w_norm ] ) and integrated , again neglecting surface terms . adding the results we obtain the energy theorem @xmath227 .         = \int d\bm{x}\ , s_{\|}.\end{aligned}\ ] ] the energy sinks due to parallel losses are contained in @xmath228 as discussed in sec .
[ sec : energyconservation ] the electron and ion energy equilibration terms @xmath120 , the frictional heat exchange terms @xmath121 , and the viscous terms exactly cancel and therefore do not enter the energy theorem .
just as in the energy theorem for the full model given in eq .
( [ eq : cons_energy ] ) , collisions do not give rise to energy sinks and sources but merely give rise to conservative energy transfer terms in the model .
a thorough discussion of energy transfer mechanisms and the impact on mean flows can be found in refs . .
in order to test the simplified model for perpendicular collisional tranport applied in the hesel model , we have solved a subset of the governing equations ( [ eq : n_norm])-([eq : ionpressure_thinlayer ] ) numerically using the bout++ framework@xcite .
specifically , we assume that : 1 ) the system is homogeneous in the y - direction , 2 ) only terms describing collisional transport are retained , and 3 ) parallel sinks are not active i.e. @xmath229 . in this particular limit
, the hesel equations reduce to the 1d coupled set of non - linear equations @xmath230    in order to validate the simplified model , we have solved the full collisional model in the same limit .
the full model is described in details in section [ sec : moments ] . in the absence of turbulence
the model equations reduce to : @xmath231 \bigg )       + \bm{u}_r \cdot { \nabla_{\perp}}p_i&= -q_{\delta } , \\      \frac{3}{2}{\frac{\partial}{\partialt}}p_i       + \nabla \cdot \bm{q}_{i\perp }      + \frac{5}{2}\nabla \cdot ( p_i \bm{u}_r )          - \bm{u}_r \cdot { \nabla_{\perp}}p_i               & = q_{\delta }       \end{aligned }       \right\ } \quad \text{full model ( fm ) }       \label{eq : fm}\end{aligned}\ ] ] the resistive drift @xmath47 is given in eqs .
( [ eq : uru])-([eq : urt ] ) , the heat exchange term @xmath120 is given in eq .
( [ eq : energy_equilibration ] ) , and the perpendicular ion heat flux @xmath232 is given in eq .
( [ eq : ion_heatflux ] ) .
we stress that no terms are linearised and that collision frequencies are evaluated as dynamical functions of temperatures and particle density .    for comparison
we have also solved simple diffusion equations with constant diffusion coefficients @xmath233 and @xmath234 as in the hesel model @xmath235 the fields are initialized as @xmath236 and the constant reference values used in the hesel and simple diffusion model equations were @xmath237 where @xmath238 denotes the proton mass . in figure
[ fig : diff_solutions ]     solving the full non - linear model ( fm ) eqs .
( [ eq : fm ] ) ( green ) , the hesel model given in eqs .
( [ eq : hesel ] ) ( blue ) , and a simple diffusion model given in eqs .
( [ eq : sd ] ) ( purple ) .
solid black lines shows the initial condition @xmath239 given in eq .
( [ eq : initialcondition ] ) .
only one third of the simulation domain is shown.,title="fig:",scaledwidth=100.0% ] [ fig : n ]     solving the full non - linear model ( fm ) eqs .
( [ eq : fm ] ) ( green ) , the hesel model given in eqs .
( [ eq : hesel ] ) ( blue ) , and a simple diffusion model given in eqs .
( [ eq : sd ] ) ( purple ) .
solid black lines shows the initial condition @xmath239 given in eq .
( [ eq : initialcondition ] ) .
only one third of the simulation domain is shown.,title="fig:",scaledwidth=100.0% ] [ fig : pe ]     solving the full non - linear model ( fm ) eqs .
( [ eq : fm ] ) ( green ) , the hesel model given in eqs .
( [ eq : hesel ] ) ( blue ) , and a simple diffusion model given in eqs .
( [ eq : sd ] ) ( purple ) .
solid black lines shows the initial condition @xmath239 given in eq .
( [ eq : initialcondition ] ) .
only one third of the simulation domain is shown.,title="fig:",scaledwidth=100.0% ] [ fig : te ]     solving the full non - linear model ( fm ) eqs .
( [ eq : fm ] ) ( green ) , the hesel model given in eqs .
( [ eq : hesel ] ) ( blue ) , and a simple diffusion model given in eqs .
( [ eq : sd ] ) ( purple ) .
solid black lines shows the initial condition @xmath239 given in eq .
( [ eq : initialcondition ] ) .
only one third of the simulation domain is shown.,title="fig:",scaledwidth=100.0% ] [ fig : pi ]     solving the full non - linear model ( fm ) eqs .
( [ eq : fm ] ) ( green ) , the hesel model given in eqs .
( [ eq : hesel ] ) ( blue ) , and a simple diffusion model given in eqs .
( [ eq : sd ] ) ( purple ) .
solid black lines shows the initial condition @xmath239 given in eq .
( [ eq : initialcondition ] ) .
only one third of the simulation domain is shown.,title="fig:",scaledwidth=100.0% ] [ fig : ti ]    we show solutions of all three set of equations at @xmath240 for @xmath241 .
@xmath242 is the time at which the particle density amplitude is half the initial value according to the analytical solution of the 1d diffusion equation .
we have investigated various initial conditions but have chosen this specific case where the initial temperatures are not identical in order to demonstrate the energy exchange mechanisms within the model .
the size @xmath243 of the simulation domain is @xmath244 . in the simulations
the time evolution of the particle density @xmath62 does not vary significantly between the three models .
the sd and hesel models give almost identical results whereas the fm model shows a slightly flatter profile .
for the electron and ion temperatures and pressures the fm model and the hesel models agree quite well .
the sd model profiles , on the other hand , significantly differs from the profiles of the hesel and fm models .
the presence of energy exchange terms in the hesel and fm models is evident . here , the electron and ion temperatures are equilibrated whereas the absence of an energy coupling in the sd model is evident , particularly in the ion temperature . in conclusion
, the hesel model matches the full model quite well .
the simulations also demonstrate that a simple diffusion model should be used with care .
the implications of the choice of perpendicular collisional model in turbulence simulations is not investigated here but is left for future work .
lastly , in the 1d simulations we have observed that solving the full model takes significantly longer time than solving the hesel model . a comparison between the computational requirements of the full model fm and hesel is likewise left for future work .
in this paper we have demonstrated how perpendicular collisional effects can be incorporated in low - frequency drift - fluid turbulence models .
based on the electrostatic limit of the braginskii fluid equations@xcite , it was shown how resistivity and viscosity give rise to perpendicular fluid drifts which when inserted into the fluid moment equations provide an energy conserving closed model governing the time - evolution of the particle density , vorticity , and electron and ion pressure .
we demonstrated that resistivity gives rise to energy exchange between the electron and ion thermal energies and that viscous effects due to ion - ion collisions transfer energy between the modified kinetic and the ion thermal energies .
it was also shown that collisions do not imply energy sinks and sources in the global energy theorem , but merely provide conservative energy transfer terms . in conclusion
, it was shown that the inclusion of perpendicular collisional transport in drift - fluid models gives rise to perpendicular diffusion in all moment equations , which obviates the need for explicit artificial diffusion .
based on the general results , we also derived a simplified 2d model aiming at describing interchange driven turbulent transport in the vicinity of the outboard midplane .
the model is named hesel and , as the predecessor esel , it is aimed at being computational efficient allowing the model equations to be time integrated for multiple turbulence de - correlation times .
the validity of the simplified collisional dynamics in hesel was investigated by means of numerical simulations without turbulence .
the simulations showed good agreement between hesel and the full model .
this work has been carried out within the framework of the eurofusion consortium and has received funding from the euratom research and training programme 2014 - 2018 under grant agreement no .
the views and opinions expressed herein do not necessarily reflect those of the european commission .
gunn , c.  boucher , m.  dionne , i.  uran , v.  fuchs , t.  loarer , i.  nanobashvili , r.  pnek , j .- y .
pascal , f.  saint - laurent , j.  stckel , t.  van rompuy , r.  zagrski , j.  admek , j.  bucalossi , r.  dejarnac , p.  devynck , p.  hertout , m.  hron , g.  lebrun , p.  moreau , f.  rimini , a.  sarkissian , and g.  van oost , , * 363 * , 484 ( 2007 ) .