relative dispersion of passive particles in turbulent flows is one of the fundamental problems in turbulence research .
it characterizes the transport and mixing properties of turbulence and is important from both theoretical and practical points of view . reflecting universal behavior of turbulent fluctuations
, relative dispersion also has some universal properties because of its locality - in - scale nature @xcite .
in particular , the dispersion process exhibits anomalous dispersion in the inertial range .
this is first observed by richardson ( 1926 ) @xcite , and since then , a number of theoretical , experimental , and numerical investigations have been devoted to understand and model relative dispersion process @xcite .
however , comprehensive understanding has not been obtained yet .
recently , a few works focusing on underlying mechanism of the anomalous dispersion were reported . in the inertial range ,
the mean free path , @xmath2 , the mean length for persistent expansion of relative separation without changing its moving direction , is an order of relative separation itself , @xmath0 : @xmath3 @xcite . here
, @xmath4 is a non - dimensional constant called the persistent parameter . in two - dimensional inverse cascade turbulence ( 2d - ic )
, @xmath4 is estimated as @xmath5 @xcite .
this means separating motions are not purely diffusive but composed of an appreciable amount of persistent motions .
in addition , it was also reported that there is a relation between stagnation - point structures and richardson s law @xcite , and that dispersion process is described by persistent streamline topology @xcite . from these results
, it is expected that coherence of turbulent field , which must share its origin with fine coherent vortical structures such as worms in three - dimensional navier - stokes ( 3d - ns ) turbulence , has a significant role in turbulent relative dispersion .
the correlations in turbulence are characterized in scale - space due to their self - similarity , and are not made disappear by coarse graining in real space . reflecting these correlations ,
relative separation moves persistently to some extent , so that , unlike the brownian motion , relative separation process should not be described only by random collision motions even as an approximation@xcite . in other words ,
the characteristic length can not be defined because the mean free path , @xmath2 , varies depending on the spatial scale . whereas there are several experiments and numerical simulations of which results are rather close to the prediction of richardson s diffusion equation @xcite that closely relates to random collision motions .
thus , these results raise a question ; how are the effects of persistent motions wiped out ? in the previous paper @xcite , we have introduced a self - similar telegraph model of turbulent relative dispersion , and showed that the separation pdf can be close to the prediction of richardson s diffusion equation for slowly separating particle pairs even in the presence of persistent motions .    in the present paper
, we check the consistency of the physical picture of the self - similar telegraph model by carrying out direct numerical simulations ( dns ) of 2d free convection ( 2d - fc ) turbulence instead of 3d - ns turbulence .
this is because ( i ) dns of 3d - ns turbulence requires extremely large computer resources , so that it is difficult to track particles for a long time , which is necessary to investigate dynamical properties of relative dispersion process , and ( ii ) 2d - fc turbulence has both statistical and dynamical characteristics similar to those of 3d - ns turbulence ( see appendix a ) @xcite . among them , the existence of coherent structures , which are approximated by the burgers t - vortex layer , is notable .
this is a crucial difference from coherent structures in 2d - ic turbulence that are nested cat s eye vortices @xcite .
therefore , comparing the results of the 2d - fc case with those of the 2d - ic case , we can investigate the effects of coherent structure on turbulent relative dispersion .
this comparison gives a physical meaning of the drift term of the self - similar telegraph model .
the inertial range achieved by our dns is not so wide that the relative motions of particle - pairs in the dissipation , the inertial , and the energy containing scales are not sufficiently resolved by usual fixed time statistics . in order to investigate the scaling natures of relative dispersion in such a narrow and limited inertial range , we utilize exit - time statistics introduced into research of turbulent relative dispersion by boffetta and sokolov ( see appendix b ) @xcite . by detailed investigation of the pdf of exit - time ,
we show that the pdf is divided into two region , the region - i and -ii , corresponding to persistent expansion and random transition between expansion and compression of relative separation , respectively .
this result agrees with the picture of the self - similar telegraph model .
in addition , we provide a method for estimation of the parameters of the self - similar telegraph model by using exit - time pdf .    in the following sections ,
first we provide a brief review of the self - similar telegraph model in  2 .
section 3 presents a summary of the numerical scheme and parameters of our dns of 2d - fc turbulence .
some of the results by the fixed - time and the exit - time statistics are provided and discussed in  4 .
concluding remarks are made in  5 . in addition , some properties of 2d - fc turbulence and exit - time statistics are presented in appendix a and b , respectively .
in the previous paper @xcite , we have introduced a self - similar telegraph model for turbulent relative dispersion , which is a model describing the evolution of the pdf of relative separation of particle pairs in the inertial range .
the relative separation of two particles , @xmath6 is defined as follows : @xmath7 where @xmath8 and @xmath9 are the lagrangian positions of the particles . the model is based on sokolov s model @xcite and consists of persistent expansion and compression of relative separation , @xmath0 , according to the relative velocity , @xmath10 , where @xmath11 and @xmath12 are a dimensional constant and a scaling exponent , respectively : @xmath13 for kolmogorov scaling and @xmath14 for bolgiano - obukhov scaling . the transition rate from expansion to compression and that from compression to expansion are given by @xmath15 and @xmath16 , respectively .
here @xmath17 is a characteristic time scale at a spatial scale @xmath0 , @xmath18 , and @xmath19 are the inverses of persistent parameters , @xmath20 for expansion and @xmath21 for compression . introducing @xmath22 and @xmath23 , the evolution equation of the separation pdf , @xmath24 ,
is derived as follows : @xmath25   \\
+ \sigma\frac{{\partial}}{{\partial}r }   \left[v(r)p\right ]   ,   \label{eq : t - model }   \end{aligned}\ ] ] where @xmath26 is richardson s diffusion coefficient , @xmath27 , and @xmath28 .
the parameters of the model are @xmath29 and @xmath30 : @xmath29 represents the strength of persistency of moving directions , and @xmath30 does the difference in persistency between expansion and compression .
therefore , @xmath29 and @xmath30 reflect the strength and structure of coherence of the flow , respectively .    for slowly - separating particle pairs , i.e. , in the case of @xmath31 , the first term of the l.h.s .
( [ eq : t - model ] ) can be neglected and the approximated equation is given by @xmath32   + \sigma\frac{{\partial}}{{\partial}r }   \left[v(r)p\right ] .
\label{eq : palm - eq}\ ] ] this form of the equation was first derived by palm @xcite and is also the same as the diffusion equation of goto - vassilicos model @xcite .
we call eq .
( [ eq : palm - eq ] ) palm s equation .
the similarity solution of eq .
( [ eq : palm - eq ] ) is @xmath33 where @xmath34 is the normalization factor .
because the tail of the exit - time pdf , @xmath35 , consists of slowly - separating particle pairs , it is calculated from palm s equation ( [ eq : palm - eq ] ) with the method used by boffetta and sokolov @xcite .
the asymptotic form of @xmath35 is given by @xmath36 where @xmath37 is the @xmath38-th zero of the @xmath39-th order bessel function and @xmath40 is the mean exit - time from @xmath0 to @xmath41 : @xmath42 note that the mean exit - time calculated from the self - similar telegraph model , eq .
( [ eq : t - model ] ) , is the same as eq .
( [ eq : avetime - palm ] ) .
this is because the mean exit - time is calculated from a steady solution of the equation @xcite , and the solution of eq .
( [ eq : t - model ] ) is the same as that of eq .
( [ eq : palm - eq ] ) . by comparing the tail of the exit - time pdf obtained by dns and eq .
( [ eq : palm - etime ] ) , we can estimate the value of @xmath30 .
the last term of the r.h.s .  of eq .
( [ eq : t - model ] ) is a drift term consistent with the scaling law , and the direction of the drift is determined by @xmath43 .
the parameter @xmath43 consists of two parts , the `` scaling - determined '' one , @xmath44 , and the `` dynamics - determined '' one , @xmath30 . in order for eq .
( [ eq : palm - eq ] ) to recover richardson s equation , the drift term has to disappear , which means @xmath45 .
we call this case the richardson case or the zero - drift case , where the parameters of the model reduce to one , @xmath46 ; @xmath47 is determined by the relation @xmath48 .
rescaled energy and entropy spectrum @xmath49 and @xmath50 obtained by dns at resolution @xmath51 ( nv10 ) and @xmath52 ( nv11 ) .
the straight lines refer to the bolgiano - obukhov scaling : @xmath53 , @xmath54 .
, width=332 ]    .[table : dnsparms ] parameters used in the present dns . [ cols="^,^,^,^,^,^",options="header " , ]     in this section , we explain the method of dns used in the present work and show basic properties of turbulent fields produced by our simulation .
we generate turbulent field by dns of the vorticity equation with a large - scale friction term @xmath55 and the temperature equation with a large - scale forcing term @xmath56 : @xmath57 where @xmath58 , @xmath59 , @xmath60 , @xmath39 , @xmath61 , @xmath62 , and @xmath63 represent the vorticity field , the temperature field , the velocity field , the kinematic viscosity , the thermal diffusivity , the thermal expansion coefficient , and the gravitational acceleration , respectively .
the large - scale forcing term used here is @xmath64 where @xmath65 is a constant .
the large - scale friction term is written in the fourier space as @xmath66 where @xmath67 , @xmath68 , @xmath69 , and @xmath70 are a constant , the wave number vector , the fourier mode of the friction term , and the fourier mode of the vorticity field , respectively .
our dns is performed on a @xmath71 domain with the doubly periodic boundary conditions at resolutions @xmath72 : @xmath51 ( nv10 ) and @xmath73 ( nv11 ) .
we employ a pseudo - spectral method for accurately calculating convolutions and spatial derivatives , and 4-th order runge - kutta method for time integration .
aliasing error is removed by adopting the phase - shift method ( nv10 ) and the @xmath74 method ( nv11 ) .
all results presented in this paper are obtained for statistically stationary and locally isotropic turbulence .
we summarize the parameters used in our simulation in table [ table : dnsparms ] and characteristic quantities of generated turbulence in table [ table : char_quantities ] .
figure [ fig : spectrum ] shows the entropy and energy spectra obtained by our dns .
the spectra are rescaled with entropy dissipation scales and the entropy dissipation rate , and there is a region consistent with the bolgiano - obukhov scaling , eqs .
( [ eq : bo - ek ] ) and ( [ eq : bo - sk ] ) , ( see appendix a ) .
we call this region the inertial range .
most investigations in the present paper are carried out in this range .    in the velocity field generated by dns , we track a number of particle pairs ( @xmath75 pairs ) according to the advection equation : @xmath76 where @xmath77 is the position of the @xmath78-th particle at time @xmath79 .
we employ the 4-th order runge - kutta method for numerical integration , and the linear interpolation to obtain the velocity of each particles .
particle pairs are distributed homogeneously with relative separation of the grid scale @xmath80 at the initial time .
temporal evolution of mean relative separation @xmath81 obtained in the case of nv11 .
solid , dotted and dashed line refer to @xmath82 and @xmath83 , respectively .
straight lines indicate richardson s law .
the inset shows the local slope of @xmath81 .
, width=332 ]         first , we briefly discuss the results obtained by standard fixed - time statistics , which concerns the distribution of relative separation at a certain ( fixed ) time .    in fig .
[ fig : ftime ] we plot temporal evolutions of the mean relative separation of particle pairs , @xmath84 , for different powers , @xmath85 , @xmath86 , and @xmath83 . in the cases of @xmath85 and @xmath86 , the generalized richardson s relations , @xmath87 , are realized , though the times at which they start and the time intervals in which they hold differ .
these differences result from the limited width of the inertial range in dns .
for example , in the case of nv11 , the inertial range is roughly estimated as @xmath88 . from fig .
[ fig : fpdf](a ) , it is clear that the separation pdf broadens rapidly , so that @xmath84 is contributed to by a quite broad range of relative separation including the dissipation , the inertial , and the energy containing ranges . as a result ,
the weight of each range varies with the power of the moment of the relative separation and thus , the range of time in which the contribution of the inertial range dominates differs accordingly .
the slope of @xmath89 is less steep than richardson s law .
this is because , in the case of @xmath90 , the contribution of particle pairs in the energy containing range to @xmath84 is much larger than that in the cases of @xmath91 and @xmath92 . in the energy containing range ,
relative separation process is described by the brownian motion , @xmath93 , which is less steep than richardson s law .
hence , the larger the contribution of the range is , the less steep the slope of @xmath84 is .    even though the richardson s law is observed for @xmath85 and @xmath86 , the fact does nt necessarily support the validity of the richardson s law because it is reported that the temporal evolution of @xmath84 strongly depends on the initial separation of particle pairs @xcite .
hence we have to check the scaling law by adopting a different method that is independent of the initial separation .
in order to check self - similarity of the pdf of relative separation , we plot the pdf rescaled with @xmath94 in fig .
[ fig : fpdf](b ) .
the values of @xmath94 at @xmath95 , @xmath96 , @xmath97 , and @xmath98 are @xmath99 , @xmath100 , @xmath101 , and @xmath102 , respectively .
all pdfs collapse well for @xmath103 and they are in good agreement with the similarity solution of palm s equation with @xmath104 .
however , we can not conclude that the temporal evolution of the separation pdf is self - similar and governed by palm s equation because fig .
[ fig : fpdf](b ) is obtained not from separations only in the inertial range but in the much wider range , @xmath105 . although the collapse in fig .
[ fig : fpdf](b ) implies existence of a self - similar stage governed by palm s equation , we have not had any reasonable explanation of the collapse .
residual ratio of particle pairs .
dotted lines with open and closed circles refer to data obtained by dns at resolution @xmath51 ( nv10 ) and @xmath52 ( nv11 ) respectively .
the inset is the log - log plot of the ratio .
, width=332 ]     scale dependence of the mean exit - time @xmath106 for @xmath107 rescaled with the dissipation time scale , @xmath108 .
open and closed circles refer to the results obtained by the present dns at resolution @xmath51 ( nv10 ) and @xmath52 ( nv11 ) , respectively .
the inset is the compensated plot with @xmath109 .
the dashed line represents estimated values of the coefficients @xmath110 .
, width=332 ]                  the scale dependence of the mean exit - time obtained by our dns is shown in fig .
[ fig : etime ] .
it is observed that , although the width is narrow , there is a region consistent with the bolgiano - obukhov scaling , @xmath111 .
the exit - time statistics is independent of initial separation of particle pairs if spatial scale @xmath0 is large enough for them to forget information of their initial conditions @xcite . in our simulation ,
the initial separation of particle pairs is @xmath80 that is much smaller than the scales in which the scaling law holds .
therefore , this result indicates that richardson s law is valid in the 2d - fc turbulence .    according to the scaling law
, the form of the mean exit - time is given by @xmath112 where @xmath113 is considered to be a universal constant .
the inset of fig .
[ fig : etime ] shows a compensated plot of the mean exit - time .
although the weak @xmath114 dependence of @xmath113 seen in fig .
[ fig : etime ] requires higher resolution of dns to more accurate estimation , @xmath113 is estimated as @xmath115 in the present dns .    in order to obtain statistically reliable data of exit - time
, it is important that the residual ratio of particle pairs must be small enough to take in slowly separating particle pairs . the residual ratio at a scale @xmath0 and a time @xmath79 is defined as the ratio of particle pairs of which first passage times at @xmath0 are less than @xmath79 .
if the residual ratio is not small enough at the time when statistics of exit - time are calculated , the statistics can not reflect slowly separating particle pairs .
figure [ fig : erate ] shows the residual ratio of particle pairs used in the present work at the termination time of dns , @xmath116 .
it is obvious that the ratio is almost zero in the inertial range , so that our exit - time data is reliable in the sense mentioned above .
figure [ fig : marginal - etime](b ) shows the scale dependence of the mean exit - time obtained from insufficient data and illustrates the importance of the residual ratio .
although the results at @xmath117 seem to have the wider inertial range than others , it is a fake . therefore , to obtain statistically reliable data , we have to track particle pairs until the residual ratio becomes at least less than @xmath118 .
in contrast to the inertial range , the mean exit - time is almost constant in the dissipation range ( @xmath119 in fig .
[ fig : etime ] ) except in the initial transient region , because the relative velocity is proportional to relative separation : @xmath120 . the exit - time , then , is given by @xmath121 note that the mean exit - times rescaled by @xmath108 collapse each other in fig .
[ fig : etime ] .
this is because there is only one characteristic time scale , @xmath108 , in the dissipation range @xcite .
we plot exit - time pdf , @xmath35 , in the inertial range in fig .
[ fig : pdf - et](a ) , and that rescaled by the mean exit - time , @xmath122 , in fig .
[ fig : pdf - et](b ) .
obviously the pdfs of different scales collapse onto one curve in fig  [ fig : pdf - et](b ) .
this means exit - time pdf is self - similar in the inertial range and indicates that relative dispersion process is self - similar .    as is clear from fig .
[ fig : pdf - et](b ) , the exit - time pdf consists of two regions : the sharp peak and the long exponential tail .
we call these two regions the region - i and the region - ii , respectively .
the qualitative difference in form between these indicates that the pdf reflects two different types of motions . in order to clarify the difference ,
snapshots of typical distribution of particle pairs in the region - i and -ii are shown in figs .
[ fig : pp - dist - i ] and [ fig : pp - dist - ii ] , respectively , with a snapshot of the magnitude of t - vorticity field @xmath123 .
temporal evolution of relative separation for several particles in the case of nv11 .
different lines represent different particle pairs . ,
width=332 ]    in fig .
[ fig : pp - dist - i](a ) , on the snapshot of @xmath123 at @xmath124 ( @xmath125 ) , we plot line segments representing particle pairs of which exit - time @xmath126 for @xmath127 and @xmath128 satisfies @xmath129 , that is , @xmath126 is in the shaded region in fig .
[ fig : pp - dist - i](b ) . in order to select such particle pairs ,
we first extract pairs satisfying the condition @xmath130 , and then pick out ones of which first passage times at @xmath131 are smaller than @xmath132 and those at @xmath133 are larger than @xmath132 .
figure [ fig : pp - dist - ii](a ) is drawn by the same procedure . to draw figs .
[ fig : pp - dist - i ] and [ fig : pp - dist - ii ] , we used the data of nv10 .    in 2d - fc turbulence ,
fine coherent structures are well approximated by the burgers t - vortex layer , so that shear layers are formed around the structures @xcite .
in addition , such structures are persistent in time to some extent .
hence , particles around the coherent structures are advected along them and relative separations of the particles expand or compress persistently .
we call this type of motion the persistent separation . as is shown in fig .
[ fig : pp - dist - i](a ) , particle pairs in the region - i appear to be along the fine coherent structures .
moreover relative separations of particle pairs in fig .
[ fig : pp - dist - i](a ) expand rapidly as can be known by their short exit - time . therefore , fig .
[ fig : pp - dist - i ] supports the picture of persistent separations .
in contrast to the region - i , the distribution of particle pairs in the region - ii is in disorder ( see fig .  [
fig : pp - dist - ii](a ) ) . because the particle pairs contained in the region - ii have long exit - time , the longitudinal relative velocity of them remains small positive during the passage from @xmath131 to @xmath133 or moving direction of relative separation of the pairs changes from expansion to compression at least once before they reach @xmath133 .
however it is quite unlikely that relative velocity remains small positive for a long time , so that the latter is probably the main mechanism of the formation of the region - ii .
hence , it is expected that both of expanding and compressing pairs are contained in fig .
[ fig : pp - dist - ii](a ) .
in fact , there are several pairs placed along coherent structures , which probably are expanding .
note that their relative separations are not confined to be larger than @xmath131 .
in addition , there are also pairs across the structures , which are presumably compressing because the structures are approximated by the burgers t - vortex layer
. moreover , some pairs are located at positions where the fine coherent structures twist .
the twisted regions are generated when well - stretched structures lose their activities and are folded , or when the structures are generated by plumes .
therefore the particle pairs in such regions are changing their moving directions from expansion to compression or the opposite .
figure [ fig : pp ] shows five typical evolutions of particle - pair separation .
sample 3 expands persistently without any strong compression but sample 4 experiences both persistent expansion and compression several times .
the former corresponds to motions along structures and belongs to the region - i ; the latter does across structures or twisted region and belongs to the region - ii .
note that even a single evolution contains motions that are categorized into the region - i or into the region - ii .
these facts are consistent with the assumptions of the self - similar telegraph model : relative separation process consists of persistent expansion and compression with some random transition mechanisms .       estimated values of @xmath30 in the case of nv11 .
dotted lines with open box , closed box , open circle , closed circle , open triangle , closed triangle , and reverse open triangle denote @xmath134 , @xmath135 , @xmath136 , @xmath137 , @xmath138 , @xmath139 , and @xmath140 , respectively . solid and dashed lines indicate the mean value , @xmath141 , and the value of the richardson s case , @xmath142 , respectively .
the shaded region represents the standard deviation .
, width=332 ]    it is expected that expansion and compression of relative separation differ in persistency ; the former is experienced mainly by particle pairs along coherent structures and the latter does by those across the structures . since the auto - correlations of strain and temperature along the structures have longer characteristic lengths than those across them in 2d - fc turbulence @xcite , expanding motions must be more persistent than compressing ones .
so as to confirm this consideration , we investigate the slope of the tail of the exit - time pdf .
then , with eq .
( [ eq : palm - etime ] ) , we estimate @xmath30 that describes the difference in persistency between expansion and compression of relative separation .    note that using the asymptotic form of the exit - time pdf calculated from palm s equation , eq .
( [ eq : palm - etime ] ) , may be justified as follows . introducing a time scale @xmath143 , the l.h.s .
of the self - similar telegraph model , eq .
( [ eq : t - model ] ) , is rewritten as follows : @xmath144 where @xmath145 . in the case that @xmath146 at a certain spatial scale @xmath0 , the first term of the l.h.s . can be neglected , and thus , the model is reduced to palm s equation .
hence , the tail of the exit - time pdf by the telegraph model , which consists of slowly - separating particle pairs , must coincide with that by palm s equation .
it is also easy to show that the mean exit - time is the same between eqs .
( [ eq : t - model ] ) and ( [ eq : palm - eq ] ) . we , therefore , adopt eq .
( [ eq : palm - etime ] ) for the estimation of @xmath30 .
figure [ fig : delta ] shows the estimated values of @xmath30 in the inertial range .
the mean value of @xmath30 is @xmath147 and the standard deviation is @xmath148 .
this negative value denotes that expanding motions are more persistent than compressing ones , which supports the above expectation .
in addition , @xmath30 is larger than that of the richardson case ( the zero - drift ) , @xmath149 . because @xmath23 , @xmath150 means @xmath151 , that is , the compressing motion of relative separation is more persistent than that of the richardson case .
this fact results in the negative drift in the self - similar telegraph model , eq .
( [ eq : t - model ] ) .
the negativity of the drift can be accepted considering that coherent structures in 2d - fc turbulence are string - like and scale - transversal ones @xcite .
that is , because the coherent structures in 2d - fc turbulence are not nested in scales , there are few obstacles to compressing motions . on the other hand ,
the typical structures in 2d - ic turbulence are nested vortices called `` cat s eye in a cat s eye '' . a separation process of particle pairs in 2d - ic turbulence ,
thus , is a step - by - step one each step of which consists of a trapping by one of nested vortices and sudden separation into a next larger vortex @xcite ; compressing motions are probably blocked by the nested vortices .
in fact , the estimated values of @xmath30 from dns of 2d - ic turbulence by boffetta & sokolov @xcite and goto & vassilicos @xcite are smaller than that of the richardson case @xcite ; the compressing motion of relative separation is less persistent than the richardson case .
this results in the positive drift in eq .
( [ eq : t - model ] ) .
we , therefore , conclude that the drift term in the self - similar telegraph model , eq .
( [ eq : t - model ] ) , reflects the characteristics of coherent structures of the flow .
rescaled pdf of exit - time ( @xmath152 ) in the case of nv11 and the prediction of palm s equation ( @xmath141 ) for several values of @xmath153 .
solid and dashed lines refer to the results of dns and that obtained from palm s equation through the relation eq .
( [ eq : palm - epdf ] ) .
, width=332 ]    in fig .
[ fig : epdf - palm ] , we compare the pdf of exit - time obtained by dns with that calculated from palm s equation by eq .
( [ eq : palm - epdf ] ) . in the region - i
, the form of the pdf by dns is totally different from that by palm s equation . in the region - ii , however , the two pdfs collapse onto a single curve even when @xmath153 is very small , and thus , the separation process of particle pairs in the region - ii can be described by palm s equation , eq .
( [ eq : palm - eq ] ) .
this indicates that the motion of particle pair is diffusive in the region - ii ; the separating process from @xmath154 to @xmath0 does nt affect that from @xmath0 to @xmath155 .
it is also observed that the form of the exit - time pdf varies depending on the value of @xmath153 in fig .
[ fig : epdf - palm ] .
the larger the value of @xmath153 is , the larger the proportion of the region - ii becomes .
if @xmath153 is sufficiently large , the whole pdf seems to be occupied by the region - ii and the relative separation process is substantially described by palm s equation , eq .
( [ eq : palm - eq ] ) .
as mentioned in the previous subsection , the tail of the exit - time pdf calculated from the self - similar telegraph model , @xmath156 , must agree with that calculated from palm s equation , @xmath157 .
in addition , as the value of @xmath153 gets larger , the region of @xmath156 overlapping with @xmath157 expands .
this is because , for large values of @xmath153 , most of exit - times of particle pairs are longer than the characteristic time @xmath158 .
hence , the self - similar telegraph model probably has the same properties as the results shown in fig .
[ fig : epdf - palm ] .
rescaled pdf of exit - time ( @xmath152 ) in the case of nv11 .
the peak value is normalized to be unity .
different lines represent different values of @xmath153 : @xmath159 , @xmath160 , @xmath161 , @xmath162 , @xmath137 , @xmath163 , @xmath164 , and @xmath165 from the curve with the lowest cut - off rescaled exit - time to that with the highest .
the inset is the same plot in linear scale .
, width=332 ]    figure [ fig : epdf - eqpeak ] shows the exit - time pdf for nv11 rescaled with @xmath40 and normalized by their peak values .
it is clear that the shape of the region - i is independent of @xmath153 .
that is , there exists a @xmath153-independent function representing the shape of the exit - time pdf in the region - i , @xmath166 . here
, @xmath167 denotes rescaling by the mean exit - time , @xmath40 .
we assume that @xmath168 is normalized as @xmath169 .
then , in the region - i , the exit - time pdf is written as @xmath170 where @xmath171 is a normalization factor depending on @xmath153 .
the exit - time pdf , @xmath172 , takes the peak value at @xmath173 .
it is estimated that @xmath174 in fig .
[ fig : epdf - eqpeak ] , which can be regarded as a characteristic time scale of the region - i . as shown in the inset of fig .
[ fig : epdf - eqpeak ] , @xmath166 distributes sharply around @xmath173 .
if the region - i is formed by particle pairs of which relative separations expand persistently according to @xmath10 , the exit - time in the region - i , @xmath175 , is given by the pass through time from @xmath0 to @xmath41 : @xmath176 then , using the general form of the mean exit - time , eq .
( [ eq : etime - general ] ) , the exit - time rescaled by @xmath40 in the region - i , @xmath177 , is @xmath178 which is independent of @xmath153 .
thus , @xmath166 can be connected to the pdf of @xmath11 , @xmath179 , which is the pdf of lagrangian relative velocity , @xmath180 , rescaled with @xmath181 : @xmath182 note that @xmath179 is expected to be independent both of @xmath0 and @xmath153 . since ,
in the self - similar telegraph model , the distribution of the coefficient , @xmath11 , of the relative velocity is not considered , @xmath179 is a @xmath30-function : @xmath183 .
accordingly , the region - i of the exit - time pdf is approximated by a @xmath30-function at @xmath184 in the model .    combining eqs .
( [ eq : avetime - palm ] ) and ( [ eq : t - model_region - i ] ) , we can estimate @xmath29 from @xmath184 . from eq .
( [ eq : avetime - palm ] ) , @xmath185 in the case of the self - similar telegraph model is given by @xmath186 then , @xmath29 is calculated as follows : @xmath187 if we assume that @xmath188 , then , @xmath29 is estimated as @xmath189 .
however , the value of @xmath184 is not necessarily well defined . roughly speaking
, it may take a value satisfying @xmath190 , if assuming @xmath184 is in the half - value width of the region - i around @xmath173 .
then , the estimated value of @xmath29 is in the range , @xmath191 .
we have investigated relative dispersion in 2d free convection turbulence by direct numerical simulation . in the inertial range ,
where the entropy cascade dominates , we have confirmed with exit - time statistics that relative dispersion satisfies the bolgiano - obukhov scaling and , therefore , is self - similar .
it was also shown that the exit - time pdf , @xmath192 , is divided into two parts , the region - i and -ii , and that both of them satisfy the scaling - law and the self - similarity .
@xmath192 is written as the following form : @xmath193 where @xmath194 is a division time - scale ( exit - time ) between the region - i and -ii , and @xmath195 is a smoothed step function such that @xmath196 if @xmath197 and @xmath198 if @xmath199 . @xmath200 and @xmath201 correspond to the exit - time pdf of the region - i and -ii , respectively .
the investigation of the distribution of particle pairs in the real space indicates that the region - i and -ii are formed , respectively , by particle pairs expanding along coherent structures and by those experiencing turns between expansion and compression ( figs .
[ fig : pp - dist - i ] and [ fig : pp - dist - ii ] ) .    figures [ fig : epdf - eqpeak ] and [ fig : epdf - palm ] show the following characteristics of @xmath200 and @xmath201 .
the form of @xmath200 is independent of @xmath153 if it is rescaled with the mean exit - time .
moreover , if we assume that the region - i is constituted by persistently - separating particle pairs , the form of @xmath200 is related to the pdf of lagrangian relative velocity , @xmath10 .
that is , @xmath202 where @xmath171 is a @xmath153-dependent normalization factor , and @xmath179 is the pdf of @xmath180 rescaled with @xmath181 .
on the other hand , @xmath201 agrees with the exit - time pdf calculated from palm s equation , @xmath157 , the form of which varies depending on @xmath153 even when it is rescaled with @xmath40 .
hence @xmath172 is written as @xmath203     similarity solution ( fixed - time pdf of relative separation ) of the self - similar telegraph model in the case of 2d - fc turbulence .
the parameters used in this figure are estimated from dns data for nv11 by using the exit - time pdf .
different lines refer to different values of @xmath29 .
@xmath204 corresponds to the similarity solution of palm s equation .
the inset is a blowup of the head region in the linear scale .
, width=332 ]    these results support the self - similar telegraph model . in the model
, @xmath205 is approximated by a @xmath30-function .
on the other hand , for slowly - separating particle pairs , the model is approximated by palm s equation .
that is , the head region of the fixed - time pdf as well as the tail of the exit - time pdf is approximately described by palm s equation . by taking advantage of these characteristics of the exit - time pdf , we can estimate the parameters of the model , @xmath29 and @xmath43 ( @xmath30 ) : @xmath29 is estimated in the region - i and @xmath30 is in the region - ii . by our dns of 2d - fc turbulence ,
these parameters are estimated as @xmath206 and @xmath207 ( @xmath208 ) .
the estimated value of @xmath43 corresponds to the negative drift in the model .
this is the opposite to the 2d - ic turbulence case @xcite .
we speculate that the difference in the sign of @xmath209 is caused by the difference in coherent structures . with these parameters
, we can numerically calculate the similarity solution of the self - similar telegraph model , which is shown in fig .
[ fig : tm - fpdf ] @xcite
. however we can not compare the solution directly with the results of dns because the inertial range achieved by the dns is too narrow .
we need dns with higher resolutions to the comparison .
although the self - similar telegraph model can capture the essential characteristics of the relative separation process , the distribution of the relative velocity is crucial to understand further details of separation processes and their relation to the coherent structures , as suggested by the exit - time pdf obtained by dns .
in fact , we have obtained preliminary results showing that the distribution of the relative velocity in the inertial range is tightly connected to that in the dissipation range @xcite .
this strongly indicates that the distribution directly relates to the coherent structures because their length are of the order of the inertial range and their width of the order of the dissipation range @xcite .
an extension of the model to include the distribution of relative velocity , as well as an extension to the dissipation range , will be done in the future work .
this work was supported by the grant - in - aid for the 21st century coe `` center for diversity and universality in physics '' from the ministry of education , culture , sports , science and technology ( mext ) of japan . numerical computation in this work
was carried out on a nec sx-5 at the yukawa institute computer facility .
the governing equations of the 2d - fc turbulence are @xmath210 where @xmath60 , @xmath59 , and @xmath211 represent velocity , temperature , and pressure field , respectively .
@xmath212 is the unit vector in the direction of the gravity .
@xmath39 , @xmath61 , @xmath213 , @xmath62 , and @xmath63 are the kinematic viscosity , the molecular diffusivity , the mean density of the fluid , the thermal expansion coefficient , and the gravitational acceleration , respectively .
the 2d - fc turbulence has two important properties : the bolgiano - obukhov scaling and the fine coherent structures .      in this system , the integral of the squared temperature , @xmath214 is a conserved quantity in the ideal case , where @xmath215 is the volume of the system .
we call this quantity entropy for convenience .
the entropy cascades from large to small scales similar to the energy cascade in the 3d - ns turbulence , which leads scaling laws of the energy and entropy spectra , @xmath216 and @xmath217 : @xmath218 where @xmath219 and @xmath220 are considered to be universal constants , and @xmath221 is the entropy dissipation rate defined as follows : @xmath222 in this paper , @xmath223 denotes an ensemble average . these scaling laws , eqs .
( [ eq : bo - ek ] ) and ( [ eq : bo - sk ] ) , are called the bolgiano - obukhov scaling @xcite .
entropy dissipation length scale @xmath224 and time scale @xmath108 is estimated by dimensional analysis with @xmath221 , @xmath225 , and @xmath61 : @xmath226 these scales correspond to the kolmogorov scales in the 3d - ns turbulence .
there are several other quantities characterizing 2d - fc turbulence : the thermal taylor microscale @xmath29 , and the rayleigh number at @xmath29 , @xmath227 .
these are defined as @xmath228 and @xmath229 , where @xmath230    according to the bolgiano - obukhov scaling , richardson s law in the 3d - ns turbulence is modified as follows : @xmath231 where @xmath6 is relative separation of passive particles at time @xmath79 and @xmath232 is considered to be a universal constant .
in addition , richardson s diffusion equation of the separation pdf @xmath24 is modified as follows : @xmath233 ,     \label{eq:2dfc - richardson }   \end{aligned}\ ] ] where @xmath24 is the probability density of relative separations at @xmath0 and @xmath79 , and @xmath234 is considered to be a universal constant .
this equation has a self similar solution if the initial condition is the delta function . for @xmath235 , @xmath236 ,     \label{eq:2dfc - richardson - solution }   \end{aligned}\ ] ] where @xmath237 is the normalization constant .
we introduce a vector quantity called t - vorticity , @xmath238 , as follows : @xmath239    figure [ fig : snap_tv_a ] shows a snapshot of the magnitude of the t - vorticity and the strain field .
it is observed that there are linearly concentrated area of t - vorticity .
its width and length are an order of the entropy dissipation scale and the integral scale , respectively @xcite .
we call these fine coherent structures . at the same location of the structures
, we can observe linearly straining regions in fig .
[ fig : snap_tv_a](b ) .
t - vorticity has similar properties to vorticity in the 3d - ns turbulence .
the evolution equation of the t - vorticity has the stretching term and is the same as that of vorticity in 3d - ns system : @xmath240 accordingly , there is a solution corresponding to the burgers vortex in the 3d - ns system , which is called the burgers t - vortex layer .
it is numerically confirmed that fine coherent structures in the 2d - fc turbulence are well approximated by the burgers t - vortex layer @xcite .
this is similar to fine coherent structures such as worms in the 3d - ns turbulence which are well approximated by the burgers vortex .
the exit - time statistics is one of the scale - fixed statistics . the exit - time , @xmath241 ,
is defined as @xmath242 where @xmath243 is the time when a relative separation @xmath6 reach a threshold @xmath244 for the first time ( first - passage time ) @xcite . according to the scaling law of the characteristic time , @xmath245 , and the additivity of the mean
, the form of the mean exit - time is expected as follows : @xmath246 where @xmath185 is a constant depending on the system . in the case of the bolgiano - obukhov scaling , @xmath247 , where @xmath113 is considered to be a non - dimensional universal constant .      1 .
the exit - time statistics can specify a spatial scale by choosing a threshold @xmath244 .
therefore we can extract information of the inertial range if both @xmath244 and @xmath248 are in the inertial range .
the interval between thresholds ( the width for averaging ) can be controlled by @xmath153 .
this means that we can control the degree of coarse graining of dispersion process .    in order to calculate exit - times ,
we prepare a set of thresholds @xmath249 , where @xmath38 is a positive integer , and record the time at which @xmath6 reach @xmath250 for the first time for every particle pairs and every @xmath38 . in the present work , we set @xmath153 to @xmath251 @xmath252 and @xmath253 to the grid size @xmath80 .    if palm s equation , eq .
( [ eq : palm - eq ] ) , describes relative dispersion well , from the equation we can calculate the probability density function ( pdf ) of exit - time @xmath254 , which represents the probability density for a particle pair of which exit - time from @xmath244 to @xmath248 is @xmath255 , @xcite . in the present case , to calculate the pdf , first we solve eq .
( [ eq : palm - eq ] ) with the initial condition @xmath256 where the boundary conditions are the reflecting condition at @xmath257 and the absorbing condition at @xmath258 .
and then , the pdf is obtained as the time derivative of the probability of @xmath259 : @xmath260 we have numerically calculated the pdf of exit - time by the above procedure and compared with the pdf obtained by our dns .